Exotic dynamics in a firing rate model of neural tissue with threshold accommodation

Many of the equations describing the dynamics of neural systems are written in terms of firing rate functions, which themselves are often taken to be threshold functions of synaptic activity. Dating back to work by Hill in 1936 it has been recognized that more realistic models of neural tissue can be obtained with the introduction of state-dependent dynamic thresholds. In this paper we treat a specific phenomenological model of threshold accommodation that mimics many of the properties originally described by Hill. Importantly we explore the consequences of this dynamic threshold at the tissue level, by modifying a standard neural field model of Wilson-Cowan type. As in the case without threshold accommodation classical Mexican-Hat connectivity is shown to allow for the existence of spatially localized states (bumps) in both one and two dimensions. Importantly an analysis of bump stability in one dimension, using recent Evans function techniques, shows that bumps may undergo instabilities leading to the emergence of both breathers and traveling waves. Moreover, a similar analysis for traveling pulses leads to the conditions necessary to observe a stable traveling breather. In the regime where a bump solution does not exist direct numerical simulations show the possibility of self-replicating bumps via a form of bump splitting. Simulations in two space dimensions show analogous localized and traveling solutions to those seen in one dimension. Indeed dynamical behavior in this neural model appears reminiscent of that seen in other dissipative systems that support localized structures, and in particular those of coupled cubic complex Ginzburg-Landau equations. Further numerical explorations illustrate that the traveling pulses in this model exhibit particle like properties, similar to those of dispersive solitons observed in some three component reaction-diffusion systems. A preliminary account of this work first appeared in S Coombes and M R Owen, Bumps, breathers, and waves in a neural network with spike frequency adaptation, Physical Review Letters 94 (2005), 148102(1-4).

A Measurement of `pp��`LL near threshold

Using antiprotons from the Low Energy Antiproton Ring at CERN, an investigation of the reaction pp -+ AA at threshold has been completed. This work includes a thorough scan of the 2 MeV region where a hint of a cross section resonance had been observed. No such resonance is evident in this work. In all, the cross section is measured in 30 energy bins. These measurements range from reaction threshold to 6.0 MeV excess energy in 0.25 MeV steps. An additional six points in the 13 MeV region have also been determined. Here, these cross sections and the differential cross sections in five energy bins will be presented.

A study at high energy-resolution and at low primary energies of the secondary emission from (100) tungsten

This work was supported by the Joint Services Electronics Program (U.S. Army, U.S. Navy, and U.S. Air Force) under Contract No. DA 28 043 AMC 00073(E).

A bidomain threshold model of propagating calcium waves

We present a bidomain fire-diffuse-fire model that facilitates mathematical analysis of propagating waves of elevated intracellular calcium (Ca) in living cells. Modelling Ca release as a threshold process allows the explicit construction of travelling wave solutions to probe the dependence of Ca wave speed on physiologically important parameters such as the threshold for Ca release from the endoplasmic reticulum (ER) to the cytosol, the rate of Ca resequestration from the cytosol to the ER, and the total [Ca] (cytosolic plus ER). Interestingly, linear stability analysis of the bidomain fire-diffuse-fire model predicts the onset of dynamic wave instabilities leading to the emergence of Ca waves that propagate in a back-and-forth manner. Numerical simulations are used to confirm the presence of these so-called "tango waves" and the dependence of Ca wave speed on the total [Ca]. The original publication is available at www.springerlink.com (Journal of Mathematical Biology)

Self-exciting threshold binomial autoregressive processes

We introduce a new class of integer-valued self-exciting threshold models, which is based on the binomial autoregressive model of order one as introduced by McKenzie (Water Resour Bull 21:645–650, 1985. doi:10.1111/j.1752-1688.1985. tb05379.x). Basic probabilistic and statistical properties of this class of models are discussed. Moreover, parameter estimation and forecasting are addressed. Finally, the performance of these models is illustrated through a simulation study and an empirical application to a set of measle cases in Germany.

Super-orbital variability of LS I+61 degrees 303 at TeV energies

Context. The gamma-ray binary LS I +61º303 is a well-established source from centimeter radio up to very high energy (VHE; E > 100 GeV). The broadband emission shows a periodicity of ∼26.5 days, coincident with the orbital period. A longer (super-orbital) period of 1667 ± 8 days was proposed from radio variability and confirmed using optical and high-energy (HE; E ¿ 100 MeV) gamma-ray observations. In this paper, we report on a four-year campaign performed by MAGIC together with archival data concentrating on a search for a long-timescale signature in the VHE emission from LS I +61º303. Aims. We focus on the search for super-orbital modulation of the VHE emission, similar to that observed at other energies, and on the search for correlations between TeV emission and an optical determination of the extension of the circumstellar disk. Methods. A four-year campaign has been carried out using the MAGIC telescopes. The source was observed during the orbital phases when the periodic VHE outbursts have occurred (φ = 0.55 – 0.75, one orbit = 26.496 days). Additionally, we included archival MAGIC observations and data published by the VERITAS collaboration in these studies. For the correlation studies, LS I +61◦303 has also been observed during the orbital phases where sporadic VHE emission had been detected in the past (φ = 0.75 – 1.0). These MAGIC observations were simultaneous with optical spectroscopy from the LIVERPOOL telescope. Results. The TeV flux of the periodical outburst in orbital phases φ = 0.5 – 0.75 was found to show yearly variability consistent with the long-term modulation of ∼4.5 years found in the radio band. This modulation of the TeV flux can be well described by a sine function with a best-fit period of 1610±58 days. The complete data, including archival observations, span two super-orbital periods. There is no evidence for a correlation between the TeV emission and the mass-loss rate of the Be star, but this may be affected by the strong, short-timescale (as short as intra-day) variation displayed by the Hα fluxes.

LHC main dipole magnet circuits: sustaining near-nominal beam energies

Crossing the Franco-Swiss border, the Large Hadron Collider (LHC), designed to collide 7 TeV proton beams, is the world's largest and most powerful particle accelerator the operation of which was originally intended to commence in 2008. Unfortunately, due to an interconnect discontinuity in one of the main dipole circuit's 13 kA superconducting busbars, a catastrophic quench event occurred during initial magnet training, causing significant physical system damage. Furthermore, investigation into the cause found that such discontinuities were not only present in the circuit in question, but throughout the entire LHC. This prevented further magnet training and ultimately resulted in the maximum sustainable beam energy being limited to approximately half that of the design nominal, 3.5-4 TeV, for the first three years of operation (Run 1, 2009-2012) and a major consolidation campaign being scheduled for the first long shutdown (LS 1, 2012-2014). Throughout Run 1, a series of studies attempted to predict the amount of post-installation training quenches still required to qualify each circuit to nominal-energy current levels. With predictions in excess of 80 quenches (each having a recovery time of 8-12+ hours) just to achieve 6.5 TeV and close to 1000 quenches for 7 TeV, it was decided that for Run 2, all systems be at least qualified for 6.5 TeV operation. However, even with all interconnect discontinuities scheduled to be repaired during LS 1, numerous other concerns regarding circuit stability arose. In particular, observations of an erratic behaviour of magnet bypass diodes and the degradation of other potentially weak busbar sections, as well as observations of seemingly random millisecond spikes in beam losses, known as unidentified falling object (UFO) events, which, if persist at 6.5 TeV, may eventually deposit sufficient energy to quench adjacent magnets. In light of the above, the thesis hypothesis states that, even with the observed issues, the LHC main dipole circuits can safely support and sustain near-nominal proton beam energies of at least 6.5 TeV. Research into minimising the risk of magnet training led to the development and implementation of a new qualification method, capable of providing conclusive evidence that all aspects of all circuits, other than the magnets and their internal joints, can safely withstand a quench event at near-nominal current levels, allowing for magnet training to be carried out both systematically and without risk. This method has become known as the Copper Stabiliser Continuity Measurement (CSCM). Results were a success, with all circuits eventually being subject to a full current decay from 6.5 TeV equivalent current levels, with no measurable damage occurring. Research into UFO events led to the development of a numerical model capable of simulating typical UFO events, reproducing entire Run 1 measured event data sets and extrapolating to 6.5 TeV, predicting the likelihood of UFO-induced magnet quenches. Results provided interesting insights into the involved phenomena as well as confirming the possibility of UFO-induced magnet quenches. The model was also capable of predicting that such events, if left unaccounted for, are likely to be commonplace or not, resulting in significant long-term issues for 6.5+ TeV operation. Addressing the thesis hypothesis, the following written works detail the development and results of all CSCM qualification tests and subsequent magnet training as well as the development and simulation results of both 4 TeV and 6.5 TeV UFO event modelling. The thesis concludes, post-LS 1, with the LHC successfully sustaining 6.5 TeV proton beams, but with UFO events, as predicted, resulting in otherwise uninitiated magnet quenches and being at the forefront of system availability issues.

Instabilities in threshold-diffusion equations with delay

The introduction of delays into ordinary or partial differential equation models is well known to facilitate the production of rich dynamics ranging from periodic solutions through to spatio-temporal chaos. In this paper we consider a class of scalar partial differential equations with a delayed threshold nonlinearity which admits exact solutions for equilibria, periodic orbits and travelling waves. Importantly we show how the spectra of periodic and travelling wave solutions can be determined in terms of the zeros of a complex analytic function. Using this as a computational tool to determine stability we show that delays can have very different effects on threshold systems with negative as opposed to positive feedback. Direct numerical simulations are used to confirm our bifurcation analysis, and to probe some of the rich behaviour possible for mixed feedback.

Room temperature photo-response of titanium supersaturated silicon at energies over the bandgap

Silicon samples were implanted with high Ti doses and subsequently processed with the pulsed-laser melting technique. The electronic transport properties in the 15–300 K range and the room temperature spectral photoresponse at energies over the bandgap were measured. Samples with Ti concentration below the insulator-metal (I-M) transition limit showed a progressive reduction of the carrier lifetime in the implanted layer as Ti dose is increased. However, when the Ti concentration exceeded this limit, an extraordinary recovery of the photoresponse was measured. This result supports the theory of intermediate band materials and is of utmost relevance for photovoltaic cells and Si-based detectors.

Iterated responsive threshold search for the quadratic multiple knapsack problem

International audience

Supra-treatment threshold neonatal jaundice: Incidence in HIV-exposed compared to non-exposed neonates at Queen Elizabeth Central Hospital in Blantyre, Malawi

Introduction Jaundice is the yellowish pigmentation of the skin, sclera, and mucous membranes resulting from bilirubin deposition. Children born to mothers with HIV are more likely to be born premature, with low birth weight, and to become septic—all risk factors for neonatal jaundice. Further, there has been a change in the prevention of mother-to-child transmission (PMTCT) of HIV guidelines from single-dose nevirapine to a six-week course, all of which theoretically put HIV-exposed newborns at greater risk of developing neonatal jaundice. Aim We carried out a study to determine the incidence of severe and clinical neonatal jaundice in HIV-exposed neonates admitted to the Chatinkha Nursery (CN) neonatal unit at Queen Elizabeth Central Hospital (QECH) in Blantyre. Methods Over a period of four weeks, the incidence among non-exposed neonates was also determined for comparison between the two groups of infants. Clinical jaundice was defined as transcutaneous bilirubin levels greater than 5 mg/dL and severe jaundice as bilirubin levels above the age-specific treatment threshold according the QECH guidelines. Case notes of babies admitted were retrieved and information on birth date, gestational age, birth weight, HIV status of mother, type of feeding, mode of delivery, VDRL status of mother, serum bilirubin, duration of stay in CN, and outcome were extracted. Results Of the 149 neonates who were recruited, 17 (11.4%) were HIV-exposed. One (5.88%) of the 17 HIV-exposed and 19 (14.4%) of 132 HIVnon- exposed infants developed severe jaundice requiring therapeutic intervention (p = 0.378). Eight (47%) of the HIV-exposed and 107 (81%) of the non-exposed neonates had clinical jaundice of bilirubin levels greater than 5 mg/dL (p < 0.001). Conclusions The study showed a significant difference in the incidence of clinical jaundice between the HIV-exposed and HIV-non-exposed neonates. Contrary to our hypothesis, however, the incidence was greater in HIVnon- exposed than in HIV-exposed infants.