Dynamical and Statistical Criticality in a Model of Neural Tissue

For the nervous system to work at all, a delicate balance of excitation and inhibition must be achieved. However, when such a balance is sought by global strategies, only few modes remain balanced close to instability, and all other modes are strongly stable. Here we present a simple model of neural tissue in which this balance is sought locally by neurons following `anti-Hebbian' behavior: {\sl all} degrees of freedom achieve a close balance of excitation and inhibition and become "critical" in the dynamical sense. At long timescales, the modes of our model oscillate around the instability line, so an extremely complex "breakout" dynamics ensues in which different modes of the system oscillate between prominence and extinction. We show the system develops various anomalous statistical behaviours and hence becomes self-organized critical in the statistical sense

Dynamics of Triangulations

We study a few problems related to Markov processes of flipping triangulations of the sphere. We show that these processes are ergodic and mixing, but find a natural example which does not satisfy detailed balance. In this example, the expected distribution of the degrees of the nodes seems to follow the power law $d^{-4}$

Noise-induced memory in extended excitable systems

We describe a form of memory exhibited by extended excitable systems driven by stochastic fluctuations. Under such conditions, the system self-organizes into a state characterized by power-law correlations thus retaining long-term memory of previous states. The exponents are robust and model-independent. We discuss novel implications of these results for the functioning of cortical neurons as well as for networks of neurons.Comment: 4 pages, latex + 5 eps figure

The Cochlear Tuning Curve

The tuning curve of the cochlea measures how large an input is required to elicit a given output level as a function of the frequency. It is a fundamental object of auditory theory, for it summarizes how to infer what a sound was on the basis of the cochlear output. A simple model is presented showing that only two elements are sufficient for establishing the cochlear tuning curve: a broadly tuned traveling wave, moving unidirectionally from high to low frequencies, and a set of mechanosensors poised at the threshold of an oscillatory (Hopf) instability. These two components suffice to generate the various frequency-response regimes which are needed for a cochlear tuning curve with a high slope

Irreversible and reversible modes of operation of deterministic ratchets

We discuss a problem of optimization of the energetic efficiency of a simple rocked ratchet. We concentrate on a low-temperature case in which the particle's motion in a ratchet potential is deterministic. We show that the energetic efficiency of a ratchet working adiabatically is bounded from above by a value depending on the form of ratchet potential. The ratchets with strongly asymmetric potentials can achieve ideal efficiency of unity without approaching reversibility. On the other hand we show that for any form of the ratchet potential a set of time-protocols of the outer force exist under which the operation is reversible and the ideal value of efficiency is also achieved. The mode of operation of the ratchet is still quasistatic but not adiabatic. The high values of efficiency can be preserved even under elevated temperatures

Noise and Inertia-Induced Inhomogeneity in the Distribution of Small Particles in Fluid Flows

The dynamics of small spherical neutrally buoyant particulate impurities immersed in a two-dimensional fluid flow are known to lead to particle accumulation in the regions of the flow in which rotation dominates over shear, provided that the Stokes number of the particles is sufficiently small. If the flow is viewed as a Hamiltonian dynamical system, it can be seen that the accumulations occur in the nonchaotic parts of the phase space: the Kolmogorov--Arnold--Moser tori. This has suggested a generalization of these dynamics to Hamiltonian maps, dubbed a bailout embedding. In this paper we use a bailout embedding of the standard map to mimic the dynamics of impurities subject not only to drag but also to fluctuating forces modelled as white noise. We find that the generation of inhomogeneities associated with the separation of particle from fluid trajectories is enhanced by the presence of noise, so that they appear in much broader ranges of the Stokes number than those allowing spontaneous separation

Scalable, Time-Responsive, Digital, Energy-Efficient Molecular Circuits using DNA Strand Displacement

We propose a novel theoretical biomolecular design to implement any Boolean circuit using the mechanism of DNA strand displacement. The design is scalable: all species of DNA strands can in principle be mixed and prepared in a single test tube, rather than requiring separate purification of each species, which is a barrier to large-scale synthesis. The design is time-responsive: the concentration of output species changes in response to the concentration of input species, so that time-varying inputs may be continuously processed. The design is digital: Boolean values of wires in the circuit are represented as high or low concentrations of certain species, and we show how to construct a single-input, single-output signal restoration gate that amplifies the difference between high and low, which can be distributed to each wire in the circuit to overcome signal degradation. This means we can achieve a digital abstraction of the analog values of concentrations. Finally, the design is energy-efficient: if input species are specified ideally (meaning absolutely 0 concentration of unwanted species), then output species converge to their ideal concentrations at steady-state, and the system at steady-state is in (dynamic) equilibrium, meaning that no energy is consumed by irreversible reactions until the input again changes. Drawbacks of our design include the following. If input is provided non-ideally (small positive concentration of unwanted species), then energy must be continually expended to maintain correct output concentrations even at steady-state. In addition, our fuel species - those species that are permanently consumed in irreversible reactions - are not "generic"; each gate in the circuit is powered by its own specific type of fuel species. Hence different circuits must be powered by different types of fuel. Finally, we require input to be given according to the dual-rail convention, so that an input of 0 is specified not only by the absence of a certain species, but by the presence of another. That is, we do not construct a "true NOT gate" that sets its output to high concentration if and only if its input's concentration is low. It remains an open problem to design scalable, time-responsive, digital, energy-efficient molecular circuits that additionally solve one of these problems, or to prove that some subset of their resolutions are mutually incompatible.Comment: version 2: the paper itself is unchanged from version 1, but the arXiv software stripped some asterisk characters out of the abstract whose purpose was to highlight words. These characters have been replaced with underscores in version 2. The arXiv software also removed the second paragraph of the abstract, which has been (attempted to be) re-inserted. Also, although the secondary subject is "Soft Condensed Matter", this classification was chosen by the arXiv moderators after submission, not chosen by the authors. The authors consider this submission to be a theoretical computer science paper

Structure, Scaling and Phase Transition in the Optimal Transport Network

We minimize the dissipation rate of an electrical network under a global constraint on the sum of powers of the conductances. We construct the explicit scaling relation between currents and conductances, and show equivalence to a a previous model [J. R. Banavar {\it et al} Phys. Rev. Lett. {\bf 84}, 004745 (2000)] optimizing a power-law cost function in an abstract network. We show the currents derive from a potential, and the scaling of the conductances depends only locally on the currents. A numerical study reveals that the transition in the topology of the optimal network corresponds to a discontinuity in the slope of the power dissipation.Comment: 4 pages, 3 figure

Computational Models of Adult Neurogenesis

Experimental results in recent years have shown that adult neurogenesis is a significant phenomenon in the mammalian brain. Little is known, however, about the functional role played by the generation and destruction of neurons in the context of and adult brain. Here we propose two models where new projection neurons are incorporated. We show that in both models, using incorporation and removal of neurons as a computational tool, it is possible to achieve a higher computational efficiency that in purely static, synapse-learning driven networks. We also discuss the implication for understanding the role of adult neurogenesis in specific brain areas.Comment: To appear Physica A, 7 page