In this paper we consider an infinite system of instantaneously coalescing rate 1 simple symmetric random walks on ℤ2, started from the initial condition with all sites in ℤ2 occupied. Two-dimensional coalescing random walks are a `critical' model of interacting particle systems: unlike coalescence models in dimension three or higher, the fluctuation effects are important for the description of large-time statistics in two dimensions, manifesting themselves through the logarithmic corrections to the `mean field' answers. Yet the fluctuation effects are not as strong as for the one-dimensional coalescence, in which case the fluctuation effects modify the large time statistics at the leading order. Unfortunately, unlike its one-dimensional counterpart, the two-dimensional model is not exactly solvable, which explains a relative scarcity of rigorous analytic answers for the statistics of fluctuations at large times. Our contribution is to find, for any N≥2, the leading asymptotics for the correlation functions ρN(x1,…,xN) as t→∞. This generalises the results for N=1 due to Bramson and Griffeath (1980) and confirms a prediction in the physics literature for N>1. An analogous statement holds for instantaneously annihilating random walks. The key tools are the known asymptotic ρ1(t)∼logt∕πt due to Bramson and Griffeath (1980), and the noncollision probability 𝒑NC(t), that no pair of a finite collection of N two-dimensional simple random walks meets by time t, whose asymptotic 𝒑NC(t)∼c0(logt)-(N2) was found by Cox et al. (2010). We re-derive the asymptotics, and establish new error bounds, both for ρ1(t) and 𝒑NC(t) by proving that these quantities satisfy effective rate equations; that is, approximate differential equations at large times. This approach can be regarded as a generalisation of the Smoluchowski theory of renormalised rate equations to multi-point statistics.

A new specially correlated partially coherent beam named nonuniform multi-Gaussian correlated (NMGC) partially coherent beam is introduced. The correlation functions of such beam in $x$ and $y$ directions are different from each other, i.e., nonuniform correlation function in one direction and multi-Gaussian correlated Schell-model function in the other direction. The propagation properties of an NMGC partially coherent beam in free pace are demonstrated, and we find that the intensity distribution of such beam exhibits self-focusing and self-shifting effect in one direction and self-shaping effect in the other direction on propagation. The correlation-induced self-focusing and self-shaping effect will be useful in some applications, where the high power and shaped laser is required, such as material thermal processing and laser carving.

Gaussian particles provide a flexible framework for modelling and simulating three-dimensional star-shaped random sets. In our framework, the radial function of the particle arises from a kernel smoothing, and is associated with an isotropic random field on the sphere. If the kernel is a von Mises-Fisher density, or uniform on a spherical cap, the correlation function of the associated random field admits a closed form expression. The Hausdorff dimension of the surface of the Gaussian particle reflects the decay of the correlation function at the origin, as quantified by the fractal index. Under power kernels we obtain particles with boundaries of any Hausdorff dimension between 2 and 3.

We derive explicit lower and upper bounds for the probability generating functional of a stationary locally stable Gibbs point process, which can be applied to summary statistics such as the F function. For pairwise interaction processes we obtain further estimates for the G and K functions, the intensity, and higher-order correlation functions. The proof of the main result is based on Stein's method for Poisson point process approximation.

In this paper we consider a single-server queue with Lévy input, and, in particular, its workload process (Qt)t≥0, focusing on its correlation structure. With the correlation function defined as r(t):= cov(Q0, Qt) / varQ0 (assuming that the workload process is in stationarity at time 0), we first study its transform ∫0∞r(t)e-ϑtdt, both for when the Lévy process has positive jumps and when it has negative jumps. These expressions allow us to prove that r(·) is positive, decreasing, and convex, relying on the machinery of completely monotone functions. For the light-tailed case, we estimate the behavior of r(t) for large t. We then focus on techniques to estimate r(t) by simulation. Naive simulation techniques require roughly (r(t))-2 runs to obtain an estimate of a given precision, but we develop a coupling technique that leads to substantial variance reduction (the required number of runs being roughly (r(t))-1). If this is augmented with importance sampling, it even leads to a logarithmically efficient algorithm.

We introduce a model which uses semi-analytic techniques to trace formation and evolution of galaxy disks in their cosmological context. For the first time we model the growth of gas and stellar disks separately. In contrast to previous work we follow in detail the angular momentum accumulation history through the gas cooling, merging and star formation processes. Our model successfully reproduces the stellar mass–radius distribution and gas-to-stellar disk size ratio distribution observed locally. We also investigate the dependence of clustering on galaxy size and find qualitative agreement with observation. There is still some discrepancy at small scale for less massive galaxies, indicating that our treatment of satellite galaxies needs to be improved.

Martin and Walker ((1997) J. Appl. Prob.34, 657–670) proposed the power-law ρ(v) = c|v|-β, |v| ≥ 1, as a correlation model for stationary time series with long-memory dependence. A straightforward proof of their conjecture on the permissible range of c is given, and various other models for long-range dependence are discussed. In particular, the Cauchy family ρ(v) = (1 + |v/c|α)-β/α allows for the simultaneous fitting of both the long-term and short-term correlation structure within a simple analytical model. The note closes with hints at the fast and exact simulation of fractional Gaussian noise and related processes.

Results and methods discussed in two previous papers are extended to other cases. Comparison is made with recent work by Pickard (1980) and an earlier conjecture is disproved.

The G/G/1 queueing system is studied by means of the regeneration point method, exploiting the concept of busy cycles. Recurrence relations are set-up for the distribution of the queue length at the arrival epochs. The same method is used to obtain the correlation structure of arrivals and departures for the G/M/1 queue.

We consider a two-dimensional stochastic process of binary variables xij defined so that the conditional distribution of xij given its predecessors depends only on xi–1j and xij–1 Properties of the equilibrium distribution, when this exists, are investigated using three different representations of the process and explicit results are given in some special cases.