We are interested in a fragmentation process. We observe fragments frozen when their sizes are less than $\varepsilon$ ($\varepsilon>0$). It is known (Bertoin and Martínez, 2005) that the empirical measure of these fragments converges in law, under some renormalization. Hoffmann and Krell (2011) showed a bound for the rate of convergence. Here, we show a central limit theorem, under some assumptions. This gives us an exact rate of convergence.

In large storage systems, files are often coded across several servers to improve reliability and retrieval speed. We study load balancing under the batch sampling routeing scheme for a network of n servers storing a set of files using the maximum distance separable (MDS) code (cf. Li (2016)). Specifically, each file is stored in equally sized pieces across L servers such that any k pieces can reconstruct the original file. When a request for a file is received, the dispatcher routes the job into the k-shortest queues among the L for which the corresponding server contains a piece of the file being requested. We establish a law of large numbers and a central limit theorem as the system becomes large (i.e. n → ∞), for the setting where all interarrival and service times are exponentially distributed. For the central limit theorem, the limit process take values in ℓ2, the space of square summable sequences. Due to the large size of such systems, a direct analysis of the n-server system is frequently intractable. The law of large numbers and diffusion approximations established in this work provide practical tools with which to perform such analysis. The power-of-d routeing scheme, also known as the supermarket model, is a special case of the model considered here.

In this paper we develop a collection of results associated to the analysis of the sequential Monte Carlo (SMC) samplers algorithm, in the context of high-dimensional independent and identically distributed target probabilities. The SMC samplers algorithm can be designed to sample from a single probability distribution, using Monte Carlo to approximate expectations with respect to this law. Given a target density in d dimensions our results are concerned with d → ∞, while the number of Monte Carlo samples, N, remains fixed. We deduce an explicit bound on the Monte-Carlo error for estimates derived using the SMC sampler and the exact asymptotic relative -error of the estimate of the normalising constant associated to the target. We also establish marginal propagation of chaos properties of the algorithm. These results are deduced when the cost of the algorithm is O(Nd2).

A discrete time stochastic model for a multiagent system given in terms of a large collection of interacting Markov chains is studied. The evolution of the interacting particles is described through a time inhomogeneous transition probability kernel that depends on the ‘gradient’ of the potential field. The particles, in turn, dynamically modify the potential field through their cumulative input. Interacting Markov processes of the above form have been suggested as models for active biological transport in response to external stimulus such as a chemical gradient. One of the basic mathematical challenges is to develop a general theory of stability for such interacting Markovian systems and for the corresponding nonlinear Markov processes that arise in the large agent limit. Such a theory would be key to a mathematical understanding of the interactive structure formation that results from the complex feedback between the agents and the potential field. It will also be a crucial ingredient in developing simulation schemes that are faithful to the underlying model over long periods of time. The goal of this work is to study qualitative properties of the above stochastic system as the number of particles (N) and the time parameter (n) approach infinity. In this regard asymptotic properties of a deterministic nonlinear dynamical system, that arises in the propagation of chaos limit of the stochastic model, play a key role. We show that under suitable conditions this dynamical system has a unique fixed point. This result allows us to study stability properties of the underlying stochastic model. We show that as N → ∞, the stochastic system is well approximated by the dynamical system, uniformly over time. As a consequence, for an arbitrarily initialized system, as N → ∞ and n → ∞, the potential field and the empirical measure of the interacting particles are shown to converge to the unique fixed point of the dynamical system. In general, simulation of such interacting Markovian systems is a computationally daunting task. We propose a particle based approximation for the dynamic potential field which allows for a numerically tractable simulation scheme. It is shown that this simulation scheme well approximates the true physical system, uniformly over an infinite time horizon.

In this paper we are interested in a nonlinear parabolic evolution equation occurring in rheology. We give a probabilistic interpretation to this equation by associating a nonlinear martingale problem with it. We prove the existence of a unique solution, P, to this martingale problem. For any t, the time marginal of P at time t admits a density ρ(t,x) with respect to the Lebesgue measure, where the function ρ is the unique weak solution to the evolution equation in a well-chosen energy space. Next we introduce a simulable system of n interacting particles and prove that the empirical measure of this system converges to P as n tends to ∞. This propagation-of-chaos result ensures that the solution to the equation of interest can be approximated using a Monte Carlo method. Finally, we illustrate the convergence in some numerical experiments.

A result for the propagation of chaos is obtained for a class of pure jump particle systems of two species with mean field interaction. This result leads to the corresponding result for particle systems with one species and the argument used is valid for particle systems with more than two species. The model is motivated by the study of the phenomenon of self-organization in biology, chemistry and physics, and the technical difficulty is the unboundedness of the jump rates.

The paper deals with large trunk line systems of the type appearing in telephone networks. There are many nodes or input sources, each pair of which is connected by a trunk line containing many individual circuits. Traffic arriving at either end of a trunk line wishes to communicate to the node at the other end. If the direct route is full, a rerouting might be attempted via an alternative route containing several trunks and connecting the same endpoints. The basic questions concern whether to reroute, and if so how to choose the alternative path. If the network is ‘large’ and fully connected, then the overflow traffic which is offered for rerouting to any trunk comes from many other trunks in the network with no one dominating. In this case one expects that some sort of averaging method can be used to approximate the rerouting requests and hence simplify the analysis. Essentially, the overflow traffic that a trunk offers the network for rerouting is in some average sense similar to the overflow traffic offered to that trunk. Indeed, a formalization of this idea involves the widely used (but generally heuristic) ‘fixed point' approximation method. One sets up the fixed point equations for appropriate rerouting strategies and then solves them to obtain an approximation to the system loss. In this paper we work in the heavy traffic regime, where the external offered traffic to any trunk is close to the service capacity of that trunk. It is shown that, as the number of links and circuits within each link go to infinity and for a variety of rerouting strategies, the system can be represented by an averaged limit. This limit is a reflected diffusion of the McKean–Vlasov (propagation of chaos) type, where the driving terms depend on the mean values of the solution of the equation. The averages occur due to the symmetry of the network and the averaging effects of the many interactions. This provides a partial justification for the fixed point method. The concrete dynamical systems flavor of the approach and the representations of the limit processes provide a useful way of visualizing the system and promise to be useful for the development of numerical methods and further analysis.

We consider the N-oscillator system of the van der Pol type, which contains a small positive parameter ε multiplying the non-linear damping and the random disturbance. For a formulation of the output we take the solution X(t) = (Xi(t))i=1···,N of the system of 2N-dimensional stochastic differential equations. Rotating each component Xi(t) about the origin of the plane by an angle t, we find that on time scales of order 1/ε together with sufficiently large N each Xi(t) behaves as the equi-ultimately bounded solution of an equation of the McKean type admitting a stationary probability distribution.