KERNEL TRUNCATION/LTRUNCATE: Difference between revisions
| Line 36: | Line 36: | ||
This setup corresponds to truncating the Coulomb interaction along the surface normal direction (say, along z) for a [[2D material]], using no vacuum padding and a truncation length of z/2. In this case, half of the simulation box is effectively unused and will produce a potential that is not desired. | This setup corresponds to truncating the Coulomb interaction along the surface normal direction (say, along z) for a [[2D material]], using no vacuum padding and a truncation length of z/2. In this case, half of the simulation box is effectively unused and will produce a potential that is not desired. | ||
However, the algorithm is much simpler. | However, the algorithm is much simpler. | ||
We recommend this configuration | We recommend this configuration for debugging purposes. | ||
== Related tags and articles == | == Related tags and articles == | ||
Latest revision as of 15:24, 18 March 2026
KERNEL_TRUNCATION/LTRUNCATE = .True. | .False.
Default: KERNEL_TRUNCATION/LTRUNCATE = .False.
Description: Truncates the Coulomb kernel to remove electrostatic interactions along non-periodic dimensions.
Setting KERNEL_TRUNCATION/LTRUNCATE = T switches on the Coulomb-kernel-truncation method[1][2][3]. It effectively removes interactions with periodic replicas in non-periodic directions. In other words, the interactions are removed along the surface normal for 2D materials, and along all directions for 0D systems, i.e. for isolated atoms and molecules.
In the simplest implementation of the Coulomb-kernel-truncation method (KERNEL_TRUNCATION/LCOARSEN = F), the computational cell provided in the POSCAR file is internally padded by an additional vacuum (see KERNEL_TRUNCATION/IPAD). This implies increasing the FFT-grid sizes by a certain factor and thus leads to a significant increase in computational cost.
Tip: Use the KERNEL_TRUNCATION/LCOARSEN = T to avoid the increased FFT-grid sizes.
|
Mind:
|
Detailed information about the setting are documented on respective related tags.
| Warning: When padding is used, the vaccum is added on the edges of the cell, therefore it is very important there are no atoms on the cell boundary in the non-periodic direction. We recommend centering the motif in the simulation box. If you encounter problems using Coulomb truncation with padding, try the same calculations without padding (see examples bellow). |
Example
KERNEL_TRUNCATION {
LTRUNCATE = T
IDIMENSIONALITY = 2
ISURFACE = 3
LCOARSEN = F
}
In this case, we pad the cell along the surface normal direction. The Coulomb interaction is truncated beyond the boundaries of the cell along this direction.
KERNEL_TRUNCATION {
LTRUNCATE = T
IDIMENSIONALITY = 2
ISURFACE = 3
IPAD = 1
FACTOR = 0.5
}
This setup corresponds to truncating the Coulomb interaction along the surface normal direction (say, along z) for a 2D material, using no vacuum padding and a truncation length of z/2. In this case, half of the simulation box is effectively unused and will produce a potential that is not desired. However, the algorithm is much simpler. We recommend this configuration for debugging purposes.
Related tags and articles
KERNEL_TRUNCATION/LCOARSEN, KERNEL_TRUNCATION/IDIMENSIONALITY, KERNEL_TRUNCATION/ISURFACE, KERNEL_TRUNCATION/FACTOR, KERNEL_TRUNCATION/IPAD
References
- ↑ S. Vijay, M. Schlipf, H. Miranda, F. Karsai, M. Kaltak, M. Marsman, and G. Kresse, Efficient periodic density functional theory calculations of charged molecules and surfaces using Coulomb kernel truncation, Phys. Rev. B 112, 045409 (2025).
- ↑ C. A. Rozzi, D. Varsano, A. Marini, E. K. Gross, A. J. Rubio, Phys. Rev. B 73, 20511 (2006).
- ↑ T. Sohier, M. Calandra, and F. Mauri, Phys. Rev. B 96, 75448 (2017).