<div dir="ltr"><div>Thanks, <span style="color:rgb(38,40,42);font-family:"Helvetica Neue",Helvetica,Arial,sans-serif;font-size:13px">Fangyong for the suggestion.</span></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Jan 15, 2021 at 10:01 PM Fangyong Yan <<a href="mailto:fyya...@gmail.com">fyya...@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi, Monu,<div><br></div><div>The difference energy is <span style="color:rgb(0,0,0);font-family:"Helvetica Neue",Helvetica,Arial,sans-serif;font-size:13px">-6235.818 - (-</span><span style="font-family:"Helvetica Neue",Helvetica,Arial,sans-serif;font-size:13px">6235.824) = 0.006 Hartree = 15.753 kj/mol. </span></div><div><span style="font-family:"Helvetica Neue",Helvetica,Arial,sans-serif;font-size:13px"><br></span></div><div><span style="font-family:"Helvetica Neue",Helvetica,Arial,sans-serif;font-size:13px">So if you use 800 K, the Boltzmann factor is: exp(-dE/R/T) = exp(-15.753*1000/8.314/800) = 0.09, which means you would have a probability of 0.09 to "jump" from the lower energy to higher energy, by the help of heat. </span></div><div><span style="font-family:"Helvetica Neue",Helvetica,Arial,sans-serif;font-size:13px"><br></span></div><div><span style="font-family:"Helvetica Neue",Helvetica,Arial,sans-serif;font-size:13px">So I think NVT at 800K, may be helpful for exploring the potential energy surface, even within 10ps. </span></div><div><span style="font-family:"Helvetica Neue",Helvetica,Arial,sans-serif;font-size:13px"><br></span></div><div><span style="font-family:"Helvetica Neue",Helvetica,Arial,sans-serif;font-size:13px">Regards,</span></div><div><span style="font-family:"Helvetica Neue",Helvetica,Arial,sans-serif;font-size:13px"><br></span></div><div><span style="font-family:"Helvetica Neue",Helvetica,Arial,sans-serif;font-size:13px">Fangyong</span></div><div><span style="font-family:"Helvetica Neue",Helvetica,Arial,sans-serif;font-size:13px"><br></span></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Jan 15, 2021 at 8:08 PM Monu Joy <<a href="mailto:monuj...@gmail.com" target="_blank">monuj...@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div style="font-family:verdana,helvetica,sans-serif;font-size:13px"><div><div style="font-size:13px" dir="ltr">Hi Lucas & <span><span style="font-family:"Helvetica Neue",Helvetica,Arial,sans-serif">Fangyong</span></span></div><div style="font-size:13px" dir="ltr"><span><span style="font-family:"Helvetica Neue",Helvetica,Arial,sans-serif"><br></span></span></div><div style="font-size:13px" dir="ltr"><span><span style="font-family:"Helvetica Neue",Helvetica,Arial,sans-serif">Thanks for the help. I have performed two cell-opt calculations of guest@host where the following total energies respectively represent the guest at the center of the MOF pore and near to the pore surface <span><span style="color:rgb(0,0,0);font-family:"Helvetica Neue",Helvetica,Arial,sans-serif">-6235.818 and </span></span></span></span><span style="font-family:"Helvetica Neue",Helvetica,Arial,sans-serif">-6235.824. Do you guys think that this energy difference is reasonable...?</span></div><div style="font-size:13px" dir="ltr"><span style="font-family:"Helvetica Neue",Helvetica,Arial,sans-serif"><br></span></div><div style="font-size:13px" dir="ltr"><span style="font-family:"Helvetica Neue",Helvetica,Arial,sans-serif">Also, do you by chance have any reference that shows local minimum geometry modulate the property of a system more effectively...? </span></div><div style="font-size:13px" dir="ltr"><span style="font-family:"Helvetica Neue",Helvetica,Arial,sans-serif"><br></span></div><div style="font-size:13px" dir="ltr">Thank you so much for your help!</div><div style="font-size:13px" dir="ltr"><br></div><div><div style="font-family:verdana,helvetica,sans-serif"><div><div dir="ltr"><div><div><div dir="ltr"><font size="2">-Monu</font></div></div></div></div></div></div></div></div>
        <div style="font-size:13px"><br></div><div style="font-size:13px"><br></div>
        
        </div><div id="gmail-m_5572969993991021658gmail-m_1209820807176214709yahoo_quoted_1700497702">
            <div style="font-family:"Helvetica Neue",Helvetica,Arial,sans-serif;font-size:13px;color:rgb(38,40,42)">
                
                <div>
                    On Wednesday, 13 January, 2021, 10:01:30 am GMT-5, Lucas Lodeiro <<a href="mailto:eluni...@gmail.com" target="_blank">eluni...@gmail.com</a>> wrote:
                </div>
                <div><br></div>
                <div><br></div>
                <div><div id="gmail-m_5572969993991021658gmail-m_1209820807176214709yiv4822642648"><div><div dir="ltr">Hi Monu and Fangyong,<div><br clear="none"></div><div>I disagree a little bit with Fangyong because an AIMD calculation is not cheaper, and if the barriers between local minimas are high, you cannot see a structural change with the MD, at least in the scale time of ps. Firstly, I would try to optimize a couple of structures and measure the energy differences, this will give you a little insight about the barriers (the minimum energy for the barriers), If they are high, the AIMD probably cannot move from one local minima to another. If they are little, you can see the structural changes, but in this case, probably the global minima is not sufficient to modulate the properties of the system, and some other local minima have to be taken into account.</div><div>Other way is to generate with a script random positions of the molecule into the MOF and minimize them with a semiempirical formulation, as XTB (<a rel="nofollow noopener noreferrer" shape="rect" href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/DFT/QS.html#METHOD" target="_blank">https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/DFT/QS.html#METHOD</a>), and rerun with the level of theory that you want just the low energy ones.</div><div><br clear="none"></div><div>Regards - Lucas Lodeiro</div><div><br clear="none"></div></div><br clear="none"><div id="gmail-m_5572969993991021658gmail-m_1209820807176214709yiv4822642648yqt65762"><div><div dir="ltr">El mié, 13 ene 2021 a las 0:07, Fangyong Yan (<<a rel="nofollow noopener noreferrer" shape="rect" href="mailto:fyya...@gmail.com" target="_blank">fyya...@gmail.com</a>>) escribió:<br clear="none"></div><blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi, Monu,<div><br clear="none"></div><div>To explore the free energy surface better, you may need to start at a high temperature, and then annealing down to the room temperature. And for each temperature, you can plot the radial distribution function, for your guest molecule and your MOF atoms, for example, at 800 K, 700 K, ...., 300K. The reason starting at high temperature is that you can use heat to get your guest molecule out of the potentially trapped free energy minima, so you can explore the free energy minima more completely. </div><div><br clear="none"></div><div>Regards,</div><div><br clear="none"></div><div>Fangyong</div></div><br clear="none"><div><div dir="ltr">On Tue, Jan 12, 2021 at 10:01 PM Fangyong Yan <<a rel="nofollow noopener noreferrer" shape="rect" href="mailto:fyya...@gmail.com" target="_blank">fyya...@gmail.com</a>> wrote:<br clear="none"></div><blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi, Monu,<div><br clear="none"></div><div>You can run a molecular dynamics simulation at room temperature, for several picosecond, and then you can plot the radial distribution function, between your guest molecule and the MOF, by doing this you would have a general picture of the free energy surface of your guest molecule in the MOF. Then based on the free energy minima structures, you can start to optimize your guest-MOF complex structure. (the free energy surface may be different from the potential energy surface, but at least by doing MD simulation and calculating the free energy surface based on radial distribution function, you can have a basic idea about the free energy surface, which should not be quite different from potential energy surface).</div><div><br clear="none"></div><div>Regards,</div><div><br clear="none"></div><div>Fangyong</div></div><br clear="none"><div><div dir="ltr">On Tue, Jan 12, 2021 at 7:20 PM Monu Joy <<a rel="nofollow noopener noreferrer" shape="rect" href="mailto:monuj...@gmail.com" target="_blank">monuj...@gmail.com</a>> wrote:<br clear="none"></div><blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div style="font-family:verdana,helvetica,sans-serif;font-size:13px"><div dir="ltr">Hi there</div><div dir="ltr"><br clear="none"></div><div dir="ltr"><div><div dir="ltr">I have a general question: I would like to find out the global minimum of a guest molecule within the pore of a MOF. I know I can place the guest molecule in the pore and get the host-guest structure optimized, but if I would place the guest molecule in another position of the same pore, there also I can get the structure optimized. In this scenario how could I find the exact most probable location of the guest molecule within the MOF pore or its global minimum...?</div><div><br clear="none"></div><div>Any help would be highly appreciated...</div></div></div><br clear="none"></div><div dir="ltr" style="font-family:verdana,helvetica,sans-serif;font-size:13px">-Monu</div></div>

<p></p>

-- <br clear="none">
You received this message because you are subscribed to the Google Groups "cp2k" group.<br clear="none">
To unsubscribe from this group and stop receiving emails from it, send an email to <a rel="nofollow noopener noreferrer" shape="rect" href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br clear="none">
To view this discussion on the web visit <a rel="nofollow noopener noreferrer" shape="rect" href="https://groups.google.com/d/msgid/cp2k/1852672747.3320.1610497192914%40mail.yahoo.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/1852672747.3320.1610497192914%40mail.yahoo.com</a>.<br clear="none">
</blockquote></div>
</blockquote></div>

<p></p>

-- <br clear="none">
You received this message because you are subscribed to the Google Groups "cp2k" group.<br clear="none">
To unsubscribe from this group and stop receiving emails from it, send an email to <a rel="nofollow noopener noreferrer" shape="rect" href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br clear="none">
To view this discussion on the web visit <a rel="nofollow noopener noreferrer" shape="rect" href="https://groups.google.com/d/msgid/cp2k/CAJjmc%2BKviWAFmj5UmYqMWL1Z5wFpQZjdKkhGS%2B7CTp7-6jaLgg%40mail.gmail.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/CAJjmc%2BKviWAFmj5UmYqMWL1Z5wFpQZjdKkhGS%2B7CTp7-6jaLgg%40mail.gmail.com</a>.<br clear="none">
</blockquote></div></div>

<p></p>

-- <br clear="none">
You received this message because you are subscribed to the Google Groups "cp2k" group.<br clear="none">
To unsubscribe from this group and stop receiving emails from it, send an email to <a rel="nofollow noopener noreferrer" shape="rect" href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br clear="none">
To view this discussion on the web visit <a rel="nofollow noopener noreferrer" shape="rect" href="https://groups.google.com/d/msgid/cp2k/CAOFT4PKCKjc0vOGZUkbjsDEEqJM_w_UxYK-PC%2Bt-vnzZzLhrgg%40mail.gmail.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/CAOFT4PKCKjc0vOGZUkbjsDEEqJM_w_UxYK-PC%2Bt-vnzZzLhrgg%40mail.gmail.com</a>.<br clear="none">
</div></div></div>
            </div>
        </div></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:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/1361749041.8962.1610759305502%40mail.yahoo.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/1361749041.8962.1610759305502%40mail.yahoo.com</a>.<br>
</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:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAJjmc%2BKyP-zV-p%2B%2BZtUkTgqoWQnovqnUvaVpG7NQby7a8KURgQ%40mail.gmail.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/CAJjmc%2BKyP-zV-p%2B%2BZtUkTgqoWQnovqnUvaVpG7NQby7a8KURgQ%40mail.gmail.com</a>.<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><a style="font-family:"comic sans ms",sans-serif" href="mailto:monupl...@yahoo.com" target="_blank"><font style="color:rgb(255,0,0)" size="4"><b>Thanks with Regards</b></font></a><br></div></div>