The free energy method and the Wright-Fisher model with 2 alleles

We systematically investigate the Wright-Fisher model of population genetics with the free energy functional formalism of statistical mechanics and in the light of recent mathematical work on the connection between Fokker-Planck equations and free energy functionals. In statistical physics, entropy increases, or equivalently, free energy decreases, and the asymptotic state is given by a Gibbs-type distribution. This also works for the Wright-Fisher model when rewritten in divergence to identify the correct free energy functional. We not only recover the known results about the stationary distribution, that is, the asymptotic equilibrium state of the model, in the presence of positive mutation rates and possibly also selection, but can also provide detailed formulae for the rate of convergence towards that stationary distribution. In the present paper, the method is illustrated for the simplest case only, that of two alleles.