Nonlinear diffusion from Einstein's master equation

We generalize Einstein's master equation for random walk processes by considering that the probability for a particle at position $r$ to make a jump of length $j$ lattice sites, $P_j(r)$ is a functional of the particle distribution function $f(r,t)$. By multiscale expansion, we obtain a generalized advection-diffusion equation. We show that the power law $P_j(r) \propto f(r)^{\alpha - 1}$ (with $\alpha > 1$) follows from the requirement that the generalized equation admits of scaling solutions ($f(r;t) = t^{-\gamma}\phi (r/t^{\gamma})$). The solutions have a $q$-exponential form and are found to be in agreement with the results of Monte-Carlo simulations, so providing a microscopic basis validating the nonlinear diffusion equation. Although its hydrodynamic limit is equivalent to the phenomenological porous media equation, there are extra terms which, in general, cannot be neglected as evidenced by the Monte-Carlo computations.}Comment: 7 pages incl. 3 fig

Lattice gas with ``interaction potential''

We present an extension of a simple automaton model to incorporate non-local interactions extending over a spatial range in lattice gases. {}From the viewpoint of Statistical Mechanics, the lattice gas with interaction range may serve as a prototype for non-ideal gas behavior. {}From the density fluctuations correlation function, we obtain a quantity which is identified as a potential of mean force. Equilibrium and transport properties are computed theoretically and by numerical simulations to establish the validity of the model at macroscopic scale.Comment: 12 pages LaTeX, figures available on demand ([email protected]

Nonextensive diffusion as nonlinear response

The porous media equation has been proposed as a phenomenological ``non-extensive'' generalization of classical diffusion. Here, we show that a very similar equation can be derived, in a systematic manner, for a classical fluid by assuming nonlinear response, i.e. that the diffusive flux depends on gradients of a power of the concentration. The present equation distinguishes from the porous media equation in that it describes \emph{% generalized classical} diffusion, i.e. with $r/\sqrt Dt$ scaling, but with a generalized Einstein relation, and with power-law probability distributions typical of nonextensive statistical mechanics

Heavy Quark Diffusion from the Lattice

We study the diffusion of heavy quarks in the Quark Gluon Plasma using the Langevin equations of motion and estimate the contribution of the transport peak to the Euclidean current-current correlator. We show that the Euclidean correlator is remarkably insensitive to the heavy quark diffusion coefficient and give a simple physical interpretation of this result using the free streaming Boltzmann equation. However if the diffusion coefficient is smaller than $\sim 1/(\pi T)$, as favored by RHIC phenomenology, the transport contribution should be visible in the Euclidean correlator. We outline a procedure to isolate this contribution.Comment: 24 pages, 5 figure

Metabolomics on integrated circuit

We have demonstrated a chip-based diagnostics tool for the quantification of metabolites, using specific enzymes, to study enzyme kinetics and calculate the Michaelis-Menten constant. An array of 256×256 ion-sensitive field effect transistors (ISFETs) fabricated in a complementary metal oxide semiconductor (CMOS) process is used for this prototype. We have used hexokinase enzyme reaction on the ISFET CMOS chip with glucose concentration in the physiological range of 0.05 mM – 231 mM and successfully studied the enzyme kinetics of hexokinase in detail. This will promote future research towards multiplexing enzyme-based metabolite quantification on a single chip, ultimately opening a pathway towards a personal metabolome machine

Macroscopic evidence of microscopic dynamics in the Fermi-Pasta-Ulam oscillator chain from nonlinear time series analysis

The problem of detecting specific features of microscopic dynamics in the macroscopic behavior of a many-degrees-of-freedom system is investigated by analyzing the position and momentum time series of a heavy impurity embedded in a chain of nearest-neighbor anharmonic Fermi-Pasta-Ulam oscillators. Results obtained in a previous work [M. Romero-Bastida, Phys. Rev. E {\bf69}, 056204 (2004)] suggest that the impurity does not contribute significantly to the dynamics of the chain and can be considered as a probe for the dynamics of the system to which the impurity is coupled. The ($r,\tau$) entropy, which measures the amount of information generated by unit time at different scales $\tau$ of time and $r$ of the observable, is numerically computed by methods of nonlinear time-series analysis using the position and momentum signals of the heavy impurity for various values of the energy density $\epsilon$ (energy per degree of freedom) of the system and some values of the impurity mass $M$. Results obtained from these two time series are compared and discussed.Comment: 7 pages, 5 figures, RevTeX4 PRE format; to be published in Phys. Rev.

Statistics of precursors to fingering processes

We present an analysis of the statistical properties of hydrodynamic field fluctuations which reveal the existence of precursors to fingering processes. These precursors are found to exhibit power law distributions, and these power laws are shown to follow from spatial $q$-Gaussian structures which are solutions to the generalized non-linear diffusion equation.Comment: 7 pages incl. 5 figs; tp appear in Europhysics Letter

Computer Simulation Study of the Phase Behavior and Structural Relaxation in a Gel-Former Modeled by Three Body Interactions

We report a computer simulation study of a model gel-former obtained by modifying the three-body interactions of the Stillinger-Weber potential for silicon. This modification reduces the average coordination number and consequently shifts the liquid-gas phase coexistence curve to low densities, thus facilitating the formation of gels without phase separation. At low temperatures and densities, the structure of the system is characterized by the presence of long linear chains interconnected by a small number of three coordinated junctions at random locations. At small wave-vectors the static structure factor shows a non-monotonic dependence on temperature, a behavior which is due to the competition between the percolation transition of the particles and the stiffening of the formed chains. We compare in detail the relaxation dynamics of the system as obtained from molecular dynamics with the one obtained from Monte Carlo dynamics. We find that the bond correlation function displays stretched exponential behavior at moderately low temperatures and densities, but exponential relaxation at low temperatures. The bond lifetime shows an Arrhenius behavior, independent of the microscopic dynamics. For the molecular dynamics at low temperatures, the mean squared displacement and the (coherent and incoherent) intermediate scattering function display at intermediate times a dynamics with ballistic character and we show that this leads to compressed exponential relaxation. For the Monte Carlo dynamics we find always an exponential or stretched exponential relaxation. Thus we conclude that the compressed exponential relaxation observed in experiments is due to the out-of-equilibrium dynamics