Direct and Adjoint Sensitivity Analysis of Chemical Kinetic Systems With

Total Page:16

File Type:pdf, Size:1020Kb

Direct and Adjoint Sensitivity Analysis of Chemical Kinetic Systems With ARTICLE IN PRESS Atmospheric Environment 37 (2003) 5083–5096 Direct and adjoint sensitivity analysis ofchemical kinetic systems with KPP: Part I—theory and software tools Adrian Sandua, Dacian N. Daescub, Gregory R. Carmichaelc,* a Department of Computer Science, Virginia Polytechnic Institute and State University, 660 McBryde Hall, Blacksburg, VA 24061, USA b Institute for Mathematics and its Applications, University of Minnesota, 400 Lind Hall, 207 Church Street S.E., Minneapolis, MN 55455, USA c Center for Global and Regional Environmental Research, University of Iowa, 204 IATL, Iowa City, IA 52242-1297, USA Received 26 January 2003; accepted 7 August 2003 Abstract The analysis ofcomprehensive chemical reactions mechanisms, parameter estimation techniques, and variational chemical data assimilation applications require the development of efficient sensitivity methods for chemical kinetics systems. The new release (KPP-1.2) of the kinetic preprocessor (KPP) contains software tools that facilitate direct and adjoint sensitivity analysis. The direct-decoupled method, built using BDF formulas, has been the method of choice for direct sensitivity studies. In this work, we extend the direct-decoupled approach to Rosenbrock stiff integration methods. The need for Jacobian derivatives prevented Rosenbrock methods to be used extensively in direct sensitivity calculations; however, the new automatic and symbolic differentiation technologies make the computation of these derivatives feasible. The direct-decoupled method is known to be efficient for computing the sensitivities of a large number ofoutput parameters with respect to a small number ofinput parameters. The adjoint modeling is presented as an efficient tool to evaluate the sensitivity of a scalar response function with respect to the initial conditions and model parameters. In addition, sensitivity with respect to time-dependent model parameters may be obtained through a single backward integration ofthe adjoint model. KPP softwaremay be used to completely generate the continuous and discrete adjoint models taking fulladvantage ofthe sparsity ofthe chemical mechanism. Flexible direct-decoupled and adjoint sensitivity code implementations are achieved with minimal user intervention. In a companion paper, we present an extensive set of numerical experiments that validate the KPP software tools for several direct/adjoint sensitivity applications, and demonstrate the efficiency of KPP-generated sensitivity code implementations. r 2003 Elsevier Ltd. All rights reserved. Keywords: Chemical kinetics; Sensitivity analysis; Direct-decoupled method; Adjoint model 1. Introduction nonlinear differential equations dy 0 0 0 F The mathematical formulation of chemical reaction ¼ f ðt; y; pÞ; yðt Þ¼y ; t ptpt : ð1Þ dt mechanisms is given by a coupled system ofstiff The solution yðtÞARn represents the time evolution of the concentrations ofthe species considered in the *Corresponding author. Tel.: +1-319-335-3333; fax: +1- chemical mechanism starting from the initial configura- 0 319-335-3337. tion y : Throughout this work vectors will be repre- T E-mail addresses: [email protected] (A. Sandu), daescu@ sented in column format and an upper script ðÁÞ will ima.umn.edu (D.N. Daescu), [email protected] denote the transposition operator. The rate ofchange in (G.R. Carmichael). the concentrations y is determined by the nonlinear 1352-2310/$ - see front matter r 2003 Elsevier Ltd. All rights reserved. doi:10.1016/j.atmosenv.2003.08.019 ARTICLE IN PRESS 5084 A. Sandu et al. / Atmospheric Environment 37 (2003) 5083–5096 T production/loss function f ¼ðf1; y; fnÞ ; which depends frequent modifications of the existing ones, there is a on a vector ofparameters pARm: In practice, the vector need for software tools that facilitate the sensitivity p may represent reaction rate coefficients, the initial state analysis ofgeneral chemical kinetic mechanisms. The ofthe model ðp ¼ y0Þ; additional source/sink terms (e.g. kinetic preprocessor (KPP) (Damian-Iordache et al., emission rates), etc. We assume that problem (1) has a 1995) has been successfully used in the forward unique solution y ¼ yðt; pÞ once the model parameters integration ofthe chemical kinetics systems ( Sandu are specified. Comprehensive atmospheric reaction et al., 1997a, b; Verwer et al., 1999). The new release mechanisms take into consideration as many as 100 (KPP-1.2) presented in this paper implements a com- chemical species involved in hundreds ofchemical prehensive set of software tools for direct and adjoint reactions (see e.g. Carter, 2000), such that to efficiently sensitivity analysis. Given a chemical mechanism de- integrate system (1) fast and reliable numerical methods scribed by a list ofchemical reactions, KPP generates a must be implemented (Byrne and Dean, 1993; Jacobson flexible code for the model, its forward integration and and Turco, 1994; Sandu et al., 1997a, b). In addition, the the direct-decoupled/adjoint sensitivity analysis. model parameters are often obtained from experimental The paper is organized as follows: a review of the data and their accuracy is hard to estimate. The direct-decoupled sensitivity analysis and extensions to development and validation ofchemical reactions Runge–Kutta and Rosenbrock stiff integration methods mechanisms require a systematic sensitivity analysis to are presented in Section 2. In Section 3, we review the evaluate the effects of parameter variations on the model continuous and discrete adjoint sensitivity methods and solution. discuss practical issues ofthe adjoint code implementa- n The sensitivities ScðtÞAR are defined as the deriva- tion for chemical kinetics systems. The KPP tools that tives ofthe solution with respect to the parameters facilitate the implementation of the sensitivity methods @yðtÞ are presented in Section 4. Tutorial examples for ScðtÞ¼ ; 1pcpm: ð2Þ building the direct-decoupled code and the adjoint code @pc with KPP are presented in Sections 5 and 6, respectively. A large sensitivity value ScðtÞ shows that the para- In Section 7, we outline the new numerical methods for meter pc plays an essential role in determining the model sensitivity calculations available in the KPP numerical forecast yðtÞ; therefore one problem of interest is to library. A summary ofthe results and concluding nÂm evaluate the sensitivities SðtÞAR for t0ptptF: Some remarks are presented in Section 8. practical applications (e.g. data assimilation) require the sensitivity ofa scalar response function Z tF 2. Direct sensitivity analysis g ¼ g#ðt; yðt; pÞÞ dt ð3Þ 0 t For the direct analysis we consider the parameters p to with respect to the model parameters be constant, i.e. they do not change in time. By Z F @g t @g# differentiating (1) with respect to the parameters one ¼ STðtÞ ðt; yðt; pÞÞ dt: ð4Þ obtains the sensitivity equations (variational equations) @p t0 @y dSc Depending on the problem at hand, an appropriate ¼ Jðt; y; pÞSc þ fpc ðt; y; pÞ; dt method for sensitivity evaluation must be selected. The 0 0 @y most popular and efficient techniques for sensitivity Scðt Þ¼ ; 1pcpm; ð5Þ pc studies are the direct-decoupled and the adjoint sensi- @ tivity methods. These approaches are complementary; where J is the Jacobian matrix ofthe derivative function for example, a single direct sensitivity calculation could and fpc are its partial derivatives with respect to the 0 0 F provide @½y1ðtÞ; y; ynðtÞ=@yjðt Þ for all t ptpt ; parameters while a single adjoint calculation provides @y ðtFÞ= j @½f1; y; fn 0 F @½y1ðtÞ; y; ynðtÞ for all t ptpt : The sensitivity infor- Jðt; y; pÞ¼ ; @½y1; y; yn mation provided by the direct method may be used to y @½f1; ; fn c asses how parametric uncertainties propagate in the fpc ðt; y; pÞ¼ ; 1p pm: ð6Þ @pc system; the adjoint method is more suitable for inverse modeling and may be used to identify the origin of Ifthe parameter pc represents the initial concentration c 0 0 uncertainty in a given model output. In this paper, ofthe th species, pc ¼ yc; then fpc ¼ 0 and Scðt Þ¼ special emphasis is given to the adjoint technique which ð0; y; 0; 1; 0; y; 0ÞT; the cth vector ofthe canonical n 0 has not been used as extensively as the direct method in basis in R ; otherwise, Scðt Þ¼0: the context ofchemical systems. The variational equations (5) are linear. The direct Given the multitude ofapplications, the continuous method solves simultaneously the model equation (1) development ofnew reaction mechanisms and the together with the variational equations (5) to obtain ARTICLE IN PRESS A. Sandu et al. / Atmospheric Environment 37 (2003) 5083–5096 5085 both concentrations and sensitivities. The combined nonlinear system ofequations which implicitly defines system (1)–(5) has the Jacobian yiþ1; and which must be solved by a Newton-type þ y iterative scheme. Typically, the solution yi 1 at @½f ; JS1 þ fp1 ; ; JSm þ fpm fqþ1g y iteration q þ 1 is computed as 0@½y; S1; ; Sm 1 J 0 ? 0 ½I À hbJðti; yi; pÞðyiþ1 À yiþ1Þ B C fqþ1g fqg B ðJS Þ þ J J ? 0 C XkÀ1 B 1 y p1 C ¼ B C; ð7Þ ¼ a yiÀj þ hbf ðtiþ1; yiþ1; pÞÀyiþ1; ð10Þ @ ^ & ^ A j fqg fqg j¼0 ðJS Þ þ J 0 ? J m y pm which requires one factorization of I À hbJ and one where the subscripts in the component sub-matrices backsubstitution per iteration. denote partial differentiation. The eigenvalues of the Discretization of the continuous sensitivity equation:In combined Jacobian (7) are the eigenvalues of J (the this approach (Dunker, 1984; Caracotsios and Stewart, Jacobian of the model equations), with different multi- 1985; Leis and Kramer, 1986) one discretizes the plicities; therefore, if model (1) is stiff the sensitivity continuous sensitivity equation (5) with the same BDF equations (5) are also stiff. To maintain stability, an method used for discretizing model (9) implicit time-stepping method is needed. Implicit meth- XkÀ1 ods solve linear systems with the ‘‘prediction’’ matrix iþ1 iÀj iþ1 iþ1 iþ1 Sc ¼ ajSc þ hbJðt ; y ; pÞSc ðI À hgJÞ; where h is the stepsize and g a parameter j¼0 determined by the method.
Recommended publications
  • Wannier90 As a Community Code: New
    IOP Journal of Physics: Condensed Matter Journal of Physics: Condensed Matter J. Phys.: Condens. Matter J. Phys.: Condens. Matter 32 (2020) 165902 (25pp) https://doi.org/10.1088/1361-648X/ab51ff 32 Wannier90 as a community code: new 2020 features and applications © 2020 IOP Publishing Ltd Giovanni Pizzi1,29,30 , Valerio Vitale2,3,29 , Ryotaro Arita4,5 , Stefan Bl gel6 , Frank Freimuth6 , Guillaume G ranton6 , JCOMEL ü é Marco Gibertini1,7 , Dominik Gresch8 , Charles Johnson9 , Takashi Koretsune10,11 , Julen Ibañez-Azpiroz12 , Hyungjun Lee13,14 , 165902 Jae-Mo Lihm15 , Daniel Marchand16 , Antimo Marrazzo1 , Yuriy Mokrousov6,17 , Jamal I Mustafa18 , Yoshiro Nohara19 , 4 20 21 G Pizzi et al Yusuke Nomura , Lorenzo Paulatto , Samuel Poncé , Thomas Ponweiser22 , Junfeng Qiao23 , Florian Thöle24 , Stepan S Tsirkin12,25 , Małgorzata Wierzbowska26 , Nicola Marzari1,29 , Wannier90 as a community code: new features and applications David Vanderbilt27,29 , Ivo Souza12,28,29 , Arash A Mostofi3,29 and Jonathan R Yates21,29 Printed in the UK 1 Theory and Simulation of Materials (THEOS) and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 CM Lausanne, Switzerland 2 Cavendish Laboratory, Department of Physics, University of Cambridge, Cambridge CB3 0HE, United Kingdom 10.1088/1361-648X/ab51ff 3 Departments of Materials and Physics, and the Thomas Young Centre for Theory and Simulation of Materials, Imperial College London, London SW7 2AZ, United Kingdom 4 RIKEN
    [Show full text]
  • Accelerating Molecular Docking by Parallelized Heterogeneous Computing – a Case Study of Performance, Quality of Results
    solis vasquez, leonardo ACCELERATINGMOLECULARDOCKINGBY PARALLELIZEDHETEROGENEOUSCOMPUTING–A CASESTUDYOFPERFORMANCE,QUALITYOF RESULTS,ANDENERGY-EFFICIENCYUSINGCPUS, GPUS,ANDFPGAS Printed and/or published with the support of the German Academic Exchange Service (DAAD). ACCELERATINGMOLECULARDOCKINGBYPARALLELIZED HETEROGENEOUSCOMPUTING–ACASESTUDYOF PERFORMANCE,QUALITYOFRESULTS,AND ENERGY-EFFICIENCYUSINGCPUS,GPUS,ANDFPGAS at the Computer Science Department of the Technische Universität Darmstadt submitted in fulfilment of the requirements for the degree of Doctor of Engineering (Dr.-Ing.) Doctoral thesis by Solis Vasquez, Leonardo from Lima, Peru Reviewers Prof. Dr.-Ing. Andreas Koch Prof. Dr. Christian Plessl Date of the oral exam October 14, 2019 Darmstadt, 2019 D 17 Solis Vasquez, Leonardo: Accelerating Molecular Docking by Parallelized Heterogeneous Computing – A Case Study of Performance, Quality of Re- sults, and Energy-Efficiency using CPUs, GPUs, and FPGAs Darmstadt, Technische Universität Darmstadt Date of the oral exam: October 14, 2019 Please cite this work as: URN: urn:nbn:de:tuda-tuprints-92886 URL: https://tuprints.ulb.tu-darmstadt.de/id/eprint/9288 This document is provided by TUprints, The Publication Service of the Technische Universität Darmstadt https://tuprints.ulb.tu-darmstadt.de This work is licensed under a Creative Commons “Attribution-ShareAlike 4.0 In- ternational” license. ERKLÄRUNGENLAUTPROMOTIONSORDNUNG §8 Abs. 1 lit. c PromO Ich versichere hiermit, dass die elektronische Version meiner Disserta- tion mit der schriftlichen Version übereinstimmt. §8 Abs. 1 lit. d PromO Ich versichere hiermit, dass zu einem vorherigen Zeitpunkt noch keine Promotion versucht wurde. In diesem Fall sind nähere Angaben über Zeitpunkt, Hochschule, Dissertationsthema und Ergebnis dieses Versuchs mitzuteilen. §9 Abs. 1 PromO Ich versichere hiermit, dass die vorliegende Dissertation selbstständig und nur unter Verwendung der angegebenen Quellen verfasst wurde.
    [Show full text]
  • Universid Ade De Sã O Pa
    A hardware/software codesign for the chemical reactivity of BRAMS Carlos Alberto Oliveira de Souza Junior Dissertação de Mestrado do Programa de Pós-Graduação em Ciências de Computação e Matemática Computacional (PPG-CCMC) UNIVERSIDADE DE SÃO PAULO DE SÃO UNIVERSIDADE Instituto de Ciências Matemáticas e de Computação Instituto Matemáticas de Ciências SERVIÇO DE PÓS-GRADUAÇÃO DO ICMC-USP Data de Depósito: Assinatura: ______________________ Carlos Alberto Oliveira de Souza Junior A hardware/software codesign for the chemical reactivity of BRAMS Master dissertation submitted to the Instituto de Ciências Matemáticas e de Computação – ICMC- USP, in partial fulfillment of the requirements for the degree of the Master Program in Computer Science and Computational Mathematics. FINAL VERSION Concentration Area: Computer Science and Computational Mathematics Advisor: Prof. Dr. Eduardo Marques USP – São Carlos August 2017 Ficha catalográfica elaborada pela Biblioteca Prof. Achille Bassi e Seção Técnica de Informática, ICMC/USP, com os dados fornecidos pelo(a) autor(a) Souza Junior, Carlos Alberto Oliveira de S684h A hardware/software codesign for the chemical reactivity of BRAMS / Carlos Alberto Oliveira de Souza Junior; orientador Eduardo Marques. -- São Carlos, 2017. 109 p. Dissertação (Mestrado - Programa de Pós-Graduação em Ciências de Computação e Matemática Computacional) -- Instituto de Ciências Matemáticas e de Computação, Universidade de São Paulo, 2017. 1. Hardware. 2. FPGA. 3. OpenCL. 4. Codesign. 5. Heterogeneous-computing. I. Marques, Eduardo, orient. II. Título. Carlos Alberto Oliveira de Souza Junior Um coprojeto de hardware/software para a reatividade química do BRAMS Dissertação apresentada ao Instituto de Ciências Matemáticas e de Computação – ICMC-USP, como parte dos requisitos para obtenção do título de Mestre em Ciências – Ciências de Computação e Matemática Computacional.
    [Show full text]
  • Development of Chemoinformatic Tools to Enumerate Functional Groups in Molecules for Organic Aerosol Characterization G
    Manuscript prepared for Atmos. Chem. Phys. with version 2015/04/24 7.83 Copernicus papers of the LATEX class copernicus.cls. Date: 1 March 2016 Technical Note: Development of chemoinformatic tools to enumerate functional groups in molecules for organic aerosol characterization G. Ruggeri1 and S. Takahama1 1ENAC/IIE Swiss Federal Institute of Technology Lausanne (EPFL), Lausanne, Switzerland Correspondence to: Satoshi Takahama (satoshi.takahama@epfl.ch) Abstract. Functional groups (FGs) can be used as a reduced representation of organic aerosol composition in both ambient and environmental controlled chamber studies, as they retain a cer- tain chemical specificity. Furthermore, FG composition has been informative for source apportion- ment, and various models based on a group contribution framework have been developed to calcu- 5 late physicochemical properties of organic compounds. In this work, we provide a set of validated chemoinformatic patterns that correspond to: 1) a complete set of functional groups that can entirely describe the molecules comprised in the α-pinene and 1,3,5-trimethylbenzene MCMv3.2 oxidation schemes, 2) FGs that are measurable by Fourier transform infrared spectroscopy (FTIR), 3) groups incorporated in the SIMPOL.1 vapor pressure estimation model, and 4) bonds necessary for the cal- 10 culation of carbon oxidation state. We also provide example applications for this set of patterns. We compare available aerosol composition reported by chemical speciation measurements and FTIR for different emission sources, and calculate the FG contribution to the O:C ratio of simulated gas phase composition generated from α-pinene photooxidation (using MCMv3.2 oxidation scheme). 1 Introduction 15 Atmospheric aerosols are complex mixtures of inorganic salts, mineral dust, sea salt, black carbon, metals, organic compounds, and water (Seinfeld and Pandis, 2006).
    [Show full text]
  • Open Chemoinformatic Resources to Explore the Structure, Properties and Chemical Space of Molecules
    RSC Advances REVIEW View Article Online View Journal | View Issue Open chemoinformatic resources to explore the structure, properties and chemical space of Cite this: RSC Adv.,2017,7,54153 molecules a ab a Mariana Gonzalez-Medina,´ J. Jesus´ Naveja, Norberto Sanchez-Cruz´ a and Jose´ L. Medina-Franco * New technologies are shaping the way drug discovery data is analyzed and shared. Open data initiatives and web servers are assisting the analysis of the large amounts of data that we are now able to produce. The final goal is to accelerate the process of moving from new data to useful information that could lead to Received 27th October 2017 treatments for human diseases. This review discusses open chemoinformatic resources to analyze the Accepted 21st November 2017 diversity and coverage of the chemical space of screening libraries and to explore structure–activity DOI: 10.1039/c7ra11831g relationships of screening data sets. Free resources to implement workflows and representative web- rsc.li/rsc-advances based applications are emphasized. Future directions in this field are also discussed. Creative Commons Attribution 3.0 Unported Licence. 1. Introduction connections between biological activities, ligands and proteins.3 During the past few years, there has been an important increase Herein we review representative chemoinformatic tools in open data initiatives to promote the availability of free essential to explore the structure, chemical space and properties research-based tools and information.1 While there is still some of molecules. The review is focused on recent and representative resistance to open data in some chemistry and drug discovery free web-based applications.
    [Show full text]
  • Development of an Atmospheric Chemistry Model Coupled to the PALM Model System 6.0: Implementation and First Applications
    https://doi.org/10.5194/gmd-2020-286 Preprint. Discussion started: 25 September 2020 c Author(s) 2020. CC BY 4.0 License. Development of an atmospheric chemistry model coupled to the PALM model system 6.0: Implementation and first applications Basit Khan1, Sabine Banzhaf2, Edward C. Chan2,3, Renate Forkel1, Farah Kanani-Sühring4,7, Klaus Ketelsen5, Mona Kurppa6, Björn Maronga4,10, Matthias Mauder1, Siegfried Raasch4, Emmanuele Russo2,8,9, Martijn Schaap2, and Matthias Sühring4 1Institute of Meteorology and Climate Research, Atmospheric Environmental Research (IMK-IFU), Karlsruhe Institute of Technology, Garmisch-Partenkirchen, 82467, Germany 2Freie Universität Berlin (FUB), Institute of Meteorology, TrUmF 3Institute for Advanced Sustainability Studies (IASS), Potsdam 4Leibniz University Hannover (LUH), Institute of Meteorology and Climatology 5Independent Software Consultant 6University of Helsinki 7Harz Energie GmbH & Co. KG, Goslar, Germany 8Climate and Environmental Physics, Physics Institute, University of Bern, Sidlerstrasse 5, 3012, Bern, Switzerland 9Oeschger Centre for Climate Change Research, University of Bern, Hochschulstrasse 4, 3012, Bern, Switzerland 10University of Bergen, Geophysical Institute, Bergen, Norway Correspondence: Basit Khan ([email protected]) Abstract. In this article we describe the implementation of an online-coupled gas-phase chemistry model in the turbulence resolving PALM model system 6.0. The new chemistry model is part of the PALM-4U components (read: PALM for you; PALM for urban applications) which are designed for application of PALM model in the urban environment (Maronga et al., 2020). The 5 latest version of the Kinetic PreProcessor (KPP, 2.2.3), has been utilised for the numerical integration of gas-phase chemical reactions. A number of tropospheric gas-phase chemistry mechanisms of different complexity have been implemented ranging from the photostationary state to more complex mechanisms such as CBM4, which includes major pollutants namely O3, NO, NO2, CO, a simplified VOC chemistry and a small number of products.
    [Show full text]
  • A Consistent Cheminformatics Framework for Automated Virtual Screening
    A consistent cheminformatics framework for automated virtual screening Cumulative Dissertation to receive the degree Dr. rer. nat. at the Faculty of Mathematics, Informatics and Natural Sciences University of Hamburg submitted to the Department of Informatics of the University of Hamburg Sascha Urbaczek born in Mannheim Hamburg, August 2014 1. Reviewer: Prof. Dr. Matthias Rarey 2. Reviewer: JProf. Dr. Tobias Schwabe 3. Reviewer: Prof. Dr. Andreas Hildebrandt Date of thesis defense: 23.01.2015 F¨urmeine Lena Acknowledgements At this point, I would like to acknowledge the people who supported and encouraged me during my doctoral thesis. First and foremost, I extend my sincere gratitude to my supervisor Prof. Matthias Rarey for entrusting me with this interesting research project and providing guidance and support every step of the way. I also thank the Beiersdorf AG for funding the project and Dr. Inken Groth and Prof. Stefan Heuser in particular for their constant support as well as their great interest in my work. Furthermore, I thank the BiosolveIT GmbH for the productive working relationship during and especially after my time at the Center for Bioinformatics. As for the current and former colleagues from the Center of Bioinformatics (which are far too many to list here); I am grateful to have had the opportunity to work in such a remarkably pleasant and productive environment. I really appreciate everyone's hard work and dedication during these (very successful) last years. In particular I owe a great deal of gratitude to Dr. Tobias Lippert from whom I learned a lot and who has become a very dear friend to me.
    [Show full text]
  • Development of Chemoinformatic Tools to Enumerate Functional Groups in Molecules for Organic Aerosol Characterization
    Atmos. Chem. Phys., 16, 4401–4422, 2016 www.atmos-chem-phys.net/16/4401/2016/ doi:10.5194/acp-16-4401-2016 © Author(s) 2016. CC Attribution 3.0 License. Technical Note: Development of chemoinformatic tools to enumerate functional groups in molecules for organic aerosol characterization Giulia Ruggeri and Satoshi Takahama ENAC/IIE Swiss Federal Institute of Technology Lausanne (EPFL), Lausanne, Switzerland Correspondence to: Satoshi Takahama (satoshi.takahama@epfl.ch) Received: 1 October 2015 – Published in Atmos. Chem. Phys. Discuss.: 27 November 2015 Revised: 4 March 2016 – Accepted: 9 March 2016 – Published: 11 April 2016 Abstract. Functional groups (FGs) can be used as a re- been many proposals for reducing representations in which a duced representation of organic aerosol composition in both mixture of 10 000C different types of molecules (Hamilton ambient and controlled chamber studies, as they retain a et al., 2004) are represented by some combination of their certain chemical specificity. Furthermore, FG composition molecular size, carbon number, polarity, or elemental ratios has been informative for source apportionment, and vari- (Pankow and Barsanti, 2009; Kroll et al., 2011; Daumit et al., ous models based on a group contribution framework have 2013; Donahue et al., 2012), many of which are associated been developed to calculate physicochemical properties of with observable quantities (e.g., by aerosol mass spectrom- organic compounds. In this work, we provide a set of val- etry (AMS; Jayne et al., 2000), gas chromatography–mass idated chemoinformatic patterns that correspond to (1) a spectrometry (GC-MS and GCxGC-MS; Rogge et al., 1993; complete set of functional groups that can entirely de- Hamilton et al., 2004)).
    [Show full text]
  • Development of an Atmospheric Chemistry Model Coupled to the PALM Model System 6.0: Implementation and first Applications Basit Khan1, Sabine Banzhaf2, Edward C
    Development of an atmospheric chemistry model coupled to the PALM model system 6.0: Implementation and first applications Basit Khan1, Sabine Banzhaf2, Edward C. Chan2,3, Renate Forkel1, Farah Kanani-Sühring4,7, Klaus Ketelsen5, Mona Kurppa6, Björn Maronga4,10, Matthias Mauder1, Siegfried Raasch4, Emmanuele Russo2,8,9, Martijn Schaap2, and Matthias Sühring4 1Institute of Meteorology and Climate Research, Atmospheric Environmental Research (IMK-IFU), Karlsruhe Institute of Technology, Garmisch-Partenkirchen, 82467, Germany 2Freie Universität Berlin (FUB), Institute of Meteorology, TrUmF, Germany 3Institute for Advanced Sustainability Studies (IASS), Potsdam, Germany 4Leibniz University Hannover (LUH), Institute of Meteorology and Climatology, Germany 5Independent Software Consultant, Hannover, Germany 6University of Helsinki, Finland 7Harz Energie GmbH & Co. KG, Goslar, Germany 8Climate and Environmental Physics, Physics Institute, University of Bern, Sidlerstrasse 5, 3012, Bern, Switzerland 9Oeschger Centre for Climate Change Research, University of Bern, Hochschulstrasse 4, 3012, Bern, Switzerland 10University of Bergen, Geophysical Institute, Bergen, Norway Correspondence: Basit Khan ([email protected]) Abstract. In this article we describe the implementation of an online-coupled gas-phase chemistry model in the turbulence resolving PALM model system 6.0 (formerly an abbreviation for Parallelized Large-eddy Simulation Model and now an independent name). The new chemistry model is implemented in the PALM model as part of the PALM-4U (PALM for urban applications) 5 components, which are designed for application of PALM model in the urban environment (Maronga et al., 2020). The latest version of the Kinetic PreProcessor (KPP, 2.2.3), has been utilised for the numerical integration of gas-phase chemical reactions. A number of tropospheric gas-phase chemistry mechanisms of different complexity have been implemented ranging from the photostationary state (PHSTAT) to mechanisms with a strongly simplified VOC chemistry (e.g.
    [Show full text]
  • Organic Aerosol from the Oxidation of Biogenic Organic Compounds: a Modelling Study
    Universiteit Antwerpen Doctoral Thesis Organic aerosol from the oxidation of biogenic organic compounds: a modelling study Promoters: Author: Dr. Jean-Fran¸cois Muller¨ , Karl Ceulemans Prof. dr. Magda Claeys A thesis submitted in fulfilment of the requirements for the degree of Doctor in de wetenschappen: chemie in the Faculty of Science, Department of Chemistry 13 October 2014 Cover design: Anita Muys (Nieuwe Mediadienst, Universiteit Antwerpen) Cover photograph: \Tall Trees in Schw¨abisch Gm¨und",Matthias Wassermann (http: //www.mawpix.com/) ISBN 9789057284649 D/2014/12.293/26 Summary Aerosols exert a strong influence on climate and air quality. The oxidation of biogenic volatile organic compounds leads to the formation of secondary organic aerosol (SOA), which makes up a large, but poorly quantified fraction of atmospheric aerosols. In this work, SOA formation from α- and β-pinene, two important biogenic species, is invest- igated through modelling. A detailed mechanism describing the complex gas phase chemistry of oxidation products is generated, based on theoretical results and structure activity relationships. SOA form- ation through the partitioning of oxidation products between the gas and aerosol phases is modelled with the help of estimated vapour pressures and activity coefficients. Given the large model uncertainties, validation against smog chamber experiments is essential. Extensive comparisons with available laboratory measurements show that the model is capable of predicting most SOA yields to within a factor of two for both α- and β-pinene. However, for dark ozonolysis experiments, the model largely overestimates the temperat- ure dependence of the yields and largely underestimates their values above 30 ◦C, by up to a factor of ten.
    [Show full text]
  • Knowledge Acquisition from Chemical Computation-An Initiative Transfer
    Knowledge Acquisition from Chemical Computation: An Initiative to Transfer Tacit Knowledge to Explicit Form Mr. Sunil Khilari 1 , Dr. Sachin Kadam2 12IMED Research Center, Bharti Vidyapeeth University, Pune Abstract - Knowledge derived from chemical computation is the initiative to transfer tacit knowledge to explicit form. Here in this paper researcher is identified computational problems from chemical domain to better understand existing chemical reaction approaches, actual experimental results and its corresponding efficient and effective solutions. During this process researcher contribute to this knowledge base by identifying the areas and locations where their knowledge Management is possible. Further this paper insight on development of chemical ontology and KM framework based on proposed ontology .Also this research paper focuses on software tools used to implement the novel solutions of computational problems in chemical reaction mechanism. Keyword - Knowledge Management, Tacit Knowledge, Explicit Knowledge, Chemical Ontology I. INTRODUCTION The science of computational chemistry carries many benefits to society. Chemistry helps to explain how the nature works, principally by exploring nature on the chemical molecule. Industry uses chemical reactions to make useful product. Chemical analysis is used to identify chemical properties of substances. The science of computational chemistry is unique in that it comprises the capability to generate new forms of compound by combining the elements and existing substances in new but controlled ways. The construction of new chemical compounds is the area from which knowledge generated and need to capture and preserve. Here, we explore the different types of chemical reactions, highlighting the reactions that have solution and implementation phases. So we consider computational chemistry to comprise the study of the structure, properties, and dynamics of chemical systems.
    [Show full text]
  • Revisiting the Proposition of Binding Pockets and Bioactive Poses for GSK-3Β Allosteric Modulators Addressed to Neurodegenerative Diseases
    International Journal of Molecular Sciences Article Revisiting the Proposition of Binding Pockets and Bioactive Poses for GSK-3β Allosteric Modulators Addressed to Neurodegenerative Diseases Guilherme M. Silva 1,* , Rosivaldo S. Borges 2,3,*, Kelton L. B. Santos 2,3 , Leonardo B. Federico 4 , Isaque A. G. Francischini 4 , Suzane Q. Gomes 1, Mariana P. Barcelos 4, Rai C. Silva 1, Cleydson B. R. Santos 2,3,4 and Carlos H. T. P. Silva 1,2,4 1 Departamento de Química, Faculdade de Filosofia, Ciências e Letras de Ribeirão Preto, Universidade de São Paulo, Ribeirão Preto 14040-901, SP, Brazil; [email protected] (S.Q.G.); [email protected] (R.C.S.); [email protected] (C.H.T.P.S.) 2 Programa de Pós-Graduação em Química Medicinal e Modelagem Molecular, Instituto de Ciências da Saúde, Universidade Federal do Pará, Belém 66075-110, PA, Brazil; [email protected] (K.L.B.S.); [email protected] (C.B.R.S.) 3 Laboratório de Modelagem em Química Computacional, Departamento de Ciências Biológicas e da Saúde, Universidade Federal do Amapá, Macapá 68902-280, AP, Brazil 4 Laboratório Computacional de Química Farmacêutica, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Universidade de São Paulo, Av. do Café, s/n, Ribeirão Preto 14040-903, SP, Brazil; [email protected] (L.B.F.); [email protected] (I.A.G.F.); [email protected] (M.P.B.) * Correspondence: [email protected] (G.M.S.); [email protected] (R.S.B.) Citation: Silva, G.M.; Borges, R.S.; Santos, K.L.B.; Federico, L.B.; Abstract: Glycogen synthase kinase-3 beta (GSK-3β) is an enzyme pertinently linked to neurode- Francischini, I.A.G.; Gomes, S.Q.; generative diseases since it is associated with the regulation of key neuropathological features in Barcelos, M.P.; Silva, R.C.; Santos, the central nervous system.
    [Show full text]