Rmu_pbq _srmk_rga iglcrga kmbcj eclcp_rgml dmp fcrcpmeclcmsq a_r_jwqgq md pclcu_`jc dccbqrmaiq
G_l A_gjjg_s
Qsncptgqmpq8 Npmd, bp, gp, Hmpgq Rfw`_sr* Npmd, bp, gp, Ictgl T_l Ecck Amslqcjjmpq8 @pgegrrc Bctmafr* Djmpclac Tcpkcgpc
K_qrcp%q bgqqcpr_rgml qs`kgrrcb gl mpbcp rm m`r_gl rfc _a_bckga bcepcc md K_qrcp md Qagclac gl Afckga_j Cleglccpgle
Bcn_prkclr Md K_rcpg_jq* Rcvrgjcq ?lb Afckga_j Cleglccpgle Af_gp8 Npmd, bp, N_sj Igciclq D_asjrw md Cleglccpgle _lb ?pafgrcarspc ?a_bckga wc_p 0./5+0./6
Rmu_pbq _srmk_rga iglcrga kmbcj eclcp_rgml dmp fcrcpmeclcmsq a_r_jwqgq md pclcu_`jc dccbqrmaiq
G_l A_gjjg_s
Qsncptgqmpq8 Npmd, bp, gp, Hmpgq Rfw`_sr* Npmd, bp, gp, Ictgl T_l Ecck Amslqcjjmpq8 @pgegrrc Bctmafr* Djmpclac Tcpkcgpc
K_qrcp%q bgqqcpr_rgml qs`kgrrcb gl mpbcp rm m`r_gl rfc _a_bckga bcepcc md K_qrcp md Qagclac gl Afckga_j Cleglccpgle
Bcn_prkclr Md K_rcpg_jq* Rcvrgjcq ?lb Afckga_j Cleglccpgle Af_gp8 Npmd, bp, N_sj Igciclq D_asjrw md Cleglccpgle _lb ?pafgrcarspc ?a_bckga wc_p 0./5+0./6
Acknowledgments
Habemus thesis! We have a thesis. As with any scientific work, this work is the product of more than one person and I’d like to acknowledge the people who helped make this possible.
First of all, I would like to thank Prof. dr. ir. Guy B. Marin for providing me with the opportunity to complete my thesis and education at LCT. I’d also like to thank both of my promotors: prof. dr. ir. Joris Thybaut & Prof. dr. ir. Kevin M. Van Geem. Both have provided me with a perspective and way of approaching problems which I will carry with me for the rest of my professional life.
I have also been blessed to have two counselors, ir. Brigitte Devocht and ir. Florence Vermeire, who have gone above and beyond to make this work the best work it could have been. I’d like to thank you both for always providing me with help and guidance when I needed it. For being patient with the many silly questions and the somewhat muddled way in which I often express myself. When you approach a new subject, it is easy to take the forest for the trees and you have consistently helped me see the bigger picture. Whenever I got lost on one of my (many) tangents, you helped point my gaze towards (more) productive avenues. I have benefited immensely from your input and feedback.
I’d also like to thank the people at LCT. I’d like to thank the Genesys group for helping me during the hackathons. I’d like to thank ir. Cato Pappijn in particular for her time and effort in the last few months.
Finally, I’d like to thank my family and friends for supporting me throughout my academic career. It’s been a thrilling ride.
FACULTY OF ENGINEERING AND ARCHITECTURE
Laboratory for Chemical Technology Director: Prof. Dr. Ir. Guy B. Marin
Laboratory for Chemical Technology
Declaration concerning the accessibility of the master thesis
Undersigned,
Ian Cailliau
Graduated from Ghent University, academic year 2017-2018 and is author of the master thesis with title:
Towards automatic kinetic model generation for heterogeneous catalysis of renewable feedstocks
The author(s) gives (give) permission to make this master dissertation available for consultation and to copy parts of this master dissertation for personal use. In the case of any other use, the copyright terms have to be respected, in particular with regard to the obligation to state expressly the source when quoting results from this master dissertation.
Laboratory for Chemical Technology • Technologiepark 914, B-9052 Gent • www.lct.ugent.be Secretariat : T +32 (0)9 33 11 756 • F +32 (0)9 33 11 759 • [email protected]
Towards automatic kinetic model generation for heterogeneous catalysis of renewable feedstocks
Ian Cailliau
Master's dissertation submitted in order to obtain the academic degree of Master of Science in Chemical Engineering
Academic year 2017-2018
Promotors: Prof. dr. ir. Joris Thybaut & Prof. dr. ir. Kevin Van Geem Coaches: ir. Brigitte Devocht & ir. Florence Vermeire
GHENT UNIVERSITY Faculty of Engineering and Architecture
Department of Chemical Engineering and Technical Chemistry Laboratory for Chemical Technology Chairman: Prof. dr. ir. Guy B. Marin
Abstract
The automatic network generation tool Genesys has been extended towards catalytic reactions. As detailed kinetic models can typically contain hundreds or thousands of reactions, the manual construction of these models becomes infeasible. Genesys, originally developed for gas-phase species and reactions, is expanded to allow for the correct representation of metal surface species. Additionally, the UBI-QEP method has been implemented to allow for the on-the-fly calculation of heats of adsorption and activation energies. As a proof of concept, a case study concerning the hydrodeoxygenation of propionic acid over a nickel catalyst is discussed.
Keywords: Automatic Network Generation, Metal Catalysis, UBI-QEP, Genesys Towards automatic kinetic model generation for heterogeneous catalysis of renewable feedstocks
Ian Cailliau
Promotor(s): Prof. dr. ir. Joris Thybaut and Prof. dr. ir. Kevin Van Geem Coaches: ir. Brigitte Devocht and ir. Florence Vermeire
Abstract: The automatic network generation tool Genesys has 3, 6, 7]. For catalytic systems RMG-Cat already has a basic been extended towards catalytic reactions. As detailed kinetic implementation which makes use of a Group Additivity scheme models can typically contain hundreds or thousands of reactions, for thermodynamic properties and Brönsted-Evans-Polanyi the manual construction of these models becomes infeasible. relationships for the determination of activation energies. In Genesys, originally developed for gas-phase species and reactions, this work, the Unity Bond Index Quadratic Exponential is expanded to allow for the correct representation of metal Potential (UBI-QEP) method, first reported by Shustorovich et surface species. Additionally, the UBI-QEP method has been implemented to allow for the on-the-fly calculation of heats of al. [8] is implemented. adsorption and activation energies. As a proof of concept, a case study concerning the hydrodeoxygenation of propionic acid over II. SPECIES REPRESENTATION a nickel catalyst is discussed. Keywords: Automatic Network Generation, Metal Catalysis, A. Graph representation of molecules UBI-QEP, Genesys Genesys makes use of a cheminformatics library, the Chemistry Development Kit (CDK) [3]. In this software I. INTRODUCTION package, a graph structure is used for the unambiguous With over 90% of all chemical processes worldwide making representation and manipulation of chemical species. The use of catalysts, there is need for a more thorough graph (i.e. a molecule) is composed of various edges (i.e. understanding of the catalytic processes [1]. Microkinetic atoms) and vertices (i.e. bonds). The edges and vertices are models are suited for this purpose as they help identify which weighted, i.e. they are given more information about lone pairs, surface species and surface reaction paths control the observed single electrons, element type, bond type, etc. product spectrum [1]. Detailed kinetic models can contain hundreds or thousands of reactions making manual Besides the graph-based structure, for input and output construction of these models tedious and error-prone [2]. In the purposes, molecule identifiers such as InChI and SMILES, are past this problem has been tackled for combustion and used. InChI’s are unique for a given species, making them an pyrolysis by use of automatic network generators. As of excellent choice for database searches. writing, most software packages deal only with gas-phase reactions and are unfit for heterogenous catalysis. B. Introducing catalytically active sites
Within this work, the goal is to implement metal catalyst Metal surface species are implemented in Genesys in a functionalities in an in-house developed automatic network similar way to gas-phase species, to allow for a smooth conversion of molecule identifier to graph structures and vice generator called Genesys [2, 3]. For this purpose, the focus will versa. Additionally, manipulations of the graph structure be with two of the main aspects of automatic network generation. The first is the correct representation of metal should result in the correct products being generated. surface species. The second is the calculation of In a first instance, use is made of a dummy atom to represent thermodynamic properties for surface species and kinetic the catalytically active site. In this method, the functionality of parameters for surface reactions. The representation of surface species has been done in literature by the use of a dummy atom an unused atom (e.g. helium) is changed to act as the catalyst to represent the active sites. ReNGeP has already correctly surface. An example of how this representation looks like for formic acid is given in represented surface species by use of this dummy atom [4, 5]. Figure 1. RMG-Cat [6] has been extended in a similar fashion and also makes use of a dummy atom. In a second stage, the decision is made to move beyond the conventional dummy atom towards a dummy catalyst block Thermodynamic properties and kinetic parameters have to be consisting of multiple, connected, dummy atoms, this determined for all species and reactions. Because of the nature of automatic network generation this has to be performed on- representation is also shown in the-fly, while the network is being generated. If available from Figure 1. extensive databases, the thermodynamic or kinetic data, A catalyst block was chosen for several reasons. On a obtained via experiment, regression or theoretical calculations, can be used. If not, fast on-the-fly calculations methods are technical level it simplifies the process of generating the correct used which are less computationally expensive, yet provide graph structure for a given species in the software. It eases the implementation of multidentate adsorption, dissociate reasonably accurate values for many species and reactions [2, adsorption, etc. (pseudo-)atom) and the gas-phase bond dissociation energies (DAB, where both A and B are (pseudo-)atoms) for a small number of covalent bonds.
The result is an incredible reduction in the number of parameters which have to be determined, whilst still providing the user with reasonably accurate predictions (15-20 kJ/mole) [8, 14, 15]. UBI-QEP has successfully been applied to several case-studies [16-18]. Because UBI-QEP can make reliable predictions under conditions of data scarcity, it is the method preferred in this work and implemented in Genesys. It is however confined to heterogeneous metal catalysis.
B. Methodology The UBI-QEP analytical equations are based on a set of mathematical assumptions [8]. The first relies on the
introduction of a new quantity, the bond index. The bond index is defined by its relationship to the distance between two atoms Figure 1 Molecule identifiers to a graph structure representation with as shown in equation (1). Here x is the bond index, r is the dummy atom and dummy block. distance between 2 atoms, r0 is the equilibrium distance and a is a scaling factor [8]. Another important reason for the usage of a dummy bock is (1) the implementation of UBI-QEP in Genesys. UBI-QEP is a fast calculation method for thermodynamics and kinetics. The UBI- It is assumed that the bond index is conserved for a given QEP method as well as its implementation will be the topic of system. Thus, every bond within a system considered has a section III. corresponding bond index (xi) and the sum of those bond indices will be equal to one in case of equilibrium (equation Along with representing the species, extra information on the (2)). Equations (1) and (2) form the Unity Bond Index (UBI) catalyst is stored. This currently includes the choice of the part of UBI-QEP [8]. dummy atom, the length of the dummy block and some ∑ 1 (2) parameters required to calculate thermodynamic and kinetic properties. In the future the information on the catalyst can be For the second assumption, the potential energy contribution expended in a flexible way. The decision was made to store this (Q) of a given bond is given as a function of the corresponding additional information independently of the representation. bond index (equation (3)). This function is the Quadratic Contrary to atoms and bonds which have their information (i.e. Exponential Potential part in UBI-QEP [8]. atomic number, symbol, atomic weight, abundance, …) all stored in the graph structure. 2 (3)
III. UBI-QEP Within the UBI-QEP framework, multiple notations are used for the heat of adsorption with subtle but important differences. A. The need for fast calculation of thermodynamic properties Within the UBI-QEP framework the coordination number, n, is and kinetic parameters the number of ligands a single atom has with the surface. QA During automatic kinetic model generation, often hundreds of denotes the heat of adsorption as observed experimentally and species and thousands of reactions are formed. All of those makes no claim as to the coordination number. Q0A is used to need to be assigned thermodynamic or kinetic data in a fast and denote the heat of adsorption for A given 1 single ligand on-the-fly manner. Some methods for surface species and (coordination number 1). QnA is used to denote the atomic heat reactions have been described in literature. The Benson Group of adsorption for A given a coordination number n. Both Q0A Additivity scheme has been extended to surface species [9, 10]. and QnA are therefore explicitly parameters as determined Additionally, Brönsted-Evans-Polanyi (BEP) relationships within the UBI-QEP scheme. This as opposed to QA which is have been successful for the determination of activation simply the observed heat of adsorption. energies in heterogeneous catalysis [11, 12]. Both are implemented in RMG-Cat [6] and for gas-phase species and C. Derivation of the analytical equations reactions in Genesys, but no extension to surface species has UBI-QEP uses equations (1)-(3) to derive analytical yet been made in the latter. These methods are able to give equations for different cases. The total energy of adsorption is accurate predictions, but the downside is the lack of databases calculated by considering the total energy of a system, E(n), for the required groups or parameters. which is obtained by summing up the energetic contributions for the relevant bonds as shown in Figure 2. If the active site is The Unity Bond Index Quadratic Exponential Potential (UBI- involved in the bond considered, the QEP contribution uses the QEP) [8, 13] method is able to determine molecular heats of atomic heat of adsorption. If it is a bond within the adsorbed adsorption, as well as activation energies for surface reactions, species, the QEP contribution uses the bond dissociation energy by using simple analytical equations and only few parameters. (BDE) of the gas-phase molecule. It should be noted that it is The independent variables that must be determined by the user are the atomic heats of adsorption (QA, with A the contact assumed that there is no interaction between atom B and the The first step (step 1 in Figure 3) is the identification of the surface. reaction, thereby classifying it into one of the implemented reactions (step 1.1): dissociative adsorption, non-dissociative adsorption, surface disproportionation and surface dissociation.
Depending on the type of reaction a different UBI-QEP equation will be used for the determination of the heat of adsorption and the activation energy. The species involved will then be mapped on to the corresponding equation (step 1.2). This mapping is required to assure that the proper thermodynamic properties (i.e. heats of adsorption for surface species and heats of formation for gas-phase analogues) are determined. Additionally, these values must be plugged in Figure 2 The total energy as the sum of the QEP contributions of the correctly into the final equation. individual bonds. As can be seen in the figure, n=3 for this specific case. In step 2.1 the heats of adsorption, denoted QA … QAB for surface species A…AB are calculated. These are also The total energy in Figure 2 is subsequently minimized under determined using the UBI-QEP methodology, the exact process the UBI constraint (equation (2)). The result is an analytical is described in the next subsection. For the gas-phase analogues equation shown in equation (4). Here Q is the heat of AB,n of A…AB the heat of formation, denoted H … H also have adsorption for a molecule AB and a given coordination number A AB to be determined. This uses the gas-phase thermodynamics n. In Figure 2, n is equal to three. D is the bond dissociation AB methods already available within Genesys. Finally, in a third energy for the AB bond. step, the activation energy is calculated.