BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

David P. Sanders

David P. Sanders appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechnlin.CDmath.DSmath-phmath.MPmath.PRphysics.chem-phcond-mat.softnlin.CGnlin.SIphysics.comp-phphysics.soc-phPopulations and Evolution

Catalog footprint

What is connected

20works
13topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2016arXiv

Ergodicity of one-dimensional systems coupled to the logistic thermostat

We analyze the ergodicity of three one-dimensional Hamiltonian systems, with harmonic, quartic and Mexican-hat potentials, coupled to the logistic thermostat. As criteria for ergodicity we employ: the independence of the Lyapunov spectrum with respect to initial conditions; the absence of visual "holes" in two-dimensional Poincaré sections; the agreement between the histograms in each variable and the theoretical marginal distributions; and the convergence of the global joint distribution to the theoretical one, as measured by the Hellinger distance. Taking a large number of random initial conditions, for certain parameter values of the thermostat we find no indication of regular trajectories and show that the time distribution converges to the ensemble one for an arbitrarily long trajectory for all the systems considered. Our results thus provide a robust numerical indication that the logistic thermostat can serve as a single one-parameter thermostat for stiff one-dimensional systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Geometric integrator for simulations in the canonical ensemble

In this work we introduce a geometric integrator for molecular dynamics simulations of physical systems in the canonical ensemble. In particular, we consider the equations arising from the so-called density dynamics algorithm with any possible type of thermostat and provide an integrator that preserves the invariant distribution. Our integrator thus constitutes a unified framework that allows the study and comparison of different thermostats and of their influence on the equilibrium and non-equilibrium (thermo-)dynamic properties of the system. To show the validity and the generality of the integrator, we implement it with a second-order, time-reversible method and apply it to the simulation of a Lennard-Jones system with three different thermostats, obtaining good conservation of the geometrical properties and recovering the expected thermodynamic results.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Efficient algorithms for general periodic Lorentz gases in two and three dimensions

We present efficient algorithms to calculate trajectories for periodic Lorentz gases consisting of square lattices of circular obstacles in two dimensions, and simple cubic lattices of spheres in three dimensions; these become increasingly efficient as the radius of the obstacles tends to 0, the so-called Boltzmann-Grad limit. The 2D algorithm applies continued fractions to obtain the exact disc with which a particle will collide at each step, instead of using periodic boundary conditions as in the classical algorithm. The 3D version incorporates the 2D algorithm by projecting to the three coordinate planes. As an application, we calculate distributions of free path lengths close to the Boltzmann-Grad limit for certain Lorentz gases. We also show how the algorithms may be applied to deal with general crystal lattices.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Horizons and free path distributions in quasiperiodic Lorentz gases

We study the structure of quasiperiodic Lorentz gases, i.e., particles bouncing elastically off fixed obstacles arranged in quasiperiodic lattices. By employing a construction to embed such structures into a higher dimensional periodic hyperlattice, we give a simple and efficient algorithm for numerical simulation of the dynamics of these systems. This same construction shows that quasiperiodic Lorentz gases generically exhibit a regime with infinite horizon, that is, empty channels through which the particles move without colliding, when the obstacles are small enough; in this case, the distribution of free paths is asymptotically a power law with exponent -3, as expected from infinite-horizon periodic Lorentz gases. For the critical radius at which these channels disappear, however, a new regime with locally-finite horizon arises, where this distribution has an unexpected exponent of -5, previously observed only in a Lorentz gas formed by superposing three incommensurable periodic lattices in the Boltzmann-Grad limit where the radius of the obstacles tends to zero.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Lévy walks on lattices as multi-state processes

Continuous-time random walks combining diffusive scattering and ballistic propagation on lattices model a class of Lévy walks. The assumption that transitions in the scattering phase occur with exponentially-distributed waiting times leads to a description of the process in terms of multiple states, whose distributions evolve according to a set of delay differential equations, amenable to analytic treatment. We obtain an exact expression of the mean squared displacement associated with such processes and discuss the emergence of asymptotic scaling laws in regimes of diffusive and superdiffusive (subballistic) transport, emphasizing, in the latter case, the effect of initial conditions on the transport coefficients. Of particular interest is the case of rare ballistic propagation, in which case a regime of superdiffusion may lurk underneath one of normal diffusion.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Emergence of a sustained current by coupling equilibrating systems: Making a NESS out of equilibrium

We show that coupling together two closed thermodynamic systems that independently attain equilibrium may give rise to a nonequilibrium stationary state (NESS) with a persistent, non-vanishing current. We study a simple example that is exactly soluble, consisting of random walkers with different biases towards a reflecting boundary, modelling, for example, Brownian particles with different charge states in an electric field. We obtain exact analytical expressions for the (generating functions of) concentrations and currents in the NESS for this model, and exhibit the main features by numerical simulation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Machta-Zwanzig regime of anomalous diffusion in infinite-horizon billiards

We study diffusion on a periodic billiard table with infinite horizon in the limit of narrow corridors. An effective trapping mechanism emerges according to which the process can be modeled by a Lévy walk combining exponentially-distributed trapping times with free propagation along paths whose precise probabilities we compute. This description yields an approximation of the mean squared displacement of infinite-horizon billiards in terms of two transport coefficients which generalizes to this anomalous regime the Machta-Zwanzig approximation of normal diffusion in finite-horizon billiards [Phys. Rev. Lett. 50, 1959 (1983)].

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Measuring logarithmic corrections to normal diffusion in infinite-horizon billiards

We perform numerical measurements of the moments of the position of a tracer particle in a two-dimensional periodic billiard model (Lorentz gas) with infinite corridors. This model is known to exhibit a weak form of super-diffusion, in the sense that there is a logarithmic correction to the linear growth in time of the mean-squared displacement. We show numerically that this expected asymptotic behavior is easily overwhelmed by the subleading linear growth throughout the time-range accessible to numerical simulations. We compare our simulations to the known analytical results for the variance of the anomalously-rescaled limiting normal distributions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Transport properties of Lévy walks: an analysis in terms of multistate processes

Continuous time random walks combining diffusive and ballistic regimes are introduced to describe a class of Lévy walks on lattices. By including exponentially-distributed waiting times separating the successive jump events of a walker, we are led to a description of such Lévy walks in terms of multistate processes whose time-evolution is shown to obey a set of coupled delay differential equations. Using simple arguments, we obtain asymptotic solutions to these equations and rederive the scaling laws for the mean squared displacement of such processes. Our calculation includes the computation of all relevant transport coefficients in terms of the parameters of the models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Encounter times in overlapping domains: application to epidemic spread in a population of territorial animals

We develop an analytical method to calculate encounter times of two random walkers in one dimension when each individual is segregated in its own spatial domain and shares with its neighbor only a fraction of the available space, finding very good agreement with numerically-exact calculations. We model a population of susceptible and infected territorial individuals with this spatial arrangement, and which may transmit an epidemic when they meet. We apply the results on encounter times to determine analytically the macroscopic propagation speed of the epidemic as a function of the microscopic characteristics: the confining geometry, the animal diffusion constant, and the infection transmission probability.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Structure and evolution of strange attractors in non-elastic triangular billiards

We study pinball billiard dynamics in an equilateral triangular table. In such dynamics, collisions with the walls are non-elastic: the outgoing angle with the normal vector to the boundary is a uniform factor $λ< 1$ smaller than the incoming angle. This leads to contraction in phase space for the discrete-time dynamics between consecutive collisions, and hence to attractors of zero Lebesgue measure, which are almost always fractal strange attractors with chaotic dynamics, due to the presence of an expansion mechanism. We study the structure of these strange attractors and their evolution as the contraction parameter $λ$ is varied. For $λ$ in the interval (0, 1/3), we prove rigorously that the attractor has the structure of a Cantor set times an interval, whereas for larger values of $λ$ the billiard dynamics gives rise to nonaccessible regions in phase space. For $λ$ close to 1, the attractor splits into three transitive components, the basins of attraction of which have fractal basin boundaries.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Chaos and stability in a two-parameter family of convex billiard tables

We study, by numerical simulations and semi-rigorous arguments, a two-parameter family of convex, two-dimensional billiard tables, generalizing the one-parameter class of oval billiards of Benettin--Strelcyn [Phys. Rev. A 17, 773 (1978)]. We observe interesting dynamical phenomena when the billiard tables are continuously deformed from the integrable circular billiard to different versions of completely-chaotic stadia. In particular, we conjecture that a new class of ergodic billiard tables is obtained in certain regions of the two-dimensional parameter space, when the billiards are close to skewed stadia. We provide heuristic arguments supporting this conjecture, and give numerical confirmation using the powerful method of Lyapunov-weighted dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Chaos in cylindrical stadium billiards via a generic nonlinear mechanism

We describe conditions under which higher-dimensional billiard models in bounded, convex regions are fully chaotic, generalizing the Bunimovich stadium to dimensions above two. An example is a three-dimensional stadium bounded by a cylinder and several planes; the combination of these elements may give rise to defocusing, allowing large chaotic regions in phase space. By studying families of marginally-stable periodic orbits that populate the residual part of phase space, we identify conditions under which a nonlinear instability mechanism arises in their vicinity. For particular geometries, this mechanism rather induces stable nonlinear oscillations, including in the form of whispering-gallery modes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stable and unstable regimes in higher-dimensional convex billiards with cylindrical shape

We introduce a class of convex, higher-dimensional billiard models which generalise stadium billiards. These models correspond to the free motion of a point-particle in a region bounded by cylinders cut by planes. They are motivated by models of particles interacting via a string-type mechanism, and confined by hard walls. The combination of these elements may give rise to a defocusing mechanism, similar to that in two dimensions, which allows large chaotic regions in phase space. The remaining part of phase space is associated with marginally stable behaviour. In fact periodic orbits in these systems generically come in continuous parametric families, sociated with a pair of parabolic eigen-directions: the periodic orbits are unstable in the presence of a defocusing mechanism, but marginally stable otherwise. By performing the stability analysis of families of periodic orbits at a nonlinear level, we establish the conditions under which families are nonlinearly stable or unstable. As a result, we identify regions in the parameter space of the models which admit non-linearly stable oscillations in the form of whispering gallery modes. Where no families of periodic orbits are stable, the billiards are completely chaotic, i.e.\ the Lyapunov exponents of the billiard map are non-zero.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Diffusion coefficients for multi-step persistent random walks on lattices

We calculate the diffusion coefficients of persistent random walks on lattices, where the direction of a walker at a given step depends on the memory of a certain number of previous steps. In particular, we describe a simple method which enables us to obtain explicit expressions for the diffusion coefficients of walks with two-step memory on different classes of one-, two- and higher-dimensional lattices.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Exact encounter times for many random walkers on regular and complex networks

The exact mean time between encounters of a given particle in a system consisting of many particles undergoing random walks in discrete time is calculated, on both regular and complex networks. Analytical results are obtained both for independent walkers, where any number of walkers can occupy the same site, and for walkers with an exclusion interaction, when no site can contain more than one walker. These analytical results are then compared with numerical simulations, showing very good agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Persistence effects in deterministic diffusion

In systems which exhibit deterministic diffusion, the gross parameter dependence of the diffusion coefficient can often be understood in terms of random walk models. Provided the decay of correlations is fast enough, one can ignore memory effects and approximate the diffusion coefficient according to dimensional arguments. By successively including the effects of one and two steps of memory on this approximation, we examine the effects of ``persistence'' on the diffusion coefficients of extended two-dimensional billiard tables and show how to properly account for these effects, using walks in which a particle undergoes jumps in different directions with probabilities that depend on where they came from.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

A Definition of Metastability for Markov Processes with Detailed Balance

A definition of metastable states applicable to arbitrary finite state Markov processes satisfying detailed balance is discussed. In particular, we identify a crucial condition that distinguishes genuine metastable states from other types of slowly decaying modes and which leads to properties similar to those postulated in the restricted ensemble approach \cite{pen71}. The intuitive physical meaning of this condition is simply that the total equilibrium probability of finding the system in the metastable state is negligible. As a concrete application of our formalism we present preliminary results on a 2D kinetic Ising model.

0 citations0 reviewsNo code linkedOpen access