Differences

This shows you the differences between two versions of the page.

--- periodicity [2020/08/21 10:15] – external edit 127.0.0.1
+++ periodicity [2021/11/05 16:14] (current) – Big Update, replace all manual entries with links, remove old info, mention more things. pmelix
@@ Line 1: / Line 1: @@
-<note warning>This text is probably out of date and needs to be revised.</note>
 ====== Periodicity ======
-The Cp2k program is optimized for periodic calculations, and those are the default, still it is possible to perform also other kinds of calculations. First an important thing to note is that by convention the cell goes from 0 to h, not from -h/2 to h/2, thus if you want to put a cluster at the center of the cell you need to put it in h/2 not in 0.
+CP2K is optimized for periodic calculations, and those are the default. Still it is possible to perform other kinds of calculations. First an important thing to note is, that by convention the cell goes from 0 to h, not from -h/2 to h/2. Thus if you want to put a cluster at the center of the cell you need to put it in h/2 not in 0. This can be done automatically using the [[inp>FORCE_EVAL/SUBSYS/TOPOLOGY/CENTER_COORDINATES]] option. For certain Poisson solvers this is recommended to ensure an electron density of zero on the non-periodic faces of the cell.
-The periodicity of the neighbor list, and thus of the gaussian collocation/integration is controlled by ''FORCE_EVAL%SUBSYS%CELL%PERIODIC''.
-The poisson solver (electrostatics) is controlled in the section
+The periodicity of the neighbor list, and thus of the gaussian collocation/integration is controlled by the keyword [[inp>FORCE_EVAL/SUBSYS/CELL#PERIODIC]].
-''FORCE_EVAL%DFT%POISSON''. The possible poisson solver are chosen with the ''POISSON_SOLVER'' keyword, that can have the values
-  * ''PERIODIC'', which can do only fully periodic solvers,
+The Poisson solver (electrostatics) is controlled in the section
-  * ''ANALYTIC'' which can do 0, 1d and 2d periodic solutions using analytical green functions in the g space (slow convergence),
+[[inp>FORCE_EVAL/DFT/POISSON]]. The possible Poisson solvers are chosen with the ''POISSON_SOLVER'' keyword. Note, that not all solvers support all periodicities.
-  * ''MT'' Martyna Tuckermann decoupling, that interacts only with the nearest neighbor, beware these are completely wrong in the cell is smaller than twice the cluster size (with electronic density),
-  * ''MULTIPOLE'', uses a schema that fits the total charge with one gaussian per atom
-The periodicity of the poisson solver (electrostatics) is set in ''FORCE_EVAL%DFT%POISSON%PERIODIC'', and the poisson solver chosen should support it.
+The periodicity of the Poisson solver (electrostatics) is set in [[inp>FORCE_EVAL/DFT/POISSON#PERIODIC]], and the Poisson solver chosen should support it. In most cases the periodicity settings of the Poisson solver and the subsystem should be identical.