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 (Poschl-Rudich-Ammann, 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 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 studies, and how the results can be extrapolated to atmospheric conditions. In particular, we show how interfacial 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 (carbon-carbon double bonds) can reach chemical lifetimes of many hours only if they are embedded in a (semi-)solid matrix with very low diffusion coefficients (< 10(-10) cm(2) 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 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.

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.

The adsorption geometry and chemical state of lysine on Cu{110}

Chemisorbed layers of lysine adsorbed on Cu{110} have been studied using X-ray photoelectron spectroscopy (XPS) and near-edge X-ray absorption fine structure (NEXAFS) spectroscopy. XPS indicates that the majority (70%) of the molecules in the saturated layer at room temperature (coverage 0.27 ML) are in their zwitterionic state with no preferential molecular orientation. After annealing to 420 K a less densely packed layer is formed (0.14 ML), which shows a strong angular dependence in the characteristic π-resonance of oxygen K edge NEXAFS and no indication of zwitterions in XPS. These experimental results are best compatible with molecules bound to the substrate through the oxygen atoms of the (deprotonated) carboxylate group and the two amino groups involving Cu atoms in three different close packed rows. This μ4 bonding arrangement with an additional bond through the !-amino group is different from geometries previously suggested for lysine on Cu{110}.

Gas-phase kinetics of chlorosilylene reactions. I. ClSiH + Me3SiH: absolute rate measurements and theoretical calculations for prototype Si−H insertion reactions

Time-resolved studies of chlorosilylene, CISiH, generated by the 193 nm laser flash photolysis of 1-chloro-1-silacyclopent-3-ene, have been carried out to obtain rate constants for its bimolecular reaction with trimethylsilane, Me3SiH, in the gas phase. The reaction was studied at total pressures up to 100 torr (with and without added SF6) over the temperature range 297-407 K. The rate constants were found to be pressure independent and gave the following Arrhenius equation: log(k/cm(3) molecule(-1) s(-1)) = (-13.97 +/- 0.25) + (12.57 +/- 1.64) kJ mol(-1)/RT In 10. The Arrhenius parameters are consistent with a mechanism involving an intermediate complex, whose rearrangement is the rate-determining step. Quantum chemical calculations of the potential energy surface for this reaction and also the reactions of CISiH with SiH4 and the other methylsilanes support this conclusion. Comparisons of both experiment and theory with the analogous Si-H insertion processes of SiH2 and SiMe2 show that the main factor causing the lower reactivity of ClSiH is the secondary energy barrier. The calculations also show the existence of a novel intramolecular H-atom exchange process in the complex of ClSiH with MeSiH3.

The gas-phase reaction between silylene and 2-butyne: kinetics, isotope studies, pressure dependence studies and quantum chemical calculations

Time-resolved kinetic studies of the reactions of silylene, SiH2, and dideutero-silylene, SiD2, generated by laser. ash photolysis of phenylsilane and phenylsilane-d(3), respectively, have been carried out to obtain rate coefficients for their bimolecular reactions with 2-butyne, CH3C CCH3. The reactions were studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas at five temperatures in the range 294-612 K. The second-order rate coefficients, obtained by extrapolation to the high pressure limits at each temperature, fitted the Arrhenius equations where the error limits are single standard deviations: log(k(H)(infinity)/cm(3) molecule(-1) s(-1)) = (-9.67 +/- 0.04) + (1.71 +/- 0.33) kJ mol(-1)/RTln10 log(k(D)(infinity)/cm(3) molecule(-1) s(-1)) = (-9.65 +/- 0.01) + (1.92 +/- 0.13) kJ mol(-1)/RTln10 Additionally, pressure-dependent rate coefficients for the reaction of SiH2 with 2-butyne in the presence of He (1-100 Torr) were obtained at 301, 429 and 613 K. Quantum chemical (ab initio) calculations of the SiC4H8 reaction system at the G3 level support the formation of 2,3-dimethylsilirene [cyclo-SiH2C(CH3)=C(CH3)-] as the sole end product. However, reversible formation of 2,3-dimethylvinylsilylene [CH3CH=C(CH3)SiH] is also an important process. The calculations also indicate the probable involvement of several other intermediates, and possible products. RRKM calculations are in reasonable agreement with the pressure dependences at an enthalpy value for 2,3-dimethylsilirene fairly close to that suggested by the ab initio calculations. The experimental isotope effects deviate significantly from those predicted by RRKM theory. The differences can be explained by an isotopic scrambling mechanism, involving H - D exchange between the hydrogens of the methyl groups and the D-atoms in the ring in 2,3-dimethylsilirene-1,1-d(2). A detailed mechanism involving several intermediate species, which is consistent with the G3 energy surface, is proposed to account for this.

Direct time-resolved study of the kinetics of the reaction of silylene with phenylsilane in the gas phase. Does SiH2 react with the aromatic ring?

Laser flash photolysis studies of silylene, SiH2, generated by the 193 nm laser flash photolysis phenylsilane, PhSiH3, have been carried out to obtain rate constants for its bimolecular reaction with PhSiH3 itself, in the gas phase. The reaction was studied in SF6 (mostly at 10 Torr total pressure) over the temperature range 298-595 K. The rate constants (also found to be pressure independent) gave the following Arrhenius equation: log(k/cm(3) molecule(-1) s(-1)) = (-9.92 +/- 0.04) + (3.31 +/- 0.27) kJ mol(-1)/RT ln 10 Similar investigations of the reaction of silylene with benzene, C6H6, (295-410 K) gave data suggestive of the fact that SiH2 might be reacting with photochemical products of C6H6 as well as with C6H6 itself. However, in the latter system, apparent rate constants were sufficiently low to indicate that in the reaction of SiH2 with PhSiH3 addition to the aromatic ring was unlikely to be in excess of 3% of the total. Quantum chemical calculations of the energy surface for SiH2 + C6H6 indicate that 7-silanorcaradiene and 7-silacycloheptatriene are possible products but that PhSiH3 formation is unlikely. RRKM calculations suggest that 7-silanorcaradiene should be the initial product but that it cannot be collisionally stabilized under experimental conditions

The influence of carbon surface oxygen groups on Dubinin–Astakhov equation parameters calculated from CO2 adsorption isotherm

We present the results of a systematic study of the influence of carbon surface oxidation on Dubinin–Astakhov isotherm parameters obtained from the fitting of CO2 adsorption data. Using GCMC simulations of adsorption on realistic VPC models differing in porosity and containing the most frequently occurring carbon surface functionalities (carboxyls, hydroxyls and carbonyls) and their mixtures, it is concluded that the maximum adsorption calculated from the DA model is not strongly affected by the presence of oxygen groups. Unfortunately, the same cannot be said of the remaining two parameters of this model i.e. the heterogeneity parameter (n) and the characteristic energy of adsorption (E0). Since from the latter the pore diameters of carbons are usually calculated, by inverse-type relationships, it is concluded that they are questionable for carbons containing surface oxides, especially carboxyls.

Kinetics of layer hopping in a diblock copolymer lamellar phase

In the ordered state, symmetric diblock copolymers self-assemble into an anisotropic lamellar morphology. The equilibrium thickness of the lamellae is the result of a delicate balance between enthalpic and entropic energies, which can be tuned by controlling the temperature. Here we devise a simple yet powerful method of detecting tiny changes in the lamellar thickness using optical microscopy. From such measurements we characterize the enthalpic interaction as well as the kinetics of molecules as they hop from one layer to the next in order to adjust the lamellar thickness in response to a temperature jump. The resolution of the measurements facilitate a direct comparison to predictions from self-consistent field theory.

Kinetics of the reaction between OH radicals and monochlorodimethylsulphide (CH3SCH2Cl)

The gas-phase rate coefficient for the reaction between OH radicals and CH3SCH2Cl (MCDMS) was determined to be (2.5±1.3)×10−12 cm3 molecule−1 s−1 using the discharge–flow kinetic technique. An estimate of ≈10−10 cm3 molecule−1 s−1 was obtained for the rate coefficient for reaction of Cl with MCDMS. It would appear that the reaction with OH is not the main loss process for CH3SCH2Cl in the marine boundary layer. The possible implications for the MBL of halogen-promoted oxidation of dimethylsulphide are considered.

A kinetics and mechanistic study of the atmospherically relevant reaction between molecular chlorine and dimethyl sulfide (DMS)

A gas-phase kinetics study of the atmospherically important reaction between Cl2 and dimethyl sulfide (DMS)Cl2 + CH3SCH3 → productshas been made using a flow-tube interfaced to a photoelectron spectrometer. The rate constant for this reaction has been measured at 1.6 and 3.0 torr at T = (294 ± 2) K as (3.4 ± 0.7) × 10−14 cm3 molecule−1 s−1. Reaction (1) has been found to proceed via an intermediate, (CH3)2SCl2, to give CH3SCH2Cl and HCl as the products. The mechanism of this reaction and the structure of the intermediate were investigated using electronic structure calculations. A comparison of the mechanisms of the reactions between Cl atoms and DMS, and Cl2 and DMS has been made and the relevance of the results to atmospheric chemistry is discussed.

The influence of the carbon surface chemical composition on Dubinin-Astakhov equation parameters calculated from SF(6) adsorption data-grand canonical Monte Carlo simulation

Using grand canonical Monte Carlo simulation we show, for the first time, the influence of the carbon porosity and surface oxidation on the parameters of the Dubinin-Astakhov (DA) adsorption isotherm equation. We conclude that upon carbon surface oxidation, the adsorption decreases for all carbons studied. Moreover, the parameters of the DA model depend on the number of surface oxygen groups. That is why in the case of carbons containing surface polar groups, SF(6) adsorption isotherm data cannot be used for characterization of the porosity.

Simple model of adsorption on external surface of carbon nanotubes: a new analytical approach basing on molecular simulation data

Nitrogen adsorption on carbon nanotubes is wide- ly studied because nitrogen adsorption isotherm measurement is a standard method applied for porosity characterization. A further reason is that carbon nanotubes are potential adsorbents for separation of nitrogen from oxygen in air. The study presented here describes the results of GCMC simulations of nitrogen (three site model) adsorption on single and multi walled closed nanotubes. The results obtained are described by a new adsorption isotherm model proposed in this study. The model can be treated as the tube analogue of the GAB isotherm taking into account the lateral adsorbate-adsorbate interactions. We show that the model describes the simulated data satisfactorily. Next this new approach is applied for a description of experimental data measured on different commercially available (and characterized using HRTEM) carbon nanotubes. We show that generally a quite good fit is observed and therefore it is suggested that the observed mechanism of adsorption in the studied materials is mainly determined by adsorption on tubes separated at large distances, so the tubes behave almost independently.

Molecular dynamics simulation insight into the mechanism of phenol adsorption at low coverages from aqueous solutions on microporous carbons

MD simulation studies showing the influence of porosity and carbon surface oxidation on phenol adsorption from aqueous solutions on carbons are reported. Based on a realistic model of activated carbon, three carbon structures with gradually changed microporosity were created. Next, a different number of surface oxygen groups was introduced. The pores with diameters around 0.6 nm are optimal for phenol adsorption and after the introduction of surface oxygen functionalities, adsorption of phenol decreases (in accordance with experimental data) for all studied models. This decrease is caused by a pore blocking effect due to the saturation of surface oxygen groups by highly hydrogen-bounded water molecules.