Bayesian hybrid algorithms and models : implementation and associated issues

This thesis addresses computational challenges arising from Bayesian analysis of complex real-world problems. Many of the models and algorithms designed for such analysis are ‘hybrid’ in nature, in that they are a composition of components for which their individual properties may be easily described but the performance of the model or algorithm as a whole is less well understood. The aim of this research project is to after a better understanding of the performance of hybrid models and algorithms. The goal of this thesis is to analyse the computational aspects of hybrid models and hybrid algorithms in the Bayesian context. The first objective of the research focuses on computational aspects of hybrid models, notably a continuous finite mixture of t-distributions. In the mixture model, an inference of interest is the number of components, as this may relate to both the quality of model fit to data and the computational workload. The analysis of t-mixtures using Markov chain Monte Carlo (MCMC) is described and the model is compared to the Normal case based on the goodness of fit. Through simulation studies, it is demonstrated that the t-mixture model can be more flexible and more parsimonious in terms of number of components, particularly for skewed and heavytailed data. The study also reveals important computational issues associated with the use of t-mixtures, which have not been adequately considered in the literature. The second objective of the research focuses on computational aspects of hybrid algorithms for Bayesian analysis. Two approaches will be considered: a formal comparison of the performance of a range of hybrid algorithms and a theoretical investigation of the performance of one of these algorithms in high dimensions. For the first approach, the delayed rejection algorithm, the pinball sampler, the Metropolis adjusted Langevin algorithm, and the hybrid version of the population Monte Carlo (PMC) algorithm are selected as a set of examples of hybrid algorithms. Statistical literature shows how statistical efficiency is often the only criteria for an efficient algorithm. In this thesis the algorithms are also considered and compared from a more practical perspective. This extends to the study of how individual algorithms contribute to the overall efficiency of hybrid algorithms, and highlights weaknesses that may be introduced by the combination process of these components in a single algorithm. The second approach to considering computational aspects of hybrid algorithms involves an investigation of the performance of the PMC in high dimensions. It is well known that as a model becomes more complex, computation may become increasingly difficult in real time. In particular the importance sampling based algorithms, including the PMC, are known to be unstable in high dimensions. This thesis examines the PMC algorithm in a simplified setting, a single step of the general sampling, and explores a fundamental problem that occurs in applying importance sampling to a high-dimensional problem. The precision of the computed estimate from the simplified setting is measured by the asymptotic variance of the estimate under conditions on the importance function. Additionally, the exponential growth of the asymptotic variance with the dimension is demonstrated and we illustrates that the optimal covariance matrix for the importance function can be estimated in a special case.

Coupling of energy from quantum emitters to the plasmonic mode of V groove waveguides: A numerical study

This work is a theoretical investigation into the coupling of a single excited quantum emitter to the plasmon mode of a V groove waveguide. The V groove waveguide consists of a triangular channel milled in gold and the emitter is modeled as a dipole emitter, and could represent a quantum dot, nitrogen vacancy in diamond, or similar. In this work the dependence of coupling efficiency of emitter to plasmon mode is determined for various geometrical parameters of the emitter-waveguide system. Using the finite element method, the effect on coupling efficiency of the emitter position and orientation, groove angle, groove depth, and tip radius, is studied in detail. We demonstrate that all parameters, with the exception of groove depth, have a significant impact on the attainable coupling efficiency. Understanding the effect of various geometrical parameters on the coupling between emitters and the plasmonic mode of the waveguide is essential for the design and optimization of quantum dot–V groove devices.

Application of theory in developing peer-support based intervention to improve self-management for cardiac patients with diabetes

The biosafety of carbon nanomaterial needs to be critically evaluated with both experimental and theoretical validations before extensive biomedical applications. In this letter, we present an analysis of the binding ability of two dimensional monolayer carbon nanomaterial on actin by molecular simulation to understand their adhesive characteristics on F-actin cytoskeleton. The modelling results indicate that the positively charged carbon nanomaterial has higher binding stability on actin. Compared to crystalline graphene, graphene oxide shows higher binding influence on actin when carrying 11 positive surface charge. This theoretical investigation provides insights into the sensitivity of actin-related cellular activities on carbon nanomaterial.

CO2 capture and gas separation on boron carbon nanotubes

Concern about the increasing atmospheric CO2 concentration and its impact on the environment has led to increasing attention directed toward finding advanced materials and technologies suited for efficient CO2 capture, storage and purification of clean-burning natural gas. In this letter, we have performed comprehensive theoretical investigation of CO2, N2, CH4 and H2 adsorption on B2CNTs. Our study shows that CO2 molecules can form strong interactions with B2CNTs with different charge states. However, N2, CH4 and H2 can only form very weak interactions with B2CNTs. Therefore, the study demonstrates B2CNTs could sever as promising materials for CO2 capture and gas separation.

Transcending the inanimate : evoking the sublime in puppet-based performance

This study investigated how contemporary puppet-based theatre can create deeply imaginative experiences for adult audiences. Designed to interrogate the potential effects of double-vision (Tillis, 1992), the theories of the sublime (Kant, 2008; 2003) and the uncanny (Jentsch, 1906; Freud, 1919) were used to create a series of creative guidelines. As practice-led research, the project embraced an iterative approach consisting of two cycles for creative experimentation, and a third for the creation of the final performance work The Harbinger, presented as a part of La Boite Theatre Company’s mainstage season. A theoretical investigation was also conducted to inform the developing practice.

Controlled synthesis of inorganic semiconductor nanocrystals and their applications

This thesis is a comprehensive study of the synthesis of nanomaterials. It explores the synthetic methods on the control of the size, shape and phase of semiconductor nanocrystals. A number of important conclusions, including the mechanism behind crystal growth and the structure-relationship, have been drawn through the experimental and theoretical investigation. The synthesized nanocrystals have been tested for applications in gas sensing, photocatalysis and solar cells, which exhibit considerable commercialization potential.

Adhesive characteristics of low dimensional carbon nanomaterial on actin

The biosafety of carbon nanomaterial needs to be critically evaluated with both experimental and theoretical validations before extensive biomedical applications. In this letter, we present an analysis of the binding ability of two dimensional monolayer carbon nanomaterial on actin by molecular simulation to understand their adhesive characteristics on F-actin cytoskeleton. The modelling results indicate that the positively charged carbon nanomaterial has higher binding stability on actin. Compared to crystalline graphene, graphene oxide shows higher binding influence on actin when carrying positive surface charge. This theoretical investigation provides insights into the sensitivity of actin-related cellular activities on carbon nanomaterial.

Theorizing why in e-learning – a frontier for cognitive engagement

Asking why is an important foundation of inquiry and fundamental to the development of reasoning skills and learning. Despite this, and despite the relentless and often disruptive nature of innovations in information and communications technology (ICT), sophisticated tools that directly support this basic act of learning appear to be undeveloped, not yet recognized, or in the very early stages of development. Why is this so? To this question, there is no single satisfactory answer; instead, numerous plausible explanations and related questions arise. After learning something, however, explaining why can be revealing of a person’s understanding (or lack of it). What then differentiates explanation from information; and, explanatory from descriptive content? What ICT scaffolding might support inquiry instigated by why-questioning? What is the role of reflective practice in inquiry-based learning? These and other questions have emerged from this investigation and underscore that why-questions often propagate further questions and are a catalyst for cognitive engagement and dialogue. This paper reports on a multi-disciplinary, theoretical investigation that informs the broad discourse on e-learning and points to a specific frontier for design and development of e-learning tools. Probing why reveals that versatile and ambiguous semantics present the core challenge – asking, learning, knowing, understanding, and explaining why.

Cumulenic and heterocumulenic anions : potential interstellar species?

A recent theoretical investigation by Terzieva & Herbst of linear carbon chains, C-n where n greater than or equal to 6, in the interstellar medium has shown that these species can undergo efficient radiative association to form the corresponding anions. An experimental study by Barckholtz, Snow & Bierbaum of these anions has demonstrated that they do not react efficiently with molecular hydrogen, leading to the possibility of detectable abundances of cumulene-type anions in dense interstellar and circumstellar environments. Here we present a series of electronic structure calculations which examine possible anionic candidates for detection in these media, namely the anion analogues of the previously identified interstellar cumulenes CnH and Cn-1CH2 and heterocumulenes CnO (where n = 2-10). The extraordinary electron affinities calculated for these molecules suggest that efficient radiative electron attachment could occur, and the large dipole moments of these simple (generally) linear molecules point to the possibility of detection by radio astronomy.

Reactions of the hydroperoxide anion with dimethyl methylphosphonate in an ion trap mass spectrometer : evidence for a gas phase alpha-effect

The gas phase degradation reactions of the chemical warfare agent (CWA) simulant, dimethyl methylphosphonate (DMMP), with the hydroperoxide anion (HOO(-)) were investigated using a modified quadrupole ion trap mass spectrometer. The HOO(-) anion reacts readily with neutral DMMP forming two significant product ions at m/z 109 and m/z 123. The major reaction pathways correspond to (i) the nucleophilic substitution at carbon to form \[CH(3)P(O)(OCH(3))O](-) (m/z 109) in a highly exothermic process and (ii) exothermic proton transfer. The branching ratios of the two reaction pathways, 89% and 11% respectively, indicate that the former reaction is significantly faster than the latter. This is in contrast to the trend for the methoxide anion with DMMP, where proton transfer dominates. The difference in the observed reactivities of the HOO(-) and CH(3)O(-) anions can be considered as evidence for an a-effect in the gas phase and is supported by electronic structure calculations at the B3LYP/aug-cc-pVTZ//B3LYP/6-31+G(d) level of theory that indicate the S(N)2(carbon) process has an activation energy 7.8 kJ mol(-1) lower for HOO(-) as compared to CH(3)O(-). A similar alpha-effect was calculated for nucleophilic addition-elimination at phosphorus, but this process an important step in the perhydrolysis degradation of CWAs in solution - was not observed to occur with DMMP in the gas phase. A theoretical investigation revealed that all processes are energetically accessible with negative activation energies. However, comparison of the relative Arrhenius pre-exponential factors indicate that substitution at phosphorus is not kinetically competitive with respect to the S(N)2(carbon) and deprotonation processes.

Non-square-well potential profile and suppression of blinking in compositionally graded Cd1−xZnxSe/CdxZn1−xSe nanocrystals

Random blinking is a major problem on the way to successful applications of semiconducting nanocrystals in optoelectronics and photonics, which until recently had neither a practical solution nor a theoretical interpretation. An experimental breakthrough has recently been made by fabricating non-blinking Cd1-xZnxSe/ZnSe graded nanocrystals [Wang et al., Nature, 2009, 459, 686]. Here, we (1) report an unequivocal and detailed theoretical investigation to understand the properties (e.g., profile) of the potential-well and the distribution of Zn content with respect to the nanocrystal radius and (2) develop a strategy to find the relationship between the photoluminescence (PL) energy peaks and the potential-well due to Zn distribution in nanocrystals. It is demonstrated that the non-square-well potential can be varied in such a way that one can indeed control the PL intensity and the energy-level difference (PL energy peaks) accurately. This implies that one can either suppress the blinking altogether, or alternatively, manipulate the PL energy peaks and intensities systematically to achieve a controlled non-random intermittent luminescence. The approach developed here is based on the ionization energy approximation and as such is generic and can be applied to any non-free-electron nanocrystals.

Architecture, the body and authority in performance

This is a theoretical investigation seeking to learn more about architecture by looking at architectural practice through another discipline. In this research architecture is investigated by examining its relationship with bodies through performance and theatre set design. This thesis aims to build on existing architectural theory, in which an absence of discourse on the body has been identified, by analysing representations of architecture and the body in performance. The research specifically examines the relationship between the body, architecture and authority in performance through the analysis of several performance works.

Catalysis of tRNA Aminoacylation: Single Turnover to Steady-State Kinetics of tRNA Synthetases

Aminoacyl-tRNA synthetases (aaRS) catalyze the bimolecular association reaction between amino acid and tRNA by specifically and unerringly choosing the cognate amino acid and tRNA. There are two classes of such synthetases that perform tRNA-aminoacylation reaction. Interestingly, these two classes of aminoacyl-tRNA synthetases differ not only in their structures but they also exhibit remarkably distinct kinetics under pre-steady-state condition. The class I synthetases show initial burst of product formation followed by a slower steady-state rate. This has been argued to represent the influence of slow product release. In contrast, there is no burst in the case of class H enzymes. The tight binding of product with enzyme for class I enzymes is correlated with the enhancement of rate in presence of elongation factor. EF-TU. In spite of extensive experimental studies, there is no detailed theoretical analysis that can provide a quantitative understanding of this important problem. In this article, we present a theoretical investigation of enzyme kinetics for both classes of aminoacyl-tRNA synthetases. We present an augmented kinetic scheme and then employ the methods of time-dependent probability statistics to obtain expressions for the first passage time distribution that gives both the time-dependent and the steady-state rates. The present study quantitatively explains all the above experimental observations. We propose an alternative path way in the case of class II enzymes showing the tRNA-dependent amino acid activation and the discrepancy between the single-turnover and steady-state rate.

Tertiary amine-based glutathione peroxidase mimics: some insights into the role of steric and electronic effects on antioxidant activity

In this work, several tertiary amine-based diaryl diselenides were synthesized and evaluated for their glutathione peroxidase (GPx)-like antioxidant activities using hydrogen peroxide, tert-butyl hydroperoxide and cumene hydroperoxide as substrates and thiophenol (PhSH) and glutathione (GSH) as co-substrates. A comparison of the GPx-like activity of 4-methoxy-substituted N,N-dialkylbenzylamine-based diselenides with that of the corresponding 6-methoxy-substituted compounds indicates that the activity highly depends on the position of the methoxy substituent. Although the methoxy group at 4- and 6-position alters the electronic properties of selenium, the substitution at the 6-position provides the required steric protection for some of the key intermediates in the catalytic cycle. A detailed experimental and theoretical investigation reveals that the 6-methoxy substituent prevents the undesired thiol exchange reactions at the selenium centers in the selenenyl sulfide intermediates. The 6-methoxy substituent also prevents the formation of seleninic and selenonic acids. When PhSH is used as the thiol co-substrate, the 4-methoxy-substituted diselenides exhibit GPx-like activity similar to that of the parent compounds as the 4-methoxy substituent does not block the selenium center in the selenenyl sulfide intermediates from thiol exchange reactions. In contrast, the 4-methoxy substituent significantly enhances the GPx-like activity of the diselenides when glutathione (GSH) is used as the co-substrate. (C) 2012 Elsevier Ltd. All rights reserved.

Kinetic proofreading at single molecular level: aminoacylation of tRNA(Ile) and the role of water as an editor

Proofreading/editing in protein synthesis is essential for accurate translation of information from the genetic code. In this article we present a theoretical investigation of efficiency of a kinetic proofreading mechanism that employs hydrolysis of the wrong substrate as the discriminatory step in enzyme catalytic reactions. We consider aminoacylation of tRNA(Ile) which is a crucial step in protein synthesis and for which experimental results are now available. We present an augmented kinetic scheme and then employ methods of stochastic simulation algorithm to obtain time dependent concentrations of different substances involved in the reaction and their rates of formation. We obtain the rates of product formation and ATP hydrolysis for both correct and wrong substrates (isoleucine and valine in our case, respectively), in single molecular enzyme as well as ensemble enzyme kinetics. The present theoretical scheme correctly reproduces (i) the amplitude of the discrimination factor in the overall rates between isoleucine and valine which is obtained as (1.8x10(2)).(4.33x10(2)) = 7.8x10(4), (ii) the rates of ATP hydrolysis for both Ile and Val at different substrate concentrations in the aminoacylation of tRNA(Ile). The present study shows a non-michaelis type dependence of rate of reaction on tRNA(Ile) concentration in case of valine. The overall editing in steady state is found to be independent of amino acid concentration. Interestingly, the computed ATP hydrolysis rate for valine at high substrate concentration is same as the rate of formation of Ile-tRNA(Ile) whereas at intermediate substrate concentration the ATP hydrolysis rate is relatively low. We find that the presence of additional editing domain in class I editing enzyme makes the kinetic proofreading more efficient through enhanced hydrolysis of wrong product at the editing CP1 domain.