Lab Home | Phone | Search
Center for Nonlinear Studies  Center for Nonlinear Studies
 Home 
 People 
 Current 
 Affiliates 
 Alumni 
 Visitors 
 Students 
 Research 
 ICAM-LANL 
 Quantum 
 Publications 
 Publications 
 2007 
 2006 
 2005 
 2004 
 2003 
 2002 
 2001 
 2000 
 <1999 
 Conferences 
 Workshops 
 Sponsorship 
 Talks 
 Colloquia 
 Colloquia Archive 
 Seminars 
 Postdoc Seminars Archive 
 Quantum Lunch 
 CMS Colloquia 
 Q-Mat Seminars 
 Q-Mat Seminars Archive 
 Archive 
 Kac Lectures 
 Dist. Quant. Lecture 
 Ulam Scholar 
 Colloquia 
 
 Jobs 
 Students 
 Summer Research 
 Student Application 
 Visitors 
 Description 
 Past Visitors 
 Services 
 General 
 PD Travel Request 
 
 History of CNLS 
 
 Maps, Directions 
 CNLS Office 
 T-Division 
 LANL 
 
Thursday, December 02, 2010
2:00 PM - 3:00 PM
CNLS Conference Room (TA-3, Bldg 1690)

Postdoc Seminar

A multiscale maximum entropy moment closure for locally regulated space–time point process models on population dynamics.

Michael Raghib
CNLS, T-5

The prevalence of structure in biological populations challenges fundamental assumptions at the heart of continuum models of population dynamics based only on mean densities (local or global). Individual-based models (IBMs) were introduced during the last decade in an attempt to overcome this limitation by following explicitly each individual in the population. Although the IBM approach has been quite useful, the capability to follow each individual usually comes at the expense of analytical tractability, which limits the generality of the statements that can be made. For the specific case of spatial structure in populations of sessile (and identical) organisms, space–time point processes with local regulation seem to cover the middle ground between analytical tractability and a higher degree of biological realism. This approach has shown that simplified representations of fecundity, local dispersal and densitydependent mortality weighted by the local competitive environment are sufficient to generate spatial patterns that mimic field observations. Continuum approximations of these stochastic processes try to distill their fundamental properties, and they keep track of not only mean densities, but also higher order spatial correlations. However, due to the non–linearities involved they result in infinite hierarchies of moment equations. This leads to the problem of finding a ‘moment closure’; that is, an appropriate order of (lower order) truncation, together with a method of expressing the highest orde density not explicitly modelled in the truncated hierarchy in terms of the lower order densities.We use the principle of constrained maximum entropy to derive a closure relationship for truncation at second order using normalisation and the product densities of first and second orders as constraints, and apply it to one such hierarchy. The resulting ‘maxent’ closure is similar to the Kirkwood superposition approximation, or ‘power-3’ closure, but it is complemented with previously unknown correction terms that depend mainly on the avoidance function of an associated Poisson point process over the region for which third order correlations are irreducible. This domain of irreducible triplet correlations is found from an integral equation associated with a normalization constraint. This also serves as a validation check, since a single, non-trivial domain can only be found if the assumptions of the closure are consistent with the predictions of the hierarchy. Comparisons between simulations of the point process, alternative heuristic closures, and the maxent closure show significant improvements in the ability of the truncated hierarchy to predict equilibrium values for mildly aggregated spatial patterns. However, the maxent closure performs comparatively poorly in segregated ones. Although the closure is applied in the context of point processes, the method does not require fixed locations to be valid, and can in principle be applied to problems where the particles move, provided that their correlation functions are stationary in space and time.

Host: Peter Loxley, loxley@lanl.gov