Source author record

Thierry Huillet

Thierry Huillet 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

20works
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

20 published item(s)

preprint2021arXiv

A generating function approach to Markov chains undergoing binomial catastrophes

In a Markov chain population model subject to catastrophes, random immigration events (birth), promoting growth, are in balance with the effect of binomial catastrophes that cause recurrent mass removal (death). Using a generating function approach, we study two versions of such population models when the binomial catastrophic events are of a slightly different random nature. In both cases, we describe the subtle balance between the two birth and death conflicting effects.

preprint2021arXiv

Scaling features of two special Markov chains involving total disasters

Catastrophe Markov chain population models have received a lot of attention in the recent past. We herewith consider two special cases of such models involving total disasters, both in discrete and in continuous-time. Depending on the parameters range, the two models can show up a recurrence/transience transition and, in the critical case, a positive/null recurrence transition. The collapse transition probabilities are chosen in such a way that the models are exactly solvable and, in case of positive recurrence, intimately related to the extended Sibuya and Pareto-Zipf distributions whose divisibility and self-decomposability properties are shown relevant. The study includes: existence and shape of the invariant measure, time-reversal, return time to the origin, contact probability at the origin, extinction probability, height and length of the excursions, a renewal approach to the fraction of time spent in the catastrophic state, scale function, first time to collapse and first-passage times, divisibility properties.

preprint2016arXiv

Additional aspects of the generalized linear-fractional branching process

We derive some additional results on the Bienyamé-Galton-Watson branching process with $θ-$linear fractional branching mechanism, as studied in \cite{Sag}. This includes: the explicit expression of the limit laws in both the sub-critical cases and the super-critical cases with finite mean, the long-run behavior of the population size in the critical case, limit laws in the super-critical cases with infinite mean when the $θ$-process is either regular or explosive, results regarding the time to absorption, an expression of the probability law of the $θ$-branching mechanism involving Bell polynomials, the explicit computation of the stochastic transition matrix of the $θ-$process, together with its powers.

preprint2016arXiv

Deterministic versus stochastic aspects of superexponential population growth models

Deterministic population growth models with power-law rates can exhibit a large variety of growth behaviors, ranging from algebraic, exponential to hyperexponential (finite time explosion). In this setup, selfsimilarity considerations play a key role, together with two time substitutions. Two stochastic versions of such models are investigated, showing a much richer variety of behaviors. One is the Lamperti construction of selfsimilar positive stochastic processes based on the exponentiation of spectrally positive processes , followed by an appropriate time change. The other one is based on stable continuous-state branching processes, given by another Lamperti time substitution applied to stable spectrally positive processes.

preprint2016arXiv

On discrete evolutionary dynamics driven by quadratic interactions

After an introduction to the general topic of models for a given locus of a diploid population whose quadratic dynamics is determined by a fitness landscape, we consider more specifically the models that can be treated using genetic (or train) algebras. In this setup, any quadratic offspring interaction can produce any type of offspring and after the use of specific changes of basis, we study the evolution and possible stability of some examples. We also consider some examples that cannot be treated using the framework of genetic algebras. Among these are bistochastic matrices.

preprint2016arXiv

Some Combinatorial Aspects of Discrete Non-linear Population Dynamics

Motivated by issues arising in population dynamics, we consider the problem of iterating a given analytic function a number of times. We use the celebrated technique known as Carleman linearization that turns (for a certain class of functions) this problem into simply taking the power of a real number. We expand this method, showing in particular that it can be used for population models with immigration, and we also apply it to the famous logistic map. We also are able to give a number of results for the invariant density of this map, some being related to the Carleman linearization.

preprint2014arXiv

Did the ever dead outnumber the living and when? A birth-and-death approach

This paper is an attempt to formalize analytically the question raised in "World Population Explained: Do Dead People Outnumber Living, Or Vice Versa?" Huffington Post, \cite{HJ}. We start developing simple deterministic Malthusian growth models of the problem (with birth and death rates either constant or time-dependent) before running into both linear birth and death Markov chain models and age-structured models.

preprint2014arXiv

On Moebius duality and Coarse-Graining

We study duality relations for zeta and Möbius matrices and monotone conditions on the kernels. We focus on the cases of family of sets and partitions. The conditions for positivity of the dual kernels are stated in terms of the positive Möbius cone of functions, which is described in terms of Sylvester formulae. We study duality under coarse-graining and show that an $h-$transform is needed to preserve stochasticity. We give conditions in order that zeta and Möbius matrices admit coarse-graining, and we prove they are satisfied for sets and partitions. This is a source of relevant examples in genetics on the haploid and multi-allelic Cannings models.

preprint2014arXiv

On the spread of a branching Brownian motion whose offspring number has infinite variance

We study the impact on shape parameters of an underlying Bienaymé-Galton-Watson branching process (height, width and first hitting time), of having a non-spatial branching mechanism with infinite variance. Aiming at providing a comparative study of the spread of an epidemics whose dynamics is given by the modulus of a branching Brownian motion (BBM) we then consider spatial branching processes in dimension d, not necessarily integer. The underlying branching mechanism is then either a binary branching model or one presenting infinite variance. In particular we evaluate the chance p(x) of being hit if the epidemics started away at distance x. We compute the large x tail probabilities of this event, both when the branching mechanism is regular and when it exhibits very large fluctuations.

preprint2013arXiv

Is superhydrophobicity robust with respect to disorder?

We consider theoretically the Cassie-Baxter and Wenzel states describing the wetting contact angles for rough substrates. More precisely, we consider different types of periodic geometries such as square protrusions and disks in 2D, grooves and nanoparticles in 3D and derive explicitly the contact angle formulas. We also show how to introduce the concept of surface disorder within the problem and, inspired by biomimetism, study its effect on superhydrophobicity. Our results, quite generally, prove that introducing disorder, at fixed given roughness, will lower the contact angle: a disordered substrate will have a lower contact angle than a corresponding periodic substrate. We also show that there are some choices of disorder for which the loss of superhydrophobicity can be made small, making superhydrophobicity robust.

preprint2013arXiv

Occupancy distributions arising in sampling from Gibbs-Poisson abundance models

Estimating the number $n$ of unseen species from a $k-$sample displaying only $p\leq k$ distinct sampled species has received attention for long. It requires a model of species abundance together with a sampling model. We start with a discrete model of iid stochastic species abundances, each with Gibbs-Poisson distribution. A $k-$sample drawn from the $n-$species abundances vector is the one obtained while conditioning it on summing to $k$% . We discuss the sampling formulae (species occupancy distributions, frequency of frequencies) in this context. We then develop some aspects of the estimation of $n$ problem from the size $k$ of the sample and the observed value of $P_{n,k}$, the number of distinct sampled species. It is shown that it always makes sense to study these occupancy problems from a Gibbs-Poisson abundance model in the context of a population with infinitely many species. From this extension, a parameter $γ$ naturally appears, which is a measure of richness or diversity of species. We rederive the sampling formulae for a population with infinitely many species, together with the distribution of the number $P_{k}$ of distinct sampled species. We investigate the estimation of $γ$ problem from the sample size $k$ and the observed value of $P_{k}$. We then exhibit a large special class of Gibbs-Poisson distributions having the property that sampling from a discrete abundance model may equivalently be viewed as a sampling problem from a random partition of unity, now in the continuum. When $n$ is finite, this partition may be built upon normalizing $% n$ infinitely divisible iid positive random variables by its partial sum. It is shown that the sampling process in the continuum should generically be biased on the total length appearing in the latter normalization. A construction with size-biased sampling from the ranked normalized jumps of a subordinator is also supplied, would the problem under study present infinitely many species. We illustrate our point of view with many examples, some of which being new ones.

preprint2013arXiv

Pareto genealogies arising from a Poisson branching evolution model with selection

We study a class of coalescents derived from a sampling procedure out of N i.i.d. Pareto(alpha) random variables, normalized by their sum, including beta-size-biasing on total length effects (beta < alpha). Depending on the range of alpha, we derive the large N limit coalescents structure, leading either to a discrete-time Poisson-Dirichlet(alpha,-beta) Xi-coalescent (alpha in [0,1)), or to a family of continuous-time Beta(2-alpha,alpha-beta) Lambda-coalescents (alpha in [1,2)), or to the Kingman coalescent (alpha geq 2). We indicate that this class of coalescent processes (and their scaling limits) may be viewed as the genealogical processes of some forward in time evolving branching population models including selection effects. In such constant-size population models, the reproduction step, which is based on a fitness-dependent Poisson Point Process with scaling power-law(alpha)intensity, is coupled to a selection step consisting of sorting out the N fittest individuals issued from the reproduction step.

preprint2011arXiv

A Bose-Einstein Approach to the Random Partitioning of an Integer

Consider N equally-spaced points on a circle of circumference N. Choose at random n points out of $N$ on this circle and append clockwise an arc of integral length k to each such point. The resulting random set is made of a random number of connected components. Questions such as the evaluation of the probability of random covering and parking configurations, number and length of the gaps are addressed. They are the discrete versions of similar problems raised in the continuum. For each value of k, asymptotic results are presented when n,N both go to infinity according to two different regimes. This model may equivalently be viewed as a random partitioning problem of N items into n recipients. A grand-canonical balls in boxes approach is also supplied, giving some insight into the multiplicities of the box filling amounts or spacings. The latter model is a k-nearest neighbor random graph with N vertices and kn edges. We shall also briefly consider the covering problem in the context of a random graph model with N vertices and n (out-degree 1) edges whose endpoints are no more bound to be neighbors.

preprint2011arXiv

A branching diffusion model of selection: from the neutral Wright-Fisher case to the one including mutations

We consider diffusion processes x_{t} on the unit interval. Doob-transformation techniques consist of a selection of x_{t}-paths procedure. The law of the transformed process is the one of a branching diffusion system of particles, each diffusing like a new process tilde{x}_{t}, superposing an additional drift to the one of x_{t}. Killing and/or branching of tilde{x}_{t}-particles occur at some space-dependent rate lambda. For this transformed process, so in the class of branching diffusions, the question arises as to whether the particle system is sub-critical, critical or super-critical. In the first two cases, extinction occurs with probability one. We apply this circle of ideas to diffusion processes arising in population genetics. In this setup, the process x_{t} is a Wright-Fisher (WF) diffusion, either neutral or with mutations. We study a particular Doob transform which is based on the exponential function in the usual fitness parameter sigma. We have in mind that this is an alternative way to introduce selection or fitness in both WF-like diffusions, leading to branching diffusion models ideas. For this Doob-transform model of fitness, the usual selection drift sigma x(1-x) should be superposed to the one of x_{t} to form tilde{x}_{t} which is the process that can branch, binarily. In the first neutral case, there is a trade-off between branching events giving birth to new particles and absorption at the boundaries, killing the particles. Under our assumptions, the branching diffusion process gets eventually globally extinct in finite time with exponential tails. In the second case with mutations, there is a trade-off between killing events removing some particles from the system and reflection at the boundaries where the particles survive. This branching diffusion process also gets eventually globally extinct but in very long finite time with power-law tails. Our approach relies on the spectral expansion of the transition probability kernels of both x_{t} and tilde{x}_{t}.

preprint2011arXiv

Metastable wetting

Consider a droplet of liquid on top of a grooved substrate. The wetting or not of a groove implies the crossing of a potential barrier as the interface has to distort, to hit the bottom of the groove. We start with computing the free energies of the dry and wet states in the context of a simple thermodynamical model before switching to a random microscopic version pertaining to the Solid-on-Solid (SOS) model. For some range in parameter space (Young angle, pressure difference, aspect ratio), the dry and wet states both share the same free energy, which means coexistence. We compute these coexistence lines together with the metastable regions. In the SOS case, we describe the dynamic transition between coexisting states in wetting. We show that the expected time to switch from one state to the other grows exponentially with the free energy barrier between the stable states and the saddle state, proportional to the groove's width. This random time appears to have an exponential-like distribution.

preprint2009arXiv

Information and (co-)variances in discrete evolutionary genetics involving solely selection

The purpose of this Note is twofold: First, we introduce the general formalism of evolutionary genetics dynamics involving fitnesses, under both the deterministic and stochastic setups, and chiefly in discrete-time. In the process, we particularize it to a one-parameter model where only a selection parameter is unknown. Then and in a parallel manner, we discuss the estimation problems of the selection parameter based on a single-generation frequency distribution shift under both deterministic and stochastic evolutionary dynamics. In the stochastics, we consider both the celebrated Wright-Fisher and Moran models.

preprint2009arXiv

Random walk versus random line

We consider random walks X_n in Z+, obeying a detailed balance condition, with a weak drift towards the origin when X_n tends to infinity. We reconsider the equivalence in law between a random walk bridge and a 1+1 dimensional Solid-On-Solid bridge with a corresponding Hamiltonian. Phase diagrams are discussed in terms of recurrence versus wetting. A drift -delta/X_n of the random walk yields a Solid-On-Solid potential with an attractive well at the origin and a repulsive tail delta(delta+2)/(8X_n^2) at infinity, showing complete wetting for delta<=1 and critical partial wetting for delta>1.

preprint2006arXiv

Correlations of a bound interface over a random substrate

The correlation function of a one-dimensional interface over a random substrate, bound to the substrate by a pressure term, is studied by Monte-Carlo simulation. It is found that the height correlation < h_i ; h_{i+j} >, averaged over the substrate disorder, fits a form exp(-(j/b)^c) to a surprising precision in the full range of j where the correlation is non-negligible. The exponent c increases from 1.0 to 1.5 when the interface tension is taken larger and larger.

preprint2005arXiv

A necklace of Wulff shapes

In a probabilistic model of a film over a disordered substrate, Monte-Carlo simulations show that the film hangs from peaks of the substrate. The film profile is well approximated by a necklace of Wulff shapes. Such a necklace can be obtained as the infimum of a collection of Wulff shapes resting on the substrate. When the random substrate is given by iid heights with exponential distribution, we prove estimates on the probability density of the resulting peaks, at small density.