Modeling electrocortical activity through improved local approximations of integral neural field equations

Neural field models of firing rate activity typically take the form of integral equations with space-dependent axonal delays. Under natural assumptions on the synaptic connectivity we show how one can derive an equivalent partial differential equation (PDE) model that properly treats the axonal delay terms of the integral formulation. Our analysis avoids the so-called long-wavelength approximation that has previously been used to formulate PDE models for neural activity in two spatial dimensions. Direct numerical simulations of this PDE model show instabilities of the homogeneous steady state that are in full agreement with a Turing instability analysis of the original integral model. We discuss the benefits of such a local model and its usefulness in modeling electrocortical activity. In particular we are able to treat "patchy'" connections, whereby a homogeneous and isotropic system is modulated in a spatially periodic fashion. In this case the emergence of a "lattice-directed" traveling wave predicted by a linear instability analysis is confirmed by the numerical simulation of an appropriate set of coupled PDEs. Article published and (c) American Physical Society 2007

The Becker-Döring equations with monomer input, competition and inhibition

We investigate the Becker-Döring model of nucleation with three generalisations; an input of monomer, an input of inhibitor and finally, we allow the monomers to form two morphologies of cluster. We assume size-independent aggregation and fragmentation rates. Initially we consider the problem of constant monomer input and determine the steady-state solution approached in the large-time limit, and the manner in which it is approached. Secondly, in addition to a constant input of monomer we allow a constant input of inhibitor, which prevents clusters growing any larger and this removes them from the kinetics of the process; the inhibitor is consumed in the action of poisoning a cluster. We determine a critical ratio of poison to monomer input below which the cluster concentrations tend to a non-zero steady-state solution and the poison concentration tends to a finite value. Above the critical input ratio, the concentrations of all cluster sizes tend to zero and the poison concentration grows without limit. In both cases the solution in the large-time limit is determined. Finally we consider a model where monomers form two morphologies, but the inhibitor only acts on one morphology. Four cases are identified, depending on the relative poison to monomer input rates and the relative thermodynamic stability. In each case we determine the final cluster distribution and poison concentration. We find that poisoning the less stable cluster type can have a significant impact on the structure of the more stable cluster distribution; a counter-intuitive result. All results are shown to agree with numerical simulation.

Nonlinear development of matrix-converter instabilities

Matrix converters convert a three-phase alternating-current power supply to a power supply of a different peak voltage and frequency, and are an emerging technology in a wide variety of applications. However, they are susceptible to an instability, whose behaviour is examined herein. The desired “steady-state” mode of operation of the matrix converter becomes unstable in a Hopf bifurcation as the output/input voltage transfer ratio, q, is increased through some threshold value, qc. Through weakly nonlinear analysis and direct numerical simulation of an averaged model, we show that this bifurcation is subcritical for typical parameter values, leading to hysteresis in the transition to the oscillatory state: there may thus be undesirable large-amplitude oscillations in the output voltages even when q is below the linear stability threshold value qc.