The Validity of the Coalescent Approximation for Large Samples

The Kingman coalescent, widely used in genetics, is known to be a good approximation when the sample size is small relative to the population size. In this article, we investigate how large the sample size can get without violating the coalescent approximation. If the haploid population size is 2N, we prove that for samples of size N1/3−ϵ, ϵ > 0, coalescence under the Wright-Fisher (WF) model converges in probability to the Kingman coalescent in the limit of large N. For samples of size N2/5−ϵ or smaller, the WF coalescent converges to a mixture of the Kingman coalescent and what we call the mod-2 coalescent. For samples of size N1/2 or larger, triple collisions in the WF genealogy of the sample become important. The sample size for which the probability of conformance with the Kingman coalescent is 95% is found to be 1.47 × N0.31 for N ∈ [103, 105], showing the pertinence of the asymptotic theory. The probability of no triple collisions is found to be 95% for sample sizes equal to 0.92 × N0.49, which too is in accord with the asymptotic theory.Varying population sizes are handled using algorithms that calculate the probability of WF coalescence agreeing with the Kingman model or taking place without triple collisions. For a sample of size 100, the probabilities of coalescence according to the Kingman model are 2%, 0%, 1%, and 0% in four models of human population with constant N, constant N except for two bottlenecks, recent exponential growth, and increasing recent exponential growth, respectively. For the same four demographic models and the same sample size, the probabilities of coalescence with no triple collision are 92%, 73%, 88%, and 87%, respectively. Visualizations of the algorithm show that even distant bottlenecks can impede agreement between the coalescent and the WF model.Finally, we prove that the WF sample frequency spectrum for samples of size N1/3−ϵ or smaller converges to the classical answer for the coalescent.

A fractional model of a dynamical Brusselator reaction-diffusion system arising in triple collision and enzymatic reactions

In this paper, we study a dynamical Brusselator reaction-diffusion system arising in triple collision and enzymatic reactions with time fractional Caputo derivative. The present article involves a more generalized effective approach, proposed for the Brusselator system say

Exact solution of the planar motion of three arbitrary point vortices

We give an exact quantitative solution for the motion of three vortices of any strength, which Poincaré showed to be integrable. The absolute motion of one vortex is generally biperiodic: in uniformly rotating axes, the motion is periodic. There are two kinds of relative equilibrium configuration: two equilateral triangles and one or three colinear configurations, their stability conditions split the strengths space into three domains in which the sets of trajectories are topologically distinct. According to the values of the strengths and the initial positions, all the possible motions are classified. Two sets of strengths lead to generic motions other than biperiodic. First, when the angular momentum vanishes, besides the biperiodic regime there exists an expansion spiral motion and even a triple collision in a finite time, but the latter motion is nongeneric. Second, when two strengths are opposite, the system also exhibits the elastic diffusion of a vortex doublet by the third vortex. For given values of the invariants, the volume of the phase space of this Hamiltonian system is proportional to the period of the reduced motion, a well known result of the theory of adiabatic invariants. We then formally examine the behaviour of the quantities that Onsager defined only for a large number of interacting vortices.