Source author record

Johan S. H. van Leeuwaarden

Johan S. H. van Leeuwaarden appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

37works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

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

37 published item(s)

preprint2023arXiv

Robust knapsack ordering for a partially-informed newsvendor with budget constraint

This paper studies the multi-item newsvendor problem with a constrained budget and information about demand limited to its range, mean and mean absolute deviation. We consider a minimax model that determines order quantities by minimizing the expected overage and underage costs for the worst-case demand distributions. The resulting optimization problem turns out to be solvable by a method reminiscent of the greedy algorithm that solves the continuous knapsack problem, purchasing items in order of marginal value. This method has lower computational complexity compared to directly solving the model and leads to a simple policy that (i) sorts items based on their marginal effect on the total cost and (ii) determines order quantities according to this ranking until the budget is spent.

preprint2022arXiv

Heavy-Traffic Universality of Redundancy Systems with Assignment Constraints

Service systems often face task-server assignment-constraints due to skill-based routing or geographical conditions. Redundancy scheduling responds to this limited flexibility by replicating tasks to specific servers in agreement with these assignment constraints. We gain insight from product-form stationary distributions and weak local stability conditions to establish a state space collapse in heavy traffic. In this limiting regime, the parallel-server system with redundancy scheduling operates as a multi-class single-server system, achieving full resource pooling and exhibiting strong insensitivity to the underlying assignment constraints. In particular, the performance of a fully flexible (unconstrained) system can be matched even with rather strict assignment constraints.

preprint2022arXiv

Power-of-two sampling in redundancy systems: the impact of assignment constraints

A classical sampling strategy for load balancing policies is power-of-two, where any server pair is sampled with equal probability. This does not cover practical settings with assignment constraints which force non-uniform sampling. While intuition suggests that non-uniform sampling adversely impacts performance, this was only supported through simulations, and rigorous statements have remained elusive. Building on product-form distributions for redundancy systems, we prove the stochastic dominance of uniform sampling for a four-server system as well as arbitrary-size systems in light traffic.

preprint2021arXiv

Scalable load balancing in networked systems: A survey of recent advances

The basic load balancing scenario involves a single dispatcher where tasks arrive that must immediately be forwarded to one of $N$ single-server queues. We discuss recent advances on scalable load balancing schemes which provide favorable delay performance when $N$ grows large, and yet only require minimal implementation overhead. Join-the-Shortest-Queue (JSQ) yields vanishing delays as $N$ grows large, as in a centralized queueing arrangement, but involves a prohibitive communication burden. In contrast, power-of-$d$ or JSQ($d$) schemes that assign an incoming task to a server with the shortest queue among $d$ servers selected uniformly at random require little communication, but lead to constant delays. In order to examine this fundamental trade-off between delay performance and implementation overhead, we consider JSQ($d(N)$) schemes where the diversity parameter $d(N)$ depends on $N$ and investigate what growth rate of $d(N)$ is required to asymptotically match the optimal JSQ performance on fluid and diffusion scale. Stochastic coupling techniques and stochastic-process limits play an instrumental role in establishing the asymptotic optimality. We demonstrate how this methodology carries over to infinite-server settings, finite buffers, multiple dispatchers, servers arranged on graph topologies, and token-based load balancing including the popular Join-the-Idle-Queue (JIQ) scheme. In this way we provide a broad overview of the many recent advances in the field. This survey extends the short review presented at ICM 2018 (arXiv:1712.08555).

preprint2020arXiv

Critical percolation on scale-free random graphs: New universality class for the configuration model

In this paper, we study the critical behavior of percolation on a configuration model with degree distribution satisfying an infinite second-moment condition, which includes power-law degrees with exponent $τ\in (2,3)$. It is well known that, in this regime, many canonical random graph models, such as the configuration model, are robust in the sense that the giant component is not destroyed when the percolation probability stays bounded away from zero. Thus, the critical behavior is observed when the percolation probability tends to zero with the network size, despite of the fact that the average degree remains bounded. In this paper, we initiate the study of critical random graphs in the infinite second-moment regime by identifying the critical window for the configuration model. We prove scaling limits for component sizes and surplus edges, and show that the maximum diameter the critical components is of order $\log n$, which contrasts with the previous universality classes arising in the literature. This introduces a third and novel universality class for the critical behavior of percolation on random networks, that is not covered by the multiplicative coalescent framework due to Aldous and Limic (1998). We also prove concentration of the component sizes outside the critical window, and that a unique, complex giant component emerges after the critical window. This completes the picture for the percolation phase transition on the configuration model.

preprint2020arXiv

Optimal subgraph structures in scale-free configuration models

Subgraphs reveal information about the geometry and functionalities of complex networks. For scale-free networks with unbounded degree fluctuations, we obtain the asymptotics of the number of times a small connected graph occurs as a subgraph or as an induced subgraph. We obtain these results by analyzing the configuration model with degree exponent $τ\in(2,3)$ and introducing a novel class of optimization problems. For any given subgraph, the unique optimizer describes the degrees of the vertices that together span the subgraph. We find that subgraphs typically occur between vertices with specific degree ranges. In this way, we can count and characterize {\it all} subgraphs. We refrain from double counting in the case of multi-edges, essentially counting the subgraphs in the {\it erased} configuration model.

preprint2016arXiv

Asymptotic Optimality of Power-of-$d$ Load Balancing in Large-Scale Systems

We consider a system of $N$ identical server pools and a single dispatcher where tasks arrive as a Poisson process of rate $λ(N)$. Arriving tasks cannot be queued, and must immediately be assigned to one of the server pools to start execution, or discarded. The execution times are assumed to be exponentially distributed with unit mean, and do not depend on the number of other tasks receiving service. However, the experienced performance (e.g. in terms of received throughput) does degrade with an increasing number of concurrent tasks at the same server pool. The dispatcher therefore aims to evenly distribute the tasks across the various server pools. Specifically, when a task arrives, the dispatcher assigns it to the server pool with the minimum number of tasks among $d(N)$ randomly selected server pools. This assignment strategy is called the JSQ$(d(N))$ scheme, as it resembles the power-of-$d$ version of the Join-the-Shortest-Queue (JSQ) policy, and will also be referred to as such in the special case $d(N) = N$. We construct a stochastic coupling to bound the difference in the system occupancy processes between the JSQ policy and a scheme with an arbitrary value of $d(N)$. We use the coupling to derive the fluid limit in case $d(N) \to \infty$ and $λ(N)/N \to λ$ as $N \to \infty$, along with the associated fixed point. The fluid limit turns out to be insensitive to the exact growth rate of $d(N)$, and coincides with that for the JSQ policy. We further leverage the coupling to establish that the diffusion limit corresponds to that for the JSQ policy as well, as long as $d(N)/\sqrt{N} \log(N) \to \infty$, and characterize the common limiting diffusion process. These results indicate that the JSQ optimality can be preserved at the fluid-level and diffusion-level while reducing the overhead by nearly a factor O($N$) and O($\sqrt{N}/\log(N)$), respectively.

preprint2016arXiv

Epidemic spreading on complex networks with community structures

Many real-world networks display a community structure. We study two random graph models that create a network with similar community structure as a given network. One model preserves the exact community structure of the original network, while the other model only preserves the set of communities and the vertex degrees. These models show that community structure is an important determinant of the behavior of percolation processes on networks, such as information diffusion or virus spreading: the community structure can both \textit{enforce} as well as \textit{inhibit} diffusion processes. Our models further show that it is the mesoscopic set of communities that matters. The exact internal structures of communities barely influence the behavior of percolation processes across networks. This insensitivity is likely due to the relative denseness of the communities.

preprint2016arXiv

Finite-pool queues with heavy-tailed services

We consider the $Δ_{(i)}/G/1$ queue, in which a a total of $n$ customers independently demand service after an exponential time. We focus on the case of heavy-tailed service times, and assume that the tail of the service time distribution decays like $x^{-α}$, with $α\in (1,2)$. We consider the asymptotic regime in which the population size grows to infinity and establish that the scaled queue length process converges to an α-stable process with a negative quadratic drift. We leverage this asymptotic result to characterize the headstart that is needed to create a long period of activity. This result should be contrasted with the case of light-tailed service times, which was shown to have a similar scaling limit, but then with a Brownian motion instead of an α-stable process.

preprint2016arXiv

Generalized random sequential adsorption on Erdős-Rényi random graphs

We investigate Random Sequential Adsorption (RSA) on a random graph via the following greedy algorithm: Order the $n$ vertices at random, and sequentially declare each vertex either active or frozen, depending on some local rule in terms of the state of the neighboring vertices. The classical RSA rule declares a vertex active if none of its neighbors is, in which case the set of active nodes forms an independent set of the graph. We generalize this nearest-neighbor blocking rule in three ways and apply it to the Erdős-Rényi random graph. We consider these generalizations in the large-graph limit $n\to\infty$ and characterize the jamming constant, the limiting proportion of active vertices in the maximal greedy set.

preprint2016arXiv

Hierarchical Configuration Model

We introduce a class of random graphs with a community structure, which we call the hierarchical configuration model. On the inter-community level, the graph is a configuration model, and on the intra-community level, every vertex in the configuration model is replaced by a community: i.e., a small graph. These communities may have any shape, as long as they are connected. For these hierarchical graphs, we find the size of the largest component, the degree distribution and the clustering coefficient. Furthermore, we determine the conditions under which a giant percolation cluster exists, and find its size.

preprint2016arXiv

Power-law relations in random networks with communities

Most random graph models are locally tree-like - do not contain short cycles - rendering them unfit for modeling networks with a community structure. We introduce the hierarchical configuration model (HCM), a generalization of the configuration model that includes community structures, while properties such as the size of the giant component, and the size of the giant percolating cluster under bond percolation can still be derived analytically. Viewing real-world networks as realizations of the HCM, we observe two previously unobserved power-law relations: between the number of edges inside a community and the community sizes, and between the number of edges going out of a community and the community sizes. We also relate the power-law exponent $τ$ of the degree distribution with the power-law exponent of the community size distribution $γ$. In the special case of extremely dense communities (e.g., complete graphs), this relation takes the simple form $τ=γ-1$.

preprint2016arXiv

The restricted Erlang-R Queue: Finite-size effects in service systems with returning customers

Motivated by health care systems with repeated services that have both personnel (nurse and physician) and space (beds) constraints, we study a restricted version of the Erlang-R model. The space restriction policies we consider are blocking or holding in a pre-entrant queue. We develop many-server approximations for the system performance measures when either policy applies, and explore the connection between them. We show that capacity allocation of both resources should be determined simultaneously, and derive the methodology to determine it explicitly. We show that the system dynamics is captured by the fraction of needy time in the network, and that returning customers should be accounted for both in steady-state and time-varying conditions. We demonstrate the application of our policies in two case studies of resource allocation in hospitals.

preprint2015arXiv

Dominant poles and tail asymptotics in the critical Gaussian many-sources regime

The dominant pole approximation (DPA) is a classical analytic method to obtain from a generating function asymptotic estimates for its underlying coefficients. We apply DPA to a discrete queue in a critical many-sources regime, in order to obtain tail asymptotics for the stationary queue length. As it turns out, this regime leads to a clustering of the poles of the generating function, which renders the classical DPA useless, since the dominant pole is not sufficiently dominant. To resolve this, we design a new DPA method, which might also find application in other areas of mathematics, like combinatorics, particularly when Gaussian scalings related to the central limit theorem are involved.

preprint2015arXiv

Giant component sizes in scale-free networks with power-law degrees and cutoffs

Scale-free networks arise from power-law degree distributions. Due to the finite size of real-world networks, the power law inevitably has a cutoff at some maximum degree $Δ$. We investigate the relative size of the giant component $S$ in the large-network limit. We show that $S$ as a function of $Δ$ increases fast when $Δ$ is just large enough for the giant component to exist, but increases ever more slowly when $Δ$ increases further. This makes that while the degree distribution converges to a pure power law when $Δ\to\infty$, $S$ approaches its limiting value at a slow pace. The convergence rate also depends on the power-law exponent $τ$ of the degree distribution. The worst rate of convergence is found to be for the case $τ\approx2$, which concerns many of the real-world networks reported in the literature.

preprint2015arXiv

Heavy-traffic analysis through uniform acceleration of queues with diminishing populations

We consider a single server queue that serves a finite population of $n$ customers that will enter the queue (require service) only once, also known as the $Δ_{(i)}/G/1$ queue. This paper presents a method for analyzing heavy-traffic behavior by using uniform acceleration, which simultaneously lets $n$ and the service rate grow large, while the initial resource utilization approaches one. A key feature of the model is that, as time progresses, more customers have joined the queue, and fewer customers can potentially join. This diminishing population gives rise to a class of reflected stochastic processes that vanish over time, and hence do not have a stationary distribution. We establish that, when the arrival times are exponentially distributed, by suitably rescaling space and time, the queue length process converges to a Brownian motion with parabolic drift, a stochastic-process limit that captures the effect of a diminishing population by a negative quadratic drift. When the arrival times are generally distributed, our techniques provide information on the typical queue length and the first busy period.

preprint2015arXiv

Product-form solutions for a class of structured multi-dimensional Markov processes

Motivated by queueing systems with heterogeneous parallel servers, we consider a class of structured multi-dimensional Markov processes whose state space can be partitioned into two parts: a finite set of boundary states and a structured multi-dimensional set of states, exactly one dimension of which is infinite. Using separation of variables, we show that the equilibrium distribution, typically of the queue length, can be represented as a linear combination of product forms. For an important subclass of queueing systems, we characterize explicitly the waiting time distribution in terms of mixtures of exponentials.

preprint2014arXiv

Cluster tails for critical power-law inhomogeneous random graphs

Recently, the scaling limit of cluster sizes for critical inhomogeneous random graphs of rank-1 type having finite variance but infinite third moment degrees was obtained (see previous work by Bhamidi, van der Hofstad and van Leeuwaarden). It was proved that when the degrees obey a power law with exponent in the interval (3,4), the sequence of clusters ordered in decreasing size and scaled appropriately converges as n goes to infinity to a sequence of decreasing non-degenerate random variables. Here, we study the tails of the limit of the rescaled largest cluster, i.e., the probability that the scaling limit of the largest cluster takes a large value u, as a function of u. This extends a related result of Pittel for the Erdős-Rényi random graph to the setting of rank-1 inhomogeneous random graphs with infinite third moment degrees. We make use of delicate large deviations and weak convergence arguments.

preprint2014arXiv

Large deviations for power-law thinned Levy processes

This paper deals with the large deviations behavior of a stochastic process called thinned Levy process. This process appeared recently as a stochastic-process limit in the context of critical inhomogeneous random graphs. The process has a strong negative drift, while we are interested in the rare event of the process being positive at large times. To characterize this rare event, we identify a tilted measure. This presents some challenges inherent to the power-law nature of the thinned Levy process. General principles prescribe that the tilt should follow from a variational problem, but in the case of the thinned Levy process this involves a Riemann sum that is hard to control. We choose to approximate the Riemann sum by its limiting integral, derive the first-order correction term, and prove that the tilt that follows from the corresponding approximate variational problem is sufficient to establish the large deviations results.

preprint2014arXiv

Novel heavy-traffic regimes for large-scale service systems

We introduce a family of heavy-traffic regimes for large scale service systems, presenting a range of scalings that include both moderate and extreme heavy traffic, as compared to classical heavy traffic. The heavy-traffic regimes can be translated into capacity sizing rules that lead to Economies-of-Scales, so that the system utilization approaches 100% while congestion remains limited. We obtain heavy-traffic approximations for stationary performance measures in terms of asymptotic expansions, using a non-standard saddle point method, tailored to the specific form of integral expressions for the performance measures, in combination with the heavy-traffic regimes.

preprint2014arXiv

Slow transitions, slow mixing and starvation in dense random-access networks

We consider dense wireless random-access networks, modeled as systems of particles with hard-core interaction. The particles represent the network users that try to become active after an exponential back-off time, and stay active for an exponential transmission time. Due to wireless interference, active users prevent other nearby users from simultaneous activity, which we describe as hard-core interaction on a conflict graph. We show that dense networks with aggressive back-off schemes lead to extremely slow transitions between dominant states, and inevitably cause long mixing times and starvation effects.

preprint2014arXiv

Staffing many-server systems with admission control and retrials

In many-server systems it is crucial to staff the right number of servers so that targeted service levels are met. These staffing problems typically lead to constraint satisfaction problems that are hard to solve. During the last decade, a powerful many-server asymptotic theory has been developed to solve such problems and optimal staffing rules are known to obey the square-root staffing principle. This paper develops many-server asymptotics in the so-called QED regime, and presents refinements to many-server asymptotics and square-root staffing for a Markovian queueing model with admission control and retrials.

preprint2014arXiv

Transient analysis of the Erlang A model

We consider the Erlang A model, or $M/M/m+M$ queue, with Poisson arrivals, exponential service times, and $m$ parallel servers, and the property that waiting customers abandon the queue after an exponential time. The queue length process is in this case a birth-death process, for which we obtain explicit expressions for the Laplace transforms of the time-dependent distribution and the first passage time. These two transient characteristics were generally presumed to be intractable. Solving for the Laplace transforms involves using Green's functions and contour integrals related to hypergeometric functions. Our results are specialized to the $M/M/\infty$ queue, the $M/M/m$ queue, and the $M/M/m/m$ loss model. We also obtain some corresponding results for diffusion approximations to these models.

preprint2013arXiv

BRAVO for many-server QED systems with finite buffers

This paper demonstrates the occurrence of the feature called BRAVO (Balancing Reduces Asymptotic Variance of Output) for the departure process of a finite-buffer Markovian many-server system in the QED (Quality and Efficiency-Driven) heavy-traffic regime. The results are based on evaluating the limit of a formula for the asymptotic variance of death counts in finite birth--death processes.

preprint2013arXiv

Delay performance in random-access grid networks

We examine the impact of torpid mixing and meta-stability issues on the delay performance in wireless random-access networks. Focusing on regular meshes as prototypical scenarios, we show that the mean delays in an $L\times L$ toric grid with normalized load $ρ$ are of the order $(\frac{1}{1-ρ})^L$. This superlinear delay scaling is to be contrasted with the usual linear growth of the order $\frac{1}{1-ρ}$ in conventional queueing networks. The intuitive explanation for the poor delay characteristics is that (i) high load requires a high activity factor, (ii) a high activity factor implies extremely slow transitions between dominant activity states, and (iii) slow transitions cause starvation and hence excessively long queues and delays. Our proof method combines both renewal and conductance arguments. A critical ingredient in quantifying the long transition times is the derivation of the communication height of the uniformized Markov chain associated with the activity process. We also discuss connections with Glauber dynamics, conductance and mixing times. Our proof framework can be applied to other topologies as well, and is also relevant for the hard-core model in statistical physics and the sampling from independent sets using single-site update Markov chains.

preprint2013arXiv

Factorization identities for reflected processes, with applications

We derive factorization identities for a class of preemptive-resume queueing systems, with batch arrivals and catastrophes that, whenever they occur, eliminate multiple customers present in the system. These processes are quite general, as they can be used to approximate Levy processes, diffusion processes, and certain types of growth-collapse processes; thus, all of the processes mentioned above also satisfy similar factorization identities. In the Levy case, our identities simplify to both the well-known Wiener-Hopf factorization, and another interesting factorization of reflected Levy processes starting at an arbitrary initial state. We also show how the ideas can be used to derive transforms for some well-known state-dependent/inhomogeneous birth-death processes and diffusion processes.

preprint2013arXiv

First passage times to congested states of many-server systems in the Halfin-Whitt regime

We consider the heavy-traffic approximation to the $GI/M/s$ queueing system in the Halfin-Whitt regime, where both the number of servers $s$ and the arrival rate $λ$ grow large (taking the service rate as unity), with $λ=s-β\sqrt{s}$ and $β$ some constant. In this asymptotic regime, the queue length process can be approximated by a diffusion process that behaves like a Brownian motion with drift above zero and like an Ornstein-Uhlenbeck process below zero. We analyze the first passage times of this hybrid diffusion process to levels in the state space that represent congested states in the original queueing system.

preprint2013arXiv

On the Sojourn Time Distribution in a Finite Population Markovian Processor Sharing Queue

We consider a finite population processor-sharing (PS) queue, with Markovian arrivals and an exponential server. Such a queue can model an interactive computer system consisting of a bank of terminals in series with a central processing unit (CPU). For systems with a large population $N$ and a commensurately rapid service rate, or infrequent arrivals, we obtain various asymptotic results. We analyze the conditional sojourn time distribution of a tagged customer, conditioned on the number $n$ of others in the system at the tagged customer's arrival instant, and also the unconditional distribution. The asymptotics are obtained by a combination of singular perturbation methods and spectral methods. We consider several space/time scales and parameter ranges, which lead to different asymptotic behaviors. We also identify precisely when the finite population model can be approximated by the standard infinite population $M/M/1$-PS queue.

preprint2013arXiv

Queues with random back-offs

We consider a broad class of queueing models with random state-dependent vacation periods, which arise in the analysis of queue-based back-off algorithms in wireless random-access networks. In contrast to conventional models, the vacation periods may be initiated after each service completion, and can be randomly terminated with certain probabilities that depend on the queue length. We examine the scaled queue length and delay in a heavy-traffic regime, and demonstrate a sharp trichotomy, depending on how the activation rate and vacation probability behave as function of the queue length. In particular, the effect of the vacation periods may either (i) completely vanish in heavy-traffic conditions, (ii) contribute an additional term to the queue lengths and delays of similar magnitude, or even (iii) give rise to an order-of-magnitude increase. The heavy-traffic asymptotics are obtained by combining stochastic lower and upper bounds with exact results for some specific cases. The heavy-traffic trichotomy provides valuable insight in the impact of the back-off algorithms on the delay performance in wireless random-access networks.

preprint2013arXiv

Spectral gap of the Erlang A model in the Halfin-Whitt regime

We consider a hybrid diffusion process that is a combination of two Ornstein-Uhlenbeck processes with different restraining forces. This process serves as the heavy-traffic approximation to the Markovian many-server queue with abandonments in the critical Halfin-Whitt regime. We obtain an expression for the Laplace transform of the time-dependent probability distribution, from which the spectral gap is explicitly characterized. The spectral gap gives the exponential rate of convergence to equilibrium. We further give various asymptotic results for the spectral gap, in the limits of small and large abandonment effects. It turns out that convergence to equilibrium becomes extremely slow for overloaded systems with small abandonment effects.

preprint2012arXiv

Achievable Performance in Product-Form Networks

We characterize the achievable range of performance measures in product-form networks where one or more system parameters can be freely set by a network operator. Given a product-form network and a set of configurable parameters, we identify which performance measures can be controlled and which target values can be attained. We also discuss an online optimization algorithm, which allows a network operator to set the system parameters so as to achieve target performance metrics. In some cases, the algorithm can be implemented in a distributed fashion, of which we give several examples. Finally, we give conditions that guarantee convergence of the algorithm, under the assumption that the target performance metrics are within the achievable range.

preprint2012arXiv

Mixing Properties of CSMA Networks on Partite Graphs

We consider a stylized stochastic model for a wireless CSMA network. Experimental results in prior studies indicate that the model provides remarkably accurate throughput estimates for IEEE 802.11 systems. In particular, the model offers an explanation for the severe spatial unfairness in throughputs observed in such networks with asymmetric interference conditions. Even in symmetric scenarios, however, it may take a long time for the activity process to move between dominant states, giving rise to potential starvation issues. In order to gain insight in the transient throughput characteristics and associated starvation effects, we examine in the present paper the behavior of the transition time between dominant activity states. We focus on partite interference graphs, and establish how the magnitude of the transition time scales with the activation rate and the sizes of the various network components. We also prove that in several cases the scaled transition time has an asymptotically exponential distribution as the activation rate grows large, and point out interesting connections with related exponentiality results for rare events and meta-stability phenomena in statistical physics. In addition, we investigate the convergence rate to equilibrium of the activity process in terms of mixing times.

preprint2012arXiv

Novel scaling limits for critical inhomogeneous random graphs

We find scaling limits for the sizes of the largest components at criticality for rank-1 inhomogeneous random graphs with power-law degrees with power-law exponent τ. We investigate the case where $τ\in(3,4)$, so that the degrees have finite variance but infinite third moment. The sizes of the largest clusters, rescaled by $n^{-(τ-2)/(τ-1)}$, converge to hitting times of a "thinned" Lévy process, a special case of the general multiplicative coalescents studied by Aldous [Ann. Probab. 25 (1997) 812-854] and Aldous and Limic [Electron. J. Probab. 3 (1998) 1-59]. Our results should be contrasted to the case τ>4, so that the third moment is finite. There, instead, the sizes of the components rescaled by $n^{-2/3}$ converge to the excursion lengths of an inhomogeneous Brownian motion, as proved in Aldous [Ann. Probab. 25 (1997) 812-854] for the Erdős-Rényi random graph and extended to the present setting in Bhamidi, van der Hofstad and van Leeuwaarden [Electron. J. Probab. 15 (2010) 1682-1703] and Turova [(2009) Preprint].

preprint2012arXiv

Online Optimization of Product-Form Networks

We develop an online gradient algorithm for optimizing the performance of product-form networks through online adjustment of control parameters. The use of standard algorithms for finding optimal parameter settings is hampered by the prohibitive computational burden of calculating the gradient in terms of the stationary probabilities. The proposed approach instead relies on measuring empirical frequencies of the various states through simulation or online operation so as to obtain estimates for the gradient. Besides the reduction in computational effort, a further benefit of the online operation lies in the natural adaptation to slow variations in ambient parameters as commonly occurring in dynamic environments. On the downside, the measurements result in inherently noisy and biased estimates. We exploit mixing time results in order to overcome the impact of the bias and establish sufficient conditions for convergence to a globally optimal solution.

preprint2012arXiv

Random walks reaching against all odds the other side of the quarter plane

For a homogeneous random walk in the quarter plane with nearest-neighbor transitions, starting from some state $(i_0,j_0)$, we study the event that the walk reaches the vertical axis, before reaching the horizontal axis. We derive an exact expression for the probability of this event, and derive an asymptotic expression for the case when $i_0$ becomes large, a situation in which the event becomes highly unlikely. The exact expression follows from the solution of a boundary value problem and is in terms of an integral that involves a conformal gluing function. The asymptotic expression follows from the asymptotic evaluation of this integral. Our results find applications in a model for nucleosome shifting, the voter model and the asymmetric exclusion process.

preprint2012arXiv

The compensation approach for walks with small steps in the quarter plane

This paper is the first application of the compensation approach to counting problems. We discuss how this method can be applied to a general class of walks in the quarter plane $Z_{+}^{2}$ with a step set that is a subset of $\{(-1,1),(-1,0),(-1,-1),(0,-1),(1,-1)\}$ in the interior of $Z_{+}^{2}$. We derive an explicit expression for the counting generating function, which turns out to be meromorphic and nonholonomic, can be easily inverted, and can be used to obtain asymptotic expressions for the counting coefficients.

preprint2009arXiv

On a Processor Sharing Queue That Models Balking

We consider the processor sharing $M/M/1$-PS queue which also models balking. A customer that arrives and sees $n$ others in the system "balks" (i.e., decides not to enter) with probability $1-b_n$. If $b_n$ is inversely proportional to $n+1$, we obtain explicit expressions for a tagged customer's sojourn time distribution. We consider both the conditional distribution, conditioned on the number of other customers present when the tagged customer arrives, as well as the unconditional distribution. We then evaluate the results in various asymptotic limits. These include large time (tail behavior) and/or large $n$, lightly loaded systems where the arrival rate $λ\to 0$, and heavily loaded systems where $λ\to\infty$. We find that the asymptotic structure for the problem with balking is much different from the standard $M/M/1$-PS queue. We also discuss a perturbation method for deriving the asymptotics, which should apply to more general balking functions.