A numerical approach to modeling the catalytic voltammetry of surface-confined redox enzymes

A finite difference method for simulating voltammograms of electrochemically driven enzyme catalysis is presented. The method enables any enzyme mechanism to be simulated. The finite difference equations can be represented as a matrix equation containing a nonlinear sparse matrix. This equation has been solved using the software package Mathematica. Our focus is on the use of cyclic voltammetry since this is the most commonly employed electrochemical method used to elucidate mechanisms. The use of cyclic voltammetry to obtain data from systems obeying Michaelis-Menten kinetics is discussed, and we then verify our observations on the Michaelis-Menten system using the finite difference simulation. Finally, we demonstrate how the method can be used to obtain mechanistic information on a real redox enzyme system, the complex bacterial molybdoenzyme xanthine dehydrogenase.

Similarity solution of axisymmetric flow in porous media

Applications of the axisymmetric Boussinesq equation to groundwater hydrology and reservoir engineering have long been recognised. An archetypal example is invasion by drilling fluid into a permeable bed where there is initially no such fluid present, a circumstance of some importance in the oil industry. It is well known that the governing Boussinesq model can be reduced to a nonlinear ordinary differential equation using a similarity variable, a transformation that is valid for a certain time-dependent flux at the origin. Here, a new analytical approximation is obtained for this case. The new solution,, which has a simple form, is demonstrated to be highly accurate. (c) 2005 Elsevier Ltd. All rights reserved.

Optimal strategies in political elections

In the Majoritarian Parliamentary System, the government has a constitutional right to call an early election. This right provides the government a control to achieve its objective to remain in power for as long as possible. We model the early election problem mathematically using opinion polls data as a stochastic process to proxy the government's probability of re-election. These data measure the difference in popularity between the government and the opposition. We fit a mean reverting Stochastic Differential Equation to describe the behaviour of the process and consider the possibility for the government to use other control tools, which are termed 'boosts' to induce shocks to the opinion polls by making timely policy announcements or economic actions. These actions improve the government's popularity and have some impact upon the early-election exercise boundary. © Austral. Mathematical Soc. 2005.

A geometric approach to quantum circuit lower bounds

What is the minimal size quantum circuit required to exactly implement a specified n-qubit unitary operation, U, without the use of ancilla qubits? We show that a lower bound on the minimal size is provided by the length of the minimal geodesic between U and the identity, I, where length is defined by a suitable Finsler metric on the manifold SU(2(n)). The geodesic curves on these manifolds have the striking property that once an initial position and velocity are set, the remainder of the geodesic is completely determined by a second order differential equation known as the geodesic equation. This is in contrast with the usual case in circuit design, either classical or quantum, where being given part of an optimal circuit does not obviously assist in the design of the rest of the circuit. Geodesic analysis thus offers a potentially powerful approach to the problem of proving quantum circuit lower bounds. In this paper we construct several Finsler metrics whose minimal length geodesics provide lower bounds on quantum circuit size. For each Finsler metric we give a procedure to compute the corresponding geodesic equation. We also construct a large class of solutions to the geodesic equation, which we call Pauli geodesics, since they arise from isometries generated by the Pauli group. For any unitary U diagonal in the computational basis, we show that: (a) provided the minimal length geodesic is unique, it must be a Pauli geodesic; (b) finding the length of the minimal Pauli geodesic passing from I to U is equivalent to solving an exponential size instance of the closest vector in a lattice problem (CVP); and (c) all but a doubly exponentially small fraction of such unitaries have minimal Pauli geodesics of exponential length.

Stochastic dynamic programming for election timing: A game theory approach

In this paper, we consider dynamic programming for the election timing in the majoritarian parliamentary system such as in Australia, where the government has a constitutional right to call an early election. This right can give the government an advantage to remain in power for as long as possible by calling an election, when its popularity is high. On the other hand, the opposition's natural objective is to gain power, and it will apply controls termed as "boosts" to reduce the chance of the government being re-elected by introducing policy and economic responses. In this paper, we explore equilibrium solutions to the government, and the opposition strategies in a political game using stochastic dynamic programming. Results are given in terms of the expected remaining life in power, call and boost probabilities at each time at any level of popularity.

Stochastic models for regulatory networks of the genetic toggle switch

Bistability arises within a wide range of biological systems from the A phage switch in bacteria to cellular signal transduction pathways in mammalian cells. Changes in regulatory mechanisms may result in genetic switching in a bistable system. Recently, more and more experimental evidence in the form of bimodal population distributions indicates that noise plays a very important role in the switching of bistable systems. Although deterministic models have been used for studying the existence of bistability properties under various system conditions, these models cannot realize cell-to-cell fluctuations in genetic switching. However, there is a lag in the development of stochastic models for studying the impact of noise in bistable systems because of the lack of detailed knowledge of biochemical reactions, kinetic rates, and molecular numbers. in this work, we develop a previously undescribed general technique for developing quantitative stochastic models for large-scale genetic regulatory networks by introducing Poisson random variables into deterministic models described by ordinary differential equations. Two stochastic models have been proposed for the genetic toggle switch interfaced with either the SOS signaling pathway or a quorum-sensing signaling pathway, and we have successfully realized experimental results showing bimodal population distributions. Because the introduced stochastic models are based on widely used ordinary differential equation models, the success of this work suggests that this approach is a very promising one for studying noise in large-scale genetic regulatory networks.

Towards a Modeling Environment for Earth Systems Simulations

The developments of models in Earth Sciences, e.g. for earthquake prediction and for the simulation of mantel convection, are fare from being finalized. Therefore there is a need for a modelling environment that allows scientist to implement and test new models in an easy but flexible way. After been verified, the models should be easy to apply within its scope, typically by setting input parameters through a GUI or web services. It should be possible to link certain parameters to external data sources, such as databases and other simulation codes. Moreover, as typically large-scale meshes have to be used to achieve appropriate resolutions, the computational efficiency of the underlying numerical methods is important. Conceptional this leads to a software system with three major layers: the application layer, the mathematical layer, and the numerical algorithm layer. The latter is implemented as a C/C++ library to solve a basic, computational intensive linear problem, such as a linear partial differential equation. The mathematical layer allows the model developer to define his model and to implement high level solution algorithms (e.g. Newton-Raphson scheme, Crank-Nicholson scheme) or choose these algorithms form an algorithm library. The kernels of the model are generic, typically linear, solvers provided through the numerical algorithm layer. Finally, to provide an easy-to-use application environment, a web interface is (semi-automatically) built to edit the XML input file for the modelling code. In the talk, we will discuss the advantages and disadvantages of this concept in more details. We will also present the modelling environment escript which is a prototype implementation toward such a software system in Python (see www.python.org). Key components of escript are the Data class and the PDE class. Objects of the Data class allow generating, holding, accessing, and manipulating data, in such a way that the actual, in the particular context best, representation is transparent to the user. They are also the key to establish connections with external data sources. PDE class objects are describing (linear) partial differential equation objects to be solved by a numerical library. The current implementation of escript has been linked to the finite element code Finley to solve general linear partial differential equations. We will give a few simple examples which will illustrate the usage escript. Moreover, we show the usage of escript together with Finley for the modelling of interacting fault systems and for the simulation of mantel convection.

Quantum nondemolition measurement of Fock states of mesoscopic mechanical oscillators

We investigate a scheme that makes a quantum nondemolition (QND) measurement of the excitation level of a mesoscopic mechanical oscillator by utilizing the anharmonic coupling between two beam bending modes. The nonlinear coupling between the two modes shifts the resonant frequency of the readout oscillator in proportion to the excitation level of the system oscillator. This frequency shift may be detected as a phase shift of the readout oscillation when driven on resonance. We derive an equation for the reduced density matrix of the system oscillator, and use this to study the conditions under which discrete jumps in the excitation level occur. The appearance of jumps in the actual quantity measured is also studied using the method of quantum trajectories. We consider the feasibility of the scheme for experimentally accessible parameters.

Gelatinisation of starch in mixtures of sugars. II. Application of differential scanning calorimetry

Differential scanning calorimetry was used to investigate the effect of mixtures of glucose and fructose, and five types of honeys on starch gelatinisation. At a 1:1 starch:water ratio, glucose generally increased the enthalpy (DeltaH(gel)) and temperatures (T-onset, T-peak and T-end) of gelatinisation more than fructose. Upon mixing, DeltaH(gel) of the low-temperature endotherm decreased in comparison to the sole sugars, but was fairly constant (7.7 +/- 0.33 J/g dry starch). DeltaH(gel) of the high-temperature endotherm increased with the fructose content. For both endotherms, the gelatinisation temperatures were unchanged (CV less than or equal to 3%) for the mixtures. With the honeys (moisture, 14.9-18.0%; fructose, 37.2-44.0%; glucose, 28.3-31.9%) added at 1.1-4.4 g per g dry starch, the enthalpy and temperatures of gelatinisation did not vary significantly (CV less than or equal to 6%). Typical thermograms are presented, and the results are interpreted in the light of the various proposed mechanisms for starch gelatinisation in sugar-water systems, total sugar content and possible sugar-sugar interactions. The thermograms were broader in the presence of the sugars and honeys, and a biphasic character was consistently exhibited. The application of an exponential equation to the gelatinisation temperatures of the starch-honey mixtures revealed an opposing influence of fructose and glucose during gelatinisation. The mechanism of starch gelatinisation may be better understood if techniques could be perfected to quantify breakage and formation of hydrogen bonds in the starch granules, and suggested techniques are discussed. (C) 2004 Elsevier Ltd. All rights reserved.

Central elements of the elliptic Zn monodromy matrix algebra at roots of unity

The central elements of the algebra of monodromy matrices associated with the Z(n) R-matrix are studied. When the crossing parameter w takes a special rational value w = n/N, where N and n are positive coprime integers, the center is substantially larger than that in the generic case for which the quantum determinant provides the center. In the trigonometric limit, the situation corresponds to the quantum group at roots of unity. This is a higher rank generalization of the recent results by Belavin and Jimbo. (c) 2004 Elsevier B.V. All rights reserved.

On the Schrodinger equation involving a critical Sobolev exponent and magnetic field

We consider the semilinear Schrodinger equation -Delta(A)u + V(x)u = Q(x)vertical bar u vertical bar(2* -2) u. Assuming that V changes sign, we establish the existence of a solution u not equal 0 in the Sobolev space H-A,V(1) + (R-N). The solution is obtained by a min-max type argument based on a topological linking. We also establish certain regularity properties of solutions for a rather general class of equations involving the operator -Delta(A).

Specific heat capacity of Australian honeys from 35 to 165C as a function of composition using differential scanning calorimetry

Modulated temperature differential scanning calorimetry was used to investigate the specific heat capacity (C-p) of 10 Australian honeys within the processing and handling temperatures. The values obtained were found to be different from the literature values at certain temperatures, and are not predictable by the additive model. The C-p of each honey exhibited a cubic relationship (P < 0.001) with the temperature (T, C). In addition, the moisture (M, %), fructose (F, %) and glucose (G, %) contents of the honeys influenced their C-p. The following equation (r(2) = 0.92) was proposed for estimating C-p of honey, and is recommended for use in the honey industry and in research: C = 996.7 + 1.4 x 10(-3)T + 5.6 x 10(-5)T(2) - 2.4 x 10(-7)T(3) - 56.5M - 25.8F - 31.0G + 1.5(M * F) + 1.8(M * G) + 0.8(F * G) - 4.6 x 10(-2) (M * F * G).

The effect of Cu addition and milling contaminations on the microstructure evolution of ball milled Al-Pb alloy during sintering

Al-10 wt.%Pb and Al-10 wt.%Pb-x wt.%Cu (x = 0-7.0) bulk alloys were prepared by sintering the mechanically alloyed powders at various temperatures. The microstructure changes of the as consolidated powders in the course of sintering were analyzed by differential scanning calorimetry, scanning electron microscopy, X-ray diffraction analysis and transmission electron microscopy. It has been found that, with respect to the Al-10 wt.%Pb-x wt.%Cu alloy, CuAl2 and Cu9Al4 phases formed in the milling process, and the amount of CuAl2 phase increased while the Cu9Al4 phase disappeared gradually in the sintering process. In both Al-10 wt.%Pb and Al-10 wt.%Pb-x wt.%Cu alloys, the sintering process results in the coarsening of Pb phase and the growth rate of Pb phase fulfills the Lifshitz-Slyozov-Wagner equation even though the size of the Pb phase was in nanometer range. The Pb particle exhibits cuboctahedral morphology and has a cubic to cubic orientation relationship with the Al matrix. The addition of Cu strongly depressed the growth rate of Pb. Contamination induced by milling has apparent influence on the microstructure of the sintered alloys. Al7Cu2Fe and aluminium oxide phases were identified in the sintered alloys. The cuboctahedral morphology of Pb particles was broken up by the presence of the oxide phase. (c) 2006 Elsevier B.V. All rights reserved.

Stochastic delay differential equations for genetic regulatory networks

Time delay is an important aspect in the modelling of genetic regulation due to slow biochemical reactions such as gene transcription and translation, and protein diffusion between the cytosol and nucleus. In this paper we introduce a general mathematical formalism via stochastic delay differential equations for describing time delays in genetic regulatory networks. Based on recent developments with the delay stochastic simulation algorithm, the delay chemical masterequation and the delay reaction rate equation are developed for describing biological reactions with time delay, which leads to stochastic delay differential equations derived from the Langevin approach. Two simple genetic regulatory networks are used to study the impact of' intrinsic noise on the system dynamics where there are delays. (c) 2006 Elsevier B.V. All rights reserved.