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 sum...@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;"><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 href data-email-masked rel="nofollow">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" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en&q=http://m3g.iqm.unicamp.br/packmol/home.shtml&source=gmail&ust=1603009100532000&usg=AFQjCNH_VUqZrEzJMCKr2sOE0yakJC6zzQ">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">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">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://ci4.googleusercontent.com/proxy/_FWANoizfRrF_IwVZ4hsWjOZYXPrnZNeC_y0vz1MXDrEor0SdjaiHAIeTEnEABfj_mjiUDCf_GrtPhGkbfn_rgtrSONJjRl6wRm_WO52zS3hDn261657Aiap9LFF-odI4D8b5t-F9yvXLn9ceMlBOdbV2ySfFoSutsD504ROkOQhNCmCKYWXvukrwr2Xq8_LmSYurhCvPSXZPSrtc_wkXkc=s0-d-e1-ft#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">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" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://link.getmailspring.com/link/FC05854F-A2E3-45BD-9677-A12F225F42A8@getmailspring.com/0?redirect%3Dhttps%253A%252F%252Fmanual.cp2k.org%252Fcp2k-7_1-branch%252FCP2K_INPUT%252FFORCE_EVAL%252FSUBSYS%252FTOPOLOGY.html%26recipient%3DY3Aya0Bnb29nbGVncm91cHMuY29t&source=gmail&ust=1603009100532000&usg=AFQjCNG6aQiPkqzQp5s-odZ1MkFE53nuog">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">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" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://link.getmailspring.com/link/FC05854F-A2E3-45BD-9677-A12F225F42A8@getmailspring.com/1?redirect%3Dmailto%253Acp2k%252Bunsubscribe%2540googlegroups.com%26recipient%3DY3Aya0Bnb29nbGVncm91cHMuY29t&source=gmail&ust=1603009100533000&usg=AFQjCNF_u9OL0yWdOv-UC1FWJcld_7MNIQ">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" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://link.getmailspring.com/link/FC05854F-A2E3-45BD-9677-A12F225F42A8@getmailspring.com/2?redirect%3Dhttps%253A%252F%252Fgroups.google.com%252Fd%252Fmsgid%252Fcp2k%252FCALSDoYad_FGm4Qh5SRvZcMun3fe%25253D5q1xpxeW3T5SoMcAoTH52A%252540mail.gmail.com%253Futm_medium%253Demail%2526utm_source%253Dfooter%26recipient%3DY3Aya0Bnb29nbGVncm91cHMuY29t&source=gmail&ust=1603009100533000&usg=AFQjCNHcC94V6MdhBa0veH1uxbFvtfBUUQ">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">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" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/FC05854F-A2E3-45BD-9677-A12F225F42A8%2540getmailspring.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1603009100533000&usg=AFQjCNGJvPbCpr2aKDOt0o9B6r0imVEU9w">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" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://link.getmailspring.com/link/E1FD0EA3-5F69-4894-B38E-7FA74A2088EE@getmailspring.com/0?redirect%3Dmailto%253Acp2k%252Bunsubscribe%2540googlegroups.com%26recipient%3DY3Aya0Bnb29nbGVncm91cHMuY29t&source=gmail&ust=1603009100533000&usg=AFQjCNGFFb5uAn5CNe8SG5-Lvm7NU4mI2A">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" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://link.getmailspring.com/link/E1FD0EA3-5F69-4894-B38E-7FA74A2088EE@getmailspring.com/1?redirect%3Dhttps%253A%252F%252Fgroups.google.com%252Fd%252Fmsgid%252Fcp2k%252FCALSDoYaSwE%25253DyxwswTkq6U5OouRZPiuev3yRJDydn%25252B37kF-G4%25252Bg%252540mail.gmail.com%253Futm_medium%253Demail%2526utm_source%253Dfooter%26recipient%3DY3Aya0Bnb29nbGVncm91cHMuY29t&source=gmail&ust=1603009100533000&usg=AFQjCNGEDxFpV39MfgzjbQ20otY_Q0i0Lw">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">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" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/E1FD0EA3-5F69-4894-B38E-7FA74A2088EE%2540getmailspring.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1603009100533000&usg=AFQjCNFIiW537LmfYhZi8ijFQZrV3CsiEQ">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 href data-email-masked rel="nofollow">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" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/ee5fea3f-2f5c-42f2-8fb3-c6ab1042e860n%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1603009100534000&usg=AFQjCNHa5bHpzOVT032lgPJ4RUWt4sGiEw">https://groups.google.com/d/msgid/cp2k/ee5fea3f-2f5c-42f2-8fb3-c6ab1042e860n%40googlegroups.com</a>.<br> </blockquote></div> </blockquote></div>