skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Perturbative field-theoretical analysis of three-species cyclic predator-prey models
We apply a perturbative Doi-Peliti field-theoretical analysis to the stochastic spatially extended symmetric Rock-paper-Scissors (RPS) and May-Leonard (ML) models, in which three species compete cyclically. Compared to the two-species Lotka-Volterra predator-prey (LV) model, according to numerical simulations, these cyclical models appear to be less affected by intrinsic stochastic fluctuations. Indeed, we demonstrate that the qualitative features of the ML model are insensitive to intrinsic reaction noise. In contrast, and although not yet observed in numerical simulations, we find that the RPS model acquires significant fluctuation-induced renormalizations in the perturbative regime, similar to the LV model. We also study the formation of spatio-temporal structures in the framework of stability analysis and provide a clearcut explanation for the absence of spatial patterns in the RPS system, whereas the spontaneous emergence of spatio-temporal structures features prominently in the LV and the ML models.  more » « less
Award ID(s):
2128587
PAR ID:
10505743
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
IOP Publ.
Date Published:
Journal Name:
Journal of Physics A: Mathematical and Theoretical
Volume:
56
Issue:
22
ISSN:
1751-8113
Page Range / eLocation ID:
225001
Subject(s) / Keyword(s):
predator-prey model cyclic competition field-theoretical analysis pattern formation fluctuation-induced behavior
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We study the stochastic spatial Lotka-Volterra model for predator-prey interaction subject to a periodically varying carrying capacity. The Lotka-Volterra model with on-site lattice occupation restrictions (i.e., finite local carrying capacity) that represent finite food resources for the prey population exhibits a continuous active-to-absorbing phase transition. The active phase is sustained by the existence of spatio-temporal patterns in the form of pursuit and evasion waves. Monte Carlo simulations on a two-dimensional lattice are utilized to investigate the effect of seasonal variations of the environment on species coexistence. The results of our simulations are also compared to a mean-field analysis in order to specifically delineate the impact of stochastic fluctuations and spatial correlations. We find that the parameter region of predator and prey coexistence is enlarged relative to the stationary situation when the carrying capacity varies periodically. The (quasi-)stationary regime of our periodically varying Lotka-Volterra predator-prey system shows qualitative agreement between the stochastic model and the mean-field approximation. However, under periodic carrying capacity-switching environments, the mean-field rate equations predict period-doubling scenarios that are washed out by internal reaction noise in the stochastic lattice model. Utilizing visual representations of the lattice simulations and dynamical correlation functions, we study how the pursuit and evasion waves are affected by ensuing resonance effects. Correlation function measurements indicate a time delay in the response of the system to sudden changes in the environment. Resonance features are observed in our simulations that cause prolonged persistent spatial correlations. Different effective static environments are explored in the extreme limits of fast and slow periodic switching. The analysis of the mean-field equations in the fast-switching regime enables a semi-quantitative description of the (quasi-)stationary state. 
    more » « less
  2. Abstract Building realistically complex models of infectious disease transmission that are relevant for informing public health is conceptually challenging and requires knowledge of coding architecture that can implement key modeling conventions. For example, many of the models built to understand COVID-19 dynamics have included stochasticity, transmission dynamics that change throughout the epidemic due to changes in host behavior or public health interventions, and spatial structures that account for important spatio-temporal heterogeneities. Here we introduce an R package, SPARSEMODr, that allows users to simulate disease models that are stochastic and spatially explicit, including a model for COVID-19 that was useful in the early phases of the epidemic. SPARSEMOD stands for SPAtial Resolution-SEnsitive Models of Outbreak Dynamics, and our goal is to demonstrate particular conventions for rapidly simulating the dynamics of more complex, spatial models of infectious disease. In this report, we outline the features and workflows of our software package that allow for user-customized simulations. We believe the example models provided in our package will be useful in educational settings, as the coding conventions are adaptable, and will help new modelers to better understand important assumptions that were built into sophisticated COVID-19 models. 
    more » « less
  3. Spiral waves are striking self-organized coherent structures that organize spatio-temporal dynamics in dissipative, spatially extended systems. In this paper, we provide a conceptual approach to various properties of spiral waves. Rather than studying existence in a specific equation, we study properties of spiral waves in general reaction-diffusion systems. We show that many features of spiral waves are robust and to some extent independent of the specific model analyzed. To accomplish this, we present a suitable analytic framework, spatial radial dynamics, that allows us to rigorously characterize features such as the shape of spiral waves and their eigenfunctions, properties of the linearization, and finite-size effects. We believe that our framework can also be used to study spiral waves further and help analyze bifurcations, as well as provide guidance and predictions for experiments and numerical simulations. From a technical point of view, we introduce non-standard function spaces for the well-posedness of the existence problem which allow us to understand properties of spiral waves using dynamical systems techniques, in particular exponential dichotomies. Using these pointwise methods, we are able to bring tools from the analysis of one-dimensional coherent structures such as fronts and pulses to bear on these inherently two-dimensional defects. 
    more » « less
  4. Abstract Computational fluid dynamics (CFD) modeling of left ventricle (LV) flow combined with patient medical imaging data has shown great potential in obtaining patient-specific hemodynamics information for functional assessment of the heart. A typical model construction pipeline usually starts with segmentation of the LV by manual delineation followed by mesh generation and registration techniques using separate software tools. However, such approaches usually require significant time and human efforts in the model generation process, limiting large-scale analysis. In this study, we propose an approach toward fully automating the model generation process for CFD simulation of LV flow to significantly reduce LV CFD model generation time. Our modeling framework leverages a novel combination of techniques including deep-learning based segmentation, geometry processing, and image registration to reliably reconstruct CFD-suitable LV models with little-to-no user intervention.1 We utilized an ensemble of two-dimensional (2D) convolutional neural networks (CNNs) for automatic segmentation of cardiac structures from three-dimensional (3D) patient images and our segmentation approach outperformed recent state-of-the-art segmentation techniques when evaluated on benchmark data containing both magnetic resonance (MR) and computed tomography(CT) cardiac scans. We demonstrate that through a combination of segmentation and geometry processing, we were able to robustly create CFD-suitable LV meshes from segmentations for 78 out of 80 test cases. Although the focus on this study is on image-to-mesh generation, we demonstrate the feasibility of this framework in supporting LV hemodynamics modeling by performing CFD simulations from two representative time-resolved patient-specific image datasets. 
    more » « less
  5. Abstract Solidification phenomenon has been an integral part of the manufacturing processes of metals, where the quantification of stochastic variations and manufacturing uncertainties is critically important. Accurate molecular dynamics (MD) simulations of metal solidification and the resulting properties require excessive computational expenses for probabilistic stochastic analyses where thousands of random realizations are necessary. The adoption of inadequate model sizes and time scales in MD simulations leads to inaccuracies in each random realization, causing a large cumulative statistical error in the probabilistic results obtained through Monte Carlo (MC) simulations. In this work, we present a machine learning (ML) approach, as a data-driven surrogate to MD simulations, which only needs a few MD simulations. This efficient yet high-fidelity ML approach enables MC simulations for full-scale probabilistic characterization of solidified metal properties considering stochasticity in influencing factors like temperature and strain rate. Unlike conventional ML models, the proposed hybrid polynomial correlated function expansion here, being a Bayesian ML approach, is data efficient. Further, it can account for the effect of uncertainty in training data by exploiting mean and standard deviation of the MD simulations, which in principle addresses the issue of repeatability in stochastic simulations with low variance. Stochastic numerical results for solidified aluminum are presented here based on complete probabilistic uncertainty quantification of mechanical properties like Young’s modulus, yield strength and ultimate strength, illustrating that the proposed error-inclusive data-driven framework can reasonably predict the properties with a significant level of computational efficiency. 
    more » « less