Overview. kmcand DFT. An example to start with. An example to start with 03/09/2014
|
|
- Curtis Fields
- 5 years ago
- Views:
Transcription
1 An Introduction of the Use of Monte Carlo Simulations in Catalysis. Tonek Jansen verview Kinetics and kinetic Monte Carlo The Master Equation Getting the Rate Constants Algorithms Examples Bangalore, 5 September 014 kmcand DFT densityfunctional theory adsorption sites other quantum chemical methods molecular mechanics experiments V ( R) reaction mechanisms rate constants molecular dynamics thermodynamics Kinetics and kinetic Monte Carlo fantasy/imagination/insight kinetics An example to start with C+ adlayer Pt(110) substrate An example to start with adlayer substrate 0.4 µm 0.4 µm time: length: seconds meter C 1x 1x1 Salazar et al., PRE, 69 (004) C 1x 1x1 Salazar et al., PRE, 69 (004)
2 Typical kmc simulations Length and time scales grid size: 18x18 104x104 system size: 5 nm 0. µm second coarsegrained kmc rate equations real time simulated: sec processor time used: 1 sec 1 hour number of reactions: number of reaction types: 5 memory used: 8 18 Mb quantum chemistry kmc MD meter Chemical kinetics an example: [ C] d C(gas) C(ads) (gas) (ads) C(ads) + (ads) C(gas) [ ] [ C][ ] = k ads * k form hat s wrong with... macroscopic rate equations, mass balance equations, reactiondiffusion equations, These equations implicitly assume that the reactants are homogeneously/randomly mixed; i.e., there is no correlation between the position of the reactants. surface science: d C = kads* k form C coverages Shortrange order The rate of bimolecular reactions depends on the number of pairs of adsorbates that react. The number of pairs are assumed to be equal to a product like C. hen there is shortrange order this is not correct. Shortrange order arises from diffusion limitations, lateral interactions, and defects. ith lateral interactions it is not even possible to talk about one rate constant for a reaction; i.e., k form C is not even correct if there is no shortrange order. Phenomenology an example: d C(gas) C(ads) (gas) (ads) C(ads) + (ads) C(gas) k C = ads* k form C Rate equations can be used to fit results of kinetic experiments. This is often successful. The drawbacks are that this only works for a limited range of reaction conditions, and the rate constants have no physical meaning. ne may also wonder why these particular equations are used.
3 The (111) surface of an fccmetal bridge/fold sites fcc/fold top/1fold sites sites hcp/fold sites The Master Equation Key concept Atoms and molecules adsorb at welldefined sites. These sites form a regular grid or lattice. This grid of sites with a description of what is adsorbed at each site is called a latticegas model. Each grid or lattice point has a label that specifies the occupation of the corresponding site and possibly other properties of the site. Phase space N particles: phase space: q, q,, q 1 K N coordinates p,, p, 1 K p N conjugate momenta ( q q, K, q, p, p, K, p ) ( q, p) 1, N 1 N = p q Evolution of a system in phase space p Time scales vibrations: subpicoseconds reactions: seconds q
4 Evolution of a system in phase space Key concept p q If you only specify the sites and their occupation, then the times when reactions take place can not be determined exactly. Instead we get probability distributions for the reaction times. If more reactions are possible, then we also can not say which reaction will take place, but can only give probabilities. Reactions in phase space dp dqdp ρ ( t) = p h t = R dqdp P ( t) = ρ p h ( q, t) D, D R dqdp h D i= 1 R ( q, t) D, ρ H ρ H pi qi qi pi The master equation dp ( TS ) kbt Q = ( IS ) h Q = [ P P ] ( IS ) dp Q = dq exp D h R B ( ) dp dp dp + dp TS 1K j 1 j 1K D Q = ds exp D 1 h S B T H k T H k The master equation dp = [ P P ] configuration probability time rate constant Key concept The evolution of a latticegas model of a surface and its adsorbates can be described by the Master Equation dp = [ P P ] 4
5 The role of the master equation ab initio kinetics (quantum chemistry) dp = [ P P ] Getting the Rate Constants numerical analytical analytical kinetic Monte Carlo rate equations Calculating rate constants dp = [ P P ] ( TS ) k = BT Q E ( ) bar exp IS h Q kbt Differences between transition probabilities and rate constants are a matter of definition. e do not distinguish between them. Rate constants from experiments Property X has value X in configuration. Statistical average: Rate equation: X d X = P X = dp X = [ P P ] X = P [ X X ] = P [ X X ] change in X average over initial configurations rate constant of change Key concept The evolution of a property can be described by the macroscopic or phenomenological equation d X = P [ X X ] Simple desorption e assume that we have just one type of adsorbate that desorbs one atom/molecule at a time, and suppose that a surface area with S sites has N adsorbates. Summation over : d N = P [ N N ] = 0, des, if is a desorption otherwise If is a desorption then N N = 1 There are N desorptions in = #nonvanishing terms. d N = des P N = des N 5
6 Simple desorption d N d = des N divide by S: = des macroscopic rate equation: d = des = k des k des The macroscopic rate equation is exact. determined from experiments coverage Bimolecular reaction e assume that we have two types of adsorbate that react with each other; A + B K ( A) d N ( A) ( A) = P [ N N ] Summation over : = 0, rx, if is a reaction otherwise ( A ) ( A ) If is a reaction then N N = 1 There are ( AB) N #AB pairs reactions in = #nonvanishing terms. ( A) d N ( AB) ( AB) = rx P N = rx N Bimolecular reaction ( A) d N ( AB) ( AB) = rx P N = rx N Two approximations need to be introduced. The adsorbatesare randomly distributed. Fluctuations can be neglected. macroscopic rate equation: rx = k rx Z d A d B = = krx A B determined from experiments da d B = = Zrx A B The macroscopic rate equations are not exact. Algorithms Kinetic Monte Carlo KMC generates an instantiation of the probabilities of the master equation. This means that it generates a configuration at time t with probability P (t), which is a solution of the master equation. t = 0 : 1 0 t 1 t L t n n Variable Step Size Method (VSSM) Initialize 1 Reaction time t = ln r, t t + t Reaction p( ) = Continuation stop, 6
7 First Reaction Method (FRM) Initialize Generate configuration, set time t, make a list of all reactions, set stop conditions, determine: 1 t = t ln r Reaction Take the first reaction to occur. Performance FRM VSSM ZGB TPD w. repulsion C+H/PtRu 60 RSM 800 FRM +RSM VSSM +RSM #sites sulphate/cu(111) C+/Pt(100) N/Rh(111) 8 7 Continuation stop NH /Pt(111) 1 ) modified model ) approximate algorithm in 10 processes per CPU second A short history of Carlos Pizzazz: a predecessor of Carlos (1994). generalpurpose code using userspecified labels firstreaction method to solve the Master Equation including linear timedependence of the temperature Carlos ( ) first version (1996): improved version of Pizzazz version.0 (1998): three algorithms, improved userfriendliness, graphical output version 4.0 (004): lateral interactions, more general timedependent rate constants version 5.0 (009): more flexible user interface Kinetix (010 ) commercial version; part of Accelrys s Materials Studio Examples Reduction of NH on Pt(111) Reduction of NH on Pt(111) NH (gas) NH (ads) NH (ads) NH (ads) + H(ads) NH (ads) NH(ads) + H(ads) NH(ads) N(ads) + H(ads) N(ads) N (gas) H(ads) H (gas) 0 58 NH (ads) 9 NH (ads) NH(ads) N(ads) 7 H (gas) 9 6 N (gas) ffermans et al., Surf. Sci., 600, 1714 (006) 7
8 kmc snapshot Site blocking N NH vacant NH Coverages Coverages Coverages (T = 1000 K, P = 1.18 atm): kmc microkinetics Boltzmann NH NH NH N H /T Sulphate on Cu(111) with steps Sulphate on Cu(111) S + 4 +** S4 e adsorption sites blocked sites 8
9 Sulphate on Cu(111) with steps Sulphate on Cu(111) with steps current current potential (V) potential (V) Sulphate on Cu(111) with steps Sulphate on Cu(111) with steps This voltammogram is called the butterfly. Surface with defects! C.G.M. Hermse et al., Surf. Sci., 57, 4760 (004). Sulphate on Cu(111) with steps Bitartrate on Cu(110) current Steps suppress the phase transition. C HCH HCH C potential (V) Lorenzo, Baddeley, Muryn, Raval, Nature,, 76 (000) 9
10 Nonlinear response of LEED STM of mixtures 50:50 60:40 Haq, Liu, Humblot, Jansen, Raval, Nature Chemistry, 1, 409 (009) kmc simulations Domain formation of bisuccinate 50:50 60:40 time Chiral induction by bitartrate TPD of N from Ru(0001) Timedependent temperature and defects! 10
11 TPD of N from Ru(0001) TPD of N from Ru(0001) increased step density Temperature (K) Temperature (K) step terrace TPD of N from Ru(0001) Tmax (K) Step density increased step density C electrooxidation on PtRu H(ads) + H H (l) + * H(ads) + H + + e H (l) + * C(ads) + H(ads) C (gas) + * + H * H H * C + H * + + e + + e Initially the electrode is at low potential and almost completely covered by C. In the experiment the potential is increased and the current is measured. All reactions take place on Pt and Ru, but with different rate constants. C electrooxidation on PtRu The rate constants depend on the electrode potential. M rx transfer coefficient M ads M des = k = k M ads M des Me0E exp kbt exp ( 1 ) k 1,M M1,M M1 = krx exp k e T M B,M B electrode potential T e E E 0 0 C electrooxidation on PtRu more reactive n Pt the C formation is fast, but the water adsorption is slow. n Ru the C formation is slow, but the water adsorption is fast. If the C formation is fast at the PtRu boundary we get the observed synergic effect. 11
12 C electrooxidation on PtRu TPD of C/Rh(100) C H Pt Ru C prefers top sites. There is a strong NN repulsion. Ru Pt 0.50 ML C Rh M.T.M. Koper et al., J. Phys. Chem. B, 10, (1999). TPD of C/Rh(100) TPD of C/Rh(100) Some C molecules move to bridge sites to reduce repulsion. Some C molecules move to bridge sites to reduce repulsion ML 0.8 ML TPD of C/Rh(100) TPD of C/Rh(100) desorption prefactor adsorption energy top sites adsorption bridge sites NN toptop interaction NNN topbridge interaction geminal bridgebridge interaction vicinal bridgebridge interaction desorption prefactor /s adsorption energy top sites 106 kj/mol adsorption bridge sites 94 kj/mol NN toptop interaction 14. kj/mol NNN topbridge interaction.6 kj/mol e used Differential Evolution to get these parameters by fitting kmc results to experimental TPD spectra. geminal bridgebridge interaction 1.8 kj/mol vicinal bridgebridge interaction 0.6 kj/mol 1
13 TPD of C/Rh(100) TPD of C/Rh(100) C on bridge site C on top site M.M.M. Jansen et al., PCCP, 1, (010). For more information see Thank you for your attention 1
arxiv:chem-ph/ v2 11 May 1995
A Monte Carlo study of temperature-programmed desorption spectra with attractive lateral interactions. A.P.J. Jansen Laboratory of Inorganic Chemistry and Catalysis arxiv:chem-ph/9502009v2 11 May 1995
More information1 Adsorption of NO 2 on Pd(100) Juan M. Lorenzi, Sebastian Matera, and Karsten Reuter,
Supporting information: Synergistic inhibition of oxide formation in oxidation catalysis: A first-principles kinetic Monte Carlo study of NO+CO oxidation at Pd(100) Juan M. Lorenzi, Sebastian Matera, and
More informationCatalysis at the Sub-Nanoscale: Complex CO Oxidation Chemistry on a Few Au Atoms
Electronic Supplementary Material (ESI) for Catalysis Science & Technology. This journal is The Royal Society of Chemistry 214 Catalysis at the Sub-Nanoscale: Complex CO Oxidation Chemistry on a Few Au
More informationAn Introduction to Kinetic Monte Carlo Simulations of Surface Reactions
Lecture Notes in Physics 856 An Introduction to Kinetic Monte Carlo Simulations of Surface Reactions Bearbeitet von A.P.J. Jansen 1. Auflage 2012. Taschenbuch. xvii, 254 S. Paperback ISBN 978 3 642 29487
More informationKINETIC MONTE CARLO MOLECULAR SIMULATIONS FOR FUEL CELL APPLICATIONS AND SURFACE REACTIONS
KINETIC MONTE CARLO MOLECULAR SIMULATIONS FOR FUEL CELL APPLICATIONS AND SURFACE REACTIONS 1 A Thesis Defense In Partial Fulfillment of the Requirements for the Degree of Master of Science in Mechanical
More informationKinetic Monte Carlo Simulation of Molecular Processes on Supported Metal Particles
INFAP CONICET UNSL Kinetic Monte Carlo Simulation of Molecular Processes on Supported Metal Particles Giorgio Zgrablich Universidad Nacional de San Luis San Luis, Argentina CONTENTS Motivation Non-equilibrium
More informationCOMPUTATIONAL PHYSICS. Research: multi-tools to study physical and chemical properties of materials at different time and length scales.
COMPUTATIONAL PHYSICS Aniket Bhattacharya 1 GS + 1UGS Archana Dubey - Abdelkader Kara 1 UGS Talat Rahman 7 GS + 1 REU + 5.2 PostDocs Patrick Schelling 2 GS + 1UGS + 3 REU Sergey Stolbov 1GS ~ 30 papers
More informationChemical reactions as network of rare events: Kinetic MonteCarlo
Chemical reactions as network of rare events: Kinetic MonteCarlo Extending the scale Length (m) 1 10 3 Potential Energy Surface: {Ri} 10 6 (3N+1) dimensional 10 9 E Thermodynamics: p, T, V, N continuum
More informationAu-C Au-Au. g(r) r/a. Supplementary Figures
g(r) Supplementary Figures 60 50 40 30 20 10 0 Au-C Au-Au 2 4 r/a 6 8 Supplementary Figure 1 Radial bond distributions for Au-C and Au-Au bond. The zero density regime between the first two peaks in g
More informationIntroduction Statistical Thermodynamics. Monday, January 6, 14
Introduction Statistical Thermodynamics 1 Molecular Simulations Molecular dynamics: solve equations of motion Monte Carlo: importance sampling r 1 r 2 r n MD MC r 1 r 2 2 r n 2 3 3 4 4 Questions How can
More informationMechanisms of H- and OH-assisted CO activation as well as C-C coupling on the flat Co(0001) surface Revisited
Electronic Supplementary Material (ESI) for Catalysis Science & Technology. This journal is The Royal Society of Chemistry 2016 Mechanisms of H- and OH-assisted CO activation as well as C-C coupling on
More informationSupporting Online Material (1)
Supporting Online Material The density functional theory (DFT) calculations were carried out using the dacapo code (http://www.fysik.dtu.dk/campos), and the RPBE (1) generalized gradient correction (GGA)
More informationKinetic Monte Carlo (KMC) Kinetic Monte Carlo (KMC)
Kinetic Monte Carlo (KMC) Molecular Dynamics (MD): high-frequency motion dictate the time-step (e.g., vibrations). Time step is short: pico-seconds. Direct Monte Carlo (MC): stochastic (non-deterministic)
More informationCover Page. The handle holds various files of this Leiden University dissertation
Cover Page The handle http://hdl.handle.net/1887/29891 holds various files of this Leiden University dissertation Author: Roobol, Sander Bas Title: The structure of a working catalyst : from flat surfaces
More informationGAS-SURFACE INTERACTIONS
Page 1 of 16 GAS-SURFACE INTERACTIONS In modern surface science, important technological processes often involve the adsorption of molecules in gaseous form onto a surface. We can treat this adsorption
More informationEnantiospecific Adsorption of (R)-3-Methylcyclohexanone on Naturally Chiral Cu(531) R&S Surfaces
Catal Lett (2008) 125:177 182 DOI 10.1007/s10562-008-9600-8 Enantiospecific Adsorption of (R)-3-Methylcyclohexanone on Naturally Chiral Cu(531) R&S Surfaces Ye Huang Æ Andrew J. Gellman Received: 22 May
More informationarxiv: v1 [cond-mat.stat-mech] 15 Feb 2008
ADSORPTION-DESORPTION KINETICS: A REVIEW OF A CLASSICAL PROBLEM arxiv:0802.2226v1 [cond-mat.stat-mech] 15 Feb 2008 S. Manzi, R.E. Belardinelli, G. Costanza V. D. Pereyra Departamento de Física, Instituto
More informationChemical Kinetics. Topic 7
Chemical Kinetics Topic 7 Corrosion of Titanic wrec Casón shipwrec 2Fe(s) + 3/2O 2 (g) + H 2 O --> Fe 2 O 3.H 2 O(s) 2Na(s) + 2H 2 O --> 2NaOH(aq) + H 2 (g) Two examples of the time needed for a chemical
More informationComputer simulation methods (1) Dr. Vania Calandrini
Computer simulation methods (1) Dr. Vania Calandrini Why computational methods To understand and predict the properties of complex systems (many degrees of freedom): liquids, solids, adsorption of molecules
More informationExtending the Tools of Chemical Reaction Engineering to the Molecular Scale
Extending the Tools of Chemical Reaction Engineering to the Molecular Scale Multiple-time-scale order reduction for stochastic kinetics James B. Rawlings Department of Chemical and Biological Engineering
More informationSurface Complexes in Catalysis
Surface Complexes in Catalysis David Karhánek Ústav organické technologie, VŠCHT Praha Institut für Materialphysik, Universität Wien XXXVII Symposium on Catalysis, Prague, October 7-8, 2005. Research Methodologies:
More informationAcidic Water Monolayer on Ruthenium(0001)
Acidic Water Monolayer on Ruthenium(0001) Youngsoon Kim, Eui-seong Moon, Sunghwan Shin, and Heon Kang Department of Chemistry, Seoul National University, 1 Gwanak-ro, Seoul 151-747, Republic of Korea.
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION Using first principles to predict bimetallic catalysts for the ammonia decomposition reaction Danielle A. Hansgen, Dionisios G. Vlachos, Jingguang G. Chen SUPPLEMENTARY INFORMATION.
More informationDynamical Monte-Carlo Simulation of Surface Kinetics
Dynamical Monte-Carlo Simulation of Surface Kinetics V. Guerra and J. Loureiro Centro de Física dos Plasmas, Instituto Superior Técnico, 1049-001 Lisboa, Portugal Abstract. This wor presents a theoretical
More information2. Thermodynamics. Introduction. Understanding Molecular Simulation
2. Thermodynamics Introduction Molecular Simulations Molecular dynamics: solve equations of motion r 1 r 2 r n Monte Carlo: importance sampling r 1 r 2 r n How do we know our simulation is correct? Molecular
More informationEnantioselective Separation on a Naturally Chiral Surface
Published on Web 10/26/2004 Enantioselective Separation on a Naturally Chiral Surface Joshua D. Horvath, Anjanette Koritnik, Preeti Kamakoti,, David S. Sholl,, and Andrew J. Gellman*, Contribution from
More informationTwo simple lattice models of the equilibrium shape and the surface morphology of supported 3D crystallites
Bull. Nov. Comp. Center, Comp. Science, 27 (2008), 63 69 c 2008 NCC Publisher Two simple lattice models of the equilibrium shape and the surface morphology of supported 3D crystallites Michael P. Krasilnikov
More informationMonte Carlo simulation of thin-film growth on a surface with a triangular lattice
Vacuum 52 (1999) 435 440 Monte Carlo simulation of thin-film growth on a surface with a triangular lattice Wei Helin*, Liu Zuli, Yao Kailun Department of Physics, Huazhong University of Science and Technology,
More informationMaster Equation for Hydrogen Recombination on Grain Surfaces
Syracuse University SURFACE Physics College of Arts and Sciences 12-12-2000 Master Equation for Hydrogen Recombination on Grain Surfaces Gianfranco Vidali Department of Physics, Syracuse University, Syracuse,
More informationEnantiospecific Desorption of Chiral Compounds from Chiral Cu(643) and Achiral Cu(111) Surfaces
Published on Web 02/13/2002 Enantiospecific Desorption of Chiral Compounds from Chiral Cu(643) and Achiral Cu(111) Surfaces Joshua D. Horvath and Andrew J. Gellman* Contribution from the Department of
More informationFirst-principles based design of Pt- and Pd-based catalysts for benzene hydrogenation
1 1 First-principles based design of Pt- and Pd-based catalysts for benzene hydrogenation Maarten K. Sabbe, Gonzalo Canduela, Marie- Françoise Reyniers, Guy B. Marin Introduction: benzene hydrogenation
More informationIntroduction. Monday, January 6, 14
Introduction 1 Introduction Why to use a simulation Some examples of questions we can address 2 Molecular Simulations Molecular dynamics: solve equations of motion Monte Carlo: importance sampling Calculate
More informationKinetic Monte Carlo simulations of the interaction of oxygen with Pt(111)
Kinetic Monte Carlo simulations of the interaction of oxygen with Pt(111) Christian Sendner 1 and Axel Groß 2 1 Physik-Department T30, Technische Universität München, D-85747 Garching, Germany 2 Institut
More informationAdsorption, desorption, and diffusion on surfaces. Joachim Schnadt Divsion of Synchrotron Radiation Research Department of Physics
Adsorption, desorption, and diffusion on surfaces Joachim Schnadt Divsion of Synchrotron Radiation Research Department of Physics Adsorption and desorption Adsorption Desorption Chemisorption: formation
More informationCatalyst structure and C-O activation during FTS: new ideas from computational catalysis. Mark Saeys
Catalyst structure and C-O activation during FTS: new ideas from computational catalysis Mark Saeys Laboratory for Chemical Technology, Ghent University Laboratory for Chemical http://www.lct.ugent.be
More informationCatalysis Today. Effects of heat and mass transfer on the kinetics of CO oxidation over RuO 2 (1 1 0) catalyst. Donghai Mei a,, Guang Lin b.
Catalysis Today 165 (2011) 56 63 Contents lists available at ScienceDirect Catalysis Today journal homepage: www.elsevier.com/locate/cattod Effects of heat and mass transfer on the kinetics of oxidation
More informationFEEDBACK CONTROL OF GROWTH RATE AND SURFACE ROUGHNESS IN THIN FILM GROWTH. Yiming Lou and Panagiotis D. Christofides
FEEDBACK CONTROL OF GROWTH RATE AND SURFACE ROUGHNESS IN THIN FILM GROWTH Yiming Lou and Panagiotis D. Christofides Department of Chemical Engineering University of California, Los Angeles IEEE 2003 Conference
More informationKinetic Monte Carlo: from transition probabilities to transition rates
Kinetic Monte Carlo: from transition probabilities to transition rates With MD we can only reproduce the dynamics of the system for 100 ns. Slow thermallyactivated processes, such as diffusion, cannot
More informationKinetic Monte Carlo (KMC)
Kinetic Monte Carlo (KMC) Molecular Dynamics (MD): high-frequency motion dictate the time-step (e.g., vibrations). Time step is short: pico-seconds. Direct Monte Carlo (MC): stochastic (non-deterministic)
More informationRate of Heating and Cooling
Rate of Heating and Cooling 35 T [ o C] Example: Heating and cooling of Water E 30 Cooling S 25 Heating exponential decay 20 0 100 200 300 400 t [sec] Newton s Law of Cooling T S > T E : System S cools
More informationTOPIC 6: Chemical kinetics
TOPIC 6: Chemical kinetics Reaction rates Reaction rate laws Integrated reaction rate laws Reaction mechanism Kinetic theories Arrhenius law Catalysis Enzimatic catalysis Fuente: Cedre http://loincognito.-iles.wordpress.com/202/04/titanic-
More informationReaction Dynamics (2) Can we predict the rate of reactions?
Reaction Dynamics (2) Can we predict the rate of reactions? Reactions in Liquid Solutions Solvent is NOT a reactant Reactive encounters in solution A reaction occurs if 1. The reactant molecules (A, B)
More informationarxiv: v1 [cond-mat.mes-hall] 8 Jan 2015
Multi-Lattice Kinetic Monte Carlo Simulations from First-Principles: Reduction of the Pd(100) Surface Oxide by CO arxiv:1501.01710v1 [cond-mat.mes-hall] 8 Jan 2015 Abstract Max J. Hoffmann,, Matthias Scheffler,
More informationSupporting Information
Supporting Information Yao et al. 10.1073/pnas.1416368111 Fig. S1. In situ LEEM imaging of graphene growth via chemical vapor deposition (CVD) on Pt(111). The growth of graphene on Pt(111) via a CVD process
More informationSelectivity in the initial C-H bond cleavage of n-butane on PdO(101)
Supporting Information for Selectivity in the initial C-H bond cleavage of n-butane on PdO(101) Can Hakanoglu (a), Feng Zhang (a), Abbin Antony (a), Aravind Asthagiri (b) and Jason F. Weaver (a) * (a)
More informationKinetic modeling in heterogeneous catalysis.
Kinetic modeling in heterogeneous catalysis. C. Franklin Goldsmith Modern Methods in Heterogeneous Catalysis Research Friday, January 27 th, 2012 Microkinetic modeling is a valuable tool in catalysis.
More informationHigh-Pressure NO-Induced Mixed Phase on Rh(111): Chemically Driven Replacement
Supporting Information for High-Pressure NO-Induced Mixed Phase on Rh(111): Chemically Driven Replacement Ryo Toyoshima, Masaaki Yoshida, Yuji Monya, Kazuma Suzuki, Kenta Amemiya, Kazuhiko Mase, Bongjin
More informationPart III: Theoretical Surface Science Adsorption at Surfaces
Technische Universität München Part III: Theoretical Surface Science Adsorption at Surfaces Karsten Reuter Lecture course: Solid State Theory Adsorption at surfaces (T,p) Phase II Phase I Corrosion Growth
More information1 Molecular Catalytic Kinetics Concepts
1 1 Molecular Catalytic Kinetics Concepts Rutger A. van Santen 1.1 Key Principles of Heterogeneous Catalysis We discuss the following topics in the subsequent sections: Sabatier principle and volcano curve;
More informationMultiscale Materials Modeling
Multiscale Materials Modeling Lecture 02 Capabilities of Classical Molecular Simulation These notes created by David Keffer, University of Tennessee, Knoxville, 2009. Outline Capabilities of Classical
More informationEfficient Synthesis of Ethanol from CH 4 and Syngas on
Efficient Synthesis of Ethanol from CH 4 and Syngas on a Cu-Co/TiO 2 Catalyst Using a Stepwise Reactor Zhi-Jun Zuo 1, Fen Peng 1,2, Wei Huang 1,* 1 Key Laboratory of Coal Science and Technology of Ministry
More informationFirst-principles based catalytic reaction engineering Matteo Maestri
CECAM International Summer School Hot topic 4 First-principles based catalytic reaction engineering Matteo Maestri July 23, 2013 Conversationshaus - Norderney, Germany Catalytic cycle Consists of the elementary
More informationUnderstanding Chemical Reactions through Computer Modeling. Tyler R. Josephson University of Delaware 4/21/16
Understanding Chemical Reactions through Computer Modeling Tyler R. Josephson University of Delaware 4/21/16 A little about me B.S. in Chem E from U of M, 2011 Currently, Ph.D. student at University of
More informationSupporting Information. DFT Study of Methane Synthesis from Syngas on Ce Doped Ni(111) Surface
Supporting Information DFT Study of Methane Synthesis from Syngas on Ce Doped Ni(111) Surface Kai Li, 1 Cong Yin, 2 Yi Zheng, 3 Feng He, 1 Ying Wang, 1 Menggai Jiao, 1 Hao Tang, 2,* Zhijian Wu 1,* 1 State
More informationAUTOMOTIVE EXHAUST AFTERTREATMENT
AUTOMOTIVE EXHAUST AFTERTREATMENT CATALYST FUNDAMENTLS Catalyst in its simplest term is a material that increase the rate (molecules converted by unit time) of a chemical reaction while itself not undergoing
More informationOxide Formation of Transition Metal Surfaces and Effect on Catalysis
Oxide Formation of Transition Metal Surfaces and Effect on Catalysis Wei-Xue Li State Key Laboratory of Catalysis, and Center for Theoretical and Computational Chemistry Dalian Institute of Chemical Physics,
More informationSupporting Information. First principles kinetic study on the effect of zeolite framework on 1-butanol dehydration
Supporting Information First principles kinetic study on the effect of zeolite framework on 1-butanol dehydration Mathew John, Konstantinos Alexopoulos, Marie-Françoise Reyniers* and Guy B. Marin Laboratory
More informationTheodore E. Madey. Department of Physics and Astronomy, and Laboratory for Surface Modification
The Science of Catalysis at the Nanometer Scale Theodore E. Madey Department of Physics and Astronomy, and Laboratory for Surface Modification http://www.physics.rutgers.edu/lsm/ Rutgers, The State University
More informationKinetic Monte Carlo: Coarsegraining
Kinetic Monte Carlo: Coarsegraining Time and Space Peter Kratzer Faculty of Physics, University Duisburg-Essen, Germany morphology Time and length scales Ga As 2D islands surface reconstruction Methods
More information3.320: Lecture 19 (4/14/05) Free Energies and physical Coarse-graining. ,T) + < σ > dµ
3.320: Lecture 19 (4/14/05) F(µ,T) = F(µ ref,t) + < σ > dµ µ µ ref Free Energies and physical Coarse-graining T S(T) = S(T ref ) + T T ref C V T dt Non-Boltzmann sampling and Umbrella sampling Simple
More informationAb initio statistical mechanics of surface adsorption and desorption. I. H 2 O on MgO 001 at low coverage
THE JOURNAL OF CHEMICAL PHYSICS 127, 114709 2007 Ab initio statistical mechanics of surface adsorption and desorption. I. H 2 O on MgO 001 at low coverage D. Alfè London Centre for Nanotechnology, University
More informationSurface Physics Surface Diffusion. Assistant: Dr. Enrico Gnecco NCCR Nanoscale Science
Surface Physics 008 8. Surface Diffusion Assistant: Dr. Enrico Gnecco NCCR Nanoscale Science Random-Walk Motion Thermal motion of an adatom on an ideal crystal surface: - Thermal excitation the adatom
More informationUnderstanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles. Srikanth Sastry
JNCASR August 20, 21 2009 Understanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles Srikanth Sastry Jawaharlal Nehru Centre for Advanced Scientific Research, Bangalore
More informationCatalytic Activity of IrO 2 (110) Surface: A DFT study
Catalytic Activity of IrO 2 (110) Surface: A DFT study Jyh-Chiang Jiang Department of Chemical Engineering, National Taiwan University of Science and Technology (NTUST) NCTS-NCKU 9/7, 2010 Computational
More informationBrief Review of Statistical Mechanics
Brief Review of Statistical Mechanics Introduction Statistical mechanics: a branch of physics which studies macroscopic systems from a microscopic or molecular point of view (McQuarrie,1976) Also see (Hill,1986;
More informationQuantum wells and Dots on surfaces
Lecture in the course Surface Physics and Nano Physics 2008 Quantum wells and Dots on surfaces Bo Hellsing Department of Physics, Göteborg University, Göteborg, S Collaborators: QW Johan Carlsson, Göteborg
More informationA Brief Introduction to Statistical Mechanics
A Brief Introduction to Statistical Mechanics E. J. Maginn, J. K. Shah Department of Chemical and Biomolecular Engineering University of Notre Dame Notre Dame, IN 46556 USA Monte Carlo Workshop Universidade
More informationMicrokinetic Modeling
Technische Universität München Microkinetic Modeling Karsten Reuter Chemistry Department and Catalysis Research Center Technische Universität München Catalysis Research Triangle Testing (Kinetics) Preparation
More informationEvidence for structure sensitivity in the high pressure CO NO reaction over Pd(111) and Pd(100)
Evidence for structure sensitivity in the high pressure CO NO reaction over Pd(111) and Pd(100) Scott M. Vesecky, Peijun Chen, Xueping Xu, and D. Wayne Goodman a) Department of Chemistry, Texas A&M University,
More informationStability, Composition and Function of Palladium Surfaces in Oxidizing Environments: A First-Principles Statistical Mechanics Approach
Stability, Composition and Function of Palladium Surfaces in Oxidizing Environments: A First-Principles Statistical Mechanics Approach von Diplom-Chemikerin Jutta Rogal im Fachbereich Physik der Freien
More informationTitle. Author(s)Yamanaka, Toshiro; Matsushima, Tatsuo. CitationPhysical Review Letters, 100(2): Issue Date Doc URL. Rights.
Title Polar- and Azimuth-Angle-Dependent Rotational and Vi Oxidation on Palladium Surfaces Author(s)Yamanaka, Toshiro; Matsushima, Tatsuo CitationPhysical Review Letters, 100(2): 026104 Issue Date 2008-01-17
More informationMaterial Surfaces, Grain Boundaries and Interfaces: Structure-Property Relationship Predictions
Material Surfaces, Grain Boundaries and Interfaces: Structure-Property Relationship Predictions Susan B. Sinnott Department of Materials Science and Engineering Penn State University September 16, 2016
More informationWhen atomic-scale resolution is not enough: Heat and mass transfer effects in in-situ model catalyst studies
Technische Universität München When atomic-scale resolution is not enough: Heat and mass transfer effects in in-situ model catalyst studies Karsten Reuter Chemistry Department and Catalysis Research Center
More informationTable 1: Residence time (τ) in seconds for adsorbed molecules
1 Surfaces We got our first hint of the importance of surface processes in the mass spectrum of a high vacuum environment. The spectrum was dominated by water and carbon monoxide, species that represent
More informationUnderstanding the enhancement effect of high-temperature steam on the carbonation reaction of CaO with CO 2
5th IEAGHG Network Meeting and Technical Workshop n High Temperature Solid Looping Cycles Cambridge University 2th-3th of September, 2013 Understanding the enhancement effect of high-temperature steam
More informationChemisorption VIII. NEVF 514 Surface Physics. Winter Term Troja, 16th December 2016
Chemisorption František Máca VIII. NEVF 514 Surface Physics Winter Term 2016-2017 Troja, 16th December 2016 Chemisorption The knowledge of chemisorption phenomena requires the determination of the geometrical
More informationDistinguish quantitatively between the adsorption isotherms of Gibbs, Freundlich and Langmuir.
Module 8 : Surface Chemistry Lecture 36 : Adsorption Objectives After studying this lecture, you will be able to Distinguish between physisorption and chemisorption. Distinguish between monolayer adsorption
More informationSupplementary Fig. 1. Progress of the surface mediated Ullmann coupling reaction using STM at 5 K. Precursor molecules
Supplementary Fig. 1. Progress of the surface mediated Ullmann coupling reaction using STM at 5 K. Precursor molecules (4-bromo-1-ethyl-2-fluorobenzene) are dosed on a Cu(111) surface and annealed to 80
More informationTheoretical Calculations of Electrochemical Ammonia Synthesis at Ambient Pressure and Temperature
Theoretical Calculations of Electrochemical Ammonia Synthesis at Ambient Pressure and Temperature Egill Skúlason 1,2, Thomas Bligaard 1,2, Jan Rossmeisl 2, Áshildur Logadóttir 2, Jens K. Nørskov 2, Hannes
More informationKinetic Monte Carlo Simulation of Two-dimensional Semiconductor Quantum Dots Growth
Kinetic Monte Carlo Simulation of Two-dimensional Semiconductor Quantum Dots Growth by Ernie Pan Richard Zhu Melissa Sun Peter Chung Computer Modeling and Simulation Group The University of Akron Outline
More information4. Systems in contact with a thermal bath
4. Systems in contact with a thermal bath So far, isolated systems microcanonical methods 4.1 Constant number of particles:kittelkroemer Chap. 3 Boltzmann factor Partition function canonical methods Ideal
More informationFormation and dynamics of water clusters on Ru 001
JOURNAL OF CHEMICAL PHYSICS VOLUME 119, NUMBER 22 8 DECEMBER 2003 Formation and dynamics of water clusters on Ru 001 Y. Lilach, a) V. Buch, b) and M. Asscher c) Department of Physical Chemistry, The Hebrew
More informationIntroduction. Statistical physics: microscopic foundation of thermodynamics degrees of freedom 2 3 state variables!
Introduction Thermodynamics: phenomenological description of equilibrium bulk properties of matter in terms of only a few state variables and thermodynamical laws. Statistical physics: microscopic foundation
More informationTheoretical Models for Chemical Kinetics
Theoretical Models for Chemical Kinetics Thus far we have calculated rate laws, rate constants, reaction orders, etc. based on observations of macroscopic properties, but what is happening at the molecular
More informationAb initio study of CNT NO 2 gas sensor
Chemical Physics Letters 387 (2004) 271 276 www.elsevier.com/locate/cplett Ab initio study of CNT NO 2 gas sensor Shu Peng a, *, Kyeongjae Cho a, Pengfei Qi b, Hongjie Dai b a Department of Mechanical
More informationTutorial on rate constants and reorganization energies
www.elsevier.nl/locate/jelechem Journal of Electroanalytical Chemistry 483 (2000) 2 6 Tutorial on rate constants reorganization energies R.A. Marcus * Noyes Laboratory of Chemical Physics, MC 127-72, California
More informationHandout 10. Applications to Solids
ME346A Introduction to Statistical Mechanics Wei Cai Stanford University Win 2011 Handout 10. Applications to Solids February 23, 2011 Contents 1 Average kinetic and potential energy 2 2 Virial theorem
More informationKinetic Monte Carlo modelling of semiconductor growth
Kinetic Monte Carlo modelling of semiconductor growth Peter Kratzer Faculty of Physics, University Duisburg-Essen, Germany Time and length scales morphology Ga As 2D islands surface reconstruction Methods
More informationIn#uence of microstructure of substrate surface on early stage of thin "lm growth
Vacuum 56 (2000) 185}190 In#uence of microstructure of substrate surface on early stage of thin "lm growth Helin Wei*, Zuli Liu, Kailun Yao Department of Physics, Huazhong University of Science and Technology.
More informationExperimental and Quantum Investigation on Ice Surface Structure and Reactivity
Experimental and Quantum Investigation on Ice Surface Structure and Reactivity A.Allouche J.P.Aycard F.Borget T.Chiavassa I.Couturier C.Manca F.Marinelli C.Martin S.Raunier P.Roubin Physique des Interactions
More informationIncluding lateral interactions into microkinetic models of catalytic reactions
Downloaded from orbit.dtu.dk on: Sep 05, 2018 Including lateral interactions into microkinetic models of catalytic reactions Hellman, Anders; Honkala, Johanna Karoliina Published in: Journal of Chemical
More informationFinite-size poisoning in heterogeneous catalysis
J. Phys. A: Math. Gen. 23 (1990) L613-J-619. Printed in the UK LE ITER TO THE EDITOR Finite-size poisoning in heterogeneous catalysis D ben-avrahamt9, S Rednert, D B Considinet and P MeakinS t Center for
More informationAdsorption and Catalysis. Radha V Jayaram Institute of Chemical Technology, Mumbai
Adsorption and Catalysis Radha V Jayaram Institute of Chemical Technology, Mumbai 1 2 Different solids would adsorb different amounts of the same gas even under similar conditions. Substances like charcoal
More informationKinetic Monte Carlo simulation of nucleation on patterned substrates
PHYSICAL REVIEW B, VOLUME 63, 035407 Kinetic Monte Carlo simulation of nucleation on patterned substrates L. Nurminen, A. Kuronen, and K. Kaski Helsinki University of Technology, Laboratory of Computational
More informationCleaner Car Exhausts Using Ion-Exchanged Zeolites: Insights From Atomistic Simulations
CLEERS conference, Ann Arbor, Michigan 20 th September 2018 Cleaner Car Exhausts Using Ion-Exchanged Zeolites: Insights From Atomistic Simulations A Combined Quasi Elastic Neutron Scattering (QENS) and
More informationMajor Concepts Lecture #11 Rigoberto Hernandez. TST & Transport 1
Major Concepts Onsager s Regression Hypothesis Relaxation of a perturbation Regression of fluctuations Fluctuation-Dissipation Theorem Proof of FDT & relation to Onsager s Regression Hypothesis Response
More informationFinal Exam Chem 260 (12/20/99) Name (printed) Signature
Final Exam Chem 260 (12/20/99) Name (printed) Signature 1. For the reaction, (14 pts total) C(graph) + CO 2 (g) 2 CO(g) (a) Write the equilibrium constant expression: (2 pts) K = P CO 2 P CO2 (b) Using
More informationDocument Version Publisher s PDF, also known as Version of Record (includes final page, issue and volume numbers)
On the temperature dependence of the Arrhenius activation energy for hydroisomerization catalyzed by Pt/Mordenite Runstraat, van de, A.; van Grondelle, J.; van Santen, R.A. Published in: Journal of Catalysis
More informationRate Law Summary. Rate Laws vary as a function of time
Rate Law Summary Measure the instantaneous rate of a reaction: this is a number with units of M/s! Measure the rate of loss of a reactant r... the rate of appearance of a product Repeat the experiment
More informationExpress the transition state equilibrium constant in terms of the partition functions of the transition state and the
Module 7 : Theories of Reaction Rates Lecture 33 : Transition State Theory Objectives After studying this Lecture you will be able to do the following. Distinguish between collision theory and transition
More information