A Model of Heat Conduction

We define a deterministic ``scattering'' model for heat conduction which is continuous in space, and which has a Boltzmann type flavor, obtained by a closure based on memory loss between collisions. We prove that this model has, for stochastic driving forces at the boundary, close to Maxwellians, a unique non-equilibrium steady state

Period Doubling Renormalization for Area-Preserving Maps and Mild Computer Assistance in Contraction Mapping Principle

It has been observed that the famous Feigenbaum-Coullet-Tresser period doubling universality has a counterpart for area-preserving maps of {\fR}^2. A renormalization approach has been used in a "hard" computer-assisted proof of existence of an area-preserving map with orbits of all binary periods in Eckmann et al (1984). As it is the case with all non-trivial universality problems in non-dissipative systems in dimensions more than one, no analytic proof of this period doubling universality exists to date. In this paper we attempt to reduce computer assistance in the argument, and present a mild computer aided proof of the analyticity and compactness of the renormalization operator in a neighborhood of a renormalization fixed point: that is a proof that does not use generalizations of interval arithmetics to functional spaces - but rather relies on interval arithmetics on real numbers only to estimate otherwise explicit expressions. The proof relies on several instance of the Contraction Mapping Principle, which is, again, verified via mild computer assistance

Memory Effects in Nonequilibrium Transport for Deterministic Hamiltonian Systems

We consider nonequilibrium transport in a simple chain of identical mechanical cells in which particles move around. In each cell, there is a rotating disc, with which these particles interact, and this is the only interaction in the model. It was shown in \cite{eckmann-young} that when the cells are weakly coupled, to a good approximation, the jump rates of particles and the energy-exchange rates from cell to cell follow linear profiles. Here, we refine that study by analyzing higher-order effects which are induced by the presence of external gradients for situations in which memory effects, typical of Hamiltonian dynamics, cannot be neglected. For the steady state we propose a set of balance equations for the particle number and energy in terms of the reflection probabilities of the cell and solve it phenomenologically. Using this approximate theory we explain how these asymmetries affect various aspects of heat and particle transport in systems of the general type described above and obtain in the infinite volume limit the deviation from the theory in \cite{eckmann-young} to first-order. We verify our assumptions with extensive numerical simulations.Comment: Several change

Leaders of neuronal cultures in a quorum percolation model

We present a theoretical framework using quorum-percolation for describing the initiation of activity in a neural culture. The cultures are modeled as random graphs, whose nodes are excitatory neurons with kin inputs and kout outputs, and whose input degrees kin = k obey given distribution functions pk. We examine the firing activity of the population of neurons according to their input degree (k) classes and calculate for each class its firing probability \Phi_k(t) as a function of t. The probability of a node to fire is found to be determined by its in-degree k, and the first-to-fire neurons are those that have a high k. A small minority of high-k classes may be called "Leaders", as they form an inter-connected subnetwork that consistently fires much before the rest of the culture. Once initiated, the activity spreads from the Leaders to the less connected majority of the culture. We then use the distribution of in-degree of the Leaders to study the growth rate of the number of neurons active in a burst, which was experimentally measured to be initially exponential. We find that this kind of growth rate is best described by a population that has an in-degree distribution that is a Gaussian centered around k = 75 with width {\sigma} = 31 for the majority of the neurons, but also has a power law tail with exponent -2 for ten percent of the population. Neurons in the tail may have as many as k = 4, 700 inputs. We explore and discuss the correspondence between the degree distribution and a dynamic neuronal threshold, showing that from the functional point of view, structure and elementary dynamics are interchangeable. We discuss possible geometric origins of this distribution, and comment on the importance of size, or of having a large number of neurons, in the culture.Comment: Keywords: Neuronal cultures, Graph theory, Activation dynamics, Percolation, Statistical mechanics of networks, Leaders of activity, Quorum. http://www.weizmann.ac.il/complex/tlusty/papers/FrontCompNeuro2010.pd

Spectral Properties of Hypoelliptic Operators

We study hypoelliptic operators with polynomially bounded coefficients that are of the form K = sum_{i=1}^m X_i^T X_i + X_0 + f, where the X_j denote first order differential operators, f is a function with at most polynomial growth, and X_i^T denotes the formal adjoint of X_i in L^2. For any e > 0 we show that an inequality of the form |u|_{delta,delta} <= C(|u|_{0,eps} + |(K+iy)u|_{0,0}) holds for suitable delta and C which are independent of y in R, in weighted Sobolev spaces (the first index is the derivative, and the second the growth). We apply this result to the Fokker-Planck operator for an anharmonic chain of oscillators coupled to two heat baths. Using a method of Herau and Nier [HN02], we conclude that its spectrum lies in a cusp {x+iy|x >= |y|^tau-c, tau in (0,1], c in R}.Comment: 3 figure

Recurrence quantification analysis as a tool for the characterization of molecular dynamics simulations

A molecular dynamics simulation of a Lennard-Jones fluid, and a trajectory of the B1 immunoglobulin G-binding domain of streptococcal protein G (B1-IgG) simulated in water are analyzed by recurrence quantification, which is noteworthy for its independence from stationarity constraints, as well as its ability to detect transients, and both linear and nonlinear state changes. The results demonstrate the sensitivity of the technique for the discrimination of phase sensitive dynamics. Physical interpretation of the recurrence measures is also discussed.Comment: 7 pages, 8 figures, revtex; revised for review for Phys. Rev. E (clarifications and expansion of discussion)-- addition of the 8 postscript figures previously omitted, but unchanged from version

Collet, Eckmann and the bifurcation measure

The moduli space $\mathcal{M}_d$ of degree $d\geq2$ rational maps can naturally be endowed with a measure $\mu_\mathrm{bif}$ detecting maximal bifurcations, called the bifurcation measure. We prove that the support of the bifurcation measure $\mu_\mathrm{bif}$ has positive Lebesgue measure. To do so, we establish a general sufficient condition for the conjugacy class of a rational map to belong to the support of $\mu_\mathrm{bif}$ and we exhibit a large set of Collet-Eckmann rational maps which satisfy this condition. As a consequence, we get a set of Collet-Eckmann rational maps of positive Lebesgue measure which are approximated by hyperbolic rational maps