<div dir="ltr"><div>Hi,</div><div><br></div><div>CP2K prints the unwrapped coordinates. Molecules diffusing out of the box are not wrapped back in. You can view the wrapped trajectory in VMD by opening the TK console,</div><div><br> <div style="background-color: rgb(250, 250, 250); border-color: rgb(187, 187, 187); border-style: solid; border-width: 1px; overflow-wrap: break-word;" class="prettyprint"><code class="prettyprint"><div class="subprettyprint"><span style="color: #000;" class="styled-by-prettify">pbc </span><span style="color: #008;" class="styled-by-prettify">set</span><span style="color: #000;" class="styled-by-prettify"> </span><span style="color: #660;" class="styled-by-prettify">{</span><span style="color: #066;" class="styled-by-prettify">15</span><span style="color: #000;" class="styled-by-prettify"> </span><span style="color: #066;" class="styled-by-prettify">15</span><span style="color: #000;" class="styled-by-prettify"> </span><span style="color: #066;" class="styled-by-prettify">20</span><span style="color: #660;" class="styled-by-prettify">}</span><span style="color: #000;" class="styled-by-prettify"> </span><span style="color: #660;" class="styled-by-prettify">-</span><span style="color: #000;" class="styled-by-prettify">all<br>pbc box<br>pbc wrap </span><span style="color: #660;" class="styled-by-prettify">-</span><span style="color: #000;" class="styled-by-prettify">all<br></span></div></code></div></div><div><br></div><div>If you're trying to perform a bulk calculation, you will need to remove the vacuum region either by CELL_OPT with fixed angles or by filling the void spaces with water.<br></div><div><br></div><div>-T</div><div><br></div><br>On Wednesday, October 16, 2019 at 5:06:25 AM UTC-3, farzaneh sarrami wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><div dir="ltr">Hi Travis,<div><br></div><div>Thanks for complete explanation. Currently I am using NVE for MD part in my system. I tried to use NPT but it doesn't work for my system. So I prefer to stick to NVE as it goes to the right direction. I am still struggling with water molecules moving out of the simulation box. I have enclosed my input file and final coordinates.</div><div>Would you please have a look and give me some advises to fix this problem?</div><div><br></div><div>Many thanks in advance</div><div>Farzaneh</div><div><br></div></div></div><br><div class="gmail_quote"><div dir="ltr">On Tue, Oct 8, 2019 at 5:22 PM Travis <<a href="javascript:" target="_blank" gdf-obfuscated-mailto="eThxRKm_DAAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">po...@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi,</div><div><br></div><div>Compile the code with your preferred Fortran compiler,<br></div><div><div style="background-color:rgb(250,250,250);border:1px solid rgb(187,187,187)"><code><div><span style="color:rgb(0,0,0)">gfortran dumpdcd</span><span style="color:rgb(102,102,0)">.</span><span style="color:rgb(0,0,0)">f90</span></div></code></div></div><div><br></div><div>Run code on a host of DCD files you want to merge, these would be coming from a long simulation that you broke up into smaller segments,<br></div><span style="color:rgb(0,0,0)"></span><div><div style="background-color:rgb(250,250,250);border:1px solid rgb(187,187,187)"><code><div><span style="color:rgb(102,102,0)">./</span><span style="color:rgb(0,0,0)">a</span><span style="color:rgb(102,102,0)">.</span><span style="color:rgb(0,0,136)">out</span><span style="color:rgb(0,0,0)"> </span><span style="color:rgb(102,102,0)">-</span><span style="color:rgb(0,0,0)">xyz reference</span><span style="color:rgb(102,102,0)">.</span><span style="color:rgb(0,0,0)">xyz </span><span style="color:rgb(102,102,0)">-</span><span style="color:rgb(0,0,0)">pbc </span><span style="color:rgb(102,102,0)">-</span><span style="color:rgb(0,0,0)">of dcd </span><span style="color:rgb(102,102,0)">-</span><span style="color:rgb(0,0,0)">o full_foo_trajectory</span><span style="color:rgb(102,102,0)">.</span><span style="color:rgb(0,0,0)">dcd foo1</span><span style="color:rgb(102,102,0)">.</span><span style="color:rgb(0,0,0)">dcd foo2</span><span style="color:rgb(102,102,0)">.</span><span style="color:rgb(0,0,0)">dcd foo3</span><span style="color:rgb(102,102,0)">.</span><span style="color:rgb(0,0,0)">dcd </span><span style="color:rgb(102,102,0)">...</span><span style="color:rgb(0,0,0)"> fooN</span><span style="color:rgb(102,102,0)">.</span><span style="color:rgb(0,0,0)">dcd<br></span></div></code></div></div><div><br></div><div>You need a reference XYZ format file for determining atom name as DCD only contains lattice parameters and coordinates. The -pbc flag says to produce a wrapped trajectory, while leaving it off produces the unwrapped trajectory.<br><br>To wrap an XYZ trajectory requires TRAVIS (<a href="https://www.travis-analyzer.de/" target="_blank" rel="nofollow" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fwww.travis-analyzer.de%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNGUfpiCJQRusRIjLI5Xz5U7ir1SnQ';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fwww.travis-analyzer.de%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNGUfpiCJQRusRIjLI5Xz5U7ir1SnQ';return true;">https://www.travis-analyzer.<wbr>de/</a>). In one of the advanced options, you can read the second line for each frame of the XYZ file as the lattice parameters in Angstroms and degrees. I give a more detailed explanation of it here, <a href="https://groups.google.com/forum/#!topic/cp2k/GNLdTcUZFUo" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/forum/#!topic/cp2k/GNLdTcUZFUo';return true;" onclick="this.href='https://groups.google.com/forum/#!topic/cp2k/GNLdTcUZFUo';return true;">https://groups.google.com/<wbr>forum/#!topic/cp2k/GNLdTcUZFUo</a> (for NPT). Alternatively, if you are comfortable with Python, you can use the Atomic Simulation Environment.</div><div><br></div><div>-T<br></div><br>On Tuesday, October 8, 2019 at 4:23:27 AM UTC-3, farzaneh sarrami wrote:<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr">Thanks Travis,<div><br></div><div>Would you please explain in more details how can I use dumpdcd.f90 script in cp2k ? as I have same problem.</div><div><br></div><div>Best,</div><div>Farzaneh</div><br>On Monday, October 7, 2019 at 5:01:01 PM UTC+2, Travis wrote:<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi,</div><div><br></div><div>CP2K dumps unwrapped coordinates, so they are not imaged back into the box if they drift beyond. A simple way to get a wrapped trajectory is to dump DCD files during MD, then use the dumpdcd.f90 script in tools/ to dump a wrapped DCD (-pbc option). I don't recommend writing XYZ format trajectories for MD. Write DCD and convert. Way easier to strip info than to add it back in.<br></div><div><br></div><div>-T</div><div><br></div><br>On Monday, October 7, 2019 at 8:09:09 AM UTC-3, Harender Singh wrote:<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Dear all, <br></div><div>I am running AIMD simulation of some electrolyte at Li electrode. While the electrolyte molecules are staying within the simulation cell, the electrode atoms are moving out of the cell. <br></div><div>Is there any possible explanation for that?</div><div>Is there something wrong with my simulation setup? If yes, then how can I correct it?</div><div><br></div><div>Thanks for the help in advance. <br></div></div></blockquote></div></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="javascript:" target="_blank" gdf-obfuscated-mailto="eThxRKm_DAAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.<wbr>com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/4b2a5d18-4e0d-43e4-8bd8-e2c906f92cfe%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/msgid/cp2k/4b2a5d18-4e0d-43e4-8bd8-e2c906f92cfe%40googlegroups.com?utm_medium\x3demail\x26utm_source\x3dfooter';return true;" onclick="this.href='https://groups.google.com/d/msgid/cp2k/4b2a5d18-4e0d-43e4-8bd8-e2c906f92cfe%40googlegroups.com?utm_medium\x3demail\x26utm_source\x3dfooter';return true;">https://groups.google.com/d/<wbr>msgid/cp2k/4b2a5d18-4e0d-43e4-<wbr>8bd8-e2c906f92cfe%<wbr>40googlegroups.com</a>.<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr"><table style="max-width:600px;direction:ltr">
<tbody><tr><td>
<table style="max-width:470px;padding-bottom:10px;margin-bottom:8px" cellspacing="0" cellpadding="0" border="0"> <tbody><tr> <td> <table style="width:470px" width="470" cellspacing="0" cellpadding="0" border="0"> <tbody><tr valign="top"> <td style="text-align:initial;font-style:normal;font-weight:normal;font-size:12px;line-height:normal;font-family:Arial;color:rgb(100,100,100);padding:0px 10px"> <div style="margin-bottom:5px;margin-top:0px"><b>Farzaneh Sarrami</b><br> <span>Post-doctoral Fellow</span> </div> <table style="width:470px;margin-top:5px" width="470" cellspacing="0" cellpadding="0" border="0"> <tbody><tr> <td style="color:rgb(141,141,141);font-size:12px"> <p style="margin:0px"> <a style="color:rgb(141,141,141);text-decoration:none;font-family:sans-serif">Laboratory of Organic Electronic</a> </p> </td> </tr> <tr> <td style="color:rgb(141,141,141);font-size:12px;font-family:sans-serif"> <span style="color:rgb(141,141,141);display:inline-block;font-family:sans-serif">Linköping University</span> </td> </tr> </tbody></table> <div style="margin-top:10px"><table cellspacing="0" cellpadding="0" border="0"><tbody><tr style="padding-top:10px"></tr></tbody></table></div> </td> </tr> </tbody></table> </td> </tr> </tbody></table> <table cellspacing="0" cellpadding="0" border="0"><tbody><tr><td></td></tr></tbody></table></td></tr></tbody></table>
<div href="http://WS_promo" style="width:auto;padding-top:2px;font-size:10px;border-top-width:1px;border-top-style:solid;border-top-color:rgb(238,238,238);margin-top:10px;display:table;direction:ltr;line-height:normal;border-spacing:initial">
<div style="padding-top:2px">
<a href="https://www.wisestamp.com/signature-in-email?utm_source=promotion&utm_medium=signature&utm_campaign=create_your_own" style="color:rgb(0,0,0);text-decoration:none" target="_blank" rel="nofollow" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fwww.wisestamp.com%2Fsignature-in-email%3Futm_source%3Dpromotion%26utm_medium%3Dsignature%26utm_campaign%3Dcreate_your_own\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHS5_thIouRUzqhfA2xi8my7FejAQ';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fwww.wisestamp.com%2Fsignature-in-email%3Futm_source%3Dpromotion%26utm_medium%3Dsignature%26utm_campaign%3Dcreate_your_own\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHS5_thIouRUzqhfA2xi8my7FejAQ';return true;">Create your own <span style="color:rgb(165,3,16)">email signature</span></a>
</div>
</div>
</div>
</blockquote></div>