Source author record

Fabio A. C. C. Chalub

Fabio A. C. C. Chalub 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

10works
2topics
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

10 published item(s)

preprint2022arXiv

Solution of the Neutral Kimura equation with two integral constraints

The Kimura equation is a degenerated partial differential equation of drift-diffusion type used in population genetics. Its solution is required to satisfy not only the equation but a series of conservation laws formulated as integral constraints. In this work, we consider a population of two types evolving without mutation or selection, the so-called neutral evolution. We obtain explicit solutions in terms of Gegenbauer polynomials. To satisfy the integral constraints it is necessary to prove new relations satisfied by the Gegenbauer polynomials. The long-term in time asymptotic is also studied.

preprint2016arXiv

Asymptotic expression for the fixation probability of a mutant in star graphs

We consider the Moran process in a graph called the "star" and obtain the asymptotic expression for the fixation probability of a single mutant when the size of the graph is large. The expression obtained corrects the previously known expression announced in reference [E Lieberman, C Hauert, and MA Nowak. Evolutionary dynamics on graphs. Nature, 433(7023):312-316, 2005] and further studied in [M. Broom and J. Rychtar. An analysis of the fixation probability of a mutant on special classes of non-directed graphs. Proc. R. Soc. A-Math. Phys. Eng. Sci., 464(2098):2609-2627, 2008]. We also show that the star graph is an accelerator of evolution, if the graph is large enough.

preprint2016arXiv

Conservative parabolic problems: non-degenerated theory and degenerated examples from population dynamics

We consider partial differential equations (PDE) of drift-diffusion type in the unit interval, supplemented by either two conservation laws or by a conservation law and a further boundary condition. We treat two different cases: (i) uniform parabolic problems; (ii) degenerated problems at the boundaries. The former can be treated in a very general and complete way, much as the traditional boundary value problems. The latter, however, bring new issues, and we restrict our study to a class of forward Kolmogorov equations that arise naturally when the corresponding stochastic process has either one or two absorbing boundaries. These equations are treated by means of a uniform parabolic regularisation, which then yields a measure solution in the vanishing regularisation limit that is unique. Two prototypical problems from population dynamics are treated in detail.

preprint2016arXiv

Optimal Vaccination Strategies and Rational Behaviour in Seasonal Epidemics

We consider a SIR model with temporary immunity and time dependent transmission rate. We assume time dependent vaccination which confers the same immunity as natural infection. We study two types of vaccination strategies: i) optimal vaccination, in the sense that it minimizes the effort of vaccination in the set of vaccination strategies for which, for any sufficiently small perturbation of the disease free state, the number of infectious individuals is monotonically decreasing; ii) Nash-equilibria strategies where all individuals simultaneously minimize the joint risk of vaccination versus the risk of the disease. The former case corresponds to an optimal solution for mandatory vaccinations, while the second correspond to the equilibrium to be expected if vaccination is fully voluntary. We are able to show the existence of both an optimal and Nash strategies in a general setting. In general, these strategies will not be functions but Radon measures. For specific forms of the transmission rate, we provide explicit formulas for the optimal and the Nash vaccination strategies.

preprint2015arXiv

Fixation in large populations: a continuous view of a discrete problem

We study fixation in large, but finite, populations with two types, and dynamics governed by birth-death processes. By considering a restricted class of such processes, we derive a continuous approximation for the probability of fixation that is valid beyond the weak-selection (WS) limit. From the continuous approximations, we then obtain asymptotic approximations for evolutionary dynamics with at most one equilibrium, in the selection-driven regime, that does not preclude a weak-selection regime. As an application, we study the fixation pattern when the infinite population limit has an interior Evolutionary Stable Strategy (ESS): (i) we show that the fixation pattern for the Hawk and Dove game satisfies what we term the one-half law: if the Evolutionary Stable Strategy (ESS) is outside a small interval around $\sfrac{1}{2}$, the fixation is of dominance type; (ii) we also show that, outside of the weak-selection regime, the long-term dynamics of large populations can have very little resemblance to the infinite population case; in addition, we also present results for the case of two equilibria. Finally, we present continuous restatements valid for large populations of two classical concepts naturally defined in the discrete case: (i) the definition of an $\textsc{ESS}_N$ strategy; (ii) the definition of a risk-dominant strategy. We then present two applications of these restatements: (i) we obtain an asymptotic definition valid in the quasi-neutral regime that recovers both the one-third law under linear fitness and the generalised one-third law for $d$-player games; (ii) we extend the ideas behind the (generalised) one-third law outside the quasi-neutral regime and, as a generalisation, we introduce the concept of critical-frequency; (iii) we recover the classification of risk-dominant strategies for $d$-player games.

preprint2013arXiv

The frequency-dependent Wright-Fisher model: diffusive and non-diffusive approximations

We study a class of processes that are akin to the Wright-Fisher model, with transition probabilities weighted in terms of the frequency-dependent fitness of the population types. By considering an approximate weak formulation of the discrete problem, we are able to derive a corresponding continuous weak formulation for the probability density. Therefore, we obtain a family of partial differential equations (PDE) for the evolution of the probability density, and which will be an approximation of the discrete process in the joint large population, small time-steps and weak selection limit. If the fitness functions are sufficiently regular, we can recast the weak formulation in a more standard formulation, without any boundary conditions, but supplemented by a number of conservation laws. The equations in this family can be purely diffusive, purely hyperbolic or of convection-diffusion type, with frequency dependent convection. The particular outcome will depend on the assumed scalings. The diffusive equations are of the degenerate type; using a duality approach, we also obtain a frequency dependent version of the Kimura equation without any further assumptions. We also show that the convective approximation is related to the replicator dynamics and provide some estimate of how accurate is the convective approximation, with respect to the convective-diffusion approximation. In particular, we show that the mode, but not the expected value, of the probability distribution is modelled by the replicator dynamics. Some numerical simulations that illustrate the results are also presented.

preprint2009arXiv

A non-standard evolution problem arising in population genetics

We study the evolution of the probability density of an asexual, one locus population under natural selection and random evolution. This evolution is governed by a Fokker-Planck equation with degenerate coefficients on the boundaries, supplemented by a pair of conservation laws. It is readily shown that no classical or standard weak solution definition yields solvability of the problem. We provide an appropriate definition of weak solution for the problem, for which we show existence and uniqueness. The solution displays a very distinctive structure and, for large time, we show convergence to a unique stationary solution that turns out to be a singular measure supported at the endpoints. An exponential rate of convergence to this steady state is also proved.

preprint2009arXiv

The SIR epidemic model from a PDE point of view

We present a derivation of the classical SIR model through a mean-field approximation from a discrete version of SIR. We then obtain a hyperbolic forward Kolmogorov equation, and show that its projected characteristics recover the standard SIR model. Moreover, we show that the long time limit of the evolution will be a Dirac measure. The exact position will depend on the well-know $R_0$ parameter, and it will be supported on the corresponding stable SIR equilibrium.

preprint2008arXiv

From discrete to continuous evolution models: a unifying approach to drift-diffusion and replicator dynamics

We study the large population limit of the Moran process, assuming weak-selection, and for different scalings. Depending on the particular choice of scalings, we obtain a continuous model that may highlight the genetic-drift (neutral evolution) or natural selection; for one precise scaling, both effects are present. For the scalings that take the genetic-drift into account, the continuous model is given by a singular diffusion equation, together with two conservation laws that are already present at the discrete level. For scalings that take into account only natural selection, we obtain a hyperbolic singular equation that embeds the Replicator Dynamics and satisfies only one conservation law. The derivation is made in two steps: a formal one, where the candidate limit model is obtained, and a rigorous one, where convergence of the probability density is proved. Additional results on the fixation probabilities are also presented.

preprint2006arXiv

Discrete versus continuous models in evolutionary dynamics: from simple to simpler -- and even simpler -- models

There are many different models--both continuous and discrete--used to describe gene mutation fixation. In particular, the Moran process, the Kimura equation and the replicator dynamics are all well known models, that might lead to different conclusions. We present a discussion of a unified framework to embrace all these models, in the large population regime.