Weighted nucleoli and dually essential coalitions

We consider linearly weighted versions of the least core and the (pre)nucleolus and investigate the reduction possibilities in their computation. We slightly extend some well-known related results and establish their counterparts by using the dual game. Our main results imply, for example, that if the core of the game is not empty, all dually inessential coalitions (which can be weakly minorized by a partition in the dual game) can be ignored when we compute the per-capita least core and the per-capita (pre)nucleolus from the dual game. This could lead to the design of polynomial time algorithms for the per-capita (and other monotone nondecreasingly weighted versions of the) least core and the (pre)nucleolus in specific classes of balanced games with polynomial many dually esential coalitions.

Gas/particle partitioning of carbonyls in the photooxidation of isoprene and 1,3,5-trimethylbenzene

A new denuder-filter sampling technique has been used to investigate the gas/particle partitioning behaviour of the carbonyl products from the photooxidation of isoprene and 1,3,5-trimethylbenzene. A series of experiments was performed in two atmospheric simulation chambers at atmospheric pressure and ambient temperature in the presence of NOx and at a relative humidity of approximately 50%. The denuder and filter were both coated with the derivatizing agent O-(2,3,4,5,6-pentafluorobenzyl)-hydroxylamine (PFBHA) to enable the efficient collection of gas- and particle-phase carbonyls respectively. The tubes and filters were extracted and carbonyls identified as their oxime derivatives by GC-MS. The carbonyl products identified in the experiments accounted for around 5% and 10% of the mass of secondary organic aerosol formed from the photooxidation of isoprene and 1,3,5-trimethylbenzene respectively. Experimental gas/particle partitioning coefficients were determined for a wide range of carbonyl products formed from the photooxidation of isoprene and 1,3,5-trimethylbenzene and compared with the theoretical values based on standard absorptive partitioning theory. Photooxidation products with a single carbonyl moiety were not observed in the particle phase, but dicarbonyls, and in particular, glyoxal and methylglyoxal, exhibited gas/particle partitioning coefficients several orders of magnitude higher than expected theoretically. These findings support the importance of heterogeneous and particle-phase chemical reactions for SOA formation and growth during the atmospheric degradation of anthropogenic and biogenic hydrocarbons.

Mechanistic Models of Anti-HIV Microbicide Drug Delivery

A new modality for preventing HIV transmission is emerging in the form of topical microbicides. Some clinical trials have shown some promising results of these methods of protection while other trials have failed to show efficacy. Due to the relatively novel nature of microbicide drug transport, a rigorous, deterministic analysis of that transport can help improve the design of microbicide vehicles and understand results from clinical trials. This type of analysis can aid microbicide product design by helping understand and organize the determinants of drug transport and the potential efficacies of candidate microbicide products.

Microbicide drug transport is modeled as a diffusion process with convection and reaction effects in appropriate compartments. This is applied here to vaginal gels and rings and a rectal enema, all delivering the microbicide drug Tenofovir. Although the focus here is on Tenofovir, the methods established in this dissertation can readily be adapted to other drugs, given knowledge of their physical and chemical properties, such as the diffusion coefficient, partition coefficient, and reaction kinetics. Other dosage forms such as tablets and fiber meshes can also be modeled using the perspective and methods developed here.

The analyses here include convective details of intravaginal flows by both ambient fluid and spreading gels with different rheological properties and applied volumes. These are input to the overall conservation equations for drug mass transport in different compartments. The results are Tenofovir concentration distributions in time and space for a variety of microbicide products and conditions. The Tenofovir concentrations in the vaginal and rectal mucosal stroma are converted, via a coupled reaction equation, to concentrations of Tenofovir diphosphate, which is the active form of the drug that functions as a reverse transcriptase inhibitor against HIV. Key model outputs are related to concentrations measured in experimental pharmacokinetic (PK) studies, e.g. concentrations in biopsies and blood. A new measure of microbicide prophylactic functionality, the Percent Protected, is calculated. This is the time dependent volume of the entire stroma (and thus fraction of host cells therein) in which Tenofovir diphosphate concentrations equal or exceed a target prophylactic value, e.g. an EC50.

Results show the prophylactic potentials of the studied microbicide vehicles against HIV infections. Key design parameters for each are addressed in application of the models. For a vaginal gel, fast spreading at small volume is more effective than slower spreading at high volume. Vaginal rings are shown to be most effective if inserted and retained as close to the fornix as possible. Because of the long half-life of Tenofovir diphosphate, temporary removal of the vaginal ring (after achieving steady state) for up to 24h does not appreciably diminish Percent Protected. However, full steady state (for the entire stromal volume) is not achieved until several days after ring insertion. Delivery of Tenofovir to the rectal mucosa by an enema is dominated by surface area of coated mucosa and whether the interiors of rectal crypts are filled with the enema fluid. For the enema 100% Percent Protected is achieved much more rapidly than for vaginal products, primarily because of the much thinner epithelial layer of the mucosa. For example, 100% Percent Protected can be achieved with a one minute enema application, and 15 minute wait time.

Results of these models have good agreement with experimental pharmacokinetic data, in animals and clinical trials. They also improve upon traditional, empirical PK modeling, and this is illustrated here. Our deterministic approach can inform design of sampling in clinical trials by indicating time periods during which significant changes in drug concentrations occur in different compartments. More fundamentally, the work here helps delineate the determinants of microbicide drug delivery. This information can be the key to improved, rational design of microbicide products and their dosage regimens.

Pyramid Multi-resolution Scanning for Two-sample Comparison

Testing for two-sample differences is challenging when the differences are local and only involve a small portion of the data. To solve this problem, we apply a multi- resolution scanning framework that performs dependent local tests on subsets of the sample space. We use a nested dyadic partition of the sample space to get a collection of windows and test for sample differences within each window. We put a joint prior on the states of local hypotheses that allows both vertical and horizontal message passing among the partition tree to reflect the spatial dependency features among windows. This information passing framework is critical to detect local sample differences. We use both the loopy belief propagation algorithm and MCMC to get the posterior null probability on each window. These probabilities are then used to report sample differences based on decision procedures. Simulation studies are conducted to illustrate the performance. Multiple testing adjustment and convergence of the algorithms are also discussed.

Distributed Feature Selection in Large n and Large p Regression Problems

Fitting statistical models is computationally challenging when the sample size or the dimension of the dataset is huge. An attractive approach for down-scaling the problem size is to first partition the dataset into subsets and then fit using distributed algorithms. The dataset can be partitioned either horizontally (in the sample space) or vertically (in the feature space), and the challenge arise in defining an algorithm with low communication, theoretical guarantees and excellent practical performance in general settings. For sample space partitioning, I propose a MEdian Selection Subset AGgregation Estimator ({\em message}) algorithm for solving these issues. The algorithm applies feature selection in parallel for each subset using regularized regression or Bayesian variable selection method, calculates the `median' feature inclusion index, estimates coefficients for the selected features in parallel for each subset, and then averages these estimates. The algorithm is simple, involves very minimal communication, scales efficiently in sample size, and has theoretical guarantees. I provide extensive experiments to show excellent performance in feature selection, estimation, prediction, and computation time relative to usual competitors.

While sample space partitioning is useful in handling datasets with large sample size, feature space partitioning is more effective when the data dimension is high. Existing methods for partitioning features, however, are either vulnerable to high correlations or inefficient in reducing the model dimension. In the thesis, I propose a new embarrassingly parallel framework named {\em DECO} for distributed variable selection and parameter estimation. In {\em DECO}, variables are first partitioned and allocated to m distributed workers. The decorrelated subset data within each worker are then fitted via any algorithm designed for high-dimensional problems. We show that by incorporating the decorrelation step, DECO can achieve consistent variable selection and parameter estimation on each subset with (almost) no assumptions. In addition, the convergence rate is nearly minimax optimal for both sparse and weakly sparse models and does NOT depend on the partition number m. Extensive numerical experiments are provided to illustrate the performance of the new framework.

For datasets with both large sample sizes and high dimensionality, I propose a new "divided-and-conquer" framework {\em DEME} (DECO-message) by leveraging both the {\em DECO} and the {\em message} algorithm. The new framework first partitions the dataset in the sample space into row cubes using {\em message} and then partition the feature space of the cubes using {\em DECO}. This procedure is equivalent to partitioning the original data matrix into multiple small blocks, each with a feasible size that can be stored and fitted in a computer in parallel. The results are then synthezied via the {\em DECO} and {\em message} algorithm in a reverse order to produce the final output. The whole framework is extremely scalable.

Boron isotope and B/Ca ratios of planktonic foraminifera

Here a new analytical methodology is described for measuring the isotopic composition of boron in foraminifera using multicollector inductively coupled plasma mass spectrometry (MC-ICPMS). This new approach is fast (~10 samples analysed in duplicate per analytical session) and accurate (to better than 0.25 per mil at 95% confidence) with acceptable sample size requirements (1-3 mg of carbonate). A core top calibration of several common planktic and two benthic species from geographically widespread localities shows a very close agreement between the isotopic composition measured by MC-ICPMS and the isotopic composition of B(OH)-4 in seawater (as predicted using the recently measured isotopic equilibrium factor of 1.0272) at the depth of habitat. A down core and core top investigation of boron concentration (B/Ca ratio) shows that the partition coefficient is influenced by [CO2-3] complicating the application of this proxy. Nevertheless, it is demonstrated that these two proxies can be used to fully constrain the carbonate system of surface water in the Caribbean Sea (ODP Site 999A) over the last 130 kyr. This reconstruction shows that during much of the Holocene and the last interglacial period surface water at Site 999A was in equilibrium with the atmosphere with respect to CO2. During the intervening colder periods although the surface water pCO2 was lower than the Holocene, it was a minor to significant source of CO2 to the atmosphere possibly due to either an expansion of the eastern equatorial Atlantic upwelling zone, or a more local expansion of coastal upwelling in the southern Caribbean. Such reorganisation of the oceanic carbonate system in favour of a larger source of CO2 to the atmosphere from the equatorial ocean may require mechanisms responsible for lowering atmospheric CO2 during glacial periods to be more efficient than previously supposed.

Meiofauna and the structure of the benthic community in the Iberian deep sea basin

1. On the cruises 3 and 15 of R.V. "Meteor" 6 grab samples, and 6 hauls with the 6 m Agassiztrawl were taken and at 2 stations the deep sea camera was lowered. This material gave quantitative results on the meiofauna and minimum counts of the macrofauna. 2. The nematodes constitute nearly 95% of the meiofauna, the copepoda only 2%. With increasing sediment depth the density of animals decrease gradually. In the uppermost centimeter of sediment 42.6% of the meiofauna are found while only 3.7% live in layer 6-7 cm. Meiofauna weight ranges from 0.6-5.7 mg/25 m**2 surface i.e. 0.24-2.8 g/m**2. 3. Mean numbers of individuals and weights show standard errors of 20-30 %. As an approximate average values for further considerations the weight of the meiofauna in the area was taken as 1 g/m**2 4. Quantitative information on the macrofauna is derived from the trawls and the photographs for the actinia Chitonanthus abyssorum only, which is found in the rate of 1 individual/36-72 m**2, but seems to be less abundant generally. 5. Animal density does not decrease steadily from nearshore to offshore biocoenoses, i.e. generally with increasing depth. The decrease is more pronounced for macro- than for meiofauna. For the deep sea the weight proportion of macrofauna : meiofauna is of the order of 1 : 1. 6. With the assumption, that adaptation of metabolism to deep sea conditions is similar in macro- and meiofauna total metabolism of invertebrates is ascribed to meiofauna to more than 80%. 7. The structure of the biocoenosis of the deep sea floor is characterized by the meiofauna living on and in the sediment and by the dominance of sediment feeders in the macrofauna. 8. Considering the large numbets and high partition rates of bacteria a comparative large part of the metabolism in the deep sea sediment must be ascribed to bacteria. This favours the hypothesis, that with increasing depth and decreasing addition of organic material to the sediment, the importance of meiofauna and microorganisms for total metabolism increases. 9. Considering the different modes of food transport to the deep sea environment, i.e. sinking of dead particles, transport by vertical migration of organisms, aggregation of organic particles, adsorption of dissoloved organic substance to inorganic particles, and heterotrophy, the sediment may be assumed to contain more food for invertebrates than the water above the bottom. 10. Suspensions feeders of macrofauna are fixed to hard substrates in the sediment surface. Some of them are shown to bend themselves down to the bottom in underwater photographs. This suggests the idea that some deep sea suspension feeders partly depend on food from the sediment surface, on which they feed directly.

Seawater carbonate chemistry and elemental composition of purple sea urchin (Strongylocentrotus purpuratus) calcite in a laboratory experiment

Ocean acidification will likely have negative impacts on invertebrates producing skeletons composed of calcium carbonate. Skeletal solubility is partly controlled by the incorporation of "foreign" ions (e.g. magnesium) into the crystal lattice of these skeletal structures, a process that is sensitive to a variety of biological and environmental factors. Here we explore effects of life stage, oceanographic region of origin, and changes in the partial pressure of carbon dioxide in seawater (pCO2) on trace elemental composition in the purple sea urchin (Strongylocentrotus purpuratus). We show that, similar to other urchin taxa, adult purple sea urchins have the ability to precipitate skeleton composed of a range of biominerals spanning low- to high-Mg calcites. Mg / Ca and Sr / Ca ratios were substantially lower in adult spines compared to adult tests. On the other hand, trace elemental composition was invariant among adults collected from four oceanographically distinct regions spanning a range of carbonate chemistry conditions (Oregon, Northern California, Central California, and Southern California). Skeletons of newly settled juvenile urchins that originated from adults from the four regions exhibited intermediate Mg / Ca and Sr / Ca between adult spine and test endmembers, indicating that skeleton precipitated during early life stages is more soluble than adult spines and less soluble than adult tests. Mean skeletal Mg / Ca or Sr / Ca of juvenile skeleton did not vary with source region when larvae were reared under present-day, global-average seawater carbonate conditions (400 µatm; pHT = 8.02 ± 0.03 1 SD; Omega calcite = 3.3 ± 0.2 1 SD). However, when reared under elevated pCO2 (900 µatm; pHT = 7.73 ± 0.03; Omega calcite = 1.8 ± 0.1), skeletal Sr / Ca in juveniles exhibited increased variance across the four regions. Although larvae from the northern populations (Oregon, Northern California, Central California) did not exhibit differences in Mg or Sr incorporation under elevated pCO2 (Sr / Ca = 2.10 ± 0.06 mmol/mol; Mg / Ca = 67.4 ± 3.9 mmol/mol), juveniles of Southern California origin partitioned ~8% more Sr into their skeletons when exposed to higher pCO2 (Sr / Ca = 2.26 ± 0.08 vs. 2.09 ± 0.005 mmol/mol 1 SD). Together these results suggest that the diversity of carbonate minerologies present across different skeletal structures and life stages in purple sea urchins does not translate into an equivalent geochemical plasticity of response associated with geographic variation or temporal shifts in seawater properties. Rather, composition of S. purpuratus skeleton precipitated during both early and adult life history stages appears relatively robust to spatial gradients and predicted future changes in carbonate chemistry. An exception to this trend may arise during early life stages, where certain populations of purple sea urchins may alter skeletal mineral precipitation rates and composition beyond a given pCO2 threshold. This potential for geochemical plasticity during early development in contrast to adult stage geochemical resilience adds to the growing body of evidence that ocean acidification can have differing effects across organismal life stages.

Sr/Ca and Cd/Ca ratios of benthic foraminifera in surface deep-sea sediments

Measurements of Sr/Ca of benthic foraminifera show a linear decrease with water depth which is superimposed upon significant variability identified by analyses of individual foraminifera. New data for Cd/Ca support previous work in defining a contrast between waters shallower and deeper than ~2500 m. Measured element partition coefficients in foraminiferal calcium carbonate relative to sea water (D) have been described by means of a one-box model in which elements are extracted by Rayleigh distillation from a biomineralization reservoir that serves for calcification with a constant fractionation factor (alpha), such that D = (1 - f**alpha)/(l - f), where f is the proportion of Ca remaining after precipitation. A modification to the model recognises differences in element speciation. The model is consistent with differences between D[Sr], D[Ba], and D[Cd] in benthic but not planktonic foraminifera. Depth variations in D for Sr and Ba are consistent with the model, as are differences in depth variation of D[Cd] in calcitic and aragonitic benthic foraminifera. The shallower depth variations may reflect increasing calcification rates with increasing water depth to an optimum of about 2500 m. Observations of unusually lower DCd for some deep waters, not accompanied by similar [Sr], or D[Ba] may be because of dissolution or a calcification response to a lower carbonate saturation state.

Distributed Extremum-Seeking Control with Applications to High-Altitude Balloons

The real-time optimization of large-scale systems is a difficult problem due to the need for complex models involving uncertain parameters and the high computational cost of solving such problems by a decentralized approach. Extremum-seeking control (ESC) is a model-free real-time optimization technique which can estimate unknown parameters and can optimize nonlinear time-varying systems using only a measurement of the cost function to be minimized. In this thesis, we develop a distributed version of extremum-seeking control which allows large-scale systems to be optimized without models and with minimal computing power. First, we develop a continuous-time distributed extremum-seeking controller. It has three main components: consensus, parameter estimation, and optimization. The consensus provides each local controller with an estimate of the cost to be minimized, allowing them to coordinate their actions. Using this cost estimate, parameters for a local input-output model are estimated, and the cost is minimized by following a gradient descent based on the estimate of the gradient. Next, a similar distributed extremum-seeking controller is developed in discrete-time. Finally, we consider an interesting application of distributed ESC: formation control of high-altitude balloons for high-speed wireless internet. These balloons must be steered into a favourable formation where they are spread out over the Earth and provide coverage to the entire planet. Distributed ESC is applied to this problem, and is shown to be effective for a system of 1200 ballons subjected to realistic wind currents. The approach does not require a wind model and uses a cost function based on a Voronoi partition of the sphere. Distributed ESC is able to steer balloons from a few initial launch sites into a formation which provides coverage to the entire Earth and can maintain a similar formation as the balloons move with the wind around the Earth.

Refractometry Studies of the Optical Properties of Polymer Films and the Development of Polymer Coated Refractive Index Sensors

Sensors for real-time monitoring of environmental contaminants are essential for protecting ecosystems and human health. Refractive index sensing is a non-selective technique that can be used to measure almost any analyte. Miniaturized refractive index sensors, such as silicon-on-insulator (SOI) microring resonators are one possible platform, but require coatings selective to the analytes of interest. A homemade prism refractometer is reported and used to characterize the interactions between polymer films and liquid or vapour-phase analytes. A camera was used to capture both Fresnel reflection and total internal reflection within the prism. For thin-films (d = 10 μm - 100 μm), interference fringes were also observed. Fourier analysis of the interferogram allowed for simultaneous extraction of the average refractive index and film thickness with accuracies of ∆n = 1-7 ×10-4 and ∆d < 3-5%. The refractive indices of 29 common organic solvents as well as aqueous solutions of sodium chloride, sucrose, ethylene glycol, glycerol, and dimethylsulfoxide were measured at λ = 1550 nm. These measurements will be useful for future calibrations of near-infrared refractive index sensors. A mathematical model is presented, where the concentration of analyte adsorbed in a film can be calculated from the refractive index and thickness changes during uptake. This model can be used with Fickian diffusion models to measure the diffusion coefficients through the bulk film and at the film-substrate interface. The diffusion of water and other organic solvents into SU-8 epoxy was explored using refractometry and the diffusion coefficient of water into SU-8 is presented. Exposure of soft baked SU-8 films to acetone, acetonitrile and methanol resulted in rapid delamination. The diffusion of volatile organic compound (VOC) vapours into polydimethylsiloxane and polydimethyl-co-polydiphenylsiloxane polymers was also studied using refractometry. Diffusion and partition coefficients are reported for several analytes. As a model system, polydimethyl-co-diphenylsiloxane films were coated onto SOI microring resonators. After the development of data acquisition software, coated devices were exposed to VOCs and the refractive index response was assessed. More studies with other polymers are required to test the viability of this platform for environmental sensing applications.

Assessing hydrothermal liquefaction for the production of bio-oil and enhanced metal recovery from microalgae cultivated on acid mine drainage

The hydrothermal liquefaction(HTL) of algal biomass is a promising route to viable second generation biofuels. In this investigation HTL was assessed for the valorisation of algae used in the remediation of acid mine drainage (AMD). Initially the HTL process was evaluated using Arthrospira platensis (Spirulina) with additional metal sulphates to simulate metal remediation. Optimised conditions were then used to process a natural algal community (predominantly Chlamydomonas sp.) cultivated under two scenarios: high uptake and low uptake of metals from AMD. High metal concentrations appear to catalyse the conversion to bio-oil, and do not significantly affect the heteroatom content or higher heating value of the bio-oil produced. The associated metals were found to partition almost exclusively into the solid residue, favourable for potential metal recovery. High metal loadings also caused partitioning of phosphates from the aqueous phase to the solid phase, potentially compromising attempts to recycle process water as a growth supplement. HTL was therefore found to be a suitable method of processing algae used in AMD remediation, producing a crude oil suitable for upgrading into hydrocarbon fuels, an aqueous and gas stream suitable for supplementing the algal growth and the partitioning of most contaminant metals to the solid residue where they would be readily amenable for recovery and/or disposal.

Assessing hydrothermal liquefaction for the production of bio-oil and enhanced metal recovery from microalgae cultivated on acid mine drainage

The hydrothermal liquefaction(HTL) of algal biomass is a promising route to viable second generation biofuels. In this investigation HTL was assessed for the valorisation of algae used in the remediation of acid mine drainage (AMD). Initially the HTL process was evaluated using Arthrospira platensis (Spirulina) with additional metal sulphates to simulate metal remediation. Optimised conditions were then used to process a natural algal community (predominantly Chlamydomonas sp.) cultivated under two scenarios: high uptake and low uptake of metals from AMD. High metal concentrations appear to catalyse the conversion to bio-oil, and do not significantly affect the heteroatom content or higher heating value of the bio-oil produced. The associated metals were found to partition almost exclusively into the solid residue, favourable for potential metal recovery. High metal loadings also caused partitioning of phosphates from the aqueous phase to the solid phase, potentially compromising attempts to recycle process water as a growth supplement. HTL was therefore found to be a suitable method of processing algae used in AMD remediation, producing a crude oil suitable for upgrading into hydrocarbon fuels, an aqueous and gas stream suitable for supplementing the algal growth and the partitioning of most contaminant metals to the solid residue where they would be readily amenable for recovery and/or disposal.

No discrimination shock avoidance with sequential presentation of stimuli but shore crabs still reduce shock exposure

Insights into the potential for pain may be obtained from examination of behavioural responses to noxious stimuli. In particular, prolonged responses coupled with long-term motivational change and avoidance learning cannot be explained by nociceptive reflex but are consistent with the idea of pain. Here, we placed shore crabs alternately in two halves of a test area divided by an opaque partition. Each area had a dark shelter and in one repeated small electric shocks were delivered in an experimental but not in a control group. Crabs showed no specific avoidance of the shock shelter either during these trials or in a subsequent test in which both were offered simultaneously; however they often emerged from the shock shelter during a trial and thus avoided further shock. More crabs emerged in later trials and took less time to emerge than in early trials. Thus, despite the lack of discrimination learning between the two shelters they used other tactics to markedly reduce the amount of shock received. We note that a previous experiment using simultaneous presentation of two shelters demonstrated rapid discrimination and avoidance learning but the paradigm of sequential presentation appears to prevent this. Nevertheless, the data show clearly that the shock is aversive and tactics, other than discrimination learning, are used to avoid it. Thus, the behaviour is only partially consistent with the idea of pain.

A new insight about pharmaceutical dosage forms for benzathine penicillin G

In this work, a micellar system of benzathine penicillin G (BPG) in sodium deoxycholate (NaDC) was developed and evaluated physicochemically. The solubility profile of the drug in water and buffer solutions at various pH was determined, as well as its n-octanol/water partition coefficient. The Critical Micellar Concentration of NaDC and its ability to incorporate BPG were also assessed. The study was carried out at low and high ionic strength which was adjusted by the addition of sodium chloride. The results demonstrated the ability of the micellar system to incorporate BPG, as well as to increase its apparent solubility in water. The enhancement of the solubility of BPG by the presence of NaDC micelles could be analyzed quantitatively within the framework of the pseudo-phase model. Concentration analysis showed that the micellar system could attain up to 90% incorporation of BPG. The incorporated drug is expected to exhibit improved stability, since the antibiotic enclosed in the hydrophobic core of micelles is rather shielded from the aqueous external environment