[QE-users] Wrong adsorption energy

Laurens Siemons laurenssiemons at hotmail.be
Thu Mar 8 16:45:54 CET 2018


Thank you all for your quick response.


I'm running a calculation now like Evren Ataman suggested. By doing this the problem of the compensating background should be resolved as well.

The reason I'm not using vdw corrections is based on an article (doi: 10.1021/jz501050f) where they did calculations with and without these corrections. There conclusion is that the corrections are negligible.


Kind regards,

Laurens Siemons

________________________________
Van: users <users-bounces at lists.quantum-espresso.org> namens Stefano de Gironcoli <degironc at sissa.it>
Verzonden: donderdag 8 maart 2018 16:19
Aan: users at lists.quantum-espresso.org
Onderwerp: Re: [QE-users] Wrong adsorption energy


two things I think might be the origin of strange results.



your ethyl-anion is charged. the code introduce a compensating background that should cancel out when you calculate the adsortion energy but you are using different cells in the two charged configurations. so something does not cancel out as it should.



there might be some vdW interaction contributing to the physisorption. I see no inclusion of vdW corrections or use of vdW-aware functionals in your calculations.


best


stefano

On 08/03/2018 15:50, Laurens Siemons wrote:

Dear all,


I'm calculating parameters like adsorption energy, geometric parameters, NMR chemical shifts of adsorbates on a surface. In this specific case I calculated the adsorption energy of an ethyl anion on an anatase 001 surface of titania. I got the following energies:


Esurface + adsorbate = -4301.05337244

Esurface = -4258.9387474

Eadsorbate = -42.37757742


So Eads = Esurface + adsorbate - (Esurface + Eadsorbate) = 0.26


Now this seems wrong. If Eads is positive, so adsorption is not favourable, then the ethyl-anion should migrate away from the surface in the relax calculation of a001cti. correct? I executed the three calculations with the exact same parameters as the other ones so I don't see how this is possible. I attached my input files of the calculations. I hope someone can guide me through this.


Thanks in advance,

Laurens Siemons

Master student chemistry

University of Antwerp




_______________________________________________
users mailing list
users at lists.quantum-espresso.org<mailto:users at lists.quantum-espresso.org>
https://lists.quantum-espresso.org/mailman/listinfo/users

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20180308/e75015fe/attachment.html>


More information about the users mailing list