howto:kg
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
howto:kg [2018/12/05 14:07] – [CP2K tutorial] 130.60.136.203 | howto:kg [2024/01/03 13:20] (current) – oschuett | ||
---|---|---|---|
Line 1: | Line 1: | ||
- | ====== How to run simulations with KG method ====== | + | This page has been moved to: https://manual.cp2k.org/trunk/methods/embedding/kim-gordon.html |
- | + | ||
- | ===== Introduction ===== | + | |
- | + | ||
- | This method is based on density embedding. Let's introduce first the subtraction scheme definition of the density embedding method: | + | |
- | + | ||
- | $\displaystyle E_{tot} = E_{HK}[\rho_{tot}] - \sum_{A}E_{HK}[\rho_{A}] + \sum_{A}E_{KS}[\rho_{A}]$. | + | |
- | + | ||
- | The total electronic density of the $\rho_{tot} = \sum_{A}\rho_{A}$ as the sum over all the subsystems $A$ of the subsystem densities $\rho_{A}$. The energy functionals $E_{HK}$ and $E_{KS}$ are the Hohenberg–Kohn and the Kohn–Sham functionals, | + | |
- | + | ||
- | $\displaystyle E_{HK}[\rho] = T_{HK}[\rho] + E_{ext}^{HK}[\rho] + \frac{1}{2} \int\int \frac{\rho(r)\rho(r' | + | |
- | $\displaystyle E_{KS}[P] = T_{S}[P] + E_{ext}[P] + \frac{1}{2} \int\int \frac{\rho(r)\rho(r' | + | |
- | + | ||
- | where $P$ is the reduced one-particle density matrix of the system. In order to arrive at the working equations, one has to introduce the restriction that the external energy functional in the Hohenberg–Kohn energy is linear in the density. | + | |
- | + | ||
- | $\displaystyle E_{ext}^{HK}[\rho_{tot}] = \sum_{A}E_{ext}^{HK}[\rho_{A}]$. | + | |
- | + | ||
- | Now, if one calls the classical Coulomb term $E_{hxc}[\rho]$ and defines the non-additive kinetic energy as $T_{nadd}[\rho, | + | |
- | + | ||
- | $\displaystyle E_{tot}[{P_{A}}] =\sum_{A}(T_{S}[P_{A}] + E_{ext}[P_{A}]) + E_{hxc}[\rho] + T_{nadd}[{P_{A}}]$. | + | |
- | + | ||
- | To avoid the integration of the kinetic energy functional for each subsystem, an atomic potential approximation can be applied. For a local potential one can write: | + | |
- | + | ||
- | $\displaystyle T_{nadd} = T_{S}[\rho]−\sum_{A}T_{S}[\rho_{A}] =$\\ | + | |
- | $\displaystyle \int\rho\mu[\rho]dr - \sum_{a}\int\rho_{A}\mu[\rho_{A}]dr =$\\ | + | |
- | $\displaystyle \sum_{a}\int\rho_{A}(\mu[\rho]-\mu[\rho_{A}])dr$ | + | |
- | + | ||
- | Making a linearization approximation for the functional $\mu[\rho]$ | + | |
- | + | ||
- | $\displaystyle \mu[\rho]-\mu[\rho_{A}] \sim \sum_{B/not = A} \frac{\partial \mu[\rho_{A}]}{\partial \rho} \rho_{B} = \mu' | + | |
- | $\displaystyle T_{nadd} = \sum_{A}T_{S}\sum_{B/not = A}\int\mu' | + | |
- | + | ||
- | A further approximation of the derivative functional in atomic contributions is: | + | |
- | + | ||
- | $\displaystyle \mu' | + | |
- | + | ||
- | and the realization that a typical kinetic energy functional is proportional to $\rho^{5/ | + | |
- | + | ||
- | $\displaystyle V_{a}^{K}(R_{a}) = N_{a}\rho_{a}^{2/ | + | |
- | + | ||
- | where $\rho_{a}$ is a model atomic density. Such local potential can help to speed up the underlying embedding calculation. | + | |
- | + | ||
- | ===== CP2K tutorial ===== | + | |
- | + | ||
- | The division of the total system into subsystems is a critical point, in order to do that properly it is important to tell the CP2K which is the ' | + | |
- | + | ||
- | < | + | |
- | & | + | |
- | &CELL | + | |
- | ABC 9.8528 9.8528 9.8528 | + | |
- | &END CELL | + | |
- | & | + | |
- | | + | |
- | | + | |
- | | + | |
- | | + | |
- | | + | |
- | | + | |
- | | + | |
- | | + | |
- | | + | |
- | ... | + | |
- | &END COORD | + | |
- | & | + | |
- | CONN_FILE_FORMAT USER | + | |
- | &END | + | |
- | </code> | + | |
- | + | ||
- | This strategy is based on the fourth column in the COORD section. At this point the code is able to find the best combination of ' | + | |
- | + | ||
- | < | + | |
- | & | + | |
- | MAX_SCF | + | |
- | EPS_FILTER | + | |
- | EPS_SCF | + | |
- | MU | + | |
- | PURIFICATION_METHOD TRS4 | + | |
- | &END | + | |
- | </code> | + | |
- | + | ||
- | This speeds up the calculation, | + | |
- | + | ||
- | < | + | |
- | &QS | + | |
- | LS_SCF | + | |
- | KG_METHOD | + | |
- | ... | + | |
- | &END QS | + | |
- | </note> | + | |
- | + | ||
- | Once all these passages are done, one has to choose the TNADD_METHOD. For the first type of calculation, | + | |
- | + | ||
- | < | + | |
- | &XC | + | |
- | & | + | |
- | & | + | |
- | FUNCTIONAL T92 # | + | |
- | &END | + | |
- | &END | + | |
- | &END | + | |
- | </code> | + | |
- | + | ||
- | And in the same section others corrections can be added (example: VDW_POTENTIAL).\\ | + | |
- | For the second type of calculation the keyword to select is ATOMIC. This method implies a supplemental atomic potential (create a file which contains all the required potentials). These templates of potentials can be found inside the tests > QS > regtest-kg folder of the CP2K and they can be generated directly from the code (look at tests > ATOM > regtest-pseudo > O_KG.inp). It's important to point out that this method is still in the experimental stage and further investigations are needed. | + | |
- | + | ||
- | < | + | |
- | </ | + |
howto/kg.1544018827.txt.gz · Last modified: 2020/08/21 10:15 (external edit)