ENCUT and k-point convergence for slab+molecule adsorption

[emilyherrmann]

Hello!

I am studying the adsorption energy of an organic semiconductor molecule (pentacene) on a metal (copper) surface and I have a few questions about ENCUT and k-point convergence.

I had previously done these convergences just for the unit Cu cell and was converging based on the E0 value. After reading further, I have learned that it is beneficial to converge based on the adsorption energy, because that is my property of interest. This means that I will need to collect the value of (E0 - KE correction term) for the surface on its own, then for the molecule on its own, then for the combined system.

1) I have read that when modelling a single molecule, you should always use just the Gamma point for your k-point mesh. I would assume, however, that you should use the same k-point mesh for all three sets of runs (surface, molecule, combined) of the process. Is this true? If that's the case, will you end up using a non-gamma point mesh for the molecule?

2) When I previously did k-point convergence for just the Cu unit cell, the intended process was: 1) Converge k-points for single unit cell 2) Do Cu layer convergence. Every time the size of the Cu cell doubles in the z direction, halve the number of k-points in the z direction. 4) Do vacuum convergence on the pentacene molecule to determine how far apart the molecules would need to be to avoid interaction. Make the Cu cell this size in the x and y directions, similarly decreasing the number of k-points as the cell increases in size. Should I instead not do any of this halving and just converge the k-points for my adsorbed system and use that value?

If it is easier for you to just explain the process of ENCUT and k-point convergence for an adsorbed system step by step instead of directly answering my questions, that is fine!

Thanks!
Emily

[christopher_sheldon1]

\(\)Dear Emily,

Thank you for your question. You are correct that, for modelling adsorption on a metal surface, it is important to converge the adsorption energy. More generally, we recommend converging with respect to your property of interest. Converge using the

Code: Select all

free  energy   TOTEN  =

found in the OUTCAR file, except for a few special cases (this is not one of them).

The adsorption energy \(\Delta E\_{ads}\) is the energy difference between the bare surface \( E\_{S}\), the molecule \( E\_{M}\), and the molecule on the surface \( E\_{M \cdot S}\):

\(\Delta E\_{ads} = E\_{M \cdot S} - E\_{S} - E\_{M} \)

1) You should use the same k-mesh for S and M · S. It is important to use only the Gamma-point on the molecule. This is because you do not want any periodic interactions. Instead, you want to eliminate these to obtain energies for a gas-phase molecule. If you use a k-mesh, then you include the periodic images of the molecule and so calculate the energy of a diffuse molecule crystal, not the gas phase.

TLDR: use only the Gamma point k-mesh for a gas-phase molecule

2) To converge with respect to k-points, use an NxNx1 mesh, where 1 is the direction with the vacuum. You need to converge the k-mesh separately for each cell that you are using. For example, you could start with a 2x2x1 k-mesh for a (2x2) cell, then increase to 4x4x1, 6x6x1, etc until you reach convergence. This needs to be done for both S and M · S. The energy is not variational with respect to k-points, so it is important to use an even set of meshes (i.e., 2x2x1, 4x4x1, 6x6x1) or and odd set (e.g., 3x3x1, 5x5x1, 7x7x1) but do not mix the two, as some of the special points in the Brillouin zone will not be in both, which introduces additional error.

TLDR:
Should I instead not do any of this halving and just converge the k-points for my adsorbed system and use that value?
Yes, converge your adsorption energy with respect to k-point mesh.

You should also perform similar convergence for the energy cutoff (ENCUT) and the vacuum. Converge each variable in turn, rather than doing all at once. I tend to do:

• Number of layers of atoms - I usually stick with 4 but this should be tested in principle

• Vacuum height - (surface - 14 A is usually sufficient but should be checked) + (molecule - i tend to use 20x20x20 A, for a molecule as large as pentacene, you may need something larger)

• ENCUT

• k-point mesh

I notice that you're also using pentacene on the metal surface. It's worth checking to see if van der Waals is important (IVDW), which can be tested by doing all the calculations with and without it. D3 is a reasonable first guess, dDsC or MBD perform better but can also come with some computational difficulties. I covered the convergence test for the adsorption energy of RPA in this paper (https://doi.org/10.1063/5.0071995) but the procedure stands true for any density functional.

I hope that this helps to answer your question.

Best wishes,

Chris

EDIT:
It's also important to check whether or not the dipole correction is important. This requires the DIPOL, IDIPOL, and LDIPOL tags; Exercise 5 of the surface tutorials goes through this (https://vasp.at/tutorials/latest/surfac ... urface-e05). Apply it to the surface and surface with molecule adsorbed. If you place the molecule at the centre of the cell