In classical MD simulations of lattice structures, whenever we are trying to incorporate vdW corrections to the atomic force calculation, we need to set up cutoff length beyond which the vdW correction is not considered to keep the computational cost manageable. As a result, the classical simulation requires reasonably large supercell(usually with its dimensions twice as long as the cutoff length) to capture vdW interactions.

When it comes to DFT calculations, however, the cutoff lengths for nonlocal vdW corrections, such as the rVV10 functional, are not explicitly set in DFT software like Quantum ESPRESSO. It is also common practice to use primitive cells in DFT calculations along with nonlocal vdW functionals. By looking at the formalism of the rVV10 vdW correction, the integration of the functional is clearly performed within a volume, $$ E_{c}^{n l}=\frac{\hbar}{2} \iint n(\mathbf{r}) \Phi\left(n, n^{\prime},|\nabla n|,\left|\nabla n^{\prime}\right|,\left|\mathbf{r}-\mathbf{r}^{\prime}\right|\right) n\left(\mathbf{r}^{\prime}\right) d \mathbf{r} d \mathbf{r}^{\prime},\tag{1} $$ where $n(\mathbf{r})$ is the charge density at location of $\mathbf{r}$. So my question is what are the lower and upper limits for Eq. (1), and how do we derive them?

To attack this question, I tried to understand Eq.(1) in terms of its Fourier transform, $$ E_{c}^{n l}=\frac{\hbar}{2} \Omega \sum_{i j} \sum_{\mathbf{G}} \tilde{\theta}_{i}^{*}(\mathbf{G}) \tilde{\Phi}^{\mathrm{rVV} 10}\left(q_{i}, q_{j},|\mathbf{G}|\right) \tilde{\theta}_{j}(\mathbf{G}),\tag{2} $$ where $\tilde{\Phi}^{\mathrm{rVV} 10}\left(q_{i}, q_{j},|\mathbf{G}|\right)$ are the Fourier transforms of the rVV10 kernel evaluated on a bidimensional grid of $q$ values, and $\tilde{\theta}_{i}$ are the Fourier transforms of $\theta_{i}$ introduced in the rVV10 paper here. But this left me with more questions than answers. I have a feeling that I might be pursuing the answer in a wrong line of thinking, so any hint will be appreciated!

  • $\begingroup$ +10. Nice first question! I see you've been here for 8 months but I'm meeting you for the first time now, so welcome to our new community and thank you very much for contributing your question here! We hope to see much more of you in the future, and thanks also for formatting yoru question so nicely, including equation number labels!!! Here's a related question which might interest you: mattermodeling.stackexchange.com/q/63/5 $\endgroup$ Nov 24, 2021 at 21:50
  • $\begingroup$ I have commented out your second question: "how is the wavevector in the reciprocal space $\mathbf{G}$ is related to the cutoff length of vdW correction?", which I encourage you to ask as a separate question if you need that! $\endgroup$ Jan 20, 2022 at 22:27
  • $\begingroup$ Have you had any luck with this since November? $\endgroup$ Feb 24, 2022 at 3:56
  • $\begingroup$ unfortunately no... this whole vdW correction, to a large extent, is still mysterious to me... Thanks for asking tho@NikeDattani $\endgroup$
    – Lonitch
    Jul 12, 2022 at 13:45
  • 1
    $\begingroup$ Good discussion here. Just a reminder that we prefer not to create new chat, so don't click the prompt to move to chat if the discussion goes long. Just ping me and I can move the comments to an existing room if you want to continue in chat. $\endgroup$
    – Tyberius
    Feb 6, 2023 at 2:33