Embedding of Markov matrices for d ⩽ 4.

The embedding problem of Markov matrices in Markov semigroups is a classic problem that regained a lot of impetus and activities through recent needs in phylogeny and population genetics. Here, we give an account for dimensions d ⩽ 4 , including a complete and simplified treatment of the case d = 3 , and derive the results in a systematic fashion, with an eye on the potential applications. Further, we reconsider the setup of the corresponding problem for time-inhomogeneous Markov chains, which is needed for real-world applications because transition rates need not be constant over time. Additional cases of this more general embedding occur for any d ⩾ 3 . We review the known case of d = 3 and describe the setting for future work on d = 4 .

Territorial behaviour of buzzards versus random matrix spacing distributions.

A deeper understanding of the processes underlying the distribution of animals in space is crucial for both basic and applied ecology. The Common buzzard (Buteo buteo) is a highly aggressive, territorial bird of prey that interacts strongly with its intra- and interspecific competitors. We propose and use random matrix theory to quantify the strength and range of repulsion as a function of the buzzard population density, thus providing a novel approach to model density dependence. As an indicator of territorial behaviour, we perform a large-scale analysis of the distribution of buzzard nests in an area of 300 square kilometres around the Teutoburger Wald, Germany, as gathered over a period of 20 years. The nearest and next-to-nearest neighbour spacing distribution between nests is compared to the two-dimensional Poisson distribution, originating from uncorrelated random variables, to the complex eigenvalues of random matrices, which are strongly correlated, and to a two-dimensional Coulomb gas interpolating between these two. A one-parameter fit to a time-moving average reveals a significant increase of repulsion between neighbouring nests, as a function of the observed increase in absolute population density over the monitored period of time, thereby proving an unexpected yet simple model for density-dependent spacing of predator territories. A similar effect is obtained for next-to-nearest neighbours, albeit with weaker repulsion, indicating a short-range interaction. Our results show that random matrix theory might be useful in the context of population ecology.

Inflation versus projection sets in aperiodic systems: the role of the window in averaging and diffraction.

Tilings based on the cut-and-project method are key model systems for the description of aperiodic solids. Typically, quantities of interest in crystallography involve averaging over large patches, and are well defined only in the infinite-volume limit. In particular, this is the case for autocorrelation and diffraction measures. For cut-and-project systems, the averaging can conveniently be transferred to internal space, which means dealing with the corresponding windows. In this topical review, this is illustrated by the example of averaged shelling numbers for the Fibonacci tiling, and the standard approach to the diffraction for this example is recapitulated. Further, recent developments are discussed for cut-and-project structures with an inflation symmetry, which are based on an internal counterpart of the renormalization cocycle. Finally, a brief review is given of the notion of hyperuniformity, which has recently gained popularity, and its application to aperiodic structures.

Mathematical diffraction of aperiodic structures.

Kinematic diffraction is well suited for a mathematical approach via measures, which has substantially been developed since the discovery of quasicrystals. The need for further insight emerged from the question of which distributions of matter, beyond perfect crystals, lead to pure point diffraction, hence to sharp Bragg peaks only. More recently, it has become apparent that one also has to study continuous diffraction in more detail, with a careful analysis of the different types of diffuse scattering involved. In this review, we summarise some key results, with particular emphasis on non-periodic structures. We choose an exposition on the basis of characteristic examples, while we refer to the existing literature for proofs and further details.

Single-crossover recombination in discrete time.

Modelling the process of recombination leads to a large coupled nonlinear dynamical system. Here, we consider a particular case of recombination in discrete time, allowing only for single crossovers. While the analogous dynamics in continuous time admits a closed solution (Baake and Baake in Can J Math 55:3-41, 2003), this no longer works for discrete time. A more general model (i.e. without the restriction to single crossovers) has been studied before (Bennett in Ann Hum Genet 18:311-317, 1954; Dawson in Theor Popul Biol 58:1-20, 2000; Linear Algebra Appl 348:115-137, 2002) and was solved algorithmically by means of Haldane linearisation. Using the special formalism introduced by Baake and Baake (Can J Math 55:3-41, 2003), we obtain further insight into the single-crossover dynamics and the particular difficulties that arise in discrete time. We then transform the equations to a solvable system in a two-step procedure: linearisation followed by diagonalisation. Still, the coefficients of the second step must be determined in a recursive manner, but once this is done for a given system, they allow for an explicit solution valid for all times.

Discrete tomography of planar model sets.

Discrete tomography is a well-established method to investigate finite point sets, in particular finite subsets of periodic systems. Here, we start to develop an efficient approach for the treatment of finite subsets of mathematical quasicrystals. To this end, the class of cyclotomic model sets is introduced, and the corresponding consistency, reconstruction and uniqueness problems of the discrete tomography of these sets are discussed.

Signal analysis of impulse response functions in MR and CT measurements of cerebral blood flow.

The impulse response function (IRF) of a localized bolus in cerebral blood flow codes important information on the tissue type. It is indirectly accessible both from MR and CT imaging methods, at least in principle. In practice, however, noise and limited signal resolution render standard deconvolution techniques almost useless. Parametric signal descriptions look more promising, and it is the aim of this contribution to develop some improvements along this line.

An asymptotic maximum principle for essentially linear evolution models.

Recent work on mutation-selection models has revealed that, under specific assumptions on the fitness function and the mutation rates, asymptotic estimates for the leading eigenvalue of the mutation-reproduction matrix may be obtained through a low-dimensional maximum principle in the limit N-->infinity (where N, or N(d) with d> or =1, is proportional to the number of types). In order to extend this variational principle to a larger class of models, we consider here a family of reversible matrices of asymptotic dimension N(d) and identify conditions under which the high-dimensional Rayleigh-Ritz variational problem may be reduced to a low-dimensional one that yields the leading eigenvalue up to an error term of order 1/N. For a large class of mutation-selection models, this implies estimates for the mean fitness, as well as a concentration result for the ancestral distribution of types.

Unequal crossover dynamics in discrete and continuous time.

We analyze a class of models for unequal crossover (UC) of sequences containing sections with repeated units that may differ in length. In these, the probability of an 'imperfect' alignment, in which the shorter sequence has d units without a partner in the longer one, scales like qd as compared to 'perfect' alignments where all these copies are paired. The class is parameterized by this penalty factor q. An effectively infinite population size and thus deterministic dynamics is assumed. For the extreme cases q = 0 and q = 1, and any initial distribution whose moments satisfy certain conditions, we prove the convergence to one of the known fixed points, uniquely determined by the mean copy number, in both discrete and continuous time. For the intermediate parameter values, the existence of fixed points is shown.