<div dir="ltr">Hi Nico, <br><br>Thanks for your patience and detailed reply! It's pretty useful for me! I'll try it.<br><br>在 2017年12月1日星期五 UTC+8下午6:31:14，Nico Holmberg写道：<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hi,<br><br>You need to supply a cube file that contains as many grid points as the grid that holds the real space electron density in your simulation system (depends on the PW_CUTOFF and cell size). You can check the size of this grid by turning on the print key FORCE_EVAL/PRINT/GRID_<wbr>INFORMATION (see example below). When you print out a cube file, unless you explicitly specify otherwise, only every second value along each axis is printed out (keyword STRIDE). You need to take this into account when you supply a file for READ_FROM_CUBE, because as far as I understand no interpolation/extrapolation from tighter/looser grids is possible.<br><br>An example output is given below. In this case, the grid contains 350^3 points at a planewave cutoff of 500 Ry. So to summarize, if your system contains 54 data points along each direction (at STRIDE 1), then you need to provide a cube with 54^3 grid points. You need a larger cube file with (2*54)^3 grid points if you printed out the cube file with STRIDE 2.<br><br> PW_GRID| Information for grid number                        <wbr>                  1<br> PW_GRID| Grid distributed over                          <wbr>         192 processors<br> PW_GRID| Real space group dimensions                    <wbr>               192    1<br> PW_GRID| the grid is blocked:                      <wbr>                          NO<br> PW_GRID| Cutoff [a.u.]                        <wbr>                            250.0<br> PW_GRID| spherical cutoff:                       <wbr>                            NO<br> PW_GRID|   Bounds   1           -175     174                Points:         350<br> PW_GRID|   Bounds   2           -175     174                Points:         350<br> PW_GRID|   Bounds   3           -175     174                Points:         350<br> PW_GRID| Volume element (a.u.^3)  0.2459E-02     Volume (a.u.^3)    105442.7279<br> PW_GRID| Grid span                          <wbr>                          FULLSPACE<br> PW_GRID|   Distribution                  <wbr>       Average         Max         Min<br> PW_GRID|   G-Vectors                     <wbr>      223307.3      223650      223300<br> PW_GRID|   G-Rays                        <wbr>         638.0         639         638<br> PW_GRID|   Real Space Points                   223307.3      245000      122500<br><br> RS_GRID| Information for grid number                        <wbr>                  1<br> RS_GRID|   Bounds   1           -175     174                Points:         350<br> RS_GRID|   Bounds   2           -175     174                Points:         350<br> RS_GRID|   Bounds   3           -175     174                Points:         350<br> RS_GRID| Real space distribution over                          <wbr>        4 groups<br> RS_GRID| Real space distribution along direction                     <wbr>         1<br> RS_GRID| Border size                          <wbr>                               29<br> RS_GRID| Real space distribution over                          <wbr>        6 groups<br> RS_GRID| Real space distribution along direction                     <wbr>         2<br> RS_GRID| Border size                          <wbr>                               29<br> RS_GRID| Real space distribution over                          <wbr>        8 groups<br> RS_GRID| Real space distribution along direction                     <wbr>         3<br> RS_GRID| Border size                          <wbr>                               29<br> RS_GRID|   Distribution                  <wbr>       Average         Max         Min<br> RS_GRID|   Planes                        <wbr>         145.5         146         145<br> RS_GRID|   Distribution                  <wbr>       Average         Max         Min<br> RS_GRID|   Planes                        <wbr>         116.3         117         116<br> RS_GRID|   Distribution                  <wbr>       Average         Max         Min<br> RS_GRID|   Planes                        <wbr>         101.8         102         101<br><br>BR, <br><br>Nico<br><br><br>perjantai 1. joulukuuta 2017 11.19.13 UTC+2 <a>jts2...@gmail.com</a> kirjoitti:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi everyone,<div><br></div><div>I have asked several questions about the same subject. Until yesterday, I have figured out the meaning of the *.cube file generated by CP2K. </div><div>What I want to do is to add potential to my simulation system, but the potential is calculated from Poisson equation using finite element method.</div><div>Hence I only have grid point values of the potential. </div><div><br></div><div>So the solution is adding the potential from pot.cube file using READ_FROM_CUBE option. I will try this in a simple simulation to figure the structure. </div><div>The simulation box in my simulation is a cube, whose length of a side is 12.4138 (<span style="background-color:rgb(252,252,254);color:rgb(67,67,67);font-family:Arial,sans-serif;font-size:14px">Angstrom, I think</span>). And from the *.cube file I know the length of the </div><div>vector along x, y and z direction is 0.434420(Angstrom, maybe). What's more, the number of data points along each direction is 54.</div><div><br></div><div>Now I need to supply the pot.cube file to CP2K. I wonder how many grid point value should be given?</div><div>(54*12.4138 / 0.434420) ^3 ?</div><div><br></div><div>I major in mathematic, so the easy thing for you may be a problem for me.</div><div><br></div><div>Thanks for your reply !</div></div></blockquote><br></div></blockquote></div>