Source author record

Sepideh Mirrahimi

Sepideh Mirrahimi 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

19works
5topics
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

19 published item(s)

preprint2022arXiv

Adaptation to a heterogeneous patchy environment with nonlocal dispersion

In this work, we provide an asymptotic analysis of the solutions to an elliptic integro-differential equation. This equation describes the evolutionary equilibria of a phenotypically structured population, subject to selection, mutation, and both local and non-local dispersion in a spatially heterogeneous, possibly patchy, environment. Considering small effects of mutations, we provide an asymptotic description of the equilibria of the phenotypic density. This asymptotic description involves a Hamilton-Jacobi equation with constraint coupled with an eigenvalue problem. Based on such analysis, we characterize some qualitative properties of the phenotypic density at equilibrium depending on the heterogeneity of the environment. In particular, we show that when the heterogeneity of the environment is low, the population concentrates around a single phenotypic trait leading to a unimodal phenotypic distribution. On the contrary, a strong fragmentation of the environment leads to multi-modal phenotypic distributions.

preprint2022arXiv

Filling the gap between individual-based evolutionary models and Hamilton-Jacobi equations

We consider a stochastic model for the evolution of a discrete population structured by a trait with values on a finite grid of the torus, and with mutation and selection. Traits are vertically inherited unless a mutation occurs, and influence the birth and death rates. We focus on a parameter scaling where population is large, individual mutations are small but not rare, and the grid mesh for the trait values is much smaller than the size of mutation steps. When considering the evolution of the population in a long time scale, the contribution of small sub-populations may strongly influence the dynamics. Our main result quantifies the asymptotic dynamics of sub-population sizes on a logarithmic scale. We establish that under the parameter scaling the logarithm of the stochastic population size process, conveniently normalized, converges to the unique viscosity solution of a Hamilton-Jacobi equation. Such Hamilton-Jacobi equations have already been derived from parabolic integro-differential equations and have been widely developed in the study of adaptation of quantitative traits. Our work provides a justification of this framework directly from a stochastic individual based model, leading to a better understanding of the results obtained within this approach. The proof makes use of almost sure maximum principles and careful controls of the martingale parts.

preprint2020arXiv

A Fokker-Planck approach to the study of robustness in gene expression

We study several Fokker-Planck equations arising from a stochastic chemical kinetic system modeling a gene regulatory network in biology. The densities solving the Fokker-Planck equations describe the joint distribution of the messenger RNA and micro RNA content in a cell. We provide theoretical and numerical evidences that the robustness of the gene expression is increased in the presence of micro RNA. At the mathematical level, increased robustness shows in a smaller coefficient of variation of the marginal density of the messenger RNA in the presence of micro RNA. These results follow from explicit formulas for solutions. Moreover, thanks to dimensional analyses and numerical simulations we provide qualitative insight into the role of each parameter in the model. As the increase of gene expression level comes from the underlying stochasticity in the models, we eventually discuss the choice of noise in our models and its influence on our results.

preprint2020arXiv

Dirac mass dynamics in multidimensional nonlocal parabolic equations

Nonlocal Lotka-Volterra models have the property that solutions concentrate as Dirac masses in the limit of small diffusion. Is it possible to describe the dynamics of the limiting concentration points and of the weights of the Dirac masses? What is the long time asymptotics of these Dirac masses? Can several Dirac masses co-exist? We will explain how these questions relate to the so-called "constrained Hamilton-Jacobi equation" and how a form of canonical equation can be established. This equation has been established assuming smoothness. Here we build a framework where smooth solutions exist and thus the full theory can be developed rigorously. We also show that our form of canonical equation comes with a structure of gradient flow. Numerical simulations show that the trajectories can exhibit unexpected dynamics well explained by this equation. Our motivation comes from population adaptive evolution a branch of mathematical ecology which models darwinian evolution.

preprint2020arXiv

Survival criterion for a population subject to selection and mutations ; Application to temporally piecewise constant environments

We study a parabolic Lotka-Volterra type equation that describes the evolution of a population structured by a phenotypic trait, under the effects of mutations and competition for resources modelled by a nonlocal feedback. The limit of small mutations is characterized by a Hamilton-Jacobi equation with constraint that describes the concentration of the population on some traits. This result was already established in Barles-Perthame 2008, Barles-Mirrahimi-Perthame 2009, Lorz-Mirrahimi-Perthame 2011 in a time-homogenous environment, when the asymptotic persistence of the population was ensured by assumptions on either the growth rate or the initial data. Here, we relax these assumptions to extend the study to situations where the population may go extinct at the limit. For that purpose, we provide conditions on the initial data for the asymptotic fate of the population. Finally, we show how this study for a time-homogenous environment allows to consider temporally piecewise constant environments.

preprint2016arXiv

A Hamilton-Jacobi approach to characterize the evolutionary equilibria in heterogeneous environments

In this work, we characterize the solution of a system of elliptic integro-differential equations describing a phenotypically structured population subject to mutation, selection and migration between two habitats. Assuming that the effects of the mutations are small but nonzero, we show that the population's distribution has at most two peaks and we give explicit conditions under which the population will be monomorphic (unimodal distribution) or dimorphic (bimodal distribution). More importantly, we provide a general method to determine the dominant terms of the population's distribution in each case. Our work, which is based on Hamilton-Jacobi equations with constraint, goes further than previous works where such tools were used, for different problems from evolutionary biology, to identify the asymptotic solutions, while the mutations vanish, as a sum of Dirac masses. In order to extend such results to the case with non-vanishing effects of mutations, the main elements are a uniqueness property and the computation of the correctors. This method allows indeed to go further than the Gaussian approximation commonly used by biologists and makes a connection between the theories of adaptive dynamics and quantitative genetics. Our work being motivated by biological questions, the objective of this article is to provide the mathematical details which are necessary for our biological results [16].

preprint2015arXiv

A class of Hamilton-Jacobi equations with constraint: uniqueness and constructive approach

We discuss a class of time-dependent Hamilton-Jacobi equations, where an unknown function of time is intended to keep the maximum of the solution to the constant value 0. Our main result is that the full problem has a unique viscosity solution, which is in fact classical. The motivation is a selection-mutation model which, in the limit of small diffusion, exhibits concentration on the zero level set of the solution of the Hamilton-Jacobi equation. Uniqueness is obtained by noticing that, as a consequence of the dynamic programming principle, the solution of the Hamilton-Jacobi equation is classical. It is then possible to write an ODE for the maximum of the solution, and treat the full problem as a nonstandard Cauchy problem.

preprint2015arXiv

Uniqueness in a class of Hamilton-Jacobi equations with constraints

In this note, we discuss a class of time-dependent Hamilton-Jacobi equations depending on a function of time, this function being chosen in order to keep the maximum of the solution to the constant value 0. The main result of the note is that the full problem has a unique classical solution. The motivation is a selection-mutation model which, in the limit of small diffusion, exhibits concentration on the zero level set of the solution of the Hamilton-Jacobi equation. The uniqueness result that we prove implies strong convergence and error estimates for the selection-mutation model.

preprint2014arXiv

Asymptotic analysis of a selection model with space

Selection of a phenotypical trait can be described in mathematical terms by 'stage structured' equations which are usually written under the form of integral equations so as to express competition for resource between individuals whatever is their trait. The solutions exhibit a concentration effect (selection of the fittest); when a small parameter is introduced they converge to a Dirac mass. An additional space variable can be considered in order to take into account local environmental conditions. Here we assume this environment is a single nutrient which diffuses in the domain. In this framework, we prove that the solution converges to a Dirac mass in the physiological trait which depends on time and on the location in space with Lipschitz continuity. The main difficulties come from the lack of compactness in time and trait variables. Strong convergence can be recovered from uniqueness in the limiting constrained equation after Hopf-Cole change of unknown. Our analysis is motivated by a model of tumor growth introduced by Lorz et al. (preprint) in order to explain emergence of resistance to therapy.

preprint2014arXiv

Direct competition results from strong competiton for limited resource

We study a model of competition for resource through a chemostat-type model where species consume the common resource that is constantly supplied. We assume that the species and resources are characterized by a continuous trait. As already proved, this model, although more complicated than the usual Lotka-Volterra direct competition model, describes competitive interactions leading to concentrated distributions of species in continuous trait space. Here we assume a very fast dynamics for the supply of the resource and a fast dynamics for death and uptake rates. In this regime we show that factors that are independent of the resource competition become as important as the competition efficiency and that the direct competition model is a good approximation of the chemostat. Assuming these two timescales allows us to establish a mathematically rigorous proof showing that our resource-competition model with continuous traits converges to a direct competition model. We also show that the two timescales assumption is required to mathematically justify the corresponding classic result on a model consisting of only finite number of species and resources (MacArthur, R. Theor. Popul. Biol. 1970:1, 1-11). This is performed through asymptotic analysis, introducing different scales for the resource renewal rate and the uptake rate. The mathematical difficulty relies in a possible initial layer for the resource dynamics. The chemostat model comes with a global convex Lyapunov functional. We show that the particular form of the competition kernel derived from the uptake kernel, satisfies a positivity property which is known to be necessary for the direct competition model to enjoy the related Lyapunov functional.

preprint2014arXiv

Influence of a spatial structure on the long time behavior of a competitive Lotka-Volterra type system

To describe population dynamics, it is crucial to take into account jointly evolution mechanisms and spatial motion. However, the models which include these both aspects, are not still well-understood. Can we extend the existing results on type structured populations, to models of populations structured by type and space, considering diffusion and nonlocal competition between individuals? We study a nonlocal competitive Lotka-Volterra type system, describing a spatially structured population which can be either monomorphic or dimorphic. Considering spatial diffusion, intrinsic death and birth rates, together with death rates due to intraspecific and interspecific competition between the individuals, leading to some integral terms, we analyze the long time behavior of the solutions. We first prove existence of steady states and next determine the long time limits, depending on the competition rates and the principal eigenvalues of some operators, corresponding somehow to the strength of traits. Numerical computations illustrate that the introduction of a new mutant population can lead to the long time evolution of the spatial niche.

preprint2014arXiv

Singular limits for reaction-diffusion equations with fractional Laplacian and local or nonlocal nonlinearity

We perform an asymptotic analysis of models of population dynamics with a fractional Laplacian and local or nonlocal reaction terms. The first part of the paper is devoted to the long time/long range rescaling of the fractional Fisher-KPP equation. This rescaling is based on the exponential speed of propagation of the population. In particular we show that the only role of the fractional Laplacian in determining this speed is at the initial layer where it determines the thickness of the tails of the solutions. Next, we show that such rescaling is also possible for models with non-local reaction terms, as selection-mutation models. However, to obtain a more relevant qualitative behavior for this second case, we introduce, in the second part of the paper, a second rescaling where we assume that the diffusion steps are small. In this way, using a WKB ansatz, we obtain a Hamilton-Jacobi equation in the limit which describes the asymptotic dynamics of the solutions, similarly to the case of selection-mutation models with a classical Laplace term or an integral kernel with thin tails. However, the rescaling introduced here is very different from the latter cases. We extend these results to the multidimensional case.

preprint2013arXiv

A Hamilton-Jacobi approach for a model of population structured by space and trait

We study a non-local parabolic Lotka-Volterra type equation describing a population structured by a space variable x 2 Rd and a phenotypical trait 2 . Considering diffusion, mutations and space-local competition between the individuals, we analyze the asymptotic (long- time/long-range in the x variable) exponential behavior of the solutions. Using some kind of real phase WKB ansatz, we prove that the propagation of the population in space can be described by a Hamilton-Jacobi equation with obstacle which is independent of . The effective Hamiltonian is derived from an eigenvalue problem. The main difficulties are the lack of regularity estimates in the space variable, and the lack of comparison principle due to the non-local term.

preprint2013arXiv

Time fluctuations in a population model of adaptive dynamics

We study the dynamics of phenotypically structured populations in environments with fluctuations. In particular, using novel arguments from the theories of Hamilton-Jacobi equations with constraints and homogenization, we obtain results about the evolution of populations in environments with time oscillations, the development of concentrations in the form of Dirac masses, the location of the dominant traits and their evolution in time. Such questions have already been studied in time homogeneous environments. More precisely we consider the dynamics of a phenotypically structured population in a changing environment under mutations and competition for a single resource. The mathematical model is a non-local parabolic equation with a periodic in time reaction term. We study the asymptotic behavior of the solutions in the limit of small diffusion and fast reaction. Under concavity assumptions on the reaction term, we prove that the solution converges to a Dirac mass whose evolution in time is driven by a Hamilton-Jacobi equation with constraint and an effective growth/death rate which is derived as a homogenization limit. We also prove that, after long-time, the population concentrates on a trait where the maximum of an effective growth rate is attained. Finally we provide an example showing that the time oscillations may lead to a strict increase of the asymptotic population size.

preprint2012arXiv

A homogenization approach for the motion of motor proteins

We consider the asymptotic behavior of an evolving weakly coupled Fokker-Planck system of two equations set in a periodic environment. The magnitudes of the diffusion and the coupling are respectively proportional and inversely proportional to the size of the period. We prove that, as the period tends to zero, the solutions of the system either propagate (concentrate) with a fixed constant velocity (determined by the data) or do not move at all. The system arises in the modeling of motor proteins which can take two different states. Our result implies that, in the limit, the molecules either move along a filament with a fixed direction and constant speed or remain immobile.

preprint2012arXiv

Adaptation and migration of a population between patches

A Hamilton-Jacobi formulation has been established previously for phenotypically structured population models where the solution concentrates as Dirac masses in the limit of small diffusion. Is it possible to extend this approach to spatial models? Are the limiting solutions still in the form of sums of Dirac masses? Does the presence of several habitats lead to polymorphic situations? We study the stationary solutions of a structured population model, while the population is structured by continuous phenotypical traits and discrete positions in space. The growth term varies from one habitable zone to another, for instance because of a change in the temperature. The individuals can migrate from one zone to another with a constant rate. The mathematical modeling of this problem, considering mutations between phenotypical traits and competitive interaction of individuals within each zone via a single resource, leads to a system of coupled parabolic integro-differential equations. We study the asymptotic behavior of the stationary solutions to this model in the limit of small mutations. The limit, which is a sum of Dirac masses, can be described with the help of an effective Hamiltonian. The presence of migration can modify the dominant traits and lead to polymorphic situations.

preprint2012arXiv

Invasion fronts with variable motility: phenotype selection, spatial sorting and wave acceleration

Invasion fronts in ecology are well studied but very few mathematical results concern the case with variable motility (possibly due to mutations). Based on an apparently simple reaction-diffusion equation, we explain the observed phenomena of front acceleration (when the motility is unbounded) as well as other quantitative results, such as the selection of the most motile individuals (when the motility is bounded). The key argument for the construction and analysis of traveling fronts is the derivation of the dispersion relation linking the speed of the wave and the spatial decay. When the motility is unbounded we show that the position of the front scales as $t^{3/2}$. When the mutation rate is low we show that the canonical equation for the dynamics of the fittest trait should be stated as a PDE in our context. It turns out to be a type of Burgers equation with source term.

preprint2011arXiv

Population structured by a space variable and a phenotypical trait

We consider populations structured by a phenotypic trait and a space variable, in a non-homogeneous environment. In the case of sex- ual populations, we are able to derive models close to existing mod- els in theoretical biology, from a structured population model. We then analyze the dynamics of the population using a simplified model, where the population either propagates through the whole space or it survives but remains confined in a limited range. For asexual pop- ulations, we show that the dynamics are simpler. In this case, the population cannot remain confined in a limited range, i.e. the popu- lation, if it does not get extinct, propagates through the whole space.

preprint2011arXiv

Singular Hamilton-Jacobi equation for the tail problem

In this paper we study the long time-long range behavior of reaction diffusion equations with negative square root -type reaction terms. In particular we investigate the exponential behavior of the solutions after a standard hyperbolic scaling. This leads to a Hamilton-Jacobi variational inequality with an obstacle that depends on the solution itself and defines the open set where the limiting solution does not vanish. Counter-examples show a nontrivial lack of uniqueness for the variational inequality depending on the conditions imposed on the boundary of this open set. Both Dirichlet and state constraints boundary conditions play a role. When the competition term does not change sign, we can identify the limit, while, in general, we find lower and upper bounds for the limit. Although models of this type are rather old and extinction phenomena are as important as blow-up, our motivation comes from the so-called "tail problem" in population biology. One way to avoid meaningless exponential tails, is to impose extra-mortality below a given survival threshold. Our study shows that the precise form of this extra-mortality term is asymptotically irrelevant and that, in the survival zone, the population profile is impacted by the survival threshold (except in the very particular case when the competition term is non-positive).