Ab initio statistical mechanics of surface adsorption and desorption.: I.: H2O on MgO (001) at low coverage

被引:18
作者
Alfe, D.
Gillan, M. J.
机构
[1] UCL, London Ctr Nanotechnol, London WC1H 0AH, England
[2] UCL, Dept Phys & Astron, London WC1E 6BT, England
[3] UCL, Dept Earth Sci, London WC1E 6BT, England
基金
英国工程与自然科学研究理事会;
关键词
1ST-ORDER PHASE-TRANSITIONS; FREE-ENERGY CALCULATIONS; AUGMENTED-WAVE METHOD; MOLECULAR-DYNAMICS; MGO(100) SURFACE; THERMAL-DESORPTION; WATER; MGO; SIMULATIONS; DISSOCIATION;
D O I
10.1063/1.2772258
中图分类号
O64 [物理化学(理论化学)、化学物理学];
学科分类号
070304 ; 081704 ;
摘要
We present a general computational scheme based on molecular dynamics (MD) simulation for calculating the chemical potential of adsorbed molecules in thermal equilibrium on the surface of a material. The scheme is based on the calculation of the mean force in MD simulations in which the height of a chosen molecule above the surface is constrained and subsequent integration of the mean force to obtain the potential of mean force and hence the chemical potential. The scheme is valid at any coverage and temperature, so that in principle it allows the calculation of the chemical potential as a function of coverage and temperature. It avoids all statistical mechanical approximations, except for the use of classical statistical mechanics for the nuclei, and assumes nothing in advance about the adsorption sites. From the chemical potential, the absolute desorption rate of the molecules can be computed, provided that the equilibration rate on the surface is faster than the desorption rate. We apply the theory by ab initio MD simulation to the case of H(2)O on MgO (001) in the low-coverage limit, using the Perdew-Burke-Ernzerhof (PBE) form of exchange correlation. The calculations yield an ab initio value of the Polanyi-Wigner frequency prefactor, which is more than two orders of magnitude greater than the value of 10(13) s(-1) often assumed in the past. Provisional comparison with experiment suggests that the PBE adsorption energy may be too low, but the extension of the calculations to higher coverages is needed before firm conclusions can be drawn. The possibility of including quantum nuclear effects by using path-integral simulations is noted. (c) 2007 American Institute of Physics.
引用
收藏
页数:12
相关论文
共 43 条
[1]   Absolute rate of thermal desorption from first-principles simulation [J].
Alfe, D. ;
Gillan, M. J. .
JOURNAL OF PHYSICS-CONDENSED MATTER, 2006, 18 (37) :L451-L457
[2]   Iron under Earth's core conditions:: Liquid-state thermodynamics and high-pressure melting curve from ab initio calculations -: art. no. 165118 [J].
Alfè, D ;
Price, GD ;
Gillan, MJ .
PHYSICAL REVIEW B, 2002, 65 (16) :1-11
[3]   The melting curve of iron at the pressures of the Earth's core from ab initio calculations [J].
Alfè, D ;
Gillan, MJ ;
Price, GD .
NATURE, 1999, 401 (6752) :462-464
[4]   Thermodynamics of hexagonal-close-packed iron under Earth's core conditions -: art. no. 045123 [J].
Alfè, D ;
Price, GD ;
Gillan, MJ .
PHYSICAL REVIEW B, 2001, 64 (04)
[5]   MOLECULAR-DYNAMICS SIMULATIONS AT CONSTANT PRESSURE AND-OR TEMPERATURE [J].
ANDERSEN, HC .
JOURNAL OF CHEMICAL PHYSICS, 1980, 72 (04) :2384-2393
[6]  
BECKER E, 2006, J CHEM PHYS, V125
[7]   PROJECTOR AUGMENTED-WAVE METHOD [J].
BLOCHL, PE .
PHYSICAL REVIEW B, 1994, 50 (24) :17953-17979
[8]   First-order phase transitions by first-principles free-energy calculations: The melting of Al [J].
de Wijs, GA ;
Kresse, G ;
Gillan, MJ .
PHYSICAL REVIEW B, 1998, 57 (14) :8223-8234
[9]   The structure and spectroscopy of monolayers of water on MgO:: An ab initio study [J].
Delle Site, L ;
Alavi, A ;
Lynden-Bell, RM .
JOURNAL OF CHEMICAL PHYSICS, 2000, 113 (08) :3344-3350
[10]   Observation of the second ordered phase of water on the MgO(100) surface: Low energy electron diffraction and helium atom scattering studies [J].
Ferry, D ;
Glebov, A ;
Senz, V ;
Suzanne, J ;
Toennies, JP ;
Weiss, H .
JOURNAL OF CHEMICAL PHYSICS, 1996, 105 (04) :1697-1701