<span style="color: rgb(55, 65, 81); font-family: Söhne, ui-sans-serif, system-ui, -apple-system, "Segoe UI", Roboto, Ubuntu, Cantarell, "Noto Sans", sans-serif, "Helvetica Neue", Arial, "Apple Color Emoji", "Segoe UI Emoji", "Segoe UI Symbol", "Noto Color Emoji"; font-size: 16px; white-space-collapse: preserve;">I think you should give it a try. There's also a rough approach - open the coordinate file in an editor, copy it to a spreadsheet program like Excel, find the most negative value in each column representing its coordinate, and add that value to all the column values. This way, your coordinates will shift in a way that allows you to accurately understand the physical dimensions of the cell. I did this when I couldn't center it through visualization programs.</span><br /><div class="gmail_quote"><div dir="auto" class="gmail_attr">пятница, 26 января 2024 г. в 22:16:40 UTC+3, Ardalan Hayatifar: <br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Great catch!<br>In VMD and using the same exact cell dimensions:<br>```<br>pbc set {18.3675 18.3675 20.1312} -all<br>pbc box<br>```<br><br><img alt="Screen Shot 2024-01-26 at 1.11.49 PM.png" width="220px" height="188px" src="https://groups.google.com/group/cp2k/attach/33577d125894d/Screen%20Shot%202024-01-26%20at%201.11.49%20PM.png?part=0.1&view=1"><div>it is off in XY plane.<br><br>It seems that in ASE-PyCP2K,<br>```<br>structure = read('thickness-1-A_1.cif')<br>structure.center(axis=2)<br>```<br>is not centering anything.<br><br>Should I try using `&CENTER_COORDINATES` to center it?<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, January 26, 2024 at 12:59:00 PM UTC-6 Станислав Бачурин 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"><span style="color:rgb(55,65,81);font-family:Söhne,ui-sans-serif,system-ui,-apple-system,"Segoe UI",Roboto,Ubuntu,Cantarell,"Noto Sans",sans-serif,"Helvetica Neue",Arial,"Apple Color Emoji","Segoe UI Emoji","Segoe UI Symbol","Noto Color Emoji";font-size:16px">Could you check if your system is aligned in the xy plane? Maybe, before applying xy periodicity, it needs to be properly oriented.</span><br>
<span style="color:rgb(55,65,81);font-family:Söhne,ui-sans-serif,system-ui,-apple-system,"Segoe UI",Roboto,Ubuntu,Cantarell,"Noto Sans",sans-serif,"Helvetica Neue",Arial,"Apple Color Emoji","Segoe UI Emoji","Segoe UI Symbol","Noto Color Emoji";font-size:16px">Also, it is necessary to set the dimensions along the Z-axis with a slight margin in the cell settings until this warning disappears. I struggled a lot with calculations under gas-phase conditions. </span>
<span style="color:rgb(55,65,81);font-family:Söhne,ui-sans-serif,system-ui,-apple-system,"Segoe UI",Roboto,Ubuntu,Cantarell,"Noto Sans",sans-serif,"Helvetica Neue",Arial,"Apple Color Emoji","Segoe UI Emoji","Segoe UI Symbol","Noto Color Emoji";font-size:16px">I had exactly the same strange behavior - the nuclei almost merged in the organic molecule. When I figured out the cell settings and Poisson, everything was optimized correctly.</span><div class="gmail_quote"><div dir="auto" class="gmail_attr">пятница, 26 января 2024 г. в 21:49:57 UTC+3, Ardalan Hayatifar: <br></div></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Thanks for your comment. I have been trying `PERIODIC XYZ` or `PERIODIC XY` and that doesn't fix it either. <br>Having more layers also shows the same behavior.<br><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, January 26, 2024 at 12:46:41 PM UTC-6 Станислав Бачурин 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"><span style="display:block;font-family:"Courier New",Courier,monospace,arial,sans-serif;margin:0px;color:rgb(0,0,0)">*** WARNING in pw_env_methods.F:738 :: The selected periodicities in the ***
*** sections &CELL and &POISSON do not match ***
You set periodicity to None, but use the periodic Poisson solver. Try wavelet or switch to Periodic in your Cell's settings.</span><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">пятница, 26 января 2024 г. в 21:38:47 UTC+3, Ardalan Hayatifar: <br></div></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<p>I am using CP2K version 2023.2. I want to perform `GEO_OPT`on FeS layers that I have made using Pymatgen. I am using .cif file from published data.</p><p>
</p><p><b>Input File:</b></p><p>Note: `MAX_SCF` set to 100 results in the same behavior.</p><p>(attached)</p><p><b>Problem:</b></p><p>This is a convergence problem. Upon visualization I can see that things get closer than they should and the system acts crazy.</p><p><b>Pictures:</b></p><p>Initial:</p><p><img alt="Screen Shot 2024-01-26 at 10.46.20 AM.png" width="213px" height="182px" src="https://groups.google.com/group/cp2k/attach/33366aeb3f013/Screen%20Shot%202024-01-26%20at%2010.46.20%20AM.png?part=0.5&view=1"><br></p><p><br></p><p>After some SCF cycels:</p><p><img alt="Screen Shot 2024-01-26 at 10.46.45 AM.png" width="213px" height="169px" src="https://groups.google.com/group/cp2k/attach/33366aeb3f013/Screen%20Shot%202024-01-26%20at%2010.46.45%20AM.png?part=0.1&view=1"><br></p><p>Another representation:</p><p><img alt="Screen Shot 2024-01-26 at 10.10.58 AM.png" width="213px" height="102px" src="https://groups.google.com/group/cp2k/attach/33366aeb3f013/Screen%20Shot%202024-01-26%20at%2010.10.58%20AM.png?part=0.2&view=1"><br></p><p>(shown atoms distance ~1A!)</p><p><b>What I’ve tried:</b></p><p><br></p><p>- I have changed the slab size in X and Y directions to remove size effects. Tried dimensions that are at least 1X1 A and 1.8 X1.8 A dimensions - all while having 15 A vacuum in Z</p><p>- I have monitored the input file for close contacts using `CHECK_INTERATOMIC_DISTANCES 1`. There are no pair-wise distances below 2.5 A.</p><p>- PBC effects: I have tried adding 1 A to cell dimensions `ABC` to again eliminate possible overlaps, but the behavior persists.</p><p>- I have tried higher cutoff values (e.g. 1200).</p><p>- I have tried writing the input so that it uses a `.xyz` file instead of PyCP2K-ASE format and that doesn’t change anything.</p><p>- Centering the solid layer using `&CENTER_COORDINATES` also did not help.</p><p><br></p><p>At this point I’d like to ask for any help or suggestions. It it possible that my parameters are wrong?</p><p><br></p><p>Thanks in advance for your time.</p><p><br></p><p>
</p><p><br></p></blockquote></div></blockquote></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="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/500b956d-00d5-48c0-8042-46420ea7040fn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/500b956d-00d5-48c0-8042-46420ea7040fn%40googlegroups.com</a>.<br />