Numerical solution of an optimal control model of dendritic cell treatment of a growing tumour

A new optimal control model of the interactions between a growing tumour and the host immune system along with an immunotherapy treatment strategy is presented. The model is based on an ordinary differential equation model of interactions between the growing tu- mour and the natural killer, cytotoxic T lymphocyte and dendritic cells of the host immune system, extended through the addition of a control function representing the application of a dendritic cell treat- ment to the system. The numerical solution of this model, obtained from a multi species Runge–Kutta forward-backward sweep scheme, is described. We investigate the effects of varying the maximum al- lowed amount of dendritic cell vaccine administered to the system and find that control of the tumour cell population is best effected via a high initial vaccine level, followed by reduced treatment and finally cessation of treatment. We also found that increasing the strength of the dendritic cell vaccine causes an increase in the number of natural killer cells and lymphocytes, which in turn reduces the growth of the tumour.

Development of optimal designs of insulated rail joints

Proper functioning of Insulated Rail Joints (IRJs) is essential for the safe operation of the railway signalling systems and broken rail identification circuitries. The Conventional IRJ (CIRJ) resembles structural butt joints consisting of two pieces of rails connected together through two joint bars on either side of their web and the assembly is held together through pre-tensioned bolts. As the IRJs should maintain electrical insulation between the two rails, a gap between the rail ends must be retained at all times and all metal contacting surfaces should be electrically isolated from each other using non-conductive material. At the gap, the rail ends lose longitudinal continuity and hence the vertical sections of the rail ends are often severely damaged, especially at the railhead, due to the passage of wheels compared to other continuously welded rail sections. Fundamentally, the reason for the severe damage can be related to the singularities of the wheel-rail contact pressure and the railhead stress. No new generation designs that have emerged in the market to date have focussed on this fundamental; they only have provided attention to either the higher strength materials or the thickness of the sections of various components of the IRJs. In this thesis a novel method of shape optimisation of the railhead is developed to eliminate the pressure and stress singularities through changes to the original sharp corner shaped railhead into an arc profile in the longitudinal direction. The optimal shape of the longitudinal railhead profile has been determined using three nongradient methods in search of accuracy and efficiency: (1) Grid Search Method; (2) Genetic Algorithm Method and (3) Hybrid Genetic Algorithm Method. All these methods have been coupled with a parametric finite element formulation for the evaluation of the objective function for each iteration or generation depending on the search algorithm employed. The optimal shape derived from these optimisation methods is termed as Stress Minimised Railhead (SMRH) in this thesis. This optimal SMRH design has exhibited significantly reduced stress concentration that remains well below the yield strength of the head hardened rail steels and has shifted the stress concentration location away from the critical zone of the railhead end. The reduction in the magnitude and the relocation of the stress concentration in the SMRH design has been validated through a full scale wheel – railhead interaction test rig; Railhead strains under the loaded wheels have been recorded using a non-contact digital image correlation method. Experimental study has confirmed the accuracy of the numerical predications. Although the SMRH shaped IRJs eliminate stress singularities, they can still fail due to joint bar or bolt hole cracking; therefore, another conceptual design, termed as Embedded IRJ (EIRJ) in this thesis, with no joint bars and pre-tensioned bolts has been developed using a multi-objective optimisation formulation based on the coupled genetic algorithm – parametric finite element method. To achieve the required structural stiffness for the safe passage of the loaded wheels, the rails were embedded into the concrete of the post tensioned sleepers; the optimal solutions for the design of the EIRJ is shown to simplify the design through the elimination of the complex interactions and failure modes of the various structural components of the CIRJ. The practical applicability of the optimal shapes SMRH and EIRJ is demonstrated through two illustrative examples, termed as improved designs (IMD1 & IMD2) in this thesis; IMD1 is a combination of the CIRJ and the SMRH designs, whilst IMD2 is a combination of the EIRJ and SMRH designs. These two improved designs have been simulated for two key operating (speed and wagon load) and design (wheel diameter) parameters that affect the wheel-rail contact; the effect of these parameters has been found to be negligible to the performance of the two improved designs and the improved designs are in turn found far superior to the current designs of the CIRJs in terms of stress singularities and deformation under the passage of the loaded wheels. Therefore, these improved designs are expected to provide longer service life in relation to the CIRJs.

Van der Waals-corrected density functional theory : benchmarking for hydrogen–nanotube and nanotube–nanotube interactions

Density functional theory (DFT) is a powerful approach to electronic structure calculations in extended systems, but suffers currently from inadequate incorporation of long-range dispersion, or Van der Waals (VdW) interactions. VdW-corrected DFT is tested for interactions involving molecular hydrogen, graphite, single-walled carbon nanotubes (SWCNTs), and SWCNT bundles. The energy correction, based on an empirical London dispersion term with a damping function at short range, allows a reasonable physisorption energy and equilibrium distance to be obtained for H2 on a model graphite surface. The VdW-corrected DFT calculation for an (8, 8) nanotube bundle reproduces accurately the experimental lattice constant. For H2 inside or outside an (8, 8) SWCNT, we find the binding energies are respectively higher and lower than that on a graphite surface, correctly predicting the well known curvature effect. We conclude that the VdW correction is a very effective method for implementing DFT calculations, allowing a reliable description of both short-range chemical bonding and long-range dispersive interactions. The method will find powerful applications in areas of SWCNT research where empirical potential functions either have not been developed, or do not capture the necessary range of both dispersion and bonding interactions.

Adatom CCV Auger rates via the local density of states

In this letter the core-core-valence Auger transitions of an atomic impurity, both in bulk or adsorbed on a jellium-like surface, are computed within a DFT framework. The Auger rates calculated by the Fermi golden rule are compared with those determined by an approximate and simpler expression. This is based on the local density of states (LDOS) with a core hole present, in a region around the impurity nucleus. Different atoms, Na and Mg, solids, Al and Ag, and several impurity locations are considered. We obtain an excellent agreement between KL1V and KL23V rates worked out with the two approaches. The radius of the sphere in which we calculate the LDOS is the relevant parameter of the simpler approach. Its value only depends on the atomic species regardless of the location of the impurity and the type of substrate. (C) 2003 Elsevier B.V. All rights reserved.

A family of Alltop functions that are EA-inequivalent to the cubic function

Sequences with optimal correlation properties are much sought after for applications in communication systems. In 1980, Alltop (\emph{IEEE Trans. Inf. Theory} 26(3):350-354, 1980) described a set of sequences based on a cubic function and showed that these sequences were optimal with respect to the known bounds on auto and crosscorrelation. Subsequently these sequences were used to construct mutually unbiased bases (MUBs), a structure of importance in quantum information theory. The key feature of this cubic function is that its difference function is a planar function. Functions with planar difference functions have been called \emph{Alltop functions}. This paper provides a new family of Alltop functions and establishes the use of Alltop functions for construction of sequence sets and MUBs.

The effect of olive leaf extract on platelet function

Background Epidemiological studies have shown a reduced incidence of cardiovascular disease in the Mediterranean population attributed to the consumption of dietary olive oil rich in antioxidants. This has lead to increased interest in the antioxidant properties of other phenolic compounds of olive tree products. It has been suggested that olive leaf extract may also have health benefits due to its antioxidant and anti-inflammatory activities. Antioxidants can prevent the effects of oxidative metabolism by scavenging free radicals and decreasing the hyperactivity of platelets associated with the development of occlusive thrombosis. No studies to date have investigated the effects of olive leaf extract on platelet function to our knowledge. Improved understanding of the antioxidant properties of olive leaf extract and its effect on platelet function could lead to improved cardiovascular health. Objective The current study used an olive leaf extract prepared from the Olea europaea L. tree. The aim was to determine if polyphenols in olive leaf extract would reduce platelet activity and, to establish an optimal dose in vitro that would reduce platelet aggregation and ATP release. Design Eleven subjects with normal platelet counts (150–400 x 109/L) were recruited for the current in vitro study. Olive leaf extract was added to citrated whole blood to obtain five concentrations ranging from 5.4 ug/mL to 54.0 ug/mL for a dose response curve. Baseline samples, without olive leaf extract were used as a negative control for each subject. After 2 hours incubation with olive leaf extract samples were analyzed for platelet aggregation and ATP release from platelets stimulated by the addition of collagen. Results Whole blood analysis (n=11) showed a clear dose-dependant reduction in platelet aggregation with the increasing olive leaf extract concentrations (p<0.0001). There was also a similar decrease in ATP release from collagen stimulated platelets (p=0.02). Conclusion In the current study the olive leaf extract obtained from Olea europaea L. inhibited platelet aggregation and ATP release from collagen stimulated platelets in vitro. This study suggests olive leaf extract may prevent occlusive thrombosis by reducing platelet hyperactivity.

Analysis and prediction of the critical current density across [001]-tilt YBa2Cu3O7- grain boundaries of arbitrary misorientation angles

A qualitative analysis of the expected dilatation strain field in the vicinity of an array of grain-boundary (GB) dislocations is presented. The analysis provides a basis for the prediction of the critical current densities (jc) across low-angle YBa2Cu3O7- (YBCO) GBs as a function of their energy. The introduction of the GB energy allows the extension of the analysis to high-angle GBs using established models which predict the GB energy as a function of misorientation angle. The results are compared to published data for jc across [001]-tilt YBCO GBs for the full range of misorientations, showing a good fit. Since the GB energy is directly related to the GB structure, the analysis may allow a generalization of the scaling behavior of jc with the GB energy. © 1995 The American Physical Society.

Effects of dissolved oxygen availability and culture biomass at induction upon the intracellular expression of monoamine oxidase by recombinant E. coli in fed batch bioprocesses

The effects of oxygen availability and induction culture biomass upon production of an industrially important monoamine oxidase (MAO) were investigated in fed-batch cultures of a recombinant E. coli. For each induction cell biomass 2 different oxygenation methods were used, aeration and oxygen enriched air. Induction at higher biomass levels increased the culture demand for oxygen, leading to fermentative metabolism and accumulation of high levels of acetate in the aerated cultures. Paradoxically, despite an almost eight fold increase in acetate accumulation to levels widely reported to be highly detrimental to protein production, when induction wet cell weight (WCW) rose from 100% to 137.5%, MAO specific activity in these aerated processes showed a 3 fold increase. By contrast, for oxygenated cultures induced at WCW's 100% and 137.5% specific activity levels were broadly similar, but fell rapidly after the maxima were reached. Induction at high biomass levels (WCW 175%) led to very low levels of specific MAO activity relative to induction at lower WCW's in both aerated and oxygenated cultures. Oxygen enrichment of these cultures was a useful strategy for boosting specific growth rates, but did not have positive effects upon specific enzyme activity. Based upon our findings, consideration of the amino acid composition of MAO and previous studies on related enzymes, we propose that this effect is due to oxidative damage to the MAO enzyme itself during these highly aerobic processes. Thus, the optimal process for MAO production is aerated, not oxygenated, and induced at moderate cell density, and clearly represents a compromise between oxygen supply effects on specific growth rate/induction cell density, acetate accumulation, and high specific MAO activity. This work shows that the negative effects of oxygen previously reported in free enzyme preparations, are not limited to these acellular environments but are also discernible in the sheltered environment of the cytosol of E. coli cells.

Output feedback optimal control with contraints

In Chapters 1 through 9 of the book (with the exception of a brief discussion on observers and integral action in Section 5.5 of Chapter 5) we considered constrained optimal control problems for systems without uncertainty, that is, with no unmodelled dynamics or disturbances, and where the full state was available for measurement. More realistically, however, it is necessary to consider control problems for systems with uncertainty. This chapter addresses some of the issues that arise in this situation. As in Chapter 9, we adopt a stochastic description of uncertainty, which associates probability distributions to the uncertain elements, that is, disturbances and initial conditions. (See Section 12.6 for references to alternative approaches to model uncertainty.) When incomplete state information exists, a popular observer-based control strategy in the presence of stochastic disturbances is to use the certainty equivalence [CE] principle, introduced in Section 5.5 of Chapter 5 for deterministic systems. In the stochastic framework, CE consists of estimating the state and then using these estimates as if they were the true state in the control law that results if the problem were formulated as a deterministic problem (that is, without uncertainty). This strategy is motivated by the unconstrained problem with a quadratic objective function, for which CE is indeed the optimal solution (˚Astr¨om 1970, Bertsekas 1976). One of the aims of this chapter is to explore the issues that arise from the use of CE in RHC in the presence of constraints. We then turn to the obvious question about the optimality of the CE principle. We show that CE is, indeed, not optimal in general. We also analyse the possibility of obtaining truly optimal solutions for single input linear systems with input constraints and uncertainty related to output feedback and stochastic disturbances.We first find the optimal solution for the case of horizon N = 1, and then we indicate the complications that arise in the case of horizon N = 2. Our conclusion is that, for the case of linear constrained systems, the extra effort involved in the optimal feedback policy is probably not justified in practice. Indeed, we show by example that CE can give near optimal performance. We thus advocate this approach in real applications.

Statistical analysis of reduction in tensile strength of cotton strips as a measure of soil microbial activity

The cotton strip assay (CSA) is an established technique for measuring soil microbial activity. The technique involves burying cotton strips and measuring their tensile strength after a certain time. This gives a measure of the rotting rate, R, of the cotton strips. R is then a measure of soil microbial activity. This paper examines properties of the technique and indicates how the assay can be optimised. Humidity conditioning of the cotton strips before measuring their tensile strength reduced the within and between day variance and enabled the distribution of the tensile strength measurements to approximate normality. The test data came from a three-way factorial experiment (two soils, two temperatures, three moisture levels). The cotton strips were buried in the soil for intervals of time ranging up to 6 weeks. This enabled the rate of loss of cotton tensile strength with time to be studied under a range of conditions. An inverse cubic model accounted for greater than 90% of the total variation within each treatment combination. This offers support for summarising the decomposition process by a single parameter R. The approximate variance of the decomposition rate was estimated from a function incorporating the variance of tensile strength and the differential of the function for the rate of decomposition, R, with respect to tensile strength. This variance function has a minimum when the measured strength is approximately 2/3 that of the original strength. The estimates of R are almost unbiased and relatively robust against the cotton strips being left in the soil for more or less than the optimal time. We conclude that the rotting rate X should be measured using the inverse cubic equation, and that the cotton strips should be left in the soil until their strength has been reduced to about 2/3.

Optimal Bayesian experimental design for models with intractable likelihoods using indirect inference applied to biological process models

This paper addresses the problem of determining optimal designs for biological process models with intractable likelihoods, with the goal of parameter inference. The Bayesian approach is to choose a design that maximises the mean of a utility, and the utility is a function of the posterior distribution. Therefore, its estimation requires likelihood evaluations. However, many problems in experimental design involve models with intractable likelihoods, that is, likelihoods that are neither analytic nor can be computed in a reasonable amount of time. We propose a novel solution using indirect inference (II), a well established method in the literature, and the Markov chain Monte Carlo (MCMC) algorithm of Müller et al. (2004). Indirect inference employs an auxiliary model with a tractable likelihood in conjunction with the generative model, the assumed true model of interest, which has an intractable likelihood. Our approach is to estimate a map between the parameters of the generative and auxiliary models, using simulations from the generative model. An II posterior distribution is formed to expedite utility estimation. We also present a modification to the utility that allows the Müller algorithm to sample from a substantially sharpened utility surface, with little computational effort. Unlike competing methods, the II approach can handle complex design problems for models with intractable likelihoods on a continuous design space, with possible extension to many observations. The methodology is demonstrated using two stochastic models; a simple tractable death process used to validate the approach, and a motivating stochastic model for the population evolution of macroparasites.

Rapid, low-temperature synthesis of nc-Si in high-density, non-equilibrium plasmas : enabling nanocrystallinity at very low hydrogen dilution

Nanocrystalline silicon thin films were deposited on single-crystal silicon and glass substrates simultaneously by inductively coupled plasma-assisted chemical vapor deposition from the reactive silane reactant gas diluted with hydrogen at a substrate temperature of 200 °C. The effect of hydrogen dilution ratio X (X is defined as the flow rate ratio of hydrogen to silane gas), ranging from 1 to 20, on the structural and optical properties of the deposited films, is extensively investigated by Raman spectroscopy, X-ray diffraction, Fourier transform infrared absorption spectroscopy, UV/VIS spectroscopy, and scanning electron microscopy. Our experimental results reveal that, with the increase of the hydrogen dilution ratio X, the deposition rate Rd and hydrogen content CH are reduced while the crystalline fraction Fc, mean grain size δ and optical bandgap ETauc are increased. In comparison with other plasma enhanced chemical vapor deposition methods of nanocrystalline silicon films where a very high hydrogen dilution ratio X is routinely required (e.g. X > 16), we have achieved nanocrystalline silicon films at a very low hydrogen dilution ratio of 1, featuring a high deposition rate of 1.57 nm/s, a high crystalline fraction of 67.1%, a very low hydrogen content of 4.4 at.%, an optical bandgap of 1.89 eV, and an almost vertically aligned columnar structure with a mean grain size of approximately 19 nm. We have also shown that a sufficient amount of atomic hydrogen on the growth surface essential for the formation of nanocrystalline silicon is obtained through highly-effective dissociation of silane and hydrogen molecules in the high-density inductively coupled plasmas. © 2009 The Royal Society of Chemistry.

Structural evolution of nanocrystalline silicon thin films synthesized in high-density, low-temperature reactive plasmas

Silicon thin films with a variable content of nanocrystalline phase were deposited on single-crystal silicon and glass substrates by inductively coupled plasma-assisted chemical vapor deposition using a silane precursor without any hydrogen dilution in the low substrate temperature range from 100 to 300 °C. The structural and optical properties of the deposited films are systematically investigated by Raman spectroscopy, x-ray diffraction, Fourier transform infrared absorption spectroscopy, UV/vis spectroscopy, scanning electron microscopy and high-resolution transmission electron microscopy. It is shown that the structure of the silicon thin films evolves from the purely amorphous phase to the nanocrystalline phase when the substrate temperature is increased from 100 to 150 °C. It is found that the variations of the crystalline fraction fc, bonded hydrogen content CH, optical bandgap ETauc, film microstructure and growth rate Rd are closely related to the substrate temperature. In particular, at a substrate temperature of 300 °C, the nanocrystalline Si thin films of our interest feature a high growth rate of 1.63nms-1, a low hydrogen content of 4.0at.%, a high crystalline fraction of 69.1%, a low optical bandgap of 1.55eV and an almost vertically aligned columnar structure with a mean grain size of approximately 10nm. It is also shown that the low-temperature synthesis of nanocrystalline Si thin films without any hydrogen dilution is attributed to the outstanding dissociation ability of the high-density inductively coupled plasmas and effective plasma-surface interactions during the growth process. Our results offer a highly effective yet simple and environmentally friendly technique to synthesize high-quality nanocrystalline Si films, vitally needed for the development of new-generation solar cells and other emerging nanotechnologies.

Electron energy distribution function in low-pressure complex plasmas

The effect of density and size of dust grains on the electron energy distribution function (EEDF) in low-temperature complex plasmas is studied. It is found that the EEDF depends strongly on the dust density and size. The behavior of the electron temperature can differ significantly from that of a pristine plasma. For low-pressure argon glow discharge, the Druyvesteyn-like EEDF often found in pristine plasmas can become nearly Maxwellian if the dust density and/or sizes are large. One can thus control the plasma parameters by the dust grains.

A review of modern computational algorithms for Bayesian optimal design

Bayesian experimental design is a fast growing area of research with many real-world applications. As computational power has increased over the years, so has the development of simulation-based design methods, which involve a number of algorithms, such as Markov chain Monte Carlo, sequential Monte Carlo and approximate Bayes methods, facilitating more complex design problems to be solved. The Bayesian framework provides a unified approach for incorporating prior information and/or uncertainties regarding the statistical model with a utility function which describes the experimental aims. In this paper, we provide a general overview on the concepts involved in Bayesian experimental design, and focus on describing some of the more commonly used Bayesian utility functions and methods for their estimation, as well as a number of algorithms that are used to search over the design space to find the Bayesian optimal design. We also discuss other computational strategies for further research in Bayesian optimal design.