ABSTRACT
Several lines of evidence point to the modification of firing patterns and of synchronization due to gap junctions (GJs) as having a role in the establishment of epileptiform activity (EA). However, previous studies consider GJs as ohmic resistors, ignoring the effects of intense variations in ionic concentration known to occur during seizures. In addition to GJs, extracellular potassium is regarded as a further important factor involved in seizure initiation and sustainment. To analyze how these two mechanisms act together to shape firing and synchronization, we use a detailed computational model for in vitro high-K(+) and low-Ca(2+) nonsynaptic EA. The model permits us to explore the modulation of electrotonic interactions under ionic concentration changes caused by electrodiffusion in the extracellular space, altered by tortuosity. In addition, we investigate the special case of null GJ current. Increased electrotonic interaction alters bursts and action potential frequencies, favoring synchronization. The particularities of pattern changes depend on the tortuosity and array size. Extracellular potassium accumulation alone modifies firing and synchronization when the GJ coupling is null.
Subject(s)
Epilepsy/metabolism , Gap Junctions/metabolism , Potassium/metabolism , Animals , Calcium/chemistry , Calcium/metabolism , Computer Simulation , Gap Junctions/chemistry , Hippocampus/chemistry , Hippocampus/metabolism , Models, Chemical , Potassium/chemistry , RatsABSTRACT
We study distributions of dissipative and nondissipative avalanches in Manna's stochastic sandpile, in one and two dimensions. Our results lead to the following conclusions: (1) avalanche distributions, in general, do not follow simple power laws, but rather have the form P(s) approximately s(-tau(s))(ln s)(gamma)f(s/s(c)), with f a cutoff function; (2) the exponents for sizes of dissipative avalanches in two dimensions differ markedly from the corresponding values for the Bak-Tang-Wiesenfeld (BTW) model, implying that the BTW and Manna models belong to distinct universality classes; (3) dissipative avalanche distributions obey finite-size scaling, unlike in the BTW model.
ABSTRACT
We study a one-dimensional fixed-energy version (that is, with no input or loss of particles) of Manna's stochastic sandpile model. The system has a continuous transition to an absorbing state at a critical value of the particle density, and exhibits the hallmarks of an absorbing-state phase transition, including finite-size scaling. Critical exponents are obtained from extensive simulations, which treat stationary and transient properties, and an associated interface representation. These exponents characterize the universality class of an absorbing-state phase transition with a static conserved density in one dimension; they differ from those expected at a linear-interface depinning transition in a medium with point disorder, and from those of directed percolation.
ABSTRACT
We study a one-dimensional lattice random walk with an absorbing boundary at the origin and a movable partial reflector. On encountering the reflector at site x, the walker is reflected (with probability r) to x-1 and the reflector is simultaneously pushed to x+1. Iteration of the transition matrix, and asymptotic analysis of the probability generating function show that the critical exponent delta governing the survival probability varies continuously between 1/2 and 1 as r varies between 0 and 1. Our study suggests a mechanism for nonuniversal kinetic critical behavior, observed in models with an infinite number of absorbing configurations.
ABSTRACT
A lattice gas with infinite repulsion between particles separated by < or = 1 lattice spacing, and nearest-neighbor hopping dynamics, is subject to a drive favoring movement along one axis of the square lattice. The equilibrium (zero drive) transition to a phase with sublattice ordering, known to be continuous, shifts to lower density, and becomes discontinuous for large bias. In the ordered nonequilibrium steady state, both the particle and order-parameter densities are nonuniform, with a large fraction of the particles occupying a jammed strip oriented along the drive. The drive thus induces separation into high- and low-density regions in a system with purely repulsive interactions. Increasing the drive can provoke a transition to the ordered phase, and thereby, a sharp reduction in current.
ABSTRACT
A correspondence between lattice models with absorbing states and models of pinned interfaces in random media can be established by defining local height variables h(x,t) as integrals of the activity at point x up to time t. Within this context we study the interface representation of a prototypical model with absorbing states, the contact process, in dimensions 1-3. Simulations confirm the scaling relation beta(W)=1-straight theta between the interface-width growth exponent beta(W) and the exponent straight theta governing the decay of the order parameter. A scaling property of the height distribution, which serves as the basis for this relation, is also verified. The height-height correlation function shows clear signs of anomalous scaling, in accord with Lopez' analysis [Phys. Rev. Lett. 83, 4594 (1999)], but no evidence of multiscaling.
ABSTRACT
We derive a pair approximation (PA) for the NO+CO model with instantaneous reactions. For both the triangular and square lattices, the PA, derived here using a simpler approach, yields a phase diagram with an active state for CO-fractions y in the interval y(1)
ABSTRACT
A simple reweighting scheme is proposed for Monte Carlo simulations of interacting particle systems, permitting one to study various parameter values in a single study, and improving efficiency by an order of magnitude. Unlike earlier reweighting schemes, the present approach does not require knowledge of the stationary probability distribution, and so is applicable out of equilibrium. The method is applied to the contact process in two and three dimensions, yielding the critical parameter and spreading exponents to unprecedented precision.
ABSTRACT
We study the stationary properties of the two-dimensional pair contact process, a nonequilibrium lattice model exhibiting a phase transition to an absorbing state with an infinite number of configurations. The critical probability and static critical exponents are determined via Monte Carlo simulations, as well as order-parameter moment ratios and the scaling of the initial density decay. The static critical properties fall in the directed percolation universality class.