On Local Bifurcations in Neural Field Models with Transmission Delays

Neural field models with transmission delay may be cast as abstract delay differential equations (DDE). The theory of dual semigroups (also called sun-star calculus) provides a natural framework for the analysis of a broad class of delay equations, among which DDE. In particular, it may be used advantageously for the investigation of stability and bifurcation of steady states. After introducing the neural field model in its basic functional analytic setting and discussing its spectral properties, we elaborate extensively an example and derive a characteristic equation. Under certain conditions the associated equilibrium may destabilise in a Hopf bifurcation. Furthermore, two Hopf curves may intersect in a double Hopf point in a two-dimensional parameter space. We provide general formulas for the corresponding critical normal form coefficients, evaluate these numerically and interpret the results

Towards a computational model for stimulation of the Pedunculopontine nucleus

The pedunculopontine nucleus (PPN) has recently been suggested as a new therapeutic target for deep brain stimulation (DBS) in patients suffering from Parkinson's disease, particularly those with severe gait and postural impairment [1]. Stimulation at this site is typically delivered at low frequencies in contrast to the high frequency stimulation required for therapeutic benefit in the subthalamic nucleus (STN) [1]. Despite real therapeutic successes, the fundamental physiological mechanisms underlying the effect of DBS are still not understood. A hypothesis is that DBS masks the pathological synchronized firing patterns of the basal ganglia that characterize the Parkinsonian state with a regularized firing pattern. It remains unclear why stimulation of PPN should be applied with low frequency in contrast to the high frequency stimulation of STN. To get a better understanding of PPN stimulation we construct a computational model for the PPN Type I neurons in a network

A Framework for Directional and Higher-Order Reconstruction in Photoacoustic Tomography

Photoacoustic tomography is a hybrid imaging technique that combines high optical tissue contrast with high ultrasound resolution. Direct reconstruction methods such as filtered backprojection, time reversal and least squares suffer from curved line artefacts and blurring, especially in case of limited angles or strong noise. In recent years, there has been great interest in regularised iterative methods. These methods employ prior knowledge on the image to provide higher quality reconstructions. However, easy comparisons between regularisers and their properties are limited, since many tomography implementations heavily rely on the specific regulariser chosen. To overcome this bottleneck, we present a modular reconstruction framework for photoacoustic tomography. It enables easy comparisons between regularisers with different properties, e.g. nonlinear, higher-order or directional. We solve the underlying minimisation problem with an efficient first-order primal-dual algorithm. Convergence rates are optimised by choosing an operator dependent preconditioning strategy. Our reconstruction methods are tested on challenging 2D synthetic and experimental data sets. They outperform direct reconstruction approaches for strong noise levels and limited angle measurements, offering immediate benefits in terms of acquisition time and quality. This work provides a basic platform for the investigation of future advanced regularisation methods in photoacoustic tomography.Comment: submitted to "Physics in Medicine and Biology". Changes from v1 to v2: regularisation with directional wavelet has been added; new experimental tests have been include

Lumping Izhikevich neurons

We present the construction of a planar vector field that yields the firing rate of a bursting Izhikevich neuron can be read out, while leaving the sub-threshold behaviour intact. This planar vector field is used to derive lumped formulations of two complex heterogeneous networks of bursting Izhikevich neurons. In both cases, the lumped model is compared with the spiking network. There is excellent agreement in terms of duration and number of action potentials within the bursts, but there is a slight mismatch of the burst frequency. The lumped model accurately accounts for both intrinsic bursting and post inhibitory rebound potentials in the neuron model, features which are absent in prevalent neural mass models

Synchronization of the parkinsonian globus pallidus by gap junctions

We introduce pallidal gap junctional coupling as a possible mechanism for synchronization of the GPe after dopamine depletion. In a confocal imaging study, we show the presence of the neural gap junction protein Cx36 in the human GPe, including a possible remodeling process in PD patients. Dopamine has been shown to down-regulate the conductance of gap junctions in different regions of the brain [2,3], making dopamine depletion a possible candidate for increased influence of gap junctional coupling in PD

Estimation and Identifiability of Model Parameters in Human Nociceptive Processing Using Yes-No Detection Responses to Electrocutaneous Stimulation

Healthy or pathological states of nociceptive subsystems determine different stimulus-response relations measured from quantitative sensory testing. In turn, stimulus-responses measurements may be used to assess these states. In a recently developed computational model, six model parameters characterize activation of nerve endings and spinal neurons. However, both model nonlinearity and limited information in yes-no detection responses to electrocutaneous stimuli challenge to estimate model parameters. Here, we address the question whether and how one can overcome these difficulties for reliable parameter estimation. First, we fit the computational model to experimental stimulus-response pairs by maximizing the likelihood. To evaluate the balance between model fit and complexity, we evaluate the Bayesian Information Criterion. We find that the computational model is better than a conventional logistic model regarding the balance. Second, our theoretical analysis suggests to vary the pulse width among applied stimuli as a necessary condition to prevent structural non-identifiability. In addition, the numerically implemented profile likelihood approach reveals structural and practical non-identifiability. Our model-based approach with integration of psychophysical measurements can be useful for a reliable assessment of states of the nociceptive system

Neural Field Models with Transmission Delays and Diffusion

A neural field models the large scale behaviour of large groups of neurons. We extend results of van Gils et al. [2013] and Dijkstra et al. [2015] by including a diffusion term into the neural field, which models direct, electrical connections. We extend known and prove new sun-star calculus results for delay equations to be able to include diffusion and explicitly characterise the essential spectrum. For a certain class of connectivity functions in the neural field model, we are able to compute its spectral properties and the first Lyapunov coefficient of a Hopf bifurcation. By examining a numerical example, we find that the addition of diffusion suppresses non-synchronised steady-states, while favouring synchronised oscillatory modes

Comparing Epileptiform Behavior of Mesoscale Detailed Models and Population Models of Neocortex

Two models of the neocortex are developed to study normal and pathologic neuronal activity. One model contains a detailed description of a neocortical microcolumn represented by 656 neurons, including superficial and deep pyramidal cells, four types of inhibitory neurons, and realistic synaptic contacts. Simulations show that neurons of a given type exhibit similar, synchronized behavior in this detailed model. This observation is captured by a population model that describes the activity of large neuronal populations with two differential equations with two delays. Both models appear to have similar sensitivity to variations of total network excitation. Analysis of the population model reveals the presence of multistability, which was also observed in various simulations of the detailed model