How the asymmetry of internal potential influences the shape of I-V characteristic of nanochannels

Ion transport in biological and synthetic nanochannels is characterized by such phenomena as ion current fluctuations, rectification, and pumping. Recently, it has been shown that the nanofabricated synthetic pores could be considered as analogous to biological channels with respect to their transport characteristics \cite{Apel, Siwy}. The ion current rectification is analyzed. Ion transport through cylindrical nanopores is described by the Smoluchowski equation. The model is considering the symmetric nanopore with asymmetric charge distribution. In this model, the current rectification in asymmetrically charged nanochannels shows a diode-like shape of $I-V$ characteristic. It is shown that this feature may be induced by the coupling between the degree of asymmetry and the depth of internal electric potential well. The role of concentration gradient is discussed

Response to sub-threshold stimulus is enhanced by spatially heterogeneous activity

Sub-threshold stimuli cannot initiate excitations in active media, but surprisingly as we show in this paper, they can alter the time-evolution of spatially heterogeneous activity by modifying the recovery dynamics. This results in significant reduction of waveback velocity which may lead to spatial coherence, terminating all activity in the medium including spatiotemporal chaos. We analytically derive model-independent conditions for which such behavior can be observed.Comment: 5 pages, 5 figure

Convective Fingering of an Autocatalytic Reaction Front

We report experimental observations of the convection-driven fingering instability of an iodate-arsenous acid chemical reaction front. The front propagated upward in a vertical slab; the thickness of the slab was varied to control the degree of instability. We observed the onset and subsequent nonlinear evolution of the fingers, which were made visible by a {\it p}H indicator. We measured the spacing of the fingers during their initial stages and compared this to the wavelength of the fastest growing linear mode predicted by the stability analysis of Huang {\it et. al.} [{\it Phys. Rev. E}, {\bf 48}, 4378 (1993), and unpublished]. We find agreement with the thickness dependence predicted by the theory.Comment: 11 pages, RevTex with 3 eps figures. To be published in Phys Rev E, [email protected], [email protected], [email protected]

Avalanche of Bifurcations and Hysteresis in a Model of Cellular Differentiation

Cellular differentiation in a developping organism is studied via a discrete bistable reaction-diffusion model. A system of undifferentiated cells is allowed to receive an inductive signal emenating from its environment. Depending on the form of the nonlinear reaction kinetics, this signal can trigger a series of bifurcations in the system. Differentiation starts at the surface where the signal is received, and cells change type up to a given distance, or under other conditions, the differentiation process propagates through the whole domain. When the signal diminishes hysteresis is observed

Dynamics of lattice spins as a model of arrhythmia

We consider evolution of initial disturbances in spatially extended systems with autonomous rhythmic activity, such as the heart. We consider the case when the activity is stable with respect to very smooth (changing little across the medium) disturbances and construct lattice models for description of not-so-smooth disturbances, in particular, topological defects; these models are modifications of the diffusive XY model. We find that when the activity on each lattice site is very rigid in maintaining its form, the topological defects - vortices or spirals - nucleate a transition to a disordered, turbulent state.Comment: 17 pages, revtex, 3 figure

Renormalization Group Theory for a Perturbed KdV Equation

We show that renormalization group(RG) theory can be used to give an analytic description of the evolution of a perturbed KdV equation. The equations describing the deformation of its shape as the effect of perturbation are RG equations. The RG approach may be simpler than inverse scattering theory(IST) and another approaches, because it dose not rely on any knowledge of IST and it is very concise and easy to understand. To the best of our knowledge, this is the first time that RG has been used in this way for the perturbed soliton dynamics.Comment: 4 pages, no figure, revte

Incorporating Inductances in Tissue-Scale Models of Cardiac Electrophysiology

In standard models of cardiac electrophysiology, including the bidomain and monodomain models, local perturbations can propagate at infinite speed. We address this unrealistic property by developing a hyperbolic bidomain model that is based on a generalization of Ohm's law with a Cattaneo-type model for the fluxes. Further, we obtain a hyperbolic monodomain model in the case that the intracellular and extracellular conductivity tensors have the same anisotropy ratio. In one spatial dimension, the hyperbolic monodomain model is equivalent to a cable model that includes axial inductances, and the relaxation times of the Cattaneo fluxes are strictly related to these inductances. A purely linear analysis shows that the inductances are negligible, but models of cardiac electrophysiology are highly nonlinear, and linear predictions may not capture the fully nonlinear dynamics. In fact, contrary to the linear analysis, we show that for simple nonlinear ionic models, an increase in conduction velocity is obtained for small and moderate values of the relaxation time. A similar behavior is also demonstrated with biophysically detailed ionic models. Using the Fenton-Karma model along with a low-order finite element spatial discretization, we numerically analyze differences between the standard monodomain model and the hyperbolic monodomain model. In a simple benchmark test, we show that the propagation of the action potential is strongly influenced by the alignment of the fibers with respect to the mesh in both the parabolic and hyperbolic models when using relatively coarse spatial discretizations. Accurate predictions of the conduction velocity require computational mesh spacings on the order of a single cardiac cell. We also compare the two formulations in the case of spiral break up and atrial fibrillation in an anatomically detailed model of the left atrium, and [...].Comment: 20 pages, 12 figure

Limit theorems for weakly subcritical branching processes in random environment

For a branching process in random environment it is assumed that the offspring distribution of the individuals varies in a random fashion, independently from one generation to the other. Interestingly there is the possibility that the process may at the same time be subcritical and, conditioned on nonextinction, 'supercritical'. This so-called weakly subcritical case is considered in this paper. We study the asymptotic survival probability and the size of the population conditioned on non-extinction. Also a functional limit theorem is proven, which makes the conditional supercriticality manifest. A main tool is a new type of functional limit theorems for conditional random walks.Comment: 35 page

ATLAS Inner Detector Event Data Model

The data model for event reconstruction (EDM) in the Inner Detector of the ATLAS experiment is presented. Different data classes represent evolving stages in the reconstruction data flow, and specific derived classes exist for the sub-detectors. The Inner Detector EDM also extends the data model for common tracking in ATLAS and is integrated into the modular design of the ATLAS high-level trigger and off-line software

Finite-gap Solutions of the Vortex Filament Equation: Isoperiodic Deformations

We study the topology of quasiperiodic solutions of the vortex filament equation in a neighborhood of multiply covered circles. We construct these solutions by means of a sequence of isoperiodic deformations, at each step of which a real double point is "unpinched" to produce a new pair of branch points and therefore a solution of higher genus. We prove that every step in this process corresponds to a cabling operation on the previous curve, and we provide a labelling scheme that matches the deformation data with the knot type of the resulting filament.Comment: 33 pages, 5 figures; submitted to Journal of Nonlinear Scienc