Diffusion on a solid surface: anomalous is normal

We present a numerical study of classical particles diffusing on a solid surface. The particles motion is modeled by an underdamped Langevin equation with ordinary thermal noise. The particle-surface interaction is described by a periodic or a random two-dimensional potential. The model leads to a rich variety of different transport regimes, some of which correspond to anomalous diffusion such as has recently been observed in experiments and Monte Carlo simulations. We show that this anomalous behavior is controlled by the friction coefficient and stress that it emerges naturally in a system described by ordinary canonical Maxwell-Boltzmann statistics.

Sancho, et al. Reply

A Reply to the Comment by Jing-Dong Bao and Yan Zhou.

Dispersionless transport in a washboard potential

We study and characterize a new dynamical regime of underdamped particles in a tilted washboard potential. We find that for small friction in a finite range of forces the particles move essentially nondispersively, that is, coherently, over long intervals of time. The associated distribution of the particle positions moves at an essentially constant velocity and is far from Gaussian-like. This new regime is complementary to, and entirely different from, well-known nonlinear response and large dispersion regimes observed for other values of the external force.

Self-sustained spatiotemporal oscillations induced by membrane-bulk coupling

We propose a novel mechanism leading to spatiotemporal oscillations in extended systems that does not rely on local bulk instabilities. Instead, oscillations arise from the interaction of two subsystems of different spatial dimensionality. Specifically, we show that coupling a passive diffusive bulk of dimension d with an excitable membrane of dimension d-1 produces a self-sustained oscillatory behavior. An analytical explanation of the phenomenon is provided for d=1. Moreover, in-phase and antiphase synchronization of oscillations are found numerically in one and two dimensions. This novel dynamic instability could be used by biological systems such as cells, where the dynamics on the cellular membrane is necessarily different from that of the cytoplasmic bulk.

Evaporation and coarsening dynamics with open boundaries

We present a study of the evaporation dynamics of a substance undergoing a coarsening process. The system is modeled by the Cahn-Hilliard equation with absorbing boundaries. We have found that the dynamics, although of a diffusive nature, is much slower than the usual one without coarsening. Analytical and simulation results are in reasonable agreement.

Symmetric Brownian motor

In this paper, we present a model of a symmetric Brownian motor which changes the sign of its velocity when the temperature gradient is inverted. The velocity, external work, and efficiency are studied as a function of the temperatures of the baths and other relevant parameters. The motor shows a current reversal when another parameter (a phase shift) is varied. Analytical predictions and results from numerical simulations are performed and agree very well. Generic properties of this type of motor are discussed.

From subdiffusion to superdiffusion of particles on solid surfaces

We present a numerical and partially analytical study of classical particles obeying a Langevin equation that describes diffusion on a surface modeled by a two-dimensional potential. The potential may be either periodic or random. Depending on the potential and the damping, we observe superdiffusion, large-step diffusion, diffusion, and subdiffusion. Superdiffusive behavior is associated with low damping and is in most cases transient, albeit often long. Subdiffusive behavior is associated with highly damped particles in random potentials. In some cases subdiffusive behavior persists over our entire simulation and may be characterized as metastable. In any case, we stress that this rich variety of behaviors emerges naturally from an ordinary Langevin equation for a system described by ordinary canonical Maxwell-Boltzmann statistics.

Competitive evaporation in arrays of droplets

We consider the evaporation of periodic arrays of initially equal droplets in two-dimensional systems with open (absorbing) boundaries. Our study is based on the numerical solution of the Cahn-Hilliard equation. We show that due to cooperative effects the droplets which are further from the boundary may evaporate earlier than those in the boundary¿s vicinity. The time evolution of the overall amount of matter in the system is also studied.

Particle dispersion in synthetic turbulent flows

We study particle dispersion advected by a synthetic turbulent flow from a Lagrangian perspective and focus on the two-particle and cluster dispersion by the flow. It has been recently reported that Richardson¿s law for the two-particle dispersion can stem from different dispersion mechanisms, and can be dominated by either diffusive or ballistic events. The nature of the Richardson dispersion depends on the parameters of our flow and is discussed in terms of the values of a persistence parameter expressing the relative importance of the two above-mentioned mechanisms. We support this analysis by studying the distribution of interparticle distances, the relative velocity correlation functions, as well as the relative trajectories.

Brownian pump powered by a white-noise flasing ratchet

A Brownian pump of particles powered by a stochastic flashing ratchet mechanism is studied. The pumping device is embedded in a finite region and bounded by particle reservoirs. In the steady state, we exactly calculate the spatial density profile, the concentration ratio between both reservoirs and the particle flux. We propose a simulation framework for the consistent evaluation of such observable quantities.

Spatial correlations and cross sections of clusters in the A + B -> 0 reaction

We consider the distribution of cross sections of clusters and the density-density correlation functions for the A+B¿0 reaction. We solve the reaction-diffusion equations numerically for random initial distributions of reactants. When both reactant species have the same diffusion coefficients the distribution of cross sections and the correlation functions scale with the diffusion length and obey superuniversal laws (independent of dimension). For different diffusion coefficients the correlation functions still scale, but the scaling functions depend on the dimension and on the diffusion coefficients. Furthermore, we display explicitly the peculiarities of the cluster-size distribution in one dimension.

Heat fluctuations in Brownian transducers

The heat fluctuation probability distribution function in Brownian transducers operating between two heat reservoirs is studied. We find, both analytically and numerically, that the recently proposed fluctuation theorem for heat exchange [C. Jarzynski and D. K. Wojcik, Phys. Rev. Lett. 92, 230602 (2004)] has to be applied carefully when the coupling mechanism between both baths is considered. We also conjecture how to extend such a relation when an external work is present.

Tight coupling in thermal Brownian motors

We study analytically a thermal Brownian motor model and calculate exactly the Onsager coefficients. We show how the reciprocity relation holds and that the determinant of the Onsager matrix vanishes. Such a condition implies that the device is built with tight coupling. This explains why Carnot¿s efficiency can be achieved in the limit of infinitely slow velocities. We also prove that the efficiency at maximum power has the maximum possible value, which corresponds to the Curzon-Alhborn bound. Finally, we discuss the model acting as a Brownian refrigerator.