<div dir="auto">Thank you... It's really helpful. </div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, 19 Oct, 2020, 9:47 pm Lucas Lodeiro, <<a href="mailto:eluni...@gmail.com">eluni...@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi,<div><br></div><div>You have to make pdb files for each molecule that you need to put in a box via packmol. If you can make a water-protein cell, the water-N2 cell is the same, just change the protein.pdb file to N2.pdb and rum packmol.</div><div>Here are some useful examples to make differents types of boxes and interfaces: <a href="http://m3g.iqm.unicamp.br/packmol/examples.shtml" target="_blank" rel="noreferrer">http://m3g.iqm.unicamp.br/packmol/examples.shtml</a></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El lun., 19 oct. 2020 a las 8:02, sumit agrawal (<<a href="mailto:sumitag...@gmail.com" target="_blank" rel="noreferrer">sumitag...@gmail.com</a>>) escribió:<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">Hii,<div><br></div><div>I tried to make a protein, containing a solvation cell via packmol. Here the solvation cell is water. Now I want to change the solvation cell. I want a N2 molecule, containing a solvation cell. How to do it via packmol? Since in packmol only water pdb file is available.</div><div><br></div><div>I am attaching a input file of a protein, containing water molecules.</div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Oct 17, 2020 at 8:57 AM sumit agrawal <<a href="mailto:sumitag...@gmail.com" target="_blank" rel="noreferrer">sumitag...@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">Okay....Thanks a lot!!!</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Oct 17, 2020 at 1:25 AM Travis <<a href="mailto:polla...@gmail.com" target="_blank" rel="noreferrer">polla...@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">Hi,<div><br></div><div>Packmol simply packs a box of some size with a randomized geometry. There will likely be overlapping atoms and other bad contacts. You'll need to <i>at least</i> optimize it but preferably equilibrate the liquid structure with a reasonable force field before using it with CP2K.</div><div><br></div><div>-T<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Saturday, October 17, 2020 at 2:57:32 AM UTC-4 <a href="mailto:sum...@gmail.com" target="_blank" rel="noreferrer">sum...@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">Dear Lucas Loderio and Travis<div><br></div><div>Thanks for the suggestions. I have tried packmol today. In packmol I can easily make my input geometry and it is very easy to use. I have given a tolerance value equal to 2 for my above system. Now I can use this input file for cp2k calculation. </div><div><br></div><div>I have one more question: Should I need to perform geo-opt before computing the vibrational spectra?</div><div><br></div><div>Thanks a lot !!!<br><br><table cellpadding="0" style="border-collapse:collapse;margin-top:0px;width:auto;font-family:Roboto,RobotoDraft,Helvetica,Arial,sans-serif;font-size:14px;letter-spacing:0.2px;display:block"><tbody style="display:block"></tbody></table></div></div><br><div class="gmail_quote"></div><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Oct 16, 2020 at 6:54 PM Travis <<a rel="nofollow noreferrer">po...@gmail.com</a>> wrote:<br></div></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi,<div><br></div><div>To add, as a general practice for condensed phase simulations: build a box with Packmol and equilibrate it with a classical force field that does reasonably well. Use several well-separated snapshots as initial configurations for CP2K. For a simulation with a simple solute in it, well-separated might mean configurations taken every 2*tau or 3*tau where tau is the solvent-solute residence time. This can be short for weakly interacting species but many nanoseconds for very strongly interacting species or viscous solutions. You may even need to simulate at elevated temperatures to reduce the residence time. If you don't know how to compute that from the output of your preferred molecular dynamics code, I would recommend LAMMPS for which there are many community-developed analysis tools. Ed Maginn's group at Notre Dame has a well-rounded analysis suite hosted on Github... PyLat or something like that. There's no unique answer for what 'well-separated' means, this is just the criteria we use.</div><div><br></div><div>For this particular system, since it sounds like this might still be an exercise, you can take one of the CP2K example files for (H2O)128, convert one of the H2O to methanol using Avogadro or GaussView, and run from there. If the concentration really matters for your case, then go the Packmol and classical force field route.</div><div><br></div><div>-T<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, October 16, 2020 at 1:59:54 PM UTC-4 Lucas Lodeiro 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">Dear Sumit Agrawal,<div><br></div><div>I suggest the Packmol software to construct initial configurations, It is a simple and tractable, where you can decide where put some particular molecules, the siz of the box, the number of different molecules, and there is a criteria to ensure your molecules are not so near between them. For little boxes is useful, but for big systems is mandatory to use a program like Packmol.</div><div><a href="http://m3g.iqm.unicamp.br/packmol/home.shtml" rel="nofollow noreferrer" target="_blank">http://m3g.iqm.unicamp.br/packmol/home.shtml</a><br></div><div><br></div><div><br></div><div>Regards</div><div><br></div></div><br><div class="gmail_quote"></div><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El vie., 16 oct. 2020 a las 14:20, Vasyl Stotskyi (<<a rel="nofollow noreferrer">vas...@gmail.com</a>>) escribió:<br></div></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div>Dear Sumit Agrawal,</div><div>Sure, if you can save it in one of avogadro readable formats (.xyz, .cif, ....). </div><div>If initial geometries of molecules in gaus view will be fine, then you don't have to optimize it in Avogadro, i.e. if you will not have O-H bonds with more than 1.5 angstrom length or overlap of molecules, or huge vacuum between molecules. Avogadro optimization is the most useful only for such bad cases.</div><div>Best regards,</div><div>Vasyl</div><br><div>On Oct 17 2020, at 12:40 am, sumit agrawal <<a rel="nofollow noreferrer">sum...@gmail.com</a>> wrote:</div><blockquote><div><div>Dear Vasyl,</div><div><br></div><div>Thanks for the answer. So I have to manually draw 100 molecules of water and methanol in Avogadro? Actually, I am more comfortable with gauss views. So can i make a coordinate file in gusse view then pre optimize in the avogadro ? </div><div><br></div></div><br><img alt="Sent from Mailspring" width="0" height="0" style="border:0px;width:0px;min-height:0px" src="https://link.getmailspring.com/open/E1FD0EA3-5F69-4894-B38E-7FA74A2088EE@getmailspring.com?me=0f52c36e&recipient=Y3Aya0Bnb29nbGVncm91cHMuY29t"><div class="gmail_quote"><div class="gmail_attr"><div>On Fri, Oct 16, 2020 at 9:04 AM Vasyl Stotskyi <<a title="mailto:vasyl....@gmail.com" rel="nofollow noreferrer">vas...@gmail.com</a>> wrote:</div></div><blockquote><div>Dear Sumit Agrawal,</div><div>You can create some model in Avogadro. There you can pre-optimize it with force fields, and save as .xyz file. In CP2K, then you can read coordinates from .xyz file. Please, read manual for this (<a href="https://link.getmailspring.com/link/FC05854F-A2E3-45BD-9677-A12F225F42A8@getmailspring.com/0?redirect=https%3A%2F%2Fmanual.cp2k.org%2Fcp2k-7_1-branch%2FCP2K_INPUT%2FFORCE_EVAL%2FSUBSYS%2FTOPOLOGY.html&recipient=Y3Aya0Bnb29nbGVncm91cHMuY29t" title="https://link.getmailspring.com/link/FC05854F-A2E3-45BD-9677-A12F225F42A8@getmailspring.com/0?redirect=https%3A%2F%2Fmanual.cp2k.org%2Fcp2k-7_1-branch%2FCP2K_INPUT%2FFORCE_EVAL%2FSUBSYS%2FTOPOLOGY.html&recipient=Y3Aya0Bnb29nbGVncm91cHMuY29t" rel="nofollow noreferrer" target="_blank">https://manual.cp2k.org/cp2k-7_1-branch/CP2K_INPUT/FORCE_EVAL/SUBSYS/TOPOLOGY.html</a>).</div><div>Best regards,</div><div>Vasyl</div><div>On Oct 16 2020, at 11:27 pm, sumit agrawal <<a title="mailto:sumitag...@gmail.com" rel="nofollow noreferrer">sum...@gmail.com</a>> wrote:</div><blockquote><div><div>Hii,</div><div><br></div><div>I have recently installed cp2k in windows. I just did some exercise which is available on cp2k website. Now I want to do some bulk phase vibrational analysis. I want to make a coordinate file in which 100 water molecules and one methanol molecule is present. This methanol molecule is present in between water molecules. So how to make a coordinate file for this system.</div><div><br></div><div>Can I use gauss view to make the coordinate file?</div></div><br><div>--</div><div>You received this message because you are subscribed to the Google Groups "cp2k" group.</div><div>To unsubscribe from this group and stop receiving emails from it, send an email to <a href="https://link.getmailspring.com/link/FC05854F-A2E3-45BD-9677-A12F225F42A8@getmailspring.com/1?redirect=mailto%3Acp2k%2Bunsubscribe%40googlegroups.com&recipient=Y3Aya0Bnb29nbGVncm91cHMuY29t" title="https://link.getmailspring.com/link/FC05854F-A2E3-45BD-9677-A12F225F42A8@getmailspring.com/1?redirect=mailto%3Acp2k%2Bunsubscribe%40googlegroups.com&recipient=Y3Aya0Bnb29nbGVncm91cHMuY29t" rel="nofollow noreferrer" target="_blank">cp...@googlegroups.com</a>.</div><div>To view this discussion on the web visit <a href="https://link.getmailspring.com/link/FC05854F-A2E3-45BD-9677-A12F225F42A8@getmailspring.com/2?redirect=https%3A%2F%2Fgroups.google.com%2Fd%2Fmsgid%2Fcp2k%2FCALSDoYad_FGm4Qh5SRvZcMun3fe%253D5q1xpxeW3T5SoMcAoTH52A%2540mail.gmail.com%3Futm_medium%3Demail%26utm_source%3Dfooter&recipient=Y3Aya0Bnb29nbGVncm91cHMuY29t" title="https://link.getmailspring.com/link/FC05854F-A2E3-45BD-9677-A12F225F42A8@getmailspring.com/2?redirect=https%3A%2F%2Fgroups.google.com%2Fd%2Fmsgid%2Fcp2k%2FCALSDoYad_FGm4Qh5SRvZcMun3fe%253D5q1xpxeW3T5SoMcAoTH52A%2540mail.gmail.com%3Futm_medium%3Demail%26utm_source%3Dfooter&recipient=Y3Aya0Bnb29nbGVncm91cHMuY29t" rel="nofollow noreferrer" target="_blank">https://groups.google.com/d/msgid/cp2k/CALSDoYad_FGm4Qh5SRvZcMun3fe%3D5q1xpxeW3T5SoMcAoTH52A%40mail.gmail.com</a>.</div></blockquote><br><div>--</div><div>You received this message because you are subscribed to the Google Groups "cp2k" group.</div><div>To unsubscribe from this group and stop receiving emails from it, send an email to <a title="mailto:cp...@googlegroups.com" rel="nofollow noreferrer">cp...@googlegroups.com</a>.</div><div>To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/FC05854F-A2E3-45BD-9677-A12F225F42A8%40getmailspring.com?utm_medium=email&utm_source=footer" title="https://groups.google.com/d/msgid/cp2k/FC05854F-A2E3-45BD-9677-A12F225F42A8%40getmailspring.com?utm_medium=email&utm_source=footer" rel="nofollow noreferrer" target="_blank">https://groups.google.com/d/msgid/cp2k/FC05854F-A2E3-45BD-9677-A12F225F42A8%40getmailspring.com</a>.</div></blockquote></div><br><div>--</div><div>You received this message because you are subscribed to the Google Groups "cp2k" group.</div><div>To unsubscribe from this group and stop receiving emails from it, send an email to <a href="https://link.getmailspring.com/link/E1FD0EA3-5F69-4894-B38E-7FA74A2088EE@getmailspring.com/0?redirect=mailto%3Acp2k%2Bunsubscribe%40googlegroups.com&recipient=Y3Aya0Bnb29nbGVncm91cHMuY29t" title="mailto:cp...@googlegroups.com" rel="nofollow noreferrer" target="_blank">cp...@googlegroups.com</a>.</div><div>To view this discussion on the web visit <a href="https://link.getmailspring.com/link/E1FD0EA3-5F69-4894-B38E-7FA74A2088EE@getmailspring.com/1?redirect=https%3A%2F%2Fgroups.google.com%2Fd%2Fmsgid%2Fcp2k%2FCALSDoYaSwE%253DyxwswTkq6U5OouRZPiuev3yRJDydn%252B37kF-G4%252Bg%2540mail.gmail.com%3Futm_medium%3Demail%26utm_source%3Dfooter&recipient=Y3Aya0Bnb29nbGVncm91cHMuY29t" title="https://groups.google.com/d/msgid/cp2k/CALSDoYaSwE%3DyxwswTkq6U5OouRZPiuev3yRJDydn%2B37kF-G4%2Bg%40mail.gmail.com?utm_medium=email&utm_source=footer" rel="nofollow noreferrer" target="_blank">https://groups.google.com/d/msgid/cp2k/CALSDoYaSwE%3DyxwswTkq6U5OouRZPiuev3yRJDydn%2B37kF-G4%2Bg%40mail.gmail.com</a>.</div></blockquote>
<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 rel="nofollow noreferrer">cp...@googlegroups.com</a>.<br></blockquote></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/E1FD0EA3-5F69-4894-B38E-7FA74A2088EE%40getmailspring.com?utm_medium=email&utm_source=footer" rel="nofollow noreferrer" target="_blank">https://groups.google.com/d/msgid/cp2k/E1FD0EA3-5F69-4894-B38E-7FA74A2088EE%40getmailspring.com</a>.<br>
</blockquote></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 rel="nofollow noreferrer">cp...@googlegroups.com</a>.<br></blockquote></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/ee5fea3f-2f5c-42f2-8fb3-c6ab1042e860n%40googlegroups.com?utm_medium=email&utm_source=footer" rel="nofollow noreferrer" target="_blank">https://groups.google.com/d/msgid/cp2k/ee5fea3f-2f5c-42f2-8fb3-c6ab1042e860n%40googlegroups.com</a>.<br>
</blockquote></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:cp...@googlegroups.com" target="_blank" rel="noreferrer">cp...@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/24cb5e6e-ba82-45ef-be8f-90b97c9f3253n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="noreferrer">https://groups.google.com/d/msgid/cp2k/24cb5e6e-ba82-45ef-be8f-90b97c9f3253n%40googlegroups.com</a>.<br>
</blockquote></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:cp...@googlegroups.com" target="_blank" rel="noreferrer">cp...@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CALSDoYY1zZUFK0OxSNi-Dqkh9wMfw-_qo1%2BQkM%3DrufvRK09e8w%40mail.gmail.com?utm_medium=email&utm_source=footer" target="_blank" rel="noreferrer">https://groups.google.com/d/msgid/cp2k/CALSDoYY1zZUFK0OxSNi-Dqkh9wMfw-_qo1%2BQkM%3DrufvRK09e8w%40mail.gmail.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" rel="noreferrer">cp...@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAOFT4PLme_F95CUF0BjDyHE%3D%3DPSkB9DKeRR0c-FgeVV-6fsi4w%40mail.gmail.com?utm_medium=email&utm_source=footer" target="_blank" rel="noreferrer">https://groups.google.com/d/msgid/cp2k/CAOFT4PLme_F95CUF0BjDyHE%3D%3DPSkB9DKeRR0c-FgeVV-6fsi4w%40mail.gmail.com</a>.<br>
</blockquote></div>