Treatment of delocalized electron transfer in periodic and embedded cluster DFT calculations: The case of Cu on ZnO (10 1¯ 0)

We assess the consequences of the interface model—embedded‐cluster or periodic‐slab model—on the ability of DFT calculations to describe charge transfer (CT) in a particularly challenging case where periodic‐slab calculations indicate a delocalized charge‐transfer state. Our example is Cu atom adsorption on ZnO(10 1¯ 0), and in fact the periodic slab calculations indicate three types of CT depending on the adsorption site: full CT, partial CT, and no CT. Interestingly, when full CT occurs in the periodic calculations, the calculated Cu atom adsorption energy depends on the underlying ZnO substrate supercell size, since when the electron enters the ZnO it delocalizes over as many atoms as possible. In the embedded‐cluster calculations, the electron transferred to the ZnO delocalizes over the entire cluster region, and as a result the calculated Cu atom adsorption energy does not agree with the value obtained using a large periodic supercell, but instead to the adsorption energy obtained for a periodic supercell of roughly the same size as the embedded cluster. Different density functionals (of GGA and hybrid types) and basis sets (local atom‐centered and plane‐waves) were assessed, and we show that embedded clusters can be used to model Cu adsorption on ZnO(10 1¯ 0), as long as care is taken to account for the effects of CT. © 2015 Wiley Periodicals, Inc.

[1]  C. Wöll,et al.  The thermally induced interaction of Cu and Au with ZnO single crystal surfaces , 2013 .

[2]  Karin Fink,et al.  Adsorption of single Cu atoms at differently stabilized polar ZnO surfaces: An ab initio study , 2008, J. Comput. Chem..

[3]  Hafner,et al.  Ab initio molecular-dynamics simulation of the liquid-metal-amorphous-semiconductor transition in germanium. , 1994, Physical review. B, Condensed matter.

[4]  Charles H. Patterson,et al.  Role of defects in ferromagnetism in Zn1- xCoxO : A+ hybrid density-functional study , 2006 .

[5]  Qi Wang,et al.  Effect of the components' interface on the synthesis of methanol over Cu/ZnO from CO2/H2: a microkinetic analysis based on DFT + U calculations. , 2015, Physical chemistry chemical physics : PCCP.

[6]  L. Pettersson,et al.  Electrostatic potential from embedded clusters , 1997 .

[7]  Harrison,et al.  Ab initio study of MnO and NiO. , 1994, Physical review. B, Condensed matter.

[8]  A. Beltrán,et al.  Theoretical study of ZnO (1010) and Cu/ZnO (1010) surfaces , 2001 .

[9]  Saunders,et al.  Magnetic interactions and the cooperative Jahn-Teller effect in KCuF3. , 1995, Physical review. B, Condensed matter.

[10]  S. F. Boys,et al.  The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors , 1970 .

[11]  G. Scuseria,et al.  Hybrid functionals based on a screened Coulomb potential , 2003 .

[12]  Robert Schlögl,et al.  The Mechanism of CO and CO2 Hydrogenation to Methanol over Cu‐Based Catalysts , 2015 .

[13]  Hess,et al.  Hartree-Fock study of phase changes in ZnO at high pressure. , 1993, Physical review. B, Condensed matter.

[14]  K. Fink,et al.  The method of local increments for the calculation of adsorption energies of atoms and small molecules on solid surfaces. Part I. A single Cu atom on the polar surfaces of ZnO. , 2009, Physical chemistry chemical physics : PCCP.

[15]  Hybrid density functional study of band alignment in ZnO-GaN and ZnO-(Ga(1-x)Zn(x))(N(1-x)O(x))-GaN heterostructures. , 2012, Physical chemistry chemical physics : PCCP.

[16]  K. Burke,et al.  Rationale for mixing exact exchange with density functional approximations , 1996 .

[17]  M. Hellström,et al.  Small Cu Clusters Adsorbed on ZnO(10(1)over-bar0) Show Even-Odd Alternations in Stability and Charge Transfer , 2014 .

[18]  W. C. Swope,et al.  A computer simulation method for the calculation of equilibrium constants for the formation of physi , 1981 .

[19]  F. Weigend,et al.  Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. , 2005, Physical chemistry chemical physics : PCCP.

[20]  K. Hermansson,et al.  A combined molecular dynamics+quantum mechanics method for investigation of dynamic effects on local surface structures. , 2004, The Journal of chemical physics.

[21]  G. Kresse,et al.  Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set , 1996 .

[22]  Blöchl,et al.  Projector augmented-wave method. , 1994, Physical review. B, Condensed matter.

[23]  O. Dulub,et al.  Growth of Copper on Single Crystalline ZnO: Surface Study of a Model Catalyst , 2005 .

[24]  E. Longo,et al.  ONIOM study of dissociated hydrogen and water on ZnO surface , 2012 .

[25]  G. Kresse,et al.  Ab initio molecular dynamics for liquid metals. , 1993 .

[26]  I. Ciofini,et al.  Modeling ZnO phases using a periodic approach: from bulk to surface and beyond. , 2009, The Journal of chemical physics.

[27]  Andreas Waag,et al.  Zinc Oxide: From Fundamental Properties Towards Novel Applications , 2010 .

[28]  S. Sinnott,et al.  Growth and Structure of Cu and Au on the Nonpolar ZnO(101̅0) Surface: STM, XPS, and DFT Studies , 2013 .

[29]  Adri C. T. van Duin,et al.  A reactive force field (ReaxFF) for zinc oxide , 2008 .

[30]  Bartolomeo Civalleri,et al.  CRYSTAL: a computational tool for the ab initio study of the electronic properties of crystals , 2005 .

[31]  J. Wagner,et al.  Water Gas Shift Catalysis , 2009 .

[32]  J. Nørskov,et al.  The Active Site of Methanol Synthesis over Cu/ZnO/Al2O3 Industrial Catalysts , 2012, Science.

[33]  A. Zunger,et al.  n -type doping of CuIn Se 2 and CuGa Se 2 , 2005 .

[34]  D. Truhlar,et al.  QM/MM: what have we learned, where are we, and where do we go from here? , 2007 .

[35]  V. Barone,et al.  Toward reliable density functional methods without adjustable parameters: The PBE0 model , 1999 .

[36]  G. Kresse,et al.  From ultrasoft pseudopotentials to the projector augmented-wave method , 1999 .

[37]  P. Sherwood,et al.  Metal Cluster Support Interactions in the Cu/ZnO System: A QM/MM Study , 2003 .

[38]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[39]  Georg Kresse,et al.  Defect energetics in ZnO: A hybrid Hartree-Fock density functional study , 2008 .

[40]  P. Sherwood,et al.  The Growth of Copper Clusters over ZnO: the Competition between Planar and Polyhedral Clusters , 2008 .

[41]  J. V. Lenthe,et al.  State of the Art in Counterpoise Theory , 1994 .

[42]  Malte Behrens,et al.  Heterogeneous catalysis of CO₂ conversion to methanol on copper surfaces. , 2014, Angewandte Chemie.

[43]  T. Liang,et al.  Cu cluster deposition on ZnO101¯0: Morphology and growth mode predicted from molecular dynamics simulations , 2014 .

[44]  Nongnuch Artrith,et al.  Neural network potentials for metals and oxides – First applications to copper clusters at zinc oxide , 2013 .

[45]  O. Dulub,et al.  STM study of Cu growth on the ZnO(1010) surface , 2002 .

[46]  Jianguo Wang,et al.  Copper Deposition and Growth over ZnO Nonpolar (101̅0) and (112̅0) Surfaces: A Density Functional Theory Study , 2009 .

[47]  Artur F Izmaylov,et al.  Influence of the exchange screening parameter on the performance of screened hybrid functionals. , 2006, The Journal of chemical physics.

[48]  M. Hellström,et al.  Band-Filling Correction Method for Accurate Adsorption Energy Calculations: A Cu/ZnO Case Study. , 2013, Journal of chemical theory and computation.

[49]  Kresse,et al.  Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. , 1996, Physical review. B, Condensed matter.

[50]  S. Abrahams,et al.  Atomic displacement, anharmonic thermal vibration, expansivity and pyroelectric coefficient thermal dependences in ZnO , 1989 .