Kinetic multi-layer model of aerosol surface and bulk chemistry (KM-SUB): the influence of interfacial transport and bulk diffusion on the oxidation of oleic acid by ozone

We present a novel kinetic multi-layer model that explicitly resolves mass transport and chemical reaction at the surface and in the bulk of aerosol particles (KM-SUB). The model is based on the PRA framework of gas–particle interactions (P¨oschl et al., 5 2007), and it includes reversible adsorption, surface reactions and surface-bulk exchange as well as bulk diffusion and reaction. Unlike earlier models, KM-SUB does not require simplifying assumptions about steady-state conditions and radial mixing. The temporal evolution and concentration profiles of volatile and non-volatile species at the gas-particle interface and in the particle bulk can be modeled along with surface 10 concentrations and gas uptake coefficients. In this study we explore and exemplify the effects of bulk diffusion on the rate of reactive gas uptake for a simple reference system, the ozonolysis of oleic acid particles, in comparison to experimental data and earlier model studies. We demonstrate how KM-SUB can be used to interpret and analyze experimental data from laboratory stud15 ies, and how the results can be extrapolated to atmospheric conditions. In particular, we show how interfacial transport and bulk transport, i.e., surface accommodation, bulk accommodation and bulk diffusion, influence the kinetics of the chemical reaction. Sensitivity studies suggest that in fine air particulate matter oleic acid and compounds with similar reactivity against ozone (C=C double bonds) can reach chemical lifetimes of 20 multiple hours only if they are embedded in a (semi-)solid matrix with very low diffusion coefficients (10−10 cm2 s−1). Depending on the complexity of the investigated system, unlimited numbers of volatile and non-volatile species and chemical reactions can be flexibly added and treated with KM-SUB. We propose and intend to pursue the application of KM-SUB 25 as a basis for the development of a detailed master mechanism of aerosol chemistry as well as for the derivation of simplified but realistic parameterizations for large-scale atmospheric and climate models.

A study of the ethene-ozone reaction with photoelectron spectroscopy: measurement of product branching ratios and atmospheric implications

The ozone-ethene reaction has been investigated at low pressure in a flow-tube interfaced to a u.v. photoelectron spectrometer. Photoelectron spectra recorded as a function of reaction time have been used to estimate partial pressures of the reagents and products, using photoionization cross-sections for selected photoelectron bands of the reagents and products, which have been measured separately. Product yields compare favourably with results of other studies, and the production of oxygen and acetaldehyde have been measured as a function of time for the first time. A reaction scheme developed for the ozone-ethene reaction has been used to simulate the reagents and products as a function of time. The results obtained are in good agreement with the experimental measurements. For each of the observed products, the simulations allow the main reaction (or reactions) for production of that product to be established. The product yields have been used in a global model to estimate their global annual emissions in the atmosphere. Of particular interest are the calculated global annual emissions of formaldehyde (0.96 ± 0.10 Tg) and formic acid, (0.05 ± 0.01 Tg) which are estimated as 0.04% and 0.7% of the total annual emission respectively.

Determination of the degradation kinetics of anthocyanins in a model juice system using isothermal and non-isothermal methods

The effect of temperature on the degradation of blackcurrant anthocyanins in a model juice system was determined over a temperature range of 4–140 °C. The thermal degradation of anthocyanins followed pseudo first-order kinetics. From 4–100 °C an isothermal method was used to determine the kinetic parameters. In order to mimic the temperature profile in retort systems, a non-isothermal method was applied to determine the kinetic parameters in the model juice over the temperature range 110–140 °C. The results from both isothermal and non-isothermal methods fit well together, indicating that the non-isothermal procedure is a reliable mathematical method to determine the kinetics of anthocyanin degradation. The reaction rate constant (k) increased from 0.16 (±0.01) × 10−3 to 9.954 (±0.004) h−1 at 4 and 140 °C, respectively. The temperature dependence of the rate of anthocyanin degradation was modelled by an extension of the Arrhenius equation, which showed a linear increase in the activation energy with temperature.

Kinetic multi-layer model of gas-particle interactions in aerosols and clouds (KM-GAP): linking condensation, evaporation and chemical reactions of organics, oxidants and water

We present a novel kinetic multi-layer model for gas-particle interactions in aerosols and clouds (KM-GAP) that treats explicitly all steps of mass transport and chemical reaction of semi-volatile species partitioning between gas phase, particle surface and particle bulk. KM-GAP is based on the PRA model framework (Pöschl-Rudich-Ammann, 2007), and it includes gas phase diffusion, reversible adsorption, surface reactions, bulk diffusion and reaction, as well as condensation, evaporation and heat transfer. The size change of atmospheric particles and the temporal evolution and spatial profile of the concentration of individual chemical species can be modelled along with gas uptake and accommodation coefficients. Depending on the complexity of the investigated system, unlimited numbers of semi-volatile species, chemical reactions, and physical processes can be treated, and the model shall help to bridge gaps in the understanding and quantification of multiphase chemistry and microphysics in atmo- spheric aerosols and clouds. In this study we demonstrate how KM-GAP can be used to analyze, interpret and design experimental investigations of changes in particle size and chemical composition in response to condensation, evaporation, and chemical reaction. For the condensational growth of water droplets, our kinetic model results provide a direct link between laboratory observations and molecular dynamic simulations, confirming that the accommodation coefficient of water at 270 K is close to unity. Literature data on the evaporation of dioctyl phthalate as a function of particle size and time can be reproduced, and the model results suggest that changes in the experimental conditions like aerosol particle concentration and chamber geometry may influence the evaporation kinetics and can be optimized for eðcient probing of specific physical effects and parameters. With regard to oxidative aging of organic aerosol particles, we illustrate how the formation and evaporation of volatile reaction products like nonanal can cause a decrease in the size of oleic acid particles exposed to ozone.

Kinetic multi-layer model of gas-particle interactions in aerosols and clouds (KM-GAP): linking condensation, evaporation and chemical reactions of organics, oxidants and water

We present a novel kinetic multi-layer model for gas-particle interactions in aerosols and clouds (KMGAP) that treats explicitly all steps of mass transport and chemical reaction of semi-volatile species partitioning between gas phase, particle surface and particle bulk. KMGAP is based on the PRA model framework (P¨oschl-Rudich- Ammann, 2007), and it includes gas phase diffusion, reversible adsorption, surface reactions, bulk diffusion and reaction, as well as condensation, evaporation and heat transfer. The size change of atmospheric particles and the temporal evolution and spatial profile of the concentration of individual chemical species can be modeled along with gas uptake and accommodation coefficients. Depending on the complexity of the investigated system and the computational constraints, unlimited numbers of semi-volatile species, chemical reactions, and physical processes can be treated, and the model shall help to bridge gaps in the understanding and quantification of multiphase chemistry and microphysics in atmospheric aerosols and clouds. In this study we demonstrate how KM-GAP can be used to analyze, interpret and design experimental investigations of changes in particle size and chemical composition in response to condensation, evaporation, and chemical reaction. For the condensational growth of water droplets, our kinetic model results provide a direct link between laboratory observations and molecular dynamic simulations, confirming that the accommodation coefficient of water at 270K is close to unity (Winkler et al., 2006). Literature data on the evaporation of dioctyl phthalate as a function of particle size and time can be reproduced, and the model results suggest that changes in the experimental conditions like aerosol particle concentration and chamber geometry may influence the evaporation kinetics and can be optimized for efficient probing of specific physical effects and parameters. With regard to oxidative aging of organic aerosol particles, we illustrate how the formation and evaporation of volatile reaction products like nonanal can cause a decrease in the size of oleic acid particles exposed to ozone.

A kinetic model for the formation of acrylamide during the finish-frying of commercial French fries

Acrylamide is formed from reducing sugars and asparagine during the preparation of French fries. The commercial preparation of French fries is a multi-stage process involving the preparation of frozen, par-fried potato strips for distribution to catering outlets where they are finish fried. The initial blanching, treatment in glucose solution and par-frying steps are crucial since they determine the levels of precursors present at the beginning of the finish frying process. In order to minimize the quantities of acrylamide in cooked fries, it is important to understand the impact of each stage on the formation of acrylamide. Acrylamide, amino acids, sugars, moisture, fat and color were monitored at time intervals during the frying of potato strips which had been dipped in varying concentrations of glucose and fructose during a typical pretreatment. A mathematical model of the finish-frying was developed based on the fundamental chemical reaction pathways, incorporating moisture and temperature gradients in the fries. This showed the contribution of both glucose and fructose to the generation of acrylamide, and accurately predicted the acrylamide content of the final fries.

Kinetic model for the formation of acrylamide during the finish-frying of commercial French fries

Acrylamide is formed from reducing sugars and asparagine during the preparation of French fries. The commercial preparation of French fries is a multistage process involving the preparation of frozen, par-fried potato strips for distribution to catering outlets, where they are finish-fried. The initial blanching, treatment in glucose solution, and par-frying steps are crucial because they determine the levels of precursors present at the beginning of the finish-frying process. To minimize the quantities of acrylamide in cooked fries, it is important to understand the impact of each stage on the formation of acrylamide. Acrylamide, amino acids, sugars, moisture, fat, and color were monitored at time intervals during the frying of potato strips that had been dipped in various concentrations of glucose and fructose during a typical pretreatment. A mathematical model based on the fundamental chemical reaction pathways of the finish-frying was developed, incorporating moisture and temperature gradients in the fries. This showed the contribution of both glucose and fructose to the generation of acrylamide and accurately predicted the acrylamide content of the final fries.

Functional model for catecholase-like activity: A mechanistic approach with manganese(III) complexes of salen type Schiff base ligands

Three new Mn(III) complexes [MnL1(OOCH)(OH2)] (1), [MnL2(OH2)(2)][Mn2L22(NO2)(3)] (2) and [Mn2L21(NO2)(2)] (3) (where H2L1 = H(2)Me(2)Salen = 2,7-bis(2-hydroxyphenyl)-2,6-diazaocta-2,6-diene and H2L2 = H(2)Salpn = 1,7-bis(2-hydroxyphenyl)-2,6-diazahepta-1,6-diene) have been synthesized. X-ray crystal structure analysis reveals that 1 is a mononuclear species whereas 2 contains a mononuclear cationic and a dinuclear nitrite bridged (mu-1 kappa O:2 kappa O') anionic unit. Complex 3 is a phenoxido bridged dimer containing terminally coordinated nitrite. Complexes 1-3 show excellent catecholase-like activity with 3,5-di-tert-butylcatechol (3,5-DTBC) as the substrate. Kinetic measurements suggest that the rate of catechol oxidation follows saturation kinetics with respect to the substrate and first order kinetics with respect to the catalyst. Formation of bis(mu-oxo)dimanganese(III,III) as an intermediate during the course of reaction is identified from ESI-MS spectra. The characteristic six line EPR spectra of complex 2 in the presence of 3,5-DTBC supports the formation of manganese(II)-semiquinonate as an intermediate species during the catalytic oxidation of 3,5-DTBC.

House price dynamics and their reaction to macroeconomic changes

This article applies a three-regime Markov switching model to investigate the impact of the macroeconomy on the dynamics of the residential real estate market in the US. Focusing on the period between 1960 and 2011, the methodology implemented allows for a clearer understanding of the drivers of the real estate market in “boom”, “steady-state” and “crash” regimes. Our results show that the sensitivity of the real estate market to economic changes is regime-dependent. The paper then proceeds to examine whether policymakers are able to influence a regime switch away from the crash regime. We find that a decrease in interest rate spreads could be an effective catalyst to precipitate such a change of state.

Metabifurcation analysis of a mean field model of the cortex

Mean field models (MFMs) of cortical tissue incorporate salient, average features of neural masses in order to model activity at the population level, thereby linking microscopic physiology to macroscopic observations, e.g., with the electroencephalogram (EEG). One of the common aspects of MFM descriptions is the presence of a high-dimensional parameter space capturing neurobiological attributes deemed relevant to the brain dynamics of interest. We study the physiological parameter space of a MFM of electrocortical activity and discover robust correlations between physiological attributes of the model cortex and its dynamical features. These correlations are revealed by the study of bifurcation plots, which show that the model responses to changes in inhibition belong to two archetypal categories or “families”. After investigating and characterizing them in depth, we discuss their essential differences in terms of four important aspects: power responses with respect to the modeled action of anesthetics, reaction to exogenous stimuli such as thalamic input, and distributions of model parameters and oscillatory repertoires when inhibition is enhanced. Furthermore, while the complexity of sustained periodic orbits differs significantly between families, we are able to show how metamorphoses between the families can be brought about by exogenous stimuli. We here unveil links between measurable physiological attributes of the brain and dynamical patterns that are not accessible by linear methods. They instead emerge when the nonlinear structure of parameter space is partitioned according to bifurcation responses. We call this general method “metabifurcation analysis”. The partitioning cannot be achieved by the investigation of only a small number of parameter sets and is instead the result of an automated bifurcation analysis of a representative sample of 73,454 physiologically admissible parameter sets. Our approach generalizes straightforwardly and is well suited to probing the dynamics of other models with large and complex parameter spaces.

Chapter 4: Simulation of reaction injection moulding

Reaction Injection Moulding is a technology that enables the rapid production of complex plastic parts directly from a mixture of two reactive materials of low viscosity. The reactants are mixed in specific quantities and injected into a mould. This process allows large complex parts to be produced without the need for high clamping pressures. This chapter explores the simulation of the complex processes involved in reaction injection moulding. The reaction processes mean that the dynamics of the material in the mould are in constant evolution and an effective model which takes full account of these changing dynamics is introduced and incorporated in to finite element procedures, which are able to provide a complete simulation of the cycle of mould filling and subsequent curing.

Ozone climatology using interactive chemistry: results from the Canadian Middle Atmosphere Model

The climatology of ozone produced by the Canadian Middle Atmosphere Model (CMAM) is presented. This three-dimensional global model incorporates the radiative feedbacks of ozone and water vapor calculated on-line with a photochemical module. This module includes a comprehensive gas-phase reaction set and a limited set of heterogeneous reactions to account for processes occurring on background sulphate aerosols. While transport is global, photochemistry is solved from about 400 hPa to the top of the model at ∼95 km. This approach provides a complete and comprehensive representation of transport, emission, and photochemistry of various constituents from the surface to the mesopause region. A comparison of model results with observations indicates that the ozone distribution and variability are in agreement with observations throughout most of the model domain. Column ozone annual variation is represented to within 5–10% of the observations except in the Southern Hemisphere for springtime high latitudes. The vertical ozone distribution is generally well represented by the model up to the mesopause region. Nevertheless, in the upper stratosphere, the model generally underestimates the amount of ozone as well as the latitudinal tilting of ozone isopleths at high latitude. Ozone variability is analyzed and compared with measurements. The comparison shows that the phase and amplitude of the seasonal variation as well as shorter timescale variations are well represented by the model at various latitudes and heights. Finally, the impact of incorporating ozone radiative feedback on the model climatology is isolated. It is found that the incorporation of ozone radiative feedback results in a cooling of ∼8 K in the summer stratopause region, which corrects a warm bias that results when climatological ozone is used.

A semi-distributed ntegrated itrogen model for multiple source assessment in tchments (INCA): Part I — model structure and process equations

A new model has been developed for assessing multiple sources of nitrogen in catchments. The model (INCA) is process based and uses reaction kinetic equations to simulate the principal mechanisms operating. The model allows for plant uptake, surface and sub-surface pathways and can simulate up to six land uses simultaneously. The model can be applied to catchment as a semi-distributed simulation and has an inbuilt multi-reach structure for river systems. Sources of nitrogen can be from atmospheric deposition, from the terrestrial environment (e.g. agriculture, leakage from forest systems etc.), from urban areas or from direct discharges via sewage or intensive farm units. The model is a daily simulation model and can provide information in the form of time series at key sites, or as profiles down river systems or as statistical distributions. The process model is described and in a companion paper the model is applied to the River Tywi catchment in South Wales and the Great Ouse in Bedfordshire.

Travelling waves in a model of species migration

A model of species migration is presented which takes the form of a reaction-diffusion system. We consider special limits of this model in which we demonstrate the existence of travelling wave solutions. These solutions can be used to describe the migration of cells, bacteria, and some organisms. © 2000 Elsevier Science Ltd. All rights reserved.

Applying robust control theory to solve problems in bio-medical sciences: study of an apoptotic model

Biological models of an apoptotic process are studied using models describing a system of differential equations derived from reaction kinetics information. The mathematical model is re-formulated in a state-space robust control theory framework where parametric and dynamic uncertainty can be modelled to account for variations naturally occurring in biological processes. We propose to handle the nonlinearities using neural networks.