Source author record

Benoît Perthame

Benoît Perthame 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

26works
9topics
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

26 published item(s)

preprint2022arXiv

Bounds and long term convergence for the voltage-conductance kinetic system arising in neuroscience

The voltage-conductance equation determines the probability distribution of a stochastic process describing a fluctuation-driven neuronal network arising in the visual cortex. Its structure and degeneracy share many common features with the kinetic Fokker-Planck equation, which has attracted much attention recently. We prove an L $\infty$ bound on the steady state solution and the long term convergence of the evolution problem towards this stationary state. Despite the hypoellipticity property, the difficulty is to treat the boundary conditions that change type along the boundary. This leads us to use specific weights in the Alikakos iterations and adapt the relative entropy method.

preprint2022arXiv

Collective motion driven by nutrient consumption

A classical problem describing the collective motion of cells, is the movement driven by consumption/depletion of a nutrient. Here we analyze one of the simplest such model written as a coupled Partial Differential Equation/Ordinary Differential Equation system which we scale so as to get a limit describing the usually observed pattern. In this limit the cell density is concentrated as a moving Dirac mass and the nutrient undergoes a discontinuity. We first carry out the analysis without diffusion, getting a complete description of the unique limit. When diffusion is included, we prove several specific a priori estimates and interpret the system as a heterogeneous monostable equation. This allow us to obtain a limiting problem which shows the concentration effect of the limiting dynamics.

preprint2022arXiv

Existence of a global weak solution for a reaction-diffusion problem with membrane conditions

Several problems, issued from physics, biology or the medical science, lead to parabolic equations set in two sub-domains separated by a membrane with selective permeability to specific molecules. The corresponding boundary conditions, describing the flow through the membrane, are compatible with mass conservation and energy dissipation, and are called the Kedem-Katchalsky conditions. Additionally, in these models, written as reaction-diffusion systems, the reaction terms have a quadratic behaviour. M. Pierre and his collaborators have developed a complete $L^1$ theory for reaction-diffusion systems with different diffusions. Here, we adapt this theory to the membrane boundary conditions and prove the existence of weak solutions when the initial data has only $L^1$ regularity using the truncation method for the nonlinearities. In particular, we establish several estimates as the $W^{1,1}$ regularity of the solutions. Also, a crucial step is to adapt the fundamental $L^2$ (space, time) integrability lemma to our situation.

preprint2022arXiv

Relaxation of the Cahn-Hilliard equation with singular single-well potential and degenerate mobility

The degenerate Cahn-Hilliard equation is a standard model to describe living tissues. It takes into account cell populations undergoing short-range attraction and long-range repulsion effects. In this framework, we consider the usual Cahn-Hilliard equation with a singular single-well potential and degenerate mobility. These degeneracy and singularity induce numerous difficulties, in particular for its numerical simulation. To overcome these issues, we propose a relaxation system formed of two second order equations which can be solved with standard packages. This system is endowed with an energy and an entropy structure compatible with the limiting equation. Here, we study the theoretical properties of this system; global existence and convergence of the relaxed system to the degenerate Cahn-Hilliard equation. We also study the long-time asymptotics which interest relies on the numerous possible steady states with given mass.

preprint2022arXiv

Selection-mutation dynamics with asymmetrical reproduction kernels

We study a family of selection-mutation models of a sexual population structured by a phenotypical trait. The main feature of these models is the asymmetric trait heredity or fecundity between the parents : we assume that each individual inherits mostly its trait from the female or that the trait acts on the female fecundity but does not affect male. Following previous works inspired from principles of adaptive dynamics, we rescale time and assume that mutations have limited effects on the phenotype. Our goal is to study the asymptotic behavior of the population distribution. We derive non-extinction conditions and BV estimates on the total population. We also obtain Lipschitz estimates on the solutions of Hamilton-Jacobi equations that arise from the study of the population distribution concentration at fittest traits. Concentration results are obtained in some special cases by using a Lyapunov functional.

preprint2021arXiv

A non-expanding transport distance for some structured equations

Structured equations are a standard modeling tool in mathematical biology. They areintegro-differential equations where the unknown depends on one or several variables, representing the state or phenotype of individuals. A large literature has been devoted to many aspects of these equations and in particular to the study of measure solutions.Here we introduce a transport distance closely related to the Monge-Kantorovich distance,which appears to be non-expanding for several (mainly linear) examples of structured equations.

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

Dynamics of concentration in a population model structured by age and a phenotypical trait

We study a mathematical model describing the growth process of a population structured by age and a phenotypical trait, subject to aging, competition between individuals and rare mutations. Our goals are to describe the asymptotic behaviour of the solution to a renewal type equation, and then to derive properties that illustrate the adaptive dynamics of such a population. We begin with a simplified model by discarding the effect of mutations, which allows us to introduce the main ideas and state the full result. Then we discuss the general model and its limitations. Our approach uses the eigenelements of a formal limiting operator, that depend on the structuring variables of the model and define an effective fitness. Then we introduce a new method which reduces the convergence proof to entropy estimates rather than estimates on the constrained Hamilton-Jacobi equation. Numerical tests illustrate the theory and show the selection of a fittest trait according to the effective fitness. For the problem with mutations, an unusual Hamiltonian arises with an exponential growth, for which we show existence of a global viscosity solution, using an uncommon a priori estimate and a new uniqueness result.

preprint2020arXiv

Dynamics of concentration in a population structured by age and a phenotypic trait with mutations. Convergence of the corrector

We study an equation structured by age and a phenotypic trait describing the growth process of a population subject to aging, competition between individuals, and mutations. This leads to a renewal equation which occurs in many evolutionary biology problems. We aim to describe precisely the asymp-totic behavior of the solution, to infer properties that illustrate the concentration and adaptive dynamics of such a population. This work is a continuation of [38] where the case without mutations is considered. When mutations are taken into account, it is necessary to control the corrector which is the main novelty of the present paper. Our approach consists in defining, by the Hopf transform, a Hamilton-Jacobi equation with an effective Hamiltonian as in homogenization problems. Its solution carries the singular part of the limiting density (typically Dirac masses) and the corrector defines the weights. The main new result of this paper is to prove that the corrector is uniformly bounded, using only the global Lipschitz and semi-convexity estimates for the viscosity solution of the Hamilton-Jacobi equation. We also establish the limiting equation satisfied by the corrector. To the best of our knowledge, this is the first example where such bounds can be proved in such a context.

preprint2020arXiv

Free boundary limit of tumor growth model with nutrient

Both compressible and incompressible porous medium models are used in the literature to describe the mechanical properties of living tissues. These two classes of models can be related using a stiff pressure law. In the incompressible limit, the compressible model generates a free boundary problem of Hele-Shaw type where incompressibility holds in the saturated phase. Here we consider the case with a nutrient. Then, a badly coupled system of equations describes the cell density number and the nutrient concentration. For that reason, the derivation of the free boundary (incompressible) limit was an open problem, in particular a difficulty is to establish the so-called complementarity relation which allows to recover the pressure using an elliptic equation. To establish the limit, we use two new ideas. The first idea, also used recently for related problems, is to extend the usual Aronson-Bénilan estimates in $L^\infty$ to an $L^2$ setting. The second idea is to derive a sharp uniform $L^4$ estimate on the pressure gradient, independently of space dimension.

preprint2020arXiv

The Aronson-Bénilan Estimate in Lebesgue Spaces

In a celebrated three-pages long paper in 1979, Aronson and Bénilan obtained a remarkable estimate on second order derivatives for the solution of the porous media equation. Since its publication, the theory of porous medium flow has expanded relentlessly with applications including thermodynamics, gas flow, ground water flow as well as ecological population dynamics. The purpose of this paper is to clarify the use of recent extensions of the Aronson and Bénilan estimate in Lp spaces, of some modifications and improvements, as well as to show certain limitations of their strategy.

preprint2019arXiv

Hele-Shaw limit for a system of two reaction-(cross-)diffusion equations for living tissues

Multiphase mechanical models are now commonly used to describe living tissues including tumour growth. The specific model we study here consists of two equations of mixed parabolic and hyperbolic type which extend the standard compressible porous medium equation, including cross-reaction terms. We study the incompressible limit, when the pressure becomes stiff, which generates a free boundary problem. We establish the complementarity relation and also a segregation result. Several major mathematical difficulties arise in the two species case. Firstly, the system structure makes comparison principles fail. Secondly, segregation and internal layers limit the regularity available on some quantities to BV. Thirdly, the Aronson-B{é}nilan estimates cannot be established in our context. We are lead, as it is classical, to add correction terms. This procedure requires technical manipulations based on BV estimates only valid in one space dimension. Another novelty is to establish an L1 version in place of the standard upper bound.

preprint2016arXiv

Toward an integrated workforce planning framework using structured equations

Strategic Workforce Planning is a company process providing best in class, economically sound, workforce management policies and goals. Despite the abundance of literature on the subject, this is a notorious challenge in terms of implementation. Reasons span from the youth of the field itself to broader data integration concerns that arise from gathering information from financial, human resource and business excellence systems. This paper aims at setting the first stones to a simple yet robust quantitative framework for Strategic Workforce Planning exercises. First a method based on structured equations is detailed. It is then used to answer two main workforce related questions: how to optimally hire to keep labor costs flat? How to build an experience constrained workforce at a minimal cost?

preprint2015arXiv

Derivation of the bacterial run-and-tumble kinetic equation from a model with biochemical pathway

Kinetic-transport equations are, by now, standard models to describe the dynamics of populations of bacteria moving by run-and-tumble. Experimental observations show that bacteria increase their run duration when encountering an increasing gradient of chemotactic molecules. This led to a first class of models which heuristically include tumbling frequencies depending on the path-wise gradient of chemotactic signal. More recently, the biochemical pathways regulating the flagellar motors were uncovered. This knowledge gave rise to a second class of kinetic-transport equations, that takes into account an intra-cellular molecular content and which relates the tumbling frequency to this information. It turns out that the tumbling frequency depends on the chemotactic signal, and not on its gradient. For these two classes of models, macroscopic equations of Keller-Segel type, have been derived using diffusion or hyperbolic rescaling. We complete this program by showing how the first class of equations can be derived from the second class with molecular content after appropriate rescaling. The main difficulty is to explain why the path-wise gradient of chemotactic signal can arise in this asymptotic process. Randomness of receptor methylation events can be included, and our approach can be used to compute the tumbling frequency in presence of such a noise.

preprint2015arXiv

Semi-discretization for stochastic scalar conservation laws with multiple rough fluxes

We develop a semi-discretization approximation for scalar conservation laws with multiple rough time dependence in inhomogeneous fluxes. The method is based on Brenier's transport-collapse algorithm and uses characteristics defined in the setting of rough paths. We prove strong $L^1$-convergence for inhomogeneous fluxes and provide a rate of convergence for homogeneous one's. The approximation scheme as well as the proofs are based on the recently developed theory of pathwise entropy solutions and uses the kinetic formulation which allows to define globally the (rough) characteristics.

preprint2014arXiv

Derivation of a Hele-Shaw type system from a cell model with active motion

We formulate a Hele-Shaw type free boundary problem for a tumor growing under the combined effects of pressure forces, cell multiplication and active motion, the latter being the novelty of the present paper. This new ingredient is considered here as a standard diffusion process. The free boundary model is derived from a description at the cell level using the asymptotic of a stiff pressure limit. Compared to the case when active motion is neglected, the pressure satisfies the same complementarity Hele-Shaw type formula. However, the cell density is smoother (Lipschitz continuous), while there is a deep change in the free boundary velocity, which is no longer given by the gradient of the pressure, because some kind of 'mushy region' prepares the tumor invasion.

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

Traveling wave solution of the Hele-Shaw model of tumor growth with nutrient

Several mathematical models of tumor growth are now commonly used to explain medical observations and predict cancer evolution based on images. These models incorporate mechanical laws for tissue compression combined with rules for nutrients availability which can differ depending on the situation under consideration, in vivo or in vitro. Numerical solutions exhibit, as expected from medical observations, a proliferative rim and a necrotic core. However, their precise profiles are rather complex, both in one and two dimensions.

preprint2013arXiv

Effects of space structure and combination therapies on phenotypic heterogeneity and drug resistance in solid tumors

Histopathological evidence supports the idea that the emergence of phenotypic heterogeneity and resistance to cytotoxic drugs can be considered as a process of adaptation, or evolution, in tumor cell populations. In this framework, can we explain intra-tumor heterogeneity in terms of cell adaptation to local conditions? How do anti-cancer therapies affect the outcome of cell competition for nutrients within solid tumors? Can we overcome the emergence of resistance and favor the eradication of cancer cells by using combination therapies? Bearing these questions in mind, we develop a model describing cell dynamics inside a tumor spheroid under the effects of cytotoxic and cytostatic drugs. Cancer cells are assumed to be structured as a population by two real variables standing for space position and the expression level of a cytotoxic resistant phenotype. The model takes explicitly into account the dynamics of resources and anti-cancer drugs as well as their interactions with the cell population under treatment. We analyze the effects of space structure and combination therapies on phenotypic heterogeneity and chemotherapeutic resistance. Furthermore, we study the efficacy of combined therapy protocols based on constant infusion and/or bang-bang delivery of cytotoxic and cytostatic drugs.

preprint2013arXiv

The Hele-Shaw asymptotics for mechanical models of tumor growth

Models of tumor growth, now commonly used, present several levels of complexity, both in terms of the biomedical ingredients and the mathematical description. The simplest ones contain competition for space using purely fluid mechanical concepts. Another possible ingredient is the supply of nutrients through vasculature. The models can describe the tissue either at the level of cell densities, or at the scale of the solid tumor, in this latter case by means of a free boundary problem. Our first goal here is to formulate a free boundary model of Hele-Shaw type, a variant including growth terms, starting from the description at the cell level and passing to a certain limit. A detailed mathematical analysis of this purely mechanical model is performed. Indeed, we are able to prove strong convergence in passing to the limit, with various uniform gradient estimates; we also prove uniqueness for the asymptotic Hele-Shaw type problem. The main tools are nonlinear regularizing effects for certain porous medium type equations, regularization techniques à la Steklov, and a Hilbert duality method for uniqueness. At variance with the classical Hele-Shaw problem, here the geometric motion governed by the pressure is not sufficient to completely describe the dynamics. A complete description requires the equation on the cell number density. Using this theory as a basis, we go on to consider the more complex model including nutrients. We obtain the equation for the limit of the coupled system; the method relies on some BV bounds and space/time a priori estimates. Here, new technical difficulties appear, and they reduce the generality of the results in terms of the initial data. Finally, we prove uniqueness for the system, a main mathematical difficulty.

preprint2012arXiv

Nonlinear stability of a Vlasov equation for magnetic plasmas

The mathematical description of laboratory fusion plasmas produced in Tokamaks is still challenging. Complete models for electrons and ions, as Vlasov-Maxwell systems, are computationally too expensive because they take into account all details and scales of magneto-hydrodynamics. In particular, for most of the relevant studies, the mass electron is negligible and the velocity of material waves is much smaller than the speed of light. Therefore it is useful to understand simplified models. Here we propose and study one of those which keeps both the complexity of the Vlasov equation for ions and the Hall effect in Maxwell's equation. Based on energy dissipation, a fundamental physical property, we show that the model is nonlinear stable and consequently prove existence.

preprint2011arXiv

Relaxation and self-sustained oscillations in the time elapsed neuron network model

The time elapsed model describes the firing activity of an homogeneous assembly of neurons thanks to the distribution of times elapsed since the last discharge. It gives a mathematical description of the probability density of neurons structured by this time. In an earlier work, based on generalized relative entropy methods, it is proved that for highly or weakly connected networks the model exhibits relaxation to the steady state and for moderately connected networks it is obtained numerical evidence of appearance of self-sustained periodic solutions. Here, we go further and, using the particular form of the model, we quantify the regime where relaxation to a stationary state occurs in terms of the network connectivity. To introduce our methodology, we first consider the case where the neurons are not connected and we give a new statement showing that total asynchronous firing of neurons appears asymptotically. In a second step, we consider the case with connections and give a low connectivity condition that still leads to asynchronous firing. Our low connectivity condition is somehow sharp because we can give an example, when this condition is not fulfilled, where synchronous rhythmic activity occurs. Indeed, we are able to build several explicit families of periodic solutions. Our construction is fully nonlinear and the resynchronization of the neural activity in the network does not follow from bifurcation analysis. It relies on an algebraically nonlinear boundary condition that occurs in the model.These analytic results are compared with numerical simulations under broader hypotheses and shown to be robust.

preprint2010arXiv

Analysis of Nonlinear Noisy Integrate\&Fire Neuron Models: blow-up and steady states

Nonlinear Noisy Leaky Integrate and Fire (NNLIF) models for neurons networks can be written as Fokker-Planck-Kolmogorov equations on the probability density of neurons, the main parameters in the model being the connectivity of the network and the noise. We analyse several aspects of the NNLIF model: the number of steady states, a priori estimates, blow-up issues and convergence toward equilibrium in the linear case. In particular, for excitatory networks, blow-up always occurs for initial data concentrated close to the firing potential. These results show how critical is the balance between noise and excitatory/inhibitory interactions to the connectivity parameter.

preprint2009arXiv

Stability Analysis of a Simplified Yet Complete Model for Chronic Myelegenous Leukemia

We analyze the asymptotic behavior of a partial differential equation (PDE) model for hematopoiesis. This PDE model is derived from the original agent-based model formulated by (Roeder et al., Nat. Med., 2006), and it describes the progression of blood cell development from the stem cell to the terminally differentiated state. To conduct our analysis, we start with the PDE model of (Kim et al, JTB, 2007), which coincides very well with the simulation results obtained by Roeder et al. We simplify the PDE model to make it amenable to analysis and justify our approximations using numerical simulations. An analysis of the simplified PDE model proves to exhibit very similar properties to those of the original agent-based model, even if for slightly different parameters. Hence, the simplified model is of value in understanding the dynamics of hematopoiesis and of chronic myelogenous leukemia, and it presents the advantage of having fewer parameters, which makes comparison with both experimental data and alternative models much easier.

preprint2008arXiv

Numerical Solution of an Inverse Problem in Size-Structured Population Dynamics

We consider a size-structured model for cell division and address the question of determining the division (birth) rate from the measured stable size distribution of the population. We propose a new regularization technique based on a filtering approach. We prove convergence of the algorithm and validate the theoretical results by implementing numerical simulations, based on classical techniques. We compare the results for direct and inverse problems, for the filtering method and for the quasi-reversibility method proposed in [Perthame-Zubelli].