Well posedness and Maximum Entropy Approximation for the Dynamics of Quantitative Traits
arXiv:1704.08757 · doi:10.1016/j.physd.2017.10.015
Abstract
We study the Fokker-Planck equation derived in the large system limit of the Markovian process describing the dynamics of quantitative traits. The Fokker-Planck equation is posed on a bounded domain and its transport and diffusion coefficients vanish on the domain's boundary. We first argue that, despite this degeneracy, the standard no-flux boundary condition is valid. We derive the weak formulation of the problem and prove the existence and uniqueness of its solutions by constructing the corresponding contraction semigroup on a suitable function space. Then, we prove that for the parameter regime with high enough mutation rate the problem exhibits a positive spectral gap, which implies exponential convergence to equilibrium. Next, we provide a simple derivation of the so-called Dynamic Maximum Entropy (DynMaxEnt) method for approximation of moments of the Fokker-Planck solution, which can be interpreted as a nonlinear Galerkin approximation. The limited applicability of the DynMaxEnt method inspires us to introduce its modified version that is valid for the whole range of admissible parameters. Finally, we present several numerical experiments to demonstrate the performance of both the original and modified DynMaxEnt methods. We observe that in the parameter regimes where both methods are valid, the modified one exhibits slightly better approximation properties compared to the original one.
28 pages, 4 tables, 5 figures
References in corpus (8)
- Identification of direct residue contacts in protein-protein interaction by message passing
- Continuous data assimilation for the three-dimensional Brinkman-Forchheimer-extended Darcy model
- A Computational Study of a Data Assimilation Algorithm for the Two-dimensional Navier--Stokes Equations
- Continuous Data Assimilation for a 2D Bénard Convection System through Horizontal Velocity Measurements Alone
- The statistical mechanics of a polygenic characterunder stabilizing selection, mutation and drift
- A general approximation for the dynamics of quantitative traits
- Schrodinger Operators with Purely Discrete Spectrum
- Postprocessing Galerkin method applied to a data assimilation algorithm: a uniform in time error estimate