File: gpw.md

package info (click to toggle)
cp2k 2025.1-1.1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 366,832 kB
  • sloc: fortran: 955,049; f90: 21,676; ansic: 18,058; python: 13,378; sh: 12,179; xml: 2,173; makefile: 964; pascal: 845; perl: 492; lisp: 272; cpp: 137; csh: 16
file content (40 lines) | stat: -rw-r--r-- 1,244 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
# Gaussian Plane Wave

The primary basis in CP2K consists of Gaussian Type Orbital (GTO) functions. To take advantage of
the efficient FFT algorithm for solving the Poisson equation, the electronic density must first be
transferred from the GTO representation to a regular grid. This is referred to as _collocation_.
After solving the Poisson equation the obtained electrostatic potential has to be transferred back
into the GTO basis, which is referred to as _integration_. This ”opportunistic” switching between a
Gaussian and a plane wave representation is the core idea of the Gaussian and Plane Waves (GPW)
method. Its high-level operations are illustrated in the following flowchart:

```mermaid
block-beta
columns 5
a["el. Density\n(Gaussian Basis)"]
space
b["el. Density\n(Regular Grid)"]
space
c["el. Density\n(Plane Waves)"]

space space space space space

f["el. Potenial\n(Gaussian Basis)"]
space
e["el. Potenial\n(Regular Grid)"]
space
d["el. Potenial\n(Plane Waves)"]

a -- "Collocate"        --> b
b -- "FFT"              --> c
c -- "Poisson Solver"   --> d
d -- "FFT<sup>-1</sup>" --> e
e -- "Integrate"        --> f
```

## See also

- <https://www.cp2k.org/gpw>
- <https://www.cp2k.org/quickstep>
- [](#VandeVondele2005)