We generalize a two-type mutation process in which particles reproduce by binary fission, inheriting the parental type, but which can mutate with small probability during their lifetimes to the opposite type. The generalization allows an arbitrary offspring distribution. The branching process structure of this scheme is exploited to obtain a variety of limit theorems, some of which extend known results for the binary case. In particular, practically usable asymptotic normality results are obtained when the initial population size is large.

We study conditions for the existence of non-trivial quasi-stationary distributions for the birth-and-death chain with 0 as absorbing state. We reduce our problem to a continued fractions one that can be solved by using extensions of classical results of this theory. We also prove that there exist normalized quasi-stationary distributions if and only if 0 is geometrically absorbing.

A unified exposition of random Johnson–Mehl tessellations in d-dimensional Euclidean space is presented. In particular, Johnson-Mehl tessellations generated by time-inhomogeneous Poisson processes and nucleation-exclusion models are studied. The ‘practical' cases d = 2 and d = 3 are discussed in detail. Several new results are established, including first- and second-order moments of various characteristics for both Johnson–Mehl tesselations and sectional Johnson–Mehl tessellations.

A three-parameter model of a random directed graph (digraph) is specified by the probability of ‘up arrows' from vertex i to vertex j where i < j, the probability of ‘down arrows' from i to j where i ≥ j, and the probability of bidirectional arrows between i and j. In this model, a phase transition—the abrupt appearance of a giant strongly connected component—takes place as the parameters cross a critical surface. The critical surface is determined explicitly. Before the giant component appears, almost surely all non-trivial components are small cycles. The asymptotic probability that the digraph contains no cycles of length 3 or more is computed explicitly. This model and its analysis are motivated by the theory of food webs in ecology.

Let X1, X2, · ·· be independent and identically distributed random variables such that ΕΧ1 < 0 and P(X1 ≥ 0) ≥ 0. Fix M ≥ 0 and let T = inf {n: X1 + X2 + · ·· + Xn ≥ M} (T = +∞, if for every n = 1,2, ···). In this paper we consider the estimation of the level-crossing probabilities P(T <∞) and , by using Monte Carlo simulation and especially importance sampling techniques. When using importance sampling, precision and efficiency of the estimation depend crucially on the choice of the simulation distribution. For this choice we introduce a new criterion which is of the type of large deviations theory; consequently, the basic large deviations theory is the main mathematical tool of this paper. We allow a wide class of possible simulation distributions and, considering the case that M →∞, we prove asymptotic optimality results for the simulation of the probabilities P(T <∞) and . The paper ends with an example.

We shall establish functional limit laws for the concentration of the various species in simple chemical reactions. These results allow us to conclude that, under quite general conditions, the concentration has an approximate normal distribution. We provide estimates for the mean and the variance which are valid at all stages of the reaction, in particular, the non-equilibrium phase. We also provide a detailed comparison of our results with the earlier work of Dunstan and Reynolds ([7], [8]).

In this paper we consider parallel and series systems, the components of which can be ‘improved'. The ‘improvement' consists of supplying the components with cold or hot standby spares or by allotting to them fixed budgets for minimal repairs. A fixed total resource of spares or minimal repairs is available. We find the optimal allocation of the resource items in several commonly encountered settings.

In flexible assembly systems, it is often necessary to coordinate jobs and materials so that specific jobs are matched with specific materials. This requires that jobs depart from upstream parallel workstations in some predetermined order. One way to satisfy this requirement is to temporarily hold the serviced jobs getting out of order at a resequencing buffer and to release them to downstream workstations as soon as all their predecessors are serviced. In this paper we consider the problem of scheduling a fixed number of non-preemptive jobs on two IHR non-identical processors with the resequencing requirement. We prove that the individually optimal policy, in which each job minimizes its own expected departure time subject to the constraint that available processors are offered to jobs in their departure order, is of a threshold type. The policy is independent of job weights and the jobs residing at the resequencing buffer and possesses the monotonicity property which states that a job will never utilize a processor in the future once it has declined the processor. Most importantly, we prove that the individually optimal policy has the stability property; namely: if at any time a job deviated from the individually optimal policy, then the departure time of every job, including its own, would be prolonged. As a direct consequence of this property, the individually optimal policy is socially optimal in the sense that it minimizes the expected total weighted departure time of the system as a whole. We identify situations under which the individually optimal policy also minimizes the expected makespan of the system.

We provide a general framework for interconnecting a collection of quasi-reversible nodes in such a way that the resulting process exhibits a product-form invariant measure. The individual nodes can be quite general, although some degree of internal balance will be assumed. Any of the nodes may possess a feedback mechanism. Indeed, we pay particular attention to a class of feedback queues, characterized by the fact that their state description allows one to maintain a record of the order in which events occur. We also examine in some detail the problem of determining for which values of the arrival rates a node does exhibit quasi-reversibility.

We consider the relationships among the stochastic ordering of random variables, of their random partial sums, and of the number of events of a point process in random intervals. Two types of result are obtained. Firstly, conditions are given under which a stochastic ordering between sequences of random variables is inherited by (vectors of) random partial sums of these variables. These results extend and generalize theorems known in the literature. Secondly, for the strong, (increasing) convex and (increasing) concave stochastic orderings, conditions are provided under which the numbers of events of a given point process in two ordered random intervals are also ordered.

These results are applied to some comparison problems in queueing systems. It is shown that if the service times in two M/GI/1 systems are compared in the sense of the strong stochastic ordering, or the (increasing) convex or (increasing) concave ordering, then the busy periods are compared for the same ordering. Stochastic bounds in the sense of increasing convex ordering on waiting times and on response times are provided for queues with bulk arrivals. The cyclic and Bernoulli policies for customer allocation to parallel queues are compared in the transient regime using the increasing convex ordering. Comparisons for the five above orderings are established for the cycle times in polling systems.

We consider the double queue arising from a system consisting of two processors serving three job streams generated by independent Poisson sources. The central job stream of rate v consists of jobs which place resource demands on both processors, which are handled separately by each processor once the request is made. In addition, the first processor receives background work at a rate of λwhile the second receives similar tasks at a rate η. Each processor has exponentially distributed service times with rates α and β respectively. A functional equation is found for P(z, w), the generating function of the joint queue-length distribution, which leads to a relation between P(z, 0) and P(0, w) in the region |z|, |w| < 1 of a complex algebraic curve associated with the problem. The curve is parametrized by elliptic functions z(ξ) and w(ξ) and the relation between Ρ (z(ξ), 0) and P(0, w(ξ)) persists on their analytic continuation as elliptic functions in the ξ-plane. This leads to their eventual determination by an appeal to the theory of elliptic functions. From this determination we obtain asymptotic limit laws for the expectations of the mean number of jobs in each queue conditioned on the other, as the number of jobs in both processors tends to∞. Transitions are observed in the asymptotic behavior of these quantities as one crosses various boundaries in the parameter space. An interpretation of these results via the theory of large deviations is presented.