when atomic-scale resolution is not enough: heat and mass...
TRANSCRIPT
![Page 1: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/1.jpg)
Technische Universität München
Karsten Reuter
Chemistry Department and Catalysis Research Center Technische Universität München
When atomic-scale resolution is not enough:
Heat and mass transfer effects in in-situ model catalyst studies
![Page 2: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/2.jpg)
Testing/Validation (Kinetics)
Preparation Characterization
Catalysis research triangle
![Page 3: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/3.jpg)
Challenges across the scales
Quantitative transient and steady-state surface kinetics
Self-consistent coupling to reactive flow field and appropriate
heat balance
Accurate (first-principles) energetics of individual elementary processes
Surface chemistry: adsorption, diffusion, reaction, desorption
![Page 4: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/4.jpg)
The Surface Science disclaimer: Mind the gap(s)!
Supported „real“ catalyst
Model catalysts
Courtesy: G. Rupprechter and Ch. Weiland, NanoToday 2, 20 (2007)
![Page 5: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/5.jpg)
Going in situ: It’s not just like opening a valve…
A. Stierle and A.M. Molenbroek (Eds.), MRS Bull. 32 (2007)
P.B. Rasmussen et al., Rev. Sci. Instrum. 69, 3879 (1998)
Reactor STM
S. Yamamoto et al., J. Phys. CM 20, 184025 (2008)
in situ XPS
![Page 6: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/6.jpg)
Part I First-principles microkinetics Part II Integrating 1p-kMC into CFD
A tour from electrons to the reactor
![Page 7: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/7.jpg)
Potential energy for fixed nuclear positions {RI}:
The “Quantum Engines”: Electronic structure methods
E0 ({RI }) = MinΨ <Ψ| He{RI} |Ψ>
⇒ E0 ({RI }) = Minn(r) E{RI}[n] Ψ = Ψ[n(r)] (Hohenberg-Kohn, 1964) 2. Density-functional theory
1. Wavefunction based methods („Quantum Chemistry“)
Ansatz for Ψ: Hartree-Fock, post-HF (MP2, CI, CC,…)
(Kohn-Sham, 1965)
)(
3 )( )( ][rnno
XCXCLDA
o
nrnrdnE=∫= ε
)(
3 ),()( ][rnnoo
XCXCGGA
o
nnrnrdnE=
∇= ∫ ε
![Page 8: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/8.jpg)
Electronic regime: Energetics of elementary processes
COcus + Ocus → CO2
- Active site model
- „Level of theory“
0.9 eV
M. Sabbe, M.F. Reyniers, and K. Reuter, Catal. Sci. Technol. 2, 2010 (2012)
![Page 9: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/9.jpg)
„The“ active site: The effective to atomistic gap
* vs.
• Generic active site independent of operation conditions
dominant over other sites (one site model) • Lately: move towards „multi-site“
models (multi = 2,…)
• Atomistic active site model
every atom counts generally insufficient experimental characterization
• Manifold of possibly active site types
on which one to focus? consider how many?
![Page 10: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/10.jpg)
A
B
Molecular Dynamics
TS
Mesoscopic regime: Tackling rare-event time scales
kA→B
kB→A
kinetic Monte Carlo
N
t
B A
∑∑ →→ +−=j
jijj
ijii tPktPkdt
tdP )()()(
ΔEA→B ΔEB→A
∆−Γ=
=
→
→→
TkEZ
ZhTkk
ji
i
jiji
B
)(TSB
exp
Transition State Theory
First-principles kinetic Monte Carlo simulations for heterogeneous catalysis: Concepts, status and frontiers K. Reuter, in “Modeling Heterogeneous Catalytic Reactions: From the Molecular Process to the Technical System”,
(Ed.) O. Deutschmann, Wiley-VCH, Weinheim (2011). http://www.fhi-berlin.mpg.de/th/paper.html
![Page 11: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/11.jpg)
Mean-field approximation: Rate equations
∑∑ →→ +−=j
jijj
ijii tPktPkdt
tdP)()(
)(
),CO( ),O(
),CO,O(
cusbr
cusbr
tt
tP
θθ ⋅
↓
θ(Ocus) θ(Obr) θ(COcus) θ(CObr)
{ }
{ }
),CO(),,CO(),,O(),,O(, ),O(
),CO(),,CO(),,O(),,O(, ),O(
brcusbrcus2
br
brcusbrcus1
cus
ttttkfdt
td
ttttkfdt
td
ji
ji
θθθθθ
θθθθθ
→
→
=
=
![Page 12: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/12.jpg)
Ru O
CO oxidation at RuO2(110): 1p-kMC model
26 elementary processes (site-specific):
- O2 adsorption/desorption (dissociative/associative) - CO adsorption/desorption (unimolecular) - O and CO diffusion - CO + O reaction
K. Reuter and M. Scheffler, Phys. Rev. B 73, 045433 (2006)
O
CO
K. Reuter, Oil&Gas Sci. Technol. 61, 471 (2006)
![Page 13: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/13.jpg)
600 K
CObr/-
pO (atm) 2
1
CObr/COcus
Obr/Ocus p C
O (
atm
)
1
10-5
105
10-5 10+5 10-15 10-10
Obr/ -
CObr/COcus
Obr/Ocus
Surface structure and composition in the reactive environment
K. Reuter, D. Frenkel and M. Scheffler, Phys. Rev. Lett. 93, 116105 (2004)
T = 600 K, pO = 1 atm, pCO = 7 atm 2
CO oxidation at RuO2(110)
![Page 14: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/14.jpg)
TPR
Steady-state and transient parameter-free turnover frequencies
x
x x
x x
x x x x x x
x
x
x x x
x x x
pCO (10-9 atm)
Exp.
Theory
0.0 1.0 2.0 3.0 6
4
2
0 TOF C
O2 (
1012
mol
/cm
2 s)
pO2 = 10-10 atm
350 K
M. Rieger, J. Rogal, and K. Reuter, Phys. Rev. Lett. 100, 016105 (2008)
K. Reuter and M. Scheffler, Phys. Rev. B 73, 045433 (2006)
![Page 15: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/15.jpg)
Key ingredients to „predictive-quality“ microkinetic modeling
Accurate rate constants: Transition state theory and beyond DFT functionals: „self-interaction“ van der Waals interactions Reaction mechanism: Process identification Lattice mapping / spatial distributions „Hot chemistry“ beyond Markov
CO + * ↔ CO* O2 + 2* ↔ O* + O* CO* + O* ↔ CO2 …
∆−Γ= →
→ TkE
k jiji
B
exp
![Page 16: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/16.jpg)
[CO]Rh + [OH]Rh+ [H]Rh
CO2
H2
[CO2]Rh +2[H]Rh[H]Rh +[OH]Rh + Rh
[H2O]Rh + 2Rh
CO
H2O
RDS
COOH*+* → CO*+OH* CO*+OH* → COOH*+* COOH*+* → CO2*+H* CO2*+H* → COOH*+* CO2* + H2O* → COOH* + OH* COOH* + OH* → CO2* + H2O* CO2*+H* → HCOO** HCOO** → CO2*+H* CO2* + OH* + *→ HCOO** + O* HCOO** + OH* → CO2* + H2O* CH* + H* → CH2* + * CH* + * → C* + H* C* + H* → CH* + * CH3* + O* → CH2* + OH* CH2* + OH* → CH3* + O* CH* + OH* → CH2* + O* CH2* + O* → CH* + OH* …
0
0.5
1
1.5
2
0 2 4 6 8 10
mm
ol/s
Axial length [mm]
H2
CO
CH4
O2
H2O CO
2
400
500
600
700
800
900
1000
1100
Tem
pera
ture
[°C
] T catalyst
T gas
Macroscale Microscale
~ 102 potential steps @ different coverages
„E pluribus unum“: Water-gas-shift at Rh(111)
M. Maestri and K. Reuter, Chem. Eng. Sci. 74, 296 (2012)
![Page 17: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/17.jpg)
Error propagation through rate-determining steps
C.T. Campbell, J. Catal. 204, 520 (2001); Nature 432, 282 (2004)
Sensitivity analysis: pO2 = 1 atm T = 600 K
COcus + Ocus
Obr/Ocus ad
Ocus/Ocus ad
CObr ad
pCO (atm)
Χrc
,i
COcus +Obr
O-poisoned CO-poisoned
H. Meskine et al., Ertl Special Issue
Surf. Sci. 603, 1724 (2009)
ij Kki
ii k
kX,
rc,TOF
TOF
∂
∂
=
CO oxidation at RuO2(110)
![Page 18: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/18.jpg)
More than Markov?
( ) 2
~ B
uc
TmkpA
TSkπ=
Z=1.5 Å
J. Meyer and K. Reuter (in preparation)
Heat dissipation during dissociative adsorption: O2/Pd(100)
![Page 19: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/19.jpg)
Part I First-principles microkinetics Part II Integrating 1p-kMC into CFD
A tour from electrons to the reactor
![Page 20: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/20.jpg)
Stagnation flow: Heat and mass transfer
T, pCO, pO2
T
pCO2
p
pO2
pCO
![Page 21: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/21.jpg)
Self-consistent coupling to the flow field
Computational Fluid Dynamics with
chemical source terms from 1p-kMC ?!
( ) 2
B
ucad
TmkpATSkπ
=
![Page 22: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/22.jpg)
Linking the scales: Quasi-instantaneous steady-state approximation
On time scales relevant to changes in flow profile surface chemistry
adapts quasi-instantaneously to new steady-state
e.g. O. Deutschmann, “Computational Fluid Dynamics Simulation of Catalytic Reactors”, In: Handbook of Heterogeneous Catalysis (2008)
Precompute steady-state 1p-kMC TOFs on dense grid in (T, pO2, pCO)-space
Interpolate using modified Shepard (ideal for grid-less scattered data)
Run CFD using continuous representation as boundary condition
S. Matera and K. Reuter, Phys. Rev. B 82, 085446 (2010)
![Page 23: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/23.jpg)
Convection + Diffusion • segregated approach
• discretization on unstructured meshes
OpenFOAM® Framework
Reactions • coupled approach • stiff ODE solvers
CatalyticFOAM Framework
Paraview ® Post-processing
Pre-Processing • OpenFOAM® BlockMesh • GAMBIT, Ansys, CFX, etc.
www.catalyticfoam.polimi.it
The CatalyticFOAM project
M. Maestri A. Cuoci
![Page 24: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/24.jpg)
In situ TOFs in the isothermal limit: Mass transfer limitations
S. Matera and K. Reuter, Catal. Lett. 133, 156 (2009)
T = const uinl = 1 cm/sec
T = 600 K pO2 = 0.3 atm
![Page 25: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/25.jpg)
When atomic scale resolution is not enough
T = 600 K, uinl(max) = 20 cm/sec pO2 = 0.3 atm, pCO = 1.9 atm
S. Matera and K. Reuter, J. Catal. 295, 261 (2012)
![Page 26: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/26.jpg)
Planar Laser Induced Fluorescence (PLIF)
Laser-sheet stimulation of known excitation (here: CO2 vibration)
→ 2D concentration profile above catalyst
Y. Zetterberg et al., Rev. Sci. Instrum. 83, 053104 (2012)
![Page 27: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/27.jpg)
Making mass transfer limitations „visible“
S. Matera et al., in preparation
in situ PLIF measurements of ambient CO oxidation at Pd(100)
together with E. Lundgren et al. (Lund University)
S. Blomberg, M.J. Hoffmann et al., Phys. Rev. Lett. 110, 117601 (2013)
CO:O2 = 1:1, ptot = 1 atm
![Page 28: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/28.jpg)
Incomplete heat dissipation in a real stagflow reactor:
Rayleigh-Bénard convection
CO:O2 = 4:1, ptot = 1 atm Tinlet = 500K, adiabatic limit
S. Döpking, S. Matera
CO oxidation at RuO2(110)
![Page 29: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/29.jpg)
First-Principles Multiscale Modeling: Where do we stand?
State-of-the-art in catalysis modeling: - Prevalence of highly coarse-grained models based on effective parameters without true microscopic meaning rate equation theory based on empirical rate constants 1D flow models - Emergence of ad-hoc 1p-microkinetic models kMC and mean-field for model catalysts & show case reactions
Steps towards a predictive character multiscale catalysis modeling: - Replace effective parameters by clean first-principles data fitted vs. DFT-based rate constants battle the curse of complexity (off-lattice, complex networks) electronic non-adiabaticity, heat dissipation - Refined modeling at each individual level reliable and efficient 1p-rate constants (where needed) necessity to resolve spatial arrangement at surface integrate 1p-surface chemistry into detailed reactor models - Robust links between theories that enable reverse-mapping sensitivity analysis to control flow of error across scales
![Page 30: When atomic-scale resolution is not enough: Heat and mass ...helper.ipam.ucla.edu/publications/msews2/msews2_10939.pdfexp Transition State Theory . First-principles kinetic Monte Carlo](https://reader036.vdocument.in/reader036/viewer/2022071408/6100094a45b0f84e0e1d54f9/html5/thumbnails/30.jpg)
Technische Universität München
www.th4.ch.tum.de
Past members: Matteo Maestri (→ U Milan, I) Hakim Meskine (→ Wiley, D) Michael Rieger (→ BASF, D) Jutta Rogal (→ RU Bochum, D)
Collaborations: Daan Frenkel (Cambridge, UK) Edvin Lundgren (U Lund, SE) Horia Metiu (UCSB, USA) Matthias Scheffler (FHI Berlin, D)
Thanks so much!!!
Present members: Max Hoffmann, Sebastian Matera, Jörg Meyer