Activities report - Basic and Applied Sciences Branch, Division of Water Supply and Pollution Control

Report covers fiscal year

Stochastic simulation of chemical reactions in spatially complex media

Discrete stochastic simulations, via techniques such as the Stochastic Simulation Algorithm (SSA) are a powerful tool for understanding the dynamics of chemical kinetics when there are low numbers of certain molecular species. However, an important constraint is the assumption of well-mixedness and homogeneity. In this paper, we show how to use Monte Carlo simulations to estimate an anomalous diffusion parameter that encapsulates the crowdedness of the spatial environment. We then use this parameter to replace the rate constants of bimolecular reactions by a time-dependent power law to produce an SSA valid in cases where anomalous diffusion occurs or the system is not well-mixed (ASSA). Simulations then show that ASSA can successfully predict the temporal dynamics of chemical kinetics in a spatially constrained environment.

Effect of chemical conjugation of L-leucine to chitosan on the dispersibility and drug release of a dry powder inhaler nanoparticle formulation

Purpose: This study investigated the effect of chemical conjugation of the amino acid L-leucine to the polysaccharide chitosan on the dispersibility and drug release pattern of a polymeric nanoparticle (NP)-based controlled release dry powder inhaler (DPI) formulation. Methods: A chemical conjugate of L-leucine with chitosan was synthesized and characterized by Infrared (IR) Spectroscopy, Nuclear Magnetic Resonance (NMR) Spectroscopy, Elemental Analysis and X-ray Photoelectron Spectroscopy (XPS). Nanoparticles of both chitosan and its conjugate were prepared by a water-in-oil emulsification – glutaraldehyde cross-linking method using the antihypertensive agent, diltiazem (Dz) hydrochloride as the model drug. The surface morphology and particle size distribution of the nanoparticles were determined by Scanning Electron Microscopy (SEM) and Dynamic Light Scattering (DLS). The dispersibility of the nanoparticle formulation was analysed by a Twin Stage Impinger (TSI) with a Rotahaler as the DPI device. Deposition of the particles in the different stages was determined by gravimetry and the amount of drug released was analysed by UV spectrophotometry. The release profile of the drug was studied in phosphate buffered saline at 37 ⁰C and analyzed by UV spectrophotometry. Results: The TSI study revealed that the fine particle fractions (FPF), as determined gravimetrically, for empty and drug-loaded conjugate nanoparticles were significantly higher than for the corresponding chitosan nanoparticles (24±1.2% and 21±0.7% vs 19±1.2% and 15±1.5% respectively; n=3, p<0.05). The FPF of drug-loaded chitosan and conjugate nanoparticles, in terms of the amount of drug determined spectrophotometrically, had similar values (21±0.7% vs 16±1.6%). After an initial burst, both chitosan and conjugate nanoparticles showed controlled release that lasted about 8 to 10 days, but conjugate nanoparticles showed twice as much total drug release compared to chitosan nanoparticles (~50% vs ~25%). Conjugate nanoparticles also showed significantly higher dug loading and entrapment efficiency than chitosan nanoparticles (conjugate: 20±1% & 46±1%, chitosan: 16±1% & 38±1%, n=3, p<0.05). Conclusion: Although L-leucine conjugation to chitosan increased dispersibility of formulated nanoparticles, the FPF values are still far from optimum. The particles showed a high level of initial burst release (chitosan, 16% and conjugate, 31%) that also will need further optimization.

Stochastic linear multistep methods for the simulation of chemical kinetics

In this paper, we introduce the Stochastic Adams-Bashforth (SAB) and Stochastic Adams-Moulton (SAM) methods as an extension of the tau-leaping framework to past information. Using the theta-trapezoidal tau-leap method of weak order two as a starting procedure, we show that the k-step SAB method with k >= 3 is order three in the mean and correlation, while a predictor-corrector implementation of the SAM method is weak order three in the mean but only order one in the correlation. These convergence results have been derived analytically for linear problems and successfully tested numerically for both linear and non-linear systems. A series of additional examples have been implemented in order to demonstrate the efficacy of this approach.

Chemical probes into the active centre of a heme thiolate monoxygenase

Linalool-8-monoxygenase, a typical bacterial P-450 heme thiolase, shows a high degree of substrate specificity towards linalool. The active site of the pure enzyme has been probed with a large number of substrate analogues with systematic alterations or conformational variations in the linalool molecule. The comparison of three parameters, the mo→mos conversion of the enzyme as a result of substrate binding monitored at 392 nm, theK D of the analogues giving information about energies of association and the relative turnover as substrate have given information about the space-filling characteristics of the substrates in the enzyme cleft, the number of contacts the molecules make with the respective domains of the enzyme and the distance of the site undergoing hydroxylation from the oxygen site, respectively. The data permit the conclusion that linalool makes contact with the enzyme by hydrogen bonding with the hydroxyl group as well through hydrophobic association with all the eight carbons carrying hydrogen in the molecules.

Understanding the chemistry of atmospheric aerosol particle formation via observations of physical properties

It is widely accepted that the global climate is heating up due to human activities, such as burning of fossil fuels. Therefore we find ourselves forced to make decisions on what measures, if any, need to be taken to decrease our warming effect on the planet before any irrevocable damage occurs. Research is being conducted in a variety of fields to better understand all relevant processes governing Earth s climate, and to assess the relative roles of anthropogenic and biogenic emissions into the atmosphere. One of the least well quantified problems is the impact of small aerosol particles (both of anthropogenic and biogenic origin) on climate, through reflecting solar radiation and their ability to act as condensation nuclei for cloud droplets. In this thesis, the compounds driving the biogenic formation of new particles in the atmosphere have been examined through detailed measurements. As directly measuring the composition of these newly formed particles is extremely difficult, the approach was to indirectly study their different characteristics by measuring the hygroscopicity (water uptake) and volatility (evaporation) of particles between 10 and 50 nm. To study the first steps of the formation process in the sub-3 nm range, the nucleation of gaseous precursors to small clusters, the chemical composition of ambient naturally charged ions were measured. The ion measurements were performed with a newly developed mass spectrometer, which was first characterized in the laboratory before being deployed at a boreal forest measurement site. It was also successfully compared to similar, low-resolution instruments. The ambient measurements showed that sulfuric acid clusters dominate the negative ion spectrum during new particle formation events. Sulfuric acid/ammonia clusters were detected in ambient air for the first time in this work. Even though sulfuric acid is believed to be the most important gas phase precursor driving the initial cluster formation, measurements of the hygroscopicity and volatility of growing 10-50 nm particles in Hyytiälä showed an increasing role of organic vapors of a variety of oxidation levels. This work has provided additional insights into the compounds participating both in the initial formation and subsequent growth of atmospheric new aerosol particles. It will hopefully prove an important step in understanding atmospheric gas-to-particle conversion, which, by influencing cloud properties, can have important climate impacts. All available knowledge needs to be constantly updated, summarized, and brought to the attention of our decision-makers. Only by increasing our understanding of all the relevant processes can we build reliable models to predict the long-term effects of decisions made today.

Quantum Chemical Studies on Tropospheric Nucleation Mechanisms Involving Sulfuric Acid

Nucleation is the first step of the process by which gas molecules in the atmosphere condense to form liquid or solid particles. Despite the importance of atmospheric new-particle formation for both climate and health-related issues, little information exists on its precise molecular-level mechanisms. In this thesis, potential nucleation mechanisms involving sulfuric acid together with either water and ammonia or reactive biogenic molecules are studied using quantum chemical methods. Quantum chemistry calculations are based on the numerical solution of Schrödinger's equation for a system of atoms and electrons subject to various sets of approximations, the precise details of which give rise to a large number of model chemistries. A comparison of several different model chemistries indicates that the computational method must be chosen with care if accurate results for sulfuric acid - water - ammonia clusters are desired. Specifically, binding energies are incorrectly predicted by some popular density functionals, and vibrational anharmonicity must be accounted for if quantitatively reliable formation free energies are desired. The calculations reported in this thesis show that a combination of different high-level energy corrections and advanced thermochemical analysis can quantitatively replicate experimental results concerning the hydration of sulfuric acid. The role of ammonia in sulfuric acid - water nucleation was revealed by a series of calculations on molecular clusters of increasing size with respect to all three co-ordinates; sulfuric acid, water and ammonia. As indicated by experimental measurements, ammonia significantly assists the growth of clusters in the sulfuric acid - co-ordinate. The calculations presented in this thesis predict that in atmospheric conditions, this effect becomes important as the number of acid molecules increases from two to three. On the other hand, small molecular clusters are unlikely to contain more than one ammonia molecule per sulfuric acid. This implies that the average NH3:H2SO4 mole ratio of small molecular clusters in atmospheric conditions is likely to be between 1:3 and 1:1. Calculations on charged clusters confirm the experimental result that the HSO4- ion is much more strongly hydrated than neutral sulfuric acid. Preliminary calculations on HSO4- NH3 clusters indicate that ammonia is likely to play at most a minor role in ion-induced nucleation in the sulfuric acid - water system. Calculations of thermodynamic and kinetic parameters for the reaction of stabilized Criegee Intermediates with sulfuric acid demonstrate that quantum chemistry is a powerful tool for investigating chemically complicated nucleation mechanisms. The calculations indicate that if the biogenic Criegee Intermediates have sufficiently long lifetimes in atmospheric conditions, the studied reaction may be an important source of nucleation precursors.

Dynamical processes in the stratosphere and upper troposphere and their influence on the distribution of trace gases in the polar atmosphere

Transport plays an important role in the distribution of long-lived gases such as ozone and water vapour in the atmosphere. Understanding of observed variability in these gases as well as prediction of the future changes depends therefore on our knowledge of the relevant atmospheric dynamics. This dissertation studies certain dynamical processes in the stratosphere and upper troposphere which influence the distribution of ozone and water vapour in the atmosphere. The planetary waves that originate in the troposphere drive the stratospheric circulation. They influence both the meridional transport of substances as well as parameters of the polar vortices. In turn, temperatures inside the polar vortices influence abundance of the Polar Stratospheric Clouds (PSC) and therefore the chemical ozone destruction. Wave forcing of the stratospheric circulation is not uniform during winter. The November-December averaged stratospheric eddy heat flux shows a significant anticorrelation with the January-February averaged eddy heat flux in the midlatitude stratosphere and troposphere. These intraseasonal variations are attributable to the internal stratospheric vacillations. In the period 1979-2002, the wave forcing exhibited a negative trend which was confined to the second half of winter only. In the period 1958-2002, area, strength and longevity of the Arctic polar vortices do not exhibit significant long-term changes while the area with temperatures lower than the threshold temperature for PSC formation shows statistically significant increase. However, the Arctic vortex parameters show significant decadal changes which are mirrored in the ozone variability. Monthly ozone tendencies in the Northern Hemisphere show significant correlations (|r|=0.7) with proxies of the stratospheric circulation. In the Antarctic, the springtime vortex in the lower stratosphere shows statistically significant trends in temperature, longevity and strength (but not in area) in the period 1979-2001. Analysis of the ozone and water vapour vertical distributions in the Arctic UTLS shows that layering below and above the tropopause is often associated with poleward Rossby wave-breaking. These observations together with calculations of cross-tropopause fluxes emphasize the importance of poleward Rossby wave breaking for the stratosphere-troposphere exchange in the Arctic.

Measurements of Volatile Organic Compounds - from Biogenic Emissions to Concentrations in Ambient Air

Volatile organic compounds (VOCs) affect atmospheric chemistry and thereafter also participate in the climate change in many ways. The long-lived greenhouse gases and tropospheric ozone are the most important radiative forcing components warming the climate, while aerosols are the most important cooling component. VOCs can have warming effects on the climate: they participate in tropospheric ozone formation and compete for oxidants with the greenhouse gases thus, for example, lengthening the atmospheric lifetime of methane. Some VOCs, on the other hand, cool the atmosphere by taking part in the formation of aerosol particles. Some VOCs, in addition, have direct health effects, such as carcinogenic benzene. VOCs are emitted into the atmosphere in various processes. Primary emissions of VOC include biogenic emissions from vegetation, biomass burning and human activities. VOCs are also produced in secondary emissions from the reactions of other organic compounds. Globally, forests are the largest source of VOC entering the atmosphere. This thesis focuses on the measurement results of emissions and concentrations of VOCs in one of the largest vegetation zones in the world, the boreal zone. An automated sampling system was designed and built for continuous VOC concentration and emission measurements with a proton transfer reaction - mass spectrometer (PTR-MS). The system measured one hour at a time in three-hourly cycles: 1) ambient volume mixing-ratios of VOCs in the Scots-pine-dominated boreal forest, 2) VOC fluxes above the canopy, and 3) VOC emissions from Scots pine shoots. In addition to the online PTR-MS measurements, we determined the composition and seasonality of the VOC emissions from a Siberian larch with adsorbent samples and GC-MS analysis. The VOC emissions from Siberian larch were reported for the fist time in the literature. The VOC emissions were 90% monoterpenes (mainly sabinene) and the rest sesquiterpenes (mainly a-farnesene). The normalized monoterpene emission potentials were highest in late summer, rising again in late autumn. The normalized sesquiterpene emission potentials were also highest in late summer, but decreased towards the autumn. The emissions of mono- and sesquiterpenes from the deciduous Siberian larch, as well as the emissions of monoterpenes measured from the evergreen Scots pine, were well described by the temperature-dependent algorithm. In the Scots-pine-dominated forest, canopy-scale emissions of monoterpenes and oxygenated VOCs (OVOCs) were of the same magnitude. Methanol and acetone were the most abundant OVOCs emitted from the forest and also in the ambient air. Annually, methanol and mixing ratios were of the order of 1 ppbv. The monoterpene and sum of isoprene 2-methyl-3-buten-2-ol (MBO) volume mixing-ratios were an order of magnitude lower. The majority of the monoterpene and methanol emissions from the Scots-pinedominated forest were explained by emissions from Scots pine shoots. The VOCs were divided into three classes based on the dynamics of the summer-time concentrations: 1) reactive compounds with local biological, anthropogenic or chemical sources (methanol, acetone, butanol and hexanal), 2) compounds whose emissions are only temperaturedependent (monoterpenes), 3) long-lived compounds (benzene, acetaldehyde). Biogenic VOC (methanol, acetone, isoprene MBO and monoterpene) volume mixing-ratios had clear diurnal patterns during summer. The ambient mixing ratios of other VOCs did not show this behaviour. During winter we did not observe systematical diurnal cycles for any of the VOCs. Different sources, removal processes and turbulent mixing explained the dynamics of the measured mixing-ratios qualitatively. However, quantitative understanding will require longterm emission measurements of the OVOCs and the use of comprehensive chemistry models. Keywords: Hydrocarbons, VOC, fluxes, volume mixing-ratio, boreal forest

Doing marriage and parenthood: The division of household labour across two life course events

The current study is a longitudinal investigation into changes in the division of household labour across transitions to marriage and parenthood in the UK. Previous research has noted a more traditional division of household labour, with women performing the majority of housework, amongst spouses and couples with children. However, the bulk of this work has been cross-sectional in nature. The few longitudinal studies that have been carried out have been rather ambiguous about the effect of marriage and parenthood on the division of housework. Theoretically, this study draws on gender construction theory. The key premise of this theory is that gender is something that is performed and created in interaction, and, as a result, something fluid and flexible rather than fixed and stable. The idea that couples ‘do gender’ through housework has been a major theoretical breakthrough. Gender-neutral explanations of the division of household labour, positing rational acting individuals, have failed to explicate why women continue to perform an unequal share of housework, regardless of socio-economic status. Contrastingly, gender construction theory situates gender as the key process in dividing household labour. By performing and avoiding certain housework chores, couples fulfill social norms of what it means to be a man and a woman although, given the emphasis on human agency in producing and contesting gender, couples are able to negotiate alternative gender roles which, in turn, feed back into the structure of social norms in an ever-changing societal landscape. This study adds extra depth to the doing gender approach by testing whether or not couples negotiate specific conjugal and parent roles in terms of the division of household labour. Both transitions hypothesise a more traditional division of household labour. Data comes from the British Household Panel Survey, a large, nationally representative quantitative survey that has been carried out annually since 1991. Here, data tracks the same 776 couples at two separate time points – 1996 and 2005. OLS regression is used to test whether or not transitions to marriage and parenthood have a significant impact on the division of household labour whilst controlling for host of relevant socio-economic factors. Results indicate that marriage has no significant effect on how couples partition housework. Those couples making the transition from cohabitation to marriage do not show significant changes in housework arrangements from those couples who remain cohabiting in both waves. On the other hand, becoming parents does lead to a more traditional division of household labour whilst controlling for socio-economic factors which accompany the move to parenthood. There is then some evidence that couples use the site of household labour to ‘do parenthood’ and generate identities which both use and inform socially prescribed notions of what it means to be a mother and a father. Support for socio-economic explanations of the division of household labour was mixed although it remains clear that they, alone, cannot explain how households divide housework.

Doing marriage and parenthood: The division of household labour across two life course events

The current study is a longitudinal investigation into changes in the division of household labour across transitions to marriage and parenthood in the UK. Previous research has noted a more traditional division of household labour, with women performing the majority of housework, amongst spouses and couples with children. However, the bulk of this work has been cross-sectional in nature. The few longitudinal studies that have been carried out have been rather ambiguous about the effect of marriage and parenthood on the division of housework. Theoretically, this study draws on gender construction theory. The key premise of this theory is that gender is something that is performed and created in interaction, and, as a result, something fluid and flexible rather than fixed and stable. The idea that couples 'do gender' through housework has been a major theoretical breakthrough. Gender-neutral explanations of the division of household labour, positing rational acting individuals, have failed to explicate why women continue to perform an unequal share of housework, regardless of socioeconomic status. Contrastingly, gender construction theory situates gender as the key process in dividing household labour. By performing and avoiding certain housework chores, couples fulfill social norms of what it means to be a man and a woman although, given the emphasis on human agency in producing and contesting gender, couples are able to negotiate alternative gender roles which, in turn, feed back into the structure of social norms in an ever-changing societal landscape. This study adds extra depth to the doing gender approach by testing whether or not couples negotiate specific conjugal and parent roles in terms of the division of household labour. Both transitions hypothesise a more traditional division of household labour. Data comes from the British Household Panel Survey, a large, nationally representative quantitative survey that has been carried out annually since 1991. Here, data tracks the same 776 couples at two separate time points - 1996 and 2005. OLS regression is used to test whether or not transitions to marriage and parenthood have a significant impact on the division of household labour whilst controlling for host of relevant socio-economic factors. Results indicate that marriage has no significant effect on how couples partition housework. Those couples making the transition from cohabitation to marriage do not show significant changes in housework arrangements from those couples who remain cohabiting in both waves. On the other hand, becoming parents does lead to a more traditional division of household labour whilst controlling for socio-economic factors which accompany the move to parenthood. There is then some evidence that couples use the site of household labour to 'do parenthood' and generate identities which both use and inform socially prescribed notions of what it means to be a mother and a father. Support for socio-economic explanations of the division of household labour was mixed although it remains clear that they, alone, cannot explain how households divide housework.

Microscopic theory of chemical-reaction rate

In this paper we deduce the formulae for rate-constant of microreaction with high resolving power of energy from the time-dependent Schrdinger equation for the general case when there is a depression on the reaetional potential surface (when the depression is zero in depth, the case is reduced to that of Eyring). Based on the assumption that Bolzmann distribution is appropriate to the description of reactants, the formula for the constant of macrorate in a form similar to Eyring's is deduced and the expression for the coefficient of transmission is given. When there is no depression on the reactional potential surface and the coefficient of transmission does not seriously depend upon temperature, it is reduced to Eyring's. Thus Eyring's is a special case of the present work.