BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
H. J. Hilhorst appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
Recently Mertens and Moore [arXiv:1909.01484v1] showed that site percolation "is odd." By this they mean that on an $M\times N$ square lattice the number of distinct site configurations that allow for vertical percolation is odd. We report here an alternative proof, based on recursive use of geometric symmetry, for both free and periodic boundary conditions.
We study a one-dimensional lattice of $N$ sites each occupied by a mathematical "polymer," that is, is a binary random sequence of arbitrary length $n$, or equivalently, a rooted path of $n$ links on an infinite binary tree. The average polymer length is controlled by the monomer fugacity $z$. A pair of polymers on adjacent sites carries a weight factor $ω$ for each link on the tree that they have in common. The phase diagram in the $zω$ plane exhibits a critical line $z=z_{\rm c}(ω)$. For $z<z_{\rm c}(ω)$ there exists an equilibrium phase with, in particular, a finite average polymer length. We investigate the equilibrium ensemble by transfer matrix and Monte Carlo methods, paying particular attention to the vicinity of the critical line. For $z>z_{\rm c}(ω)$ the equilibrium is unstable and Monte Carlo time evolution brings about a dynamical symmetry breaking which favors the evolution of a small selection of polymers to ever greater length. While of interest for its own sake, this model may also be relevant to the prelife-to-life transition that has occurred during biological evolution. We compare it to existing models of similar simplicity due to Wu and Higgs (2009, 2012) and to Chen and Nowak (2012).
We consider the 3D Poisson-Voronoi tessellation. We investigate the joint probability distribution pi_n(L) for an arbitrarily selected cell face to be n-edged and for the distance between the seeds of its adjacent cells to be equal to 2L. We derive an exact expression for this quantity, valid in the limit n->infty with n^{1/6}L fixed. The leading order correction term is determined. Good agreement with earlier Monte Carlo data is obtained. The cell face is surrounded by a three-dimensional excluded domain that is the union of n balls; it is pumpkin-shaped and analogous to the flower of the 2D Voronoi cell. For n->infty this domain tends towards a torus of equal major and minor radii. The radii scale as n^{1/3}, in agreement with earlier heuristic work. We achieve a detailed understanding of several other statistical properties of the n-edged cell face.
The master equation for the reversible reaction A+A <--> 0 is considered in Poisson representation, where it is equivalent to a Langevin equation with imaginary noise for a complex stochastic variable ϕ. Such Langevin equations appear quite generally in field-theoretic treatments of reaction-diffusion problems. For this example we study the probability flow in the complex ϕplane both analytically and by simulation. We show that this flow has various curious features that must be expected to occur similarly in other Langevin equations associated with reaction-diffusion problems.
We consider the translationally invariant but parity (left-right symmetry) breaking Ising chain Hamiltonian \begin{equation} {\cal H} = -U_2\sum_{k} s_{k}s_{k+1} - U_3\sum_{k} s_{k}s_{k+1}s_{k+3} \nonumber \end{equation} and let this system evolve by Kawasaki spin exchange dynamics. Monte Carlo simulations show that perturbations forcing this system off equilibrium make it act as a Brownian molecular motor which, in the lattice gas interpretation, transports particles along the chain. We determine the particle current under various different circumstances, in particular as a function of the ratio $U_3/U_2$ and of the conserved magnetization $M=\sum_k s_k$. The symmetry of the $U_3$ term in the Hamiltonian is discussed
We consider the three-dimensional Poisson-Voronoi tessellation and study the average facedness m_n of a cell known to neighbor an n-faced cell. Whereas Aboav's law states that m_n=A+B/n, theoretical arguments indicate an asymptotic expansion m_n = 8 + k_1 n^{-1/6} +.... Recent new Monte Carlo data due to Lazar et al., based on a very large data set, now clearly rule out Aboav's law. In this work we determine the numerical value of k_1 and compare the expansion to the Monte Carlo data. The calculation of k_1 involves an auxiliary planar cellular structure composed of circular arcs, that we will call the Poisson-Moebius diagram. It is a special case of more general Moebius diagrams (or multiplicatively weighted power diagrams) and is of interest for its own sake. We obtain exact results for the total edge length per unit area, which is a prerequisite for the coefficient k_1, and a few other quantities in this diagram.
Motivated by recent new Monte Carlo data we investigate a heuristic asymptotic theory that applies to n-faced 3D Poisson-Voronoi cells in the limit of large n. We show how this theory may be extended to n-edged cell faces. It predicts the leading order large-n behavior of the average volume and surface area of the n-faced cell, and of the average area and perimeter of the n-edged face. Such a face is shown to be surrounded by a toroidal region of volume n/lambda (with lambda the seed density) that is void of seeds. Two neighboring cells sharing an n-edged face are found to have their seeds at a typical distance that scales as n^{-1/6} and whose probability law we determine. We present a new data set of 4*10^9 Monte Carlo generated 3D Poisson-Voronoi cells, larger than any before. Full compatibility is found between the Monte Carlo data and the theory. Deviations from the asymptotic predictions are explained in terms of subleading corrections whose powers in n we estimate from the data.
At the intersection of two unidirectional traffic flows a stripe formation instability is known to occur. In this paper we consider coupled time evolution equations for the densities of the two flows in their intersection area. We show analytically how the instability arises from the randomness of the traffic entering the area. The Green function of the linearized equations is shown to form a Gaussian wave packet whose oscillations correspond to the stripes. Explicit formulas are obtained for various characteristic quantities in terms of the traffic density and comparison is made with the much simpler calculation on a torus and with numerical solution of the evolution equations.
We study two perpendicular intersecting flows of pedestrians. The latter are represented either by moving hard core particles of two types, eastbound ($\symbp$) and northbound ($\symbm$), or by two density fields, $\rhop_t(\brr)$ and $\rhom_t(\brr)$. Each flow takes place on a lattice strip of width $M$ so that the intersection is an $M\times M$ square. We investigate the spontaneous formation, observed experimentally and in simulations, of a diagonal pattern of stripes in which alternatingly one of the two particle types dominates. By a linear stability analysis of the field equations we show how this pattern formation comes about. We focus on the observation, reported recently, that the striped pattern actually consists of chevrons rather than straight lines. We demonstrate that this `chevron effect' occurs both in particle simulations with various different update schemes and in field simulations. We quantify the effect in terms of the chevron angle $Δθ_0$ and determine its dependency on the parameters governing the boundary conditions.
Domain wall theory (DWT) has proved to be a powerful tool for the analysis of one-dimensional transport processes. A simple version of it was found very accurate for the Totally Asymmetric Simple Exclusion Process (TASEP) with random sequential update. However, a general implementation of DWT is still missing in the case of updates with less fluctuations, which are often more relevant for applications. Here we develop an exact DWT for TASEP with parallel update and deterministic (p=1) bulk motion. Remarkably, the dynamics of this system can be described by the motion of a domain wall not only on the coarse-grained level but also exactly on the microscopic scale for arbitrary system size. All properties of this TASEP, time-dependent and stationary, are shown to follow from the solution of a bivariate master equation whose variables are not only the position but also the velocity of the domain wall. In the continuum limit this exactly soluble model then allows us to perform a first principle derivation of a Fokker-Planck equation for the position of the wall. The diffusion constant appearing in this equation differs from the one obtained with the traditional `simple' DWT.
We consider two parallel cyclic Ising chains counter-rotating at a relative velocity v, the motion actually being a succession of discrete steps. There is an in-chain interaction between nearest-neighbor spins and a cross-chain interaction between instantaneously opposite spins. For velocities v>0 the system, subject to a suitable markovian dynamics at a temperature T, can reach only a nonequilibrium steady state (NESS). This system was introduced by Hucht et al., who showed that for v=\infty it undergoes a para- to ferromagnetic transition, essentially due to the fact that each chain exerts an effective field on the other one. The present study of the v=\infty case determines the consequences of the fluctuations of this effective field when the system size N is finite. We show that whereas to leading order the system obeys detailed balancing with respect to an effective time-independent Hamiltonian, the higher order finite-size corrections violate detailed balancing. Expressions are given to various orders in 1/N for the interaction free energy between the chains, the spontaneous magnetization, the in-chain and cross-chain spin-spin correlations, and the spontaneus magnetization. It is shown how finite-size scaling functions may be derived explicitly. This study was motivated by recent work on a two-lane traffic problem in which a similar phase transition was found.
A q-modified version of the central limit theorem due to Umarov et al. affirms that q-Gaussians are attractors under addition and rescaling of certain classes of strongly correlated random variables. The proof of this theorem rests on a nonlinear q-modified Fourier transform. By exhibiting an invariance property we show that this Fourier transform does not have an inverse. As a consequence, the theorem falls short of achieving its stated goal.
First we consider a unidirectional flux ω_bar of vehicles each of which is characterized by its `natural' velocity v drawn from a distribution P(v). The traffic flow is modeled as a collection of straight `world lines' in the time-space plane, with overtaking events represented by a fixed queuing time tau imposed on the overtaking vehicle. This geometrical model exhibits platoon formation and allows, among many other things, for the calculation of the effective average velocity w=ϕ(v) of a vehicle of natural velocity v. Secondly, we extend the model to two opposite lanes, A and B. We argue that the queuing time τin one lane is determined by the traffic density in the opposite lane. On the basis of reasonable additional assumptions we establish a set of equations that couple the two lanes and can be solved numerically. It appears that above a critical value ω_bar_c of the control parameter ω_bar the symmetry between the lanes is spontaneously broken: there is a slow lane where long platoons form behind the slowest vehicles, and a fast lane where overtaking is easy due to the wide spacing between the platoons in the opposite direction. A variant of the model is studied in which the spatial vehicle density ρ_bar rather than the flux ω_bar is the control parameter. Unequal fluxes ω_bar_A and ω_bar_B in the two lanes are also considered. The symmetry breaking phenomenon exhibited by this model, even though no doubt hard to observe in pure form in real-life traffic, nevertheless indicates a tendency of such traffic.
We develop a set of heuristic arguments to explain several results on planar Poisson-Voronoi tessellations that were derived earlier at the cost of considerable mathematical effort. The results concern Voronoi cells having a large number n of sides. The arguments start from an entropy balance applied to the arrangement of n neighbors around a central cell. It is followed by a simplified evaluation of the phase space integral for the probability p_n that an arbitrary cell be n-sided. The limitations of the arguments are indicated. As a new application we calculate the expected number of Gabriel (or full) neighbors of an n-sided cell in the large-n limit.
In this talk I first review at an elementary level a selection of central limit theorems, including some lesser known cases, for sums and maxima of uncorrelated and correlated random variables. I recall why several of them appear in physics. Next, I show that there is room for new versions of central limit theorems applicable to specific classes of problems. Finally, I argue that we have insufficient evidence that, as a consequence of such a theorem, q-Gaussians occupy a special place in statistical physics.
We consider the d-dimensional Poisson-Voronoi tessellation and investigate the applicability of heuristic methods developed recently for two dimensions. Let p_n(d) be the probability that a cell have n neighbors (be `n-faced') and m_n(d) the average facedness of a cell adjacent to an n-faced cell. We obtain the leading order terms of the asymptotic large-n expansions for p_n(d) and m_n(3). It appears that, just as in dimension two, the Poisson-Voronoi tessellation violates Aboav's `linear law' also in dimension three. A confrontation of this statement with existing Monte Carlo work remains inconclusive. However, simulations upgraded to the level of present-day computer capacity will in principle be able to confirm (or invalidate) our theory.
We consider a family of random line tessellations of the Euclidean plane introduced in a much more formal context by Hug and Schneider [Geom. Funct. Anal. 17, 156 (2007)] and described by a parameter α\geq 1. For α=1 the zero-cell (that is, the cell containing the origin) coincides with the Crofton cell of a Poisson line tessellation, and for α=2 it coincides with the typical Poisson-Voronoi cell. Let p_n(α) be the probability for the zero-cell to have n sides. By the methods of statistical mechanics we construct the asymptotic expansion of \log p_n(α) up to terms that vanish as n\to\infty. In the large-n limit the cell is shown to become circular. The circle is centered at the origin when α>1, but gets delocalized for the Crofton cell, α=1, which is a singular point of the parameter range. The large-n expansion of \log p_n(1) is therefore different from that of the general case and we show how to carry it out. As a corollary we obtain the analogous expansion for the {\it typical} n-sided cell of a Poisson line tessellation.
I present a concise review of advances realized over the past three years on planar Poisson-Voronoi tessellations. These encompass new analytic results, a new Monte Carlo method, and application to experimental data.
Let n points be chosen randomly and independently in the unit disk. "Sylvester's question" concerns the probability p_n that they are the vertices of a convex n-sided polygon. Here we establish the link with another problem. We show that for large n this polygon, when suitably parametrized by a function r(phi) of the polar angle phi, satisfies the equation of the random acceleration process (RAP), d^2 r/d phi^2 = f(phi), where f is Gaussian noise. On the basis of this relation we derive the asymptotic expansion log p_n = -2n log n + n log(2 pi^2 e^2) - c_0 n^{1/5} + ..., of which the first two terms agree with a rigorous result due to Barany. The nonanalyticity in n of the third term is a new result. The value 1/5 of the exponent follows from recent work on the RAP due to Gyorgyi et al. [Phys. Rev. E 75, 021123 (2007)]. We show that the n-sided polygon is effectively contained in an annulus of width \sim n^{-4/5} along the edge of the disk. The distance delta_n of closest approach to the edge is exponentially distributed with average 1/(2n).
By a new Monte Carlo algorithm we evaluate the sidedness probability p_n of a planar Poisson-Voronoi cell in the range 3 \leq n \leq 1600. The algorithm is developed on the basis of earlier theoretical work; it exploits, in particular, the known asymptotic behavior of p_n as n\to\infty. Our p_n values all have between four and six significant digits. Accurate n dependent averages, second moments, and variances are obtained for the cell area and the cell perimeter. The numerical large n behavior of these quantities is analyzed in terms of asymptotic power series in 1/n. Snapshots are shown of typical occurrences of extremely rare events implicating cells of up to n=1600 sides embedded in an ordinary Poisson-Voronoi diagram. We reveal and discuss the characteristic features of such many-sided cells and their immediate environment. Their relevance for observable properties is stressed.
We show that the contact process in an ensemble with conserved total particle number, as simulated recently by Tome and de Oliveira [Phys. Rev. Lett. 86 (2001) 5463], is equivalent to the ordinary contact process, in agreement with what the authors assumed and believed. Similar conserved ensembles and equivalence proofs are easily constructed for other models.
Our recent Monte Carlo results [Phys. Rev. E 61 (2000) 6330] for the one-dimensional reaction-diffusion process A+B->2B and B->A lead to the correlation length exponent estimate nu=2.21 +/- 0.05. In a comment on our work [cond-mat/0007366] the exact value is claimed to be nu=2. We reply that the arguments advanced in the comment fail to prove the claim.
The Riemann walk is the lattice version of the Levy flight. For the one-dimensional Riemann walk of Levy exponent 0<α<2 we study the statistics of the support, i.e. the set of visited sites, after t steps. We consider a wide class of support related observables M(t), including the number S(t) of visited sites and the number I(t) of sequences of visited sites. For t->\infty we obtain the asymptotic power laws for the averages, variances, and correlations of these observables. Logarithmic correction factors appear for α=2/3 and α=1. Bulk and surface observables have different power laws for 1\leqα<2. Fluctuations are shown to be universal for 2/3\leqα<2. This means that in the limit t->\infty the deviations from average \DeltaM(t) are fully described (i) either by a single M independent stochastic process (when 2/3\leqα\leq 1) (ii) or by two such processes, one for the bulk and one for the surface observables (when 1<α<2).
We consider the d-dimensional diffusion equation for a field phi(x,t) with random initial condition, and observe that, when appropriately scaled, phi(0,t) is Gaussian and Markovian in the limit d->0. This leads via the Majumdar-Sire perturbation theory to a small-d expansion for the persistence exponent theta(d). We find theta(d) = d/4 - 0.12065...d^{3/2} + ...
Healthy and sick individuals (A and B particles) diffuse independently with diffusion constants D_A and D_B. Sick individuals upon encounter infect healthy ones (at rate k), but may also spontaneously recover (at rate 1/τ). The propagation of the epidemic therefore couples to the fluctuations in the total population density. Global extinction occurs below a critical value ρ_{c} of the spatially averaged total density. The epidemic evolves as the diffusion--reaction--decay process A + B --> 2B, B --> A , for which we write down the field theory. The stationary state properties of this theory when D_A=D_B were obtained by Kree et al. The critical behavior for D_A<D_B is governed by a new fixed point. We calculate the critical exponents of the stationary state in an $\eps$ expansion, carried out by Wilson renormalization, below the critical dimension d_{c}=4. We then go on to to obtain the critical initial time behavior at the extinction threshold, both for D_A=D_B and D_A<D_B. There is nonuniversal dependence on the initial particle distribution. The case D_A>D_B remains unsolved.
We analyse the chiral symmetry in the random $\pm J$ $XY$ model on a $N\times 2$ square lattice with periodic boundary conditions in the transverse direction. This ``tube" lattice may be seen as a two-dimensional lattice of which one dimension has been compactified. In the Villain formulation the discrete-valued {\em chiralities}\/ or {\em charges}\/ associated with the plaquettes of the lattice decouple from the continuous degrees of freedom. The difficulty of the problem lies in the fact that the chiralities interact through the long range ``strong" one-dimensional Coulomb potential - which increases linearly with distance - as well as through an exponentially decaying ``weak" interaction. By comparing the ground state energies for periodic, antiperiodic, and reflecting boundary conditions in the longitudinal direction, we show that the chiralities and the $XY$ spins have the {\em same}\/ zero-$T$ correlation length exponent, whose exact value $ν_c = 0.5564\ldots$ we determine. The equality of these correlation lengths even in the presence of long range chirality-chirality interactions lends support to the view that chiral glass order cannot be sustained without simultaneous spin glass order.