Poster

Is the van der Waals interaction relevant to the dissociative dynamics of N₂ on W(110)?

L. Martin-Gondre^1,2, J. I. Juaristi^2,3,4, M. Blanco-Rey^2,3, R. Díez Muiño^2,4, and M. Alducin^2,4

¹Institut UTINAM - CNRS UMR 6213 Université de Franche-Comté, Besançon, France

²Donostia International Physics Center (DIPC), Donostia-San Sebastián, Spain

³Dep. de Física de Materiales, UPV/EHU, Donostia-San Sebastián, Spain

⁴Centro de Física de Materiales, CFM/MPC (CSIC-UPV/EHU), Donostia-San Sebastián, Spain

Using ab initio molecular dynamics (AIMD) calculations we investigate the role of the van der Waals interaction in the dissociative adsorption of N₂ on W(110) [1]. Hitherto, existing classical dynamics calculations performed on six-dimensional potential energy surfaces based on density functional theory (DFT) and semi-local exchange-correlation functionals were unable to fully describe the dependence of the initial sticking coefficient on the molecular beam incidence conditions as found in experiments. N₂ dissociation on W(110) was shown to be very sensitive not only to close molecule-surface distances, but also to the large distances region where the van der Waals interaction, not included in semilocal-DFT, should dominate. In this work, we perform a systematic study on the dissociative adsorption using a selected representation of existing non-local functionals that include the van der Waals interaction. Starting with a pure static analysis of the potential energy surfaces, we find that the original vdW-DF [2] and vdW-DF2 [3] functionals provide a rather good description of the experimental observations regarding the adsorption well and the balance between desorption and dissociation of the chemisorbed N₂. In contrast, the optimized vdW-functionals proposed in Refs. [4,5] share the limitations of the (semi-local) PW91 potential energy surface for this system, i.e., an excessively deep adsorption well from which the dissociation path compared with the desorption one is energetically favored. Unexpectedly, the results of the AIMD performed at different incidence conditions are not as satisfactory as the preliminary static analysis suggested. We have to conclude that none of the used vdW-functionals seem to provide altogether an adequate description of the N₂/W(110) interaction at short and large distances.

[1] L. Martin-Gondre et al., J. Chem. Phys. 142, 074704 (2015)

[2] M. Dion, H. Rydberg, E. Schröder, D. C. Langreth, and B. I. Lundqvist, Phys. Rev. Lett. 92, 246401 (2004)

[3] K. Lee, E. D. Murray, L. Kong, B. I. Lundqvist, and D. C. Langreth, Phys. Rev. B 82, 081101 (2010)

[4] J. Klimeš, D. R. Bowler, and A. Michaelides, J. Phys.: Condens. Matter 22, 022201 (2010)

[5] J. Klimeš, D. R. Bowler, and A. Michaelides, Phys. Rev. B 83, 195131 (2011)