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.

The kinetics of the gas-phase reactions of selected monoterpenes and cyclo-alkenes with ozone and the NO3 radical

The relative rate method has been used to measure the room-temperature rate constants for the gasphase reactions of ozone and NO3 with selected monoterpenes and cyclo-alkenes with structural similarities to monoterpenes. Measurements were carried out at 298 ! 2 K and 760 ! 10 Torr. The following rate constants (in units of 10"18 cm3 molecule"1 s"1) were obtained for the reaction with ozone: methyl cyclohexene (132 ! 17), terpinolene (1290 ! 360), ethylidene cyclohexane (223 ! 57), norbornene (860 ! 240), t-butyl isopropylidene cyclohexane (1500 ! 460), cyclopentene (543 ! 94), cyclohexene (81 ! 18), cyclooctene (451 ! 66), dicyclopentadiene (1460 ! 170) and a-pinene (107 ! 13). For the reaction with NO3 the rate constants obtained (in units of 10"12 cm3 molecule"1 s"1) were: methyl cyclohexene (7.92 ! 0.95), terpinolene (47.9 ! 4.0), ethylidene cyclohexane (4.30 ! 0.24), norbornene (0.266 ! 0.029), cyclohexene (0.540 ! 0.017), cyclooctene (0.513 ! 0.029), dicyclopentadiene (1.20 ! 0.10) and a-pinene (5.17 ! 0.62). Errors are quoted as the root mean square of the statistical error (95% con!dence) and the quoted error in the rate constant for the reference compound. Combining these results with previous studies, new recommendations for the rate constants are presented. Molecular orbital energies were calculated for each alkene and the kinetic data are discussed in terms of the deviation from the structureeactivity relationship obtained from the rate constants for a series of simple alkenes. Lifetimes with respect to key initiators of atmospheric oxidation have been calculated suggesting that the studied reactions play dominant roles in the night-time removal of these compounds from the atmosphere.

State-resolved gas-surface reactivity of methane in the symmetric C-H stretch vibration on Ni(100)

The state-resolved reactivity of CH4 in its totally symmetric C-H stretch vibration (�1) has been measured on a Ni(100) surface. Methane molecules were accelerated to kinetic energies of 49 and 63:5 kJ=mol in a molecular beam and vibrationally excited to �1 by stimulated Raman pumping before surface impact at normal incidence. The reactivity of the symmetric-stretch excited CH4 is about an order of magnitude higher than that of methane excited to the antisymmetric stretch (�3) reported by Juurlink et al. [Phys. Rev. Lett. 83, 868 (1999)] and is similar to that we have previously observed for the excitation of the first overtone (2�3). The difference between the state-resolved reactivity for �1 and �3 is consistent with predictions of a vibrationally adiabatic model of the methane reaction dynamics and indicates that statistical models cannot correctly describe the chemisorption of CH4 on nickel.

Modelling the warm H(2) infrared emission of the Helix nebula cometary knots

Molecular hydrogen emission is commonly observed in planetary nebulae. Images taken in infrared H(2) emission lines show that at least part of the molecular emission is produced inside the ionized region. In the best studied case, the Helix nebula, the H(2) emission is produced inside cometary knots (CKs), comet-shaped structures believed to be clumps of dense neutral gas embedded within the ionized gas. Most of the H(2) emission of the CKs seems to be produced in a thin layer between the ionized diffuse gas and the neutral material of the knot, in a mini-photodissociation region (mini-PDR). However, PDR models published so far cannot fully explain all the characteristics of the H(2) emission of the CKs. In this work, we use the photoionization code AANGABA to study the H(2) emission of the CKs, particularly that produced in the interface H(+)/H(0) of the knot, where a significant fraction of the H(2) 1-0 S(1) emission seems to be produced. Our results show that the production of molecular hydrogen in such a region may explain several characteristics of the observed emission, particularly the high excitation temperature of the H(2) infrared lines. We find that the temperature derived from H(2) observations, even of a single knot, will depend very strongly on the observed transitions, with much higher temperatures derived from excited levels. We also proposed that the separation between the H alpha and [N II] peak emission observed in the images of CKs may be an effect of the distance of the knot from the star, since for knots farther from the central star the [N II] line is produced closer to the border of the CK than H alpha.

Morphometrical, biochemical and molecular tools for assessing biodiversity. An example in Plebeia remota (Holmberg, 1903) (Apidae, Meliponini)

We see today many efforts to quantify biodiversity in different biomes. It is very important then to develop and to apply other methodologies that allow us to assess biodiversity. Here we present an example of application of three tools with this goal. We analyzed two populations of Plebeia remota from two distinct biomes that already showed several differences in morphology and behavior. Based on these differences, it has been suggested that the populations of Cunha and Prudentopolis do not represent a single species. In order to verify the existence or absence of gene flow between these two groups, we characterized the patterns of mtDNA through RFLP, the patterns of wing venation through geometric morphometry, and the cuticular hydrocarbons through gas chromatography-mass spectrometry. We used bees collected in these two locations and also from colonies which have being kept for around 9 years at Sao Paulo University. We found six different haplotypes in these specimens, of which three of them occurred exclusively in the population of Cunha and three only in the Prudentopolis population. The fact that the populations do not share haplotypes suggests no maternal gene flow between them. The two populations were differentiated by the pattern of the wing veins. They also had different mixtures of cuticle hydrocarbons. Furthermore it was shown that the colonies kept at the university did not hybridize. These two groups may constitute different species. We also show here the importance of using other methodologies than traditional taxonomy to assess and understand biodiversity, especially in bees.

Electronic properties of liquid hydrogen fluoride: A sequential quantum mechanical/Born-Oppenheimer molecular dynamics approach

The electronic properties of liquid hydrogen fluoride (HF) were investigated by carrying out sequential quantum mechanics/Born-Oppenheimer molecular dynamics. The structure of the liquid is in good agreement with recent experimental information. Emphasis was placed on the analysis of polarisation effects, dynamic polarisability and electronic excitations in liquid HF. Our results indicate an increase in liquid phase of the dipole moment (similar to 0.5 D) and isotropic polarisability (5%) relative to their gas-phase values. Our best estimate for the first vertical excitation energy in liquid HF indicates a blue-shift of 0.4 +/- 0.2 eV relative to that of the gas-phase monomer (10.4 eV). (C) 2010 Elsevier B.V. All rights reserved.

Many-body energy decomposition of hydrogen-bonded glycine clusters in gas-phase

A detailed analysis of the many-body contribution to the interaction energies of the gas-phase hydrogen-bonded glycine clusters, (Gly)(N), N = 1-4 is presented. The energetics of the hydrogen-bonded dimer, trimer and tetramer complexes have been analyzed using density-functional theory. The magnitude of the two-through four-body energy terms have been calculated and compared. The relaxation energy and the two-body energy terms are the principal contributors to the total binding energy. Four-body contribution is negligible. However, the three-body contribution is found to be sizable and the formation of the cyclic glycine trimer presents geometric strains that make it less favorable. (C) 2010 Elsevier B.V. All rights reserved.

A further exploration of a nucleophilicity index based on the gas-phase ionization potentials

An empirical nucleophilicity index based on the gas-phase ionization potentials has been recently shown to be useful categorizing and settling the nucleophilicity power of a series of captodative ethylenes reacting in cycloaddition reactions (L.R. Domingo, E. Chamorro, P. Perez, Journal of Organic Chemistry 73 (2008) 4615-4624). In the present work, the applicability of such model is tested within a broader series of substituted alkenes, substituted aromatic compounds and simple nucleophilic molecules. This index obtained within a Koopman`s theorem framework has been evaluated here in both gas and solution phases for several well-known nucleophiles. These results are found to be linearly correlated. Finally, the feasibility of the predictive character of this index has been discussed in comparison to the available experimental nucleophilicities of some amines in water. These results further support and validate the usefulness of such approximation in the modeling of the global nucleophilicity. (C) 2008 Elsevier B.V. All rights reserved.

Engineering of 3D self-directed quantum dot ordering in multilayer InGaAs/GaAs nanostructures by means of flux gas composition

Lateral ordering of InGaAs quantum dots on the GaAs (001) surface has been achieved in earlier reports, resembling an anisotropic pattern. In this work, we present a method of breaking the anisotropy of ordered quantum dots (QDs) by changing the growth environment. We show experimentally that using As(2) molecules instead of As(4) as a background flux is efficient in controlling the diffusion of distant Ga adatoms to make it possible to produce isotropic ordering of InGaAs QDs over GaAs (001). The control of the lateral ordering of QDs under As(2) flux has enabled us to improve their optical properties. Our results are consistent with reported experimental and theoretical data for structure and diffusion on the GaAs surface.

Only subtle protein conformational adaptations are required for ligand binding to thyroid hormone receptors: Simulations using a novel multipoint steered molecular dynamics approach

Thyroid hormone receptors (TR) are hormone-dependent transcription regulators that play a major role in human health, development, and metabolic functions. The thyroid hormone resistance syndrome, diabetes, obesity, and some types of cancer are just a few examples of important diseases that are related to TR malfunctioning, particularly impaired hormone binding. Ligand binding to and dissociation from the receptor ultimately control gene transcription and, thus, detailed knowledge of binding and release mechanisms are fundamental for the comprehension of the receptor`s biological function and development of pharmaceuticals. In this work, we present the first computational study of ligand entry into the ligand binding domain (LBD) of a nuclear receptor. We report molecular dynamics simulations of ligand binding to TRs using a generalization of the steered molecular dynamics technique designed to perform single-molecule pulling simulations along arbitrarily nonlinear driving pathways. We show that only gentle protein movements and conformational adaptations are required for ligand entry into the LBDs and that the magnitude of the forces applied to assist ligand binding are of the order of the forces involved in ligand dissociation. Our simulations suggest an alternative view for the mechanisms ligand binding and dissociation of ligands from nuclear receptors in which ligands can simply diffuse through the protein surface to reach proper positioning within the binding pocket. The proposed picture indicates that the large-amplitude protein motions suggested by the apo- and holo-RXR alpha crystallographic structures are not required, reconciling conformational changes of LBDs required for ligand entry with other nuclear receptors apo-structures that resemble the ligand-bound LBDs.

The effects of temperature on the molecular orientation of zinc phthalocyanine films

Phthalocyanine compounds have been widely investigated as candidate materials for technological applications, which is mainly due to their thermal stability and possibility of processing in the form of thin films. In most applications, the controlled growth of thin films with high crystalline quality is essential. In this study, zinc phthalocyanine (ZnPc) thin films were prepared by evaporation on glass and Au-coated glass substrates with subsequent annealing at different temperatures in ambient atmosphere. The morphological and structural features of 80 nm thick zinc phthalocyanine films were investigated, evidencing an alpha -> beta phase transformation after annealing the films at 200 A degrees C, as indicated by UV-Vis spectroscopy and FTIR analyses. A better uniformity of the annealed films was also evidenced via AFM analysis, which may be of importance for applications where film homogeneity and excellent optical quality are required.

Gas-Phase Nucleophilic and Elimination Reactions in Simple Alkyl Nitrates

There has been increasing interest in the gas-phase reactivity of alkyl nitrates because of their well-known applications as explosives and because of then role in atmospheric and in marine processes This manuscript describes an experimental study by FT-ICR techniques of the gas-phase reactions of OH(-) and F(-) with methyl and ethyl Innate For methyl nitrate, the main reaction channel is found to be an elimination process promoted by abstraction of an a proton from the methyl group. Nucleophilic displacement of nitrate anion through an S(N)2 process at the carbon center Is also found to he an important reaction channel with methyl nitrate In ethyl nitrate, Ruination of NO(3)(-) is greatly enhanced and this is attributed to the ease of an E2-type elimination process promoted by proton abstraction at the beta position of the ethyl group. Theoretical calculations at the MP2/6-311+G(3df,2p)//MP2/6-31+G(d) level of theory ale consistent with the relative importance of the reaction channels and suggest that these reactions proceed through a double well potential The calculations also predict that nucleophilic attack by OH(-) at the nitrogen center (Sn2@N) is energetically the rueful ad pathway but experiments with (18)OH(-) showed no evidence for this channel. Single-point calculations reveal a strong preference for approach to the emboli center and may explain the lack of reactivity at the nitrogen center. Calculations were also carried out or NH(2)(-) and SH(-) to establish the reactivity pattern to provide a better understanding of environmentally relevant nitrate esters.

Flow injection analysis of ethyl xanthate by gas diffusion and UV detection as CS(2) for process monitoring of sulfide ore flotation

A sensitive and robust analytical method for spectrophotometric determination of ethyl xanthate, CH(3)CH(2)OCS(2)(-) at trace concentrations in pulp solutions from froth flotation process is proposed. The analytical method is based on the decomposition of ethyl xanthate. EtX(-), with 2.0 mol L(-1) HCl generating ethanol and carbon disulfide. CS(2). A gas diffusion cell assures that only the volatile compounds diffuse through a PTFE membrane towards an acceptor stream of deionized water, thus avoiding the interferences of non-volatile compounds and suspended particles. The CS(2) is selectively detected by UV absorbance at 206 nm (epsilon = 65,000 L mol(-1) cm(-1)). The measured absorbance is directly proportional to EtX(-) concentration present in the sample solutions. The Beer`s law is obeyed in a 1 x 10(-6) to 2 x 10(-4) mol L(-1) concentration range of ethyl xanthate in the pulp with an excellent correlation coefficient (r = 0.999) and a detection limit of 3.1 x 10(-7) mol L(-1), corresponding to 38 mu g L. At flow rates of 200 mu L min(-1) of the donor stream and 100 mu L min(-1) of the acceptor channel a sampling rate of 15 injections per hour could be achieved with RSD < 2.3% (n = 10, 300 mu L injections of 1 x 10(-5) mol L(-1) EtX(-)). Two practical applications demonstrate the versatility of the FIA method: (i) evaluation the free EtX(-) concentration during a laboratory study of the EtX(-) adsorption capacity on pulverized sulfide ore (pyrite) and (ii) monitoring of EtX(-) at different stages (from starting load to washing effluents) of a flotation pilot plant processing a Cu-Zn sulfide ore. (C) 2010 Elsevier By. All rights reserved.

Gas-Phase Acylium Ion Transfer Reactions Mediated by a Proton Shuttle Mechanism

The mechanism and the energy profile of the gas-phase reaction that mimics esterification under acidic conditions have been investigated at different levels of theory. These reactions are known to proceed with rate constants close to the collision limit in the gas-phase and questions have been raised as to whether the typical addition-elimination mechanism via a tetrahedral intermediate can explain the ease of these processes. Because these reactions are common to many organic and biochemical processes it is important to understand the intrinsic reactivity of these systems. Our calculations at different levels of theory reveal that a stepwise mechanism via a tetrahedral species is characterized by energy barriers that are inconsistent with the experimental results. For the thermoneutral exchange between protonated acetic acid and water and the exothermic reaction of protonated acetic acid and methanol our calculations show that these reactions proceed initially by a proton shuttle between the carbonyl oxygen and the hydroxy oxygen of acetic acid mediated by water, or methanol, followed by displacement at the acylium ion center. These findings suggest that the reactions in the gas-phase should be viewed as an acylium ion transfer reaction. (C) 2010 Wiley Periodicals, Inc. Int J Quantum Chem 111: 1596-1606, 2011