68 research outputs found

    Self-optimization, community stability, and fluctuations in two individual-based models of biological coevolution

    Full text link
    We compare and contrast the long-time dynamical properties of two individual-based models of biological coevolution. Selection occurs via multispecies, stochastic population dynamics with reproduction probabilities that depend nonlinearly on the population densities of all species resident in the community. New species are introduced through mutation. Both models are amenable to exact linear stability analysis, and we compare the analytic results with large-scale kinetic Monte Carlo simulations, obtaining the population size as a function of an average interspecies interaction strength. Over time, the models self-optimize through mutation and selection to approximately maximize a community fitness function, subject only to constraints internal to the particular model. If the interspecies interactions are randomly distributed on an interval including positive values, the system evolves toward self-sustaining, mutualistic communities. In contrast, for the predator-prey case the matrix of interactions is antisymmetric, and a nonzero population size must be sustained by an external resource. Time series of the diversity and population size for both models show approximate 1/f noise and power-law distributions for the lifetimes of communities and species. For the mutualistic model, these two lifetime distributions have the same exponent, while their exponents are different for the predator-prey model. The difference is probably due to greater resilience toward mass extinctions in the food-web like communities produced by the predator-prey model.Comment: 26 pages, 12 figures. Discussion of early-time dynamics added. J. Math. Biol., in pres

    Point sets on the sphere S2\mathbb{S}^2 with small spherical cap discrepancy

    Full text link
    In this paper we study the geometric discrepancy of explicit constructions of uniformly distributed points on the two-dimensional unit sphere. We show that the spherical cap discrepancy of random point sets, of spherical digital nets and of spherical Fibonacci lattices converges with order N1/2N^{-1/2}. Such point sets are therefore useful for numerical integration and other computational simulations. The proof uses an area-preserving Lambert map. A detailed analysis of the level curves and sets of the pre-images of spherical caps under this map is given

    Quasi-Monte Carlo rules for numerical integration over the unit sphere S2\mathbb{S}^2

    Full text link
    We study numerical integration on the unit sphere S2R3\mathbb{S}^2 \subset \mathbb{R}^3 using equal weight quadrature rules, where the weights are such that constant functions are integrated exactly. The quadrature points are constructed by lifting a (0,m,2)(0,m,2)-net given in the unit square [0,1]2[0,1]^2 to the sphere S2\mathbb{S}^2 by means of an area preserving map. A similar approach has previously been suggested by Cui and Freeden [SIAM J. Sci. Comput. 18 (1997), no. 2]. We prove three results. The first one is that the construction is (almost) optimal with respect to discrepancies based on spherical rectangles. Further we prove that the point set is asymptotically uniformly distributed on S2\mathbb{S}^2. And finally, we prove an upper bound on the spherical cap L2L_2-discrepancy of order N1/2(logN)1/2N^{-1/2} (\log N)^{1/2} (where NN denotes the number of points). This slightly improves upon the bound on the spherical cap L2L_2-discrepancy of the construction by Lubotzky, Phillips and Sarnak [Comm. Pure Appl. Math. 39 (1986), 149--186]. Numerical results suggest that the (0,m,2)(0,m,2)-nets lifted to the sphere S2\mathbb{S}^2 have spherical cap L2L_2-discrepancy converging with the optimal order of N3/4N^{-3/4}

    Broadband Noise Generation by a Vortex Model of Cavity Flow

    No full text
    corecore