Centrality measures aim to indicate who is important in a network. Various notions of ‘being important’ give rise to different centrality measures. In this paper, we study how important the central vertices are for the connectivity structure of the network, by investigating how the removal of the most central vertices affects the number of connected components and the size of the giant component. We use local convergence techniques to identify the limiting number of connected components for locally converging graphs and centrality measures that depend on the vertex’s neighbourhood. For the size of the giant, we prove a general upper bound. For the matching lower bound, we specialise to the case of degree centrality on one of the most popular models in network science, the configuration model, for which we show that removal of the highest-degree vertices destroys the giant most.

Let $f(z)=\sum\limits_{j=0}^{\infty} a_j z^j$ be a transcendental entire function and let $f_\omega(z)=\sum\limits_{j=0}^{\infty}\chi_j(\omega) a_j z^j$ be a random entire function, where $\chi_j(\omega)$ are independent and identically distributed random variables defined on a probability space $(\Omega, \mathcal{F}, \mu)$. In this paper, we first define a family of random entire functions, which includes Gaussian, Rademacher and Steinhaus entire functions. We prove that, for almost all functions in the family and for any constant C > 1, there exist a constant $r_0=r_0(\omega)$ and a set $E\subset [e, \infty)$ of finite logarithmic measure such that, for $r \gt r_0$ and $r\notin E$,

\begin{equation*}|\log M(r, f)- N(r,0, f_\omega)|\le (C/A)^{\frac1{B}}\,\log^{\frac1{B}}\,\log M(r,f) +\log\,\log M(r, f), \qquad a.s.\end{equation*}

where $A, B$ are constants, $M(r, f)$ is the maximum modulus and $N(r, 0, f)$ is the integrated zero-counting function of f. As a by-product of our main results, we prove Nevanlinna’s second main theorem for random entire functions. Thus, the characteristic function of almost all functions in the family is bounded above by an integrated counting function, rather than by two integrated counting functions as in the classical Nevanlinna theory. For instance, we show that, for almost all Gaussian entire functions fω and for any ϵ > 0, there is r0 such that, for $r \gt r_0$,

\begin{equation*}T(r, f) \le N(r,0, f_\omega)+\left(\tfrac12+\epsilon\right) \log T(r, f).\end{equation*}

Taylor’s power law (or fluctuation scaling) states that on comparable populations, the variance of each sample is approximately proportional to a power of the mean of the population. The law has been shown to hold by empirical observations in a broad class of disciplines including demography, biology, economics, physics, and mathematics. In particular, it has been observed in problems involving population dynamics, market trading, thermodynamics, and number theory. In applications, many authors consider panel data in order to obtain laws of large numbers. Essentially, we aim to consider ergodic behaviors without independence. We restrict our study to stationary time series, and develop different Taylor exponents in this setting. From a theoretical point of view, there has been a growing interest in the study of the behavior of such a phenomenon. Most of these works focused on the so-called static Taylor’s law related to independent samples. In this paper we introduce a dynamic Taylor’s law for dependent samples using self-normalized expressions involving Bernstein blocks. A central limit theorem (CLT) is proved under either weak dependence or strong mixing assumptions for the marginal process. The limit behavior of the estimation involves a series of covariances, unlike the classic framework where the limit behavior involves the marginal variance. We also provide an asymptotic result for a goodness-of-fit procedure suitable for checking whether the corresponding dynamic Taylor’s law holds in empirical studies.

The logistic birth and death process is perhaps the simplest stochastic population model that has both density-dependent reproduction and a phase transition, and a lot can be learned about the process by studying its extinction time, $\tau_n$ , as a function of system size n. A number of existing results describe the scaling of $\tau_n$ as $n\to\infty$ for various choices of reproductive rate $r_n$ and initial population $X_n(0)$ as a function of n. We collect and complete this picture, obtaining a complete classification of all sequences $(r_n)$ and $(X_n(0))$ for which there exist rescaling parameters $(s_n)$ and $(t_n)$ such that $(\tau_n-t_n)/s_n$ converges in distribution as $n\to\infty$ , and identifying the limits in each case.

Let A be a finite set with , let n be a positive integer, and let $A^n$ denote the discrete $n\text {-dimensional}$ hypercube (that is, $A^n$ is the Cartesian product of n many copies of A). Given a family $\langle D_t:t\in A^n\rangle $ of measurable events in a probability space (a stochastic process), what structural information can be obtained assuming that the events $\langle D_t:t\in A^n\rangle $ are not behaving as if they were independent? We obtain an answer to this problem (in a strong quantitative sense) subject to a mild ‘stationarity’ condition. Our result has a number of combinatorial consequences, including a new (and the most informative so far) proof of the density Hales-Jewett theorem.

It is well understood that a supercritical superprocess is equal in law to a discrete Markov branching process whose genealogy is dressed in a Poissonian way with immigration which initiates subcritical superprocesses. The Markov branching process corresponds to the genealogical description of prolific individuals, that is, individuals who produce eternal genealogical lines of descent, and is often referred to as the skeleton or backbone of the original superprocess. The Poissonian dressing along the skeleton may be considered to be the remaining non-prolific genealogical mass in the superprocess. Such skeletal decompositions are equally well understood for continuous-state branching processes (CSBP).

In a previous article [16] we developed an SDE approach to study the skeletal representation of CSBPs, which provided a common framework for the skeletal decompositions of supercritical and (sub)critical CSBPs. It also helped us to understand how the skeleton thins down onto one infinite line of descent when conditioning on survival until larger and larger times, and eventually forever.

Here our main motivation is to show the robustness of the SDE approach by expanding it to the spatial setting of superprocesses. The current article only considers supercritical superprocesses, leaving the subcritical case open.

We construct a multitype constant-size population model allowing for general selective interactions as well as extreme reproductive events. Our multidimensional model aims for the generality of adaptive dynamics and the tractability of population genetics. It generalises the idea of Krone and Neuhauser [39] and González Casanova and Spanò [29], who represented the selection by allowing individuals to sample several potential parents in the previous generation before choosing the ‘strongest’ one, by allowing individuals to use any rule to choose their parent. The type of the newborn can even not be one of the types of the potential parents, which allows modelling mutations. Via a large population limit, we obtain a generalisation of $\Lambda$ -Fleming–Viot processes, with a diffusion term and a general frequency-dependent selection, which allows for non-transitive interactions between the different types present in the population. We provide some properties of these processes related to extinction and fixation events, and give conditions for them to be realised as unique strong solutions of multidimensional stochastic differential equations with jumps. Finally, we illustrate the generality of our model with applications to some classical biological interactions. This framework provides a natural bridge between two of the most prominent modelling frameworks of biological evolution: population genetics and eco-evolutionary models.

Let V be an n-set, and let X be a random variable taking values in the power-set of V. Suppose we are given a sequence of random coupons $X_1, X_2, \ldots $ , where the $X_i$ are independent random variables with distribution given by X. The covering time T is the smallest integer $t\geq 0$ such that $\bigcup_{i=1}^t X_i=V$ . The distribution of T is important in many applications in combinatorial probability, and has been extensively studied. However the literature has focused almost exclusively on the case where X is assumed to be symmetric and/or uniform in some way.

In this paper we study the covering time for much more general random variables X; we give general criteria for T being sharply concentrated around its mean, precise tools to estimate that mean, as well as examples where T fails to be concentrated and when structural properties in the distribution of X allow for a very different behaviour of T relative to the symmetric/uniform case.

In this paper, we study the nonlinear diffusion equation associated with a particle system where the common drift depends on the rate of absorption of particles at a boundary. We provide an interpretation of this equation, which is also related to the supercooled Stefan problem, as a structural credit risk model with default contagion in a large interconnected banking system. Using the method of heat potentials, we derive a coupled system of Volterra integral equations for the transition density and for the loss through absorption. An approximation by expansion is given for a small interaction parameter. We also present a numerical solution algorithm and conduct computational tests.

We consider the behaviour of minimax recursions defined on random trees. Such recursions give the value of a general class of two-player combinatorial games. We examine in particular the case where the tree is given by a Galton–Watson branching process, truncated at some depth 2n, and the terminal values of the level 2n nodes are drawn independently from some common distribution. The case of a regular tree was previously considered by Pearl, who showed that as n → ∞ the value of the game converges to a constant, and by Ali Khan, Devroye and Neininger, who obtained a distributional limit under a suitable rescaling.

For a general offspring distribution, there is a surprisingly rich variety of behaviour: the (unrescaled) value of the game may converge to a constant, or to a discrete limit with several atoms, or to a continuous distribution. We also give distributional limits under suitable rescalings in various cases.

We also address questions of endogeny. Suppose the game is played on a tree with many levels, so that the terminal values are far from the root. To be confident of playing a good first move, do we need to see the whole tree and its terminal values, or can we play close to optimally by inspecting just the first few levels of the tree? The answers again depend in an interesting way on the offspring distribution.

We also mention several open questions.

In this paper, we study limit behavior for a Markov-modulated binomial counting process, also called a binomial counting process under regime switching. Such a process naturally appears in the context of credit risk when multiple obligors are present. Markov-modulation takes place when the failure/default rate of each individual obligor depends on an underlying Markov chain. The limit behavior under consideration occurs when the number of obligors increases unboundedly, and/or by accelerating the modulating Markov process, called rapid switching. We establish diffusion approximations, obtained by application of (semi)martingale central limit theorems. Depending on the specific circumstances, different approximations are found.

We extend the work of Antunović et al. (2016) on competing types in preferential attachment models to include cases where the types have different fitnesses, which may be either multiplicative or additive. We show that, depending on the values of the parameters of the models, there are different possible limiting behaviours depending on the zeros of a certain function. In particular, we show the existence of choices of the parameters where one type is favoured both by having higher fitness and by the type of attachment mechanism, but the other type has a positive probability of dominating the network in the limit.

In the paper we consider the following modification of a discrete-time branching process with stationary immigration. In each generation a binomially distributed subset of the population will be observed. The number of observed individuals constitute a partially observed branching process. After inspection both observed and unobserved individuals may change their offspring distributions. In the subcritical case we investigate the possibility of using the known estimators for the offspring mean and for the mean of the stationary-limiting distribution of the process when the observation of the population sizes is restricted. We prove that, if both the population and the number of immigrants are partially observed, the estimators are still strongly consistent. We also prove that the `skipped' version of the estimator for the offspring mean is asymptotically normal and the estimator of the stationary distribution's mean is asymptotically normal under additional assumptions.

We study the exciton diffusion in organic semiconductors from a macroscopic viewpoint. In a unified way, we conduct the equivalence analysis between Monte-Carlo method and diffusion equation model for photoluminescence quenching and photocurrent spectrum measurements, in both the presence and the absence of Förster energy transfer effect. Connections of these two models to Stern-Volmer method and exciton-exciton annihilation method are also specified for the photoluminescence quenching measurement.

We introduce a new model of competition on growing networks. This extends the preferential attachment model, with the key property that node choices evolve simultaneously with the network. When a new node joins the network, it chooses neighbours by preferential attachment, and selects its type based on the number of initial neighbours of each type. The model is analysed in detail, and in particular, we determine the possible proportions of the various types in the limit of large networks. An important qualitative feature we find is that, in contrast to many current theoretical models, often several competitors will coexist. This matches empirical observations in many real-world networks.

Retransmission-based failure recovery represents a primary approach in existing communication networks that guarantees data delivery in the presence of channel failures. Recent work has shown that, when data sizes have infinite support, retransmissions can cause long (-tailed) delays even if all traffic and network characteristics are light-tailed. In this paper we investigate the practically important case of bounded data units 0 ≤ Lb ≤ b under the condition that the hazard functions of the distributions of data sizes and channel statistics are proportional. To this end, we provide an explicit and uniform characterization of the entire body of the retransmission distribution ℙ[Nb > n] in both n and b. Our main discovery is that this distribution can be represented as the product of a power law and gamma distribution. This rigorous approximation clearly demonstrates the coupling of a power law distribution, dominating the main body, and the gamma distribution, determining the exponential tail. Our results are validated via simulation experiments and can be useful for designing retransmission-based systems with the required performance characteristics. From a broader perspective, this study applies to any other system, e.g. computing, where restart mechanisms are employed after a job processing failure.

Identifiability of evolutionary tree models has been a recent topic of discussion and some models have been shown to be nonidentifiable. A coalescent-based rooted population tree model, originally proposed by Nielsen et al. (1998), has been used by many authors in the last few years and is a simple tool to accurately model the changes in allele frequencies in the tree. However, the identifiability of this model has never been proven. Here we prove this model to be identifiable by showing that the model parameters can be expressed as functions of the probability distributions of subsamples, assuming that there are at least two (haploid) individuals sampled from each population. This a step toward proving the consistency of the maximum likelihood estimator of the population tree based on this model.

The d-dimensional Λ-Fleming-Viot generator acting on functions g(x), with x being a vector of d allele frequencies, can be written as a Wright-Fisher generator acting on functions g with a modified random linear argument of x induced by partitioning occurring in the Λ-Fleming-Viot process. The eigenvalues and right polynomial eigenvectors are easy to see from this representation. The two-dimensional process, which has a one-dimensional generator, is considered in detail. A nonlinear equation is found for the Green's function. In a model with genic selection a proof is given that there is a critical selection value such that if the selection coefficient is greater than or equal to the critical value then fixation, when the boundary 1 is hit, has probability 1 beginning from any nonzero frequency. This is an analytic proof different from the proofs of Der, Epstein and Plotkin (2011) and Foucart (2013). An application in the infinitely-many-alleles Λ-Fleming-Viot process is finding an interesting identity for the frequency spectrum of alleles that is based on size biasing. The moment dual process in the Fleming-Viot process is the usual Λ-coalescent tree back in time. The Wright-Fisher representation using a different set of polynomials gn(x) as test functions produces a dual death process which has a similarity to the Kingman coalescent and decreases by units of one. The eigenvalues of the process are analogous to the Jacobi polynomials when expressed in terms of gn(x), playing the role of xn. Under the stationary distribution when there is mutation, is analogous to the nth moment in a beta distribution. There is a d-dimensional version gn(X), and even an intriguing Ewens' sampling formula analogy when d → ∞.

The yeast Saccharomyces cerevisiae has emerged as an ideal model system to study the dynamics of prion proteins which are responsible for a number of fatal neurodegenerative diseases in humans. Within an infected cell, prion proteins aggregate in complexes which may increase in size or be fragmented and are transmitted upon cell division. Recent work in yeast suggests that only aggregates below a critical size are transmitted efficiently. We formulate a continuous-time branching process model of a yeast colony under conditions of prion curing. We generalize previous approaches by providing an explicit formula approximating prion loss as influenced by both aggregate growth and size-dependent transmission.

The critical paths of a max-plus linear system with noise are random variables. In this paper we introduce the edge criticalities which measure how often the critical paths traverse each edge in the precedence graph. We also present the parallel path approximation, a novel method for approximating these new statistics as well as the previously studied max-plus exponent. We show that, for low amplitude noise, the critical paths spend most of their time traversing the deterministic maximally weighted cycle and that, as the noise amplitude is increased, the critical paths become more random and their distribution over the edges in the precedence graph approaches a highly uniform measure of maximal entropy.