howto:lrigpw
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revisionLast revisionBoth sides next revision | ||
howto:lrigpw [2017/05/08 08:33] – [When to use it] dgolze | howto:lrigpw [2023/11/22 10:59] – oschuett | ||
---|---|---|---|
Line 1: | Line 1: | ||
- | ====== How to run DFT with local resolution-of-the-identity in CP2K ====== | + | This page has been moved to: https://manual.cp2k.org/trunk/methods/lrigpw.html |
- | + | ||
- | ===== Introduction ===== | + | |
- | Density functional theory (DFT) calculations in CP2K employ the Gaussian and plane waves (GPW) method. In GPW, the description of the total density on realspace grids is typically the computationally most expensive part. By introducing a local resolution-of-the-identity (LRI) approach, the linear scaling of the GPW approach can be retained, while reducing the prefactor for the grid operations. The combined approach, LRIGPW, is comprehensively described in [[doi> | + | |
- | In LRIGPW, the atomic pair densities $\rho_{\mathrm{AB}}$ are approximated by an expansion in a set of fit functions centered at atom A $\{f_i^{\mathrm{A}}(\mathbf{r})\}$ and atom B $\{f_j^{\mathrm{B}}(\mathbf{r})\}$, | + | |
- | + | ||
- | \begin{equation} | + | |
- | \rho_{\mathrm{AB}}\approx \sum_i{a_i^{\mathrm{A}, | + | |
- | \sum_j{a_j^{\mathrm{B}, | + | |
- | \end{equation} | + | |
- | The fit functions are also Gaussian-type functions and provided as auxiliary basis set. | + | |
- | + | ||
- | ===== How to use it ===== | + | |
- | + | ||
- | LRIGPW is specified in the [[inp> | + | |
- | + | ||
- | < | + | |
- | &QS | + | |
- | METHOD LRIGPW | + | |
- | & | + | |
- | | + | |
- | | + | |
- | &END | + | |
- | &END QS | + | |
- | </ | + | |
- | + | ||
- | Further specifications can be given in the [[inp> | + | |
- | + | ||
- | <code cp2k> | + | |
- | & | + | |
- | BASIS_SET_FILE_NAME BASIS_LRIGPW_AUXMOLOPT | + | |
- | BASIS_SET_FILE_NAME BASIS_MOLOPT | + | |
- | ... | + | |
- | & | + | |
- | & | + | |
- | &KIND O | + | |
- | BASIS_SET DZVP-MOLOPT-GTH | + | |
- | POTENTIAL GTH-PBE-q6 | + | |
- | LRI_BASIS_SET LRI-DZVP-MOLOPT-GTH-MEDIUM | + | |
- | &END KIND | + | |
- | ... | + | |
- | & | + | |
- | </code> | + | |
- | Auxiliary basis sets are available for the MOLOPT basis sets. All auxiliary basis sets have been generated by simple geometric progression without any need for further optimization. These basis sets are available in different sizes: MEDIUM and LARGE. Using the large auxiliary basis sets, the accuracy is improved, but the computational overhead increases.\\ | + | |
- | The LRI auxiliary basis sets are generally quite large leading to a potentially ill-conditioned overlap matrix, Equation (10) in [[doi> | + | |
- | The LRI integrals, Equations (31)-(34) in [[doi> | + | |
- | + | ||
- | ===== When to use it ===== | + | |
- | LRIGPW is only beneficial when the operations on the realspace grids, i.e. collocation of the density and integration of the potential, are dominating the timings. This is typically not the case for metallic systems, where the diagonalization of the Kohn-Sham matrix contributes strongly to the computational cost. LRIGPW is efficient for condensed phase systems such as liquids, molecular crystals etc. Particularly large speed-ups can be obtained for | + | |
- | * non-orthorhombic cells | + | |
- | * large grid cutoffs | + | |
- | * many SCF steps | + | |
- | Using LRI, the SCF step is accelerated and therefore single point calculations profit most. For molecular dynamics, where the wave function can be extrapolated from the previous step, the SCF converges quickly. Also in this case speed-ups can be obtained depending on the grid cutoff and system. | + | |
- | Note that LRIGPW comes with higher memory requirements than the standard GPW scheme. However, this is typically not a problem on HPC platforms, but might limit the usage on smaller clusters. | + | |
- | + | ||
- | ===== Example input files ===== | + | |
- | + | ||
- | * Ice XV: {{: | + |
howto/lrigpw.txt · Last modified: 2024/01/03 13:17 by oschuett