Numerical computation of an Evans function for travelling waves

We demonstrate a geometrically inspired technique for computing Evans functions for the linearised operators about travelling waves. Using the examples of the F-KPP equation and a Keller-Segel model of bacterial chemotaxis, we produce an Evans function which is computable through several orders of magnitude in the spectral parameter and show how such a function can naturally be extended into the continuous spectrum. In both examples, we use this function to numerically verify the absence of eigenvalues in a large region of the right half of the spectral plane. We also include a new proof of spectral stability in the appropriate weighted space of travelling waves of speed $c \geq 2 \sqrt{\delta}$ in the F-KPP equation.Comment: 37 pages, 11 figure

Application of the control volume method to a mathematical model of cell migration

In recent years, mathematical models of cell migration have become increasingly complex. These models have evolved from simple diffusion models to computationally troublesome reaction-diffusion-advection models. As such, the use of ``black box'' numerical solvers has become less appropriate. We discuss the application of the control volume technique for resolving a complicated nonlinear cell migration model. The nonlinearity is treated using an inexact Newton solver and flux limiting ensures that the cell migration fronts are captured adequately. Specifically, we analyse the model due to Perumpanani et al. (1999), comparing the numerical results of the proposed computational model developed in this research to previous results published by other researchers. We show that the finite volume computational model captures the physics of the processes with good accuracy using coarse meshes

Quiescience as a mechanism for cyclical hypoxia and acidosis

Tumour tissue characteristically experiences fluctuations in substrate supply. This unstable microenvironment drives constitutive metabolic changes within cellular populations and, ultimately, leads to a more aggressive phenotype. Previously, variations in substrate levels were assumed to occur through oscillations in the hæmodynamics of nearby and distant blood vessels. In this paper we examine an alternative hypothesis, that cycles of metabolite concentrations are also driven by cycles of cellular quiescence and proliferation. Using a mathematical modelling approach, we show that the interdependence between cell cycle and the microenvironment will induce typical cycles with the period of order hours in tumour acidity and oxygenation. As a corollary, this means that the standard assumption of metabolites entering diffusive equilibrium around the tumour is not valid; instead temporal dynamics must be considered

A reaction-diffusion model for the growth of avascular tumor

A nutrient-limited model for avascular cancer growth including cell proliferation, motility and death is presented. The model qualitatively reproduces commonly observed morphologies for primary tumors, and the simulated patterns are characterized by its gyration radius, total number of cancer cells, and number of cells on tumor periphery. These very distinct morphological patterns follow Gompertz growth curves, but exhibit different scaling laws for their surfaces. Also, the simulated tumors incorporate a spatial structure composed of a central necrotic core, an inner rim of quiescent cells and a narrow outer shell of proliferating cells in agreement with biological data. Finally, our results indicate that the competition for nutrients among normal and cancer cells may be a determinant factor in generating papillary tumor morphology.Comment: 9 pages, 6 figures, to appear in PR

A Three Species Model to Simulate Application of Hyperbaric Oxygen Therapy to Chronic Wounds

Chronic wounds are a significant socioeconomic problem for governments worldwide. Approximately 15% of people who suffer from diabetes will experience a lower-limb ulcer at some stage of their lives, and 24% of these wounds will ultimately result in amputation of the lower limb. Hyperbaric Oxygen Therapy (HBOT) has been shown to aid the healing of chronic wounds; however, the causal reasons for the improved healing remain unclear and hence current HBOT protocols remain empirical. Here we develop a three-species mathematical model of wound healing that is used to simulate the application of hyperbaric oxygen therapy in the treatment of wounds. Based on our modelling, we predict that intermittent HBOT will assist chronic wound healing while normobaric oxygen is ineffective in treating such wounds. Furthermore, treatment should continue until healing is complete, and HBOT will not stimulate healing under all circumstances, leading us to conclude that finding the right protocol for an individual patient is crucial if HBOT is to be effective. We provide constraints that depend on the model parameters for the range of HBOT protocols that will stimulate healing. More specifically, we predict that patients with a poor arterial supply of oxygen, high consumption of oxygen by the wound tissue, chronically hypoxic wounds, and/or a dysfunctional endothelial cell response to oxygen are at risk of nonresponsiveness to HBOT. The work of this paper can, in some way, highlight which patients are most likely to respond well to HBOT (for example, those with a good arterial supply), and thus has the potential to assist in improving both the success rate and hence the cost-effectiveness of this therapy

Population Receptive Field Dynamics in Human Visual Cortex

Seminal work in the early nineties revealed that the visual receptive field of neurons in cat primary visual cortex can change in location and size when artificial scotomas are applied. Recent work now suggests that these single neuron receptive field dynamics also pertain to the neuronal population receptive field (pRF) that can be measured in humans with functional magnetic resonance imaging (fMRI). To examine this further, we estimated the pRF in twelve healthy participants while masking the central portion of the visual field. We found that the pRF changes in location and size for two differently sized artificial scotomas, and that these pRF dynamics are most likely due to a combination of the neuronal receptive field position and size scatter as well as modulatory feedback signals from extrastriate visual areas

Cortical Surround Interactions and Perceptual Salience via Natural Scene Statistics

Spatial context in images induces perceptual phenomena associated with salience and modulates the responses of neurons in primary visual cortex (V1). However, the computational and ecological principles underlying contextual effects are incompletely understood. We introduce a model of natural images that includes grouping and segmentation of neighboring features based on their joint statistics, and we interpret the firing rates of V1 neurons as performing optimal recognition in this model. We show that this leads to a substantial generalization of divisive normalization, a computation that is ubiquitous in many neural areas and systems. A main novelty in our model is that the influence of the context on a target stimulus is determined by their degree of statistical dependence. We optimized the parameters of the model on natural image patches, and then simulated neural and perceptual responses on stimuli used in classical experiments. The model reproduces some rich and complex response patterns observed in V1, such as the contrast dependence, orientation tuning and spatial asymmetry of surround suppression, while also allowing for surround facilitation under conditions of weak stimulation. It also mimics the perceptual salience produced by simple displays, and leads to readily testable predictions. Our results provide a principled account of orientation-based contextual modulation in early vision and its sensitivity to the homogeneity and spatial arrangement of inputs, and lends statistical support to the theory that V1 computes visual salience

Model Cortical Association Fields Account for the Time Course and Dependence on Target Complexity of Human Contour Perception

Can lateral connectivity in the primary visual cortex account for the time dependence and intrinsic task difficulty of human contour detection? To answer this question, we created a synthetic image set that prevents sole reliance on either low-level visual features or high-level context for the detection of target objects. Rendered images consist of smoothly varying, globally aligned contour fragments (amoebas) distributed among groups of randomly rotated fragments (clutter). The time course and accuracy of amoeba detection by humans was measured using a two-alternative forced choice protocol with self-reported confidence and variable image presentation time (20-200 ms), followed by an image mask optimized so as to interrupt visual processing. Measured psychometric functions were well fit by sigmoidal functions with exponential time constants of 30-91 ms, depending on amoeba complexity. Key aspects of the psychophysical experiments were accounted for by a computational network model, in which simulated responses across retinotopic arrays of orientation-selective elements were modulated by cortical association fields, represented as multiplicative kernels computed from the differences in pairwise edge statistics between target and distractor images. Comparing the experimental and the computational results suggests that each iteration of the lateral interactions takes at least ms of cortical processing time. Our results provide evidence that cortical association fields between orientation selective elements in early visual areas can account for important temporal and task-dependent aspects of the psychometric curves characterizing human contour perception, with the remaining discrepancies postulated to arise from the influence of higher cortical areas