

Project Head: Sebastian Matera
Research Staff: Sandra Döpking
Funding: This research is carried out in the framework of Matheon supported by Einstein Foundation Berlin
Internal cooperation: F. Schmidt (ZIB), Ch. Schütte (FU, ZIB), R. Schneider (TU)
External cooperation: Ch. Scheurer, K. Reuter (TUM), M. J. Hoffmann (Stanford), Aditya (Ashi) Savara (Oak Ridge National Laboratory)
On the route to a more efficient exploitation of energy and materials, the design of new heterogenous catalysts is a central aspect, be it for the production of fine chemicals or the conversion of solar energy to fuels. Ideally, such a design is based on an atomistic understanding of the origin of catalytic activity. In order to enable this detailed undestanding, firstprinciples kinetic Monte Carlo approaches have been established during the last decade. Despite their success, these still have some some limitations. On the one hand, the electronic structure methods, employed to determine required rate parameters, introduce a nonnegligible error into the later. On the other hand, the need for stochastic simulation and the typically large dimension of the parameter space hampers the determination of the rate determining steps by sensitivity analysis, i.e. the most interesting input for a rational catalyst design.
The purpose of this project is to address the aforementioned problems by development of tools for local and global sensitivity analyses. These are applied to models from classical catalysis, but also to models for photocatalytic processes.
Finite difference approximations for gradient based local sensitivity analysis suffer from the statistical noise in the stochastic simulation results, but also from the typically large number of parameters. In order to accelerate gradient estimation, we have developed a three stage strategy, which significantly reduces the number of necessary finite difference approximations. In the first stage, bounds for the gradients are sampled from a single trajectory. In the second stage, those gradients with nonvanishing bounds are estimated using a direct sampling approach based on linear response theory. Only those sensitivities, which can not be estimated with sufficient accuracy by one of the previous approaches, are numerically estimated. The above image shows the socalled Degree of Rate Control for the CO oxidation on ruthenium dioxide. The model has 22 parameters, of which only the nonvanishing are shown. Those sensitivities, which needed to be estimated by finite differences, are marked with symbols.
With significant errors in the rate parameters, local sensitivity analysis can only be the first step. For a more complete picture and to properly quantify the error propagation, global sensitivity analyses are needed. These generally require the numerical evalutaion of integrals over the parameter space. For this, we have implemented dimension and locally adaptive sparse grids based on locally supported basis functions. We have modified the original approach such that multiprocessor enviroments can be better exploited for the simulation of the data. The above image shows the plane spanned by the first two parameters of the 22parameter model for the CO oxidation on ruthenium dioxide. Left are the grid points in the plane and right the reaction rate.
Besides established models from classical catalysis for testing purposes, we investigate the photocatalytic CO_{2} reduction on titania anatase. Together with colleagues from the Zuse Institute Berlin (CSE6), we also consider this system to elucidate the effect of particle scale field heterogeneities on the intrinsic and particle averaged activity.
Döpking, S., Plaisance,C.P., Strobusch, D., Reuter, K. Scheurer, Ch. & Matera,S. Addressing global uncertainty and sensitivity in firstprinciples based microkinetic models by an adaptive sparse grid approach (submitted)
Döpking, S., & Matera, S. (2017). Error propagation in firstprinciples kinetic Monte Carlo simulation. Chemical Physics Letters, 674, 2832.
Hoffmann, M. J., Engelmann, F., & Matera, S. (2017). A practical approach to the sensitivity analysis for kinetic Monte Carlo simulation of heterogeneous catalysis. The Journal of Chemical Physics, 146(4), 044118
Gelß, P., Matera, S., Schütte, Ch. (2016). Solving the master equation without kinetic Monte Carlo: Tensor train approximations for a CO oxidation model. Journal of Computational Physics 314(1), 489502
Matera, S., Maestri, M., Cuoci, A., & Reuter, K. (2014). Predictivequality surface reaction chemistry in real reactor models: Integrating firstprinciples kinetic Monte Carlo simulations into computational fluid dynamics. ACS Catalysis, 4(11), 40814092.
Döpking, S. (2014). Sensitivity analysis of the error propagation of first principles kinetic methods with a sparse grid approach. Master’s thesis, TU Munich
Engelmann, F. (2016). Linear Response and Sensitivity Analysis for kinetic Monte Carlo, Project report, TU Munich