Source author record

Ira B. Schwartz

Ira B. Schwartz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

36works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

36 published item(s)

preprint2022arXiv

Extreme outbreak dynamics in epidemic models

Motivated by recent epidemic outbreaks, including those of COVID-19, we solve the canonical problem of calculating the dynamics and likelihood of extensive outbreaks in a population within a large class of stochastic epidemic models with demographic noise, including the Susceptible-Infected-Recovered (SIR) model and its general extensions. In the limit of large populations, we compute the probability distribution for all extensive outbreaks, including those that entail unusually large or small (extreme) proportions of the population infected. Our approach reveals that, unlike other well-known examples of rare events occurring in discrete-state stochastic systems, the statistics of extreme outbreaks emanate from a full continuum of Hamiltonian paths, each satisfying unique boundary conditions with a conserved probability flux.

preprint2021arXiv

Optimal periodic closure for minimizing risk in emerging disease outbreaks

Without vaccines and treatments, societies must rely on non-pharmaceutical intervention strategies to control the spread of emerging diseases such as COVID-19. Though complete lockdown is epidemiologically effective, because it eliminates infectious contacts, it comes with significant costs. Several recent studies have suggested that a plausible compromise strategy for minimizing epidemic risk is periodic closure, in which populations oscillate between wide-spread social restrictions and relaxation. However, no underlying theory has been proposed to predict and explain optimal closure periods as a function of epidemiological and social parameters. In this work we develop such an analytical theory for SEIR-like model diseases, showing how characteristic closure periods emerge that minimize the total outbreak, and increase predictably with the reproductive number and incubation periods of a disease, as long as both are within predictable limits. Using our approach we demonstrate a sweet-spot effect in which optimal periodic closure is maximally effective for diseases with similar incubation and recovery periods. Our results compare well to numerical simulations, including in COVID-19 models where infectivity and recovery show significant variability.

preprint2020arXiv

Stability of milling patterns in self-propelled swarms on surfaces

In some physical and biological swarms, agents effectively move and interact along curved surfaces. The associated constraints and symmetries can affect collective-motion patterns, but little is known about pattern stability in the presence of surface curvature. To make progress, we construct a general model for self-propelled swarms moving on surfaces using Lagrangian mechanics. We find that the combination of self-propulsion, friction, mutual attraction, and surface curvature produce milling patterns where each agent in a swarm oscillates on a limit cycle, with different agents splayed along the cycle such that the swarm's center-of-mass remains stationary. In general, such patterns loose stability when mutual attraction is insufficient to overcome the constraint of curvature, and we uncover two broad classes of stationary milling-state bifurcations. In the first, a spatially periodic mode undergoes a Hopf bifurcation as curvature is increased which results in unstable spatiotemporal oscillations. This generic bifurcation is analyzed for the sphere and demonstrated numerically for several surfaces. In the second, a saddle-node-of-periodic-orbits occurs in which stable and unstable milling states collide and annihilate. The latter is analyzed for milling states on cylindrical surfaces. Our results contribute to the general understanding of swarm pattern-formation and stability in the presence of surface curvature, and may aid in designing robotic swarms that can be controlled to move over complex surfaces.

preprint2020arXiv

Unstable oscillations and bistability in delay-coupled swarms

It is known from both theory and experiments that introducing time delays into the communication network of mobile-agent swarms produces coherent rotational patterns. Often such spatio-temporal rotations can be bistable with other swarming patterns, such as milling and flocking. Yet, most known bifurcation results related to delay-coupled swarms rely on inaccurate mean-field techniques. As a consequence, the utility of applying macroscopic theory as a guide for predicting and controlling swarms of mobile robots has been limited. To overcome this limitation, we perform an exact stability analysis of two primary swarming patterns in a general model with time-delayed interactions. By correctly identifying the relevant spatio-temporal modes that determine stability in the presence of time delay, we are able to accurately predict bistability and unstable oscillations in large swarm simulations-- laying the groundwork for comparisons to robotics experiments.

preprint2016arXiv

Collective motion patterns of swarms with delay coupling: theory and experiment

The formation of coherent patterns in swarms of interacting self-propelled autonomous agents is a subject of great interest in a wide range of application areas, ranging from engineering and physics to biology. In this paper, we model and experimentally realize a mixed-reality large-scale swarm of delay-coupled agents. The coupling term is modeled as a delayed communication relay of position. Our analyses, assuming agents communicating over an Erdos-Renyi network, demonstrate the existence of stable coherent patterns that can only be achieved with delay coupling and that are robust to decreasing network connectivity and heterogeneity in agent dynamics. We also show how the bifurcation structure for emergence of different patterns changes with heterogeneity in agent acceleration capabilities and limited connectivity in the network as a function of coupling strength and delay. Our results are verified through simulation as well as preliminary experimental results of delay-induced pattern formation in a mixed-reality swarm.

preprint2016arXiv

Epidemic extinction and control in heterogeneous networks

We consider epidemic extinction in finite networks with broad variation in local connectivity. Generalizing the theory of large fluctuations to random networks with a given degree distribution, we are able to predict the most probable, or optimal, paths to extinction in various configurations, including truncated power-laws. We find that paths for heterogeneous networks follow a limiting form in which infection first decreases in low-degree nodes, which triggers a rapid extinction in high- degree nodes, and finishes with a residual low-degree extinction. The usefulness of the approach is further demonstrated through optimal control strategies that leverage finite-size fluctuations. Interestingly, we find that the optimal control is a mix of treating both high and low-degree nodes based on large-fluctuation theoretical predictions.

preprint2016arXiv

Hybrid dynamics in delay-coupled swarms with "mothership" networks

Swarming behavior continues to be a subject of immense interest because of its centrality in many naturally occurring systems in physics and biology, as well as its importance in applications such as robotics. Here we examine the effects on swarm pattern formation from delayed communication and topological heterogeneity, and in particular, the inclusion of a relatively small number of highly connected nodes, or "motherships", in a swarm's communication network. We find generalized forms of basic patterns for networks with general degree distributions, and a variety of new behaviors including new parameter regions with multi-stability and hybrid motions in bimodal networks. The latter is an interesting example of how heterogeneous networks can have dynamics that is a mix of different states in homogeneous networks, where high and low-degree nodes have simultaneously distinct behavior.

preprint2015arXiv

Pattern-Acquisition in Finite, Heterogenous, Delay-Coupled Swarms

The self-organizing behavior of swarms of inter- acting particles or agents is a topic of intense research in fields extending from biology to physics and robotics. In this paper, we carry out a systematic study of how the stable spatio-temporal patterns of a swarm depend on the number of agents, in the presence of time-delayed interaction and agent heterogeneity. We show how a coherent pattern is modified as the number of agents varies, including the time required to converge to it starting from random initial conditions. We discuss the implications of our results for curve-tracking using autonomous robotic systems.

preprint2014arXiv

Capturing pattern bi-stability dynamics in delay-coupled swarms

Swarms of large numbers of agents appear in many biological and engineering fields. Dynamic bi-stability of co-existing spatio-temporal patterns has been observed in many models of large population swarms. However, many reduced models for analysis, such as mean-field (MF), do not capture the bifurcation structure of bi-stable behavior. Here, we develop a new model for the dynamics of a large population swarm with delayed coupling. The additional physics predicts how individual particle dynamics affects the motion of the entire swarm. Specifically, (1) we correct the center of mass propulsion physics accounting for the particles velocity distribution; (2) we show that the model we develop is able to capture the pattern bi-stability displayed by the full swarm model.

preprint2014arXiv

Collective motions of heterogeneous swarms

The emerging collective motions of swarms of interacting agents are a subject of great interest in application areas ranging from biology to physics and robotics. In this paper, we conduct a careful analysis of the collective dynamics of a swarm of self-propelled heterogeneous, delay-coupled agents. We show the emergence of collective motion patterns and segregation of populations of agents with different dynamic properties; both of these behaviors (pattern formation and segregation) emerge naturally in our model, which is based on self-propulsion and attractive pairwise interactions between agents. We derive the bifurcation structure for emergence of different swarming behaviors in the mean field as a function of physical parameters and verify these results through simulation.

preprint2014arXiv

Going with the flow: enhancing stochastic switching rates in multi-gyre systems

A control strategy is employed that modifies the stochastic escape times from one basin of attraction to another in a model of a double-gyre flow. The system studied captures the behavior of a large class of fluid flows that circulate and have multiple almost invariant sets. In the presence of noise, a particle in one gyre may randomly switch to an adjacent gyre due to a rare large fluctuation. We show that large fluctuation theory may be applied for controlling autonomous agents in a stochastic environment, in fact leveraging the stochastic- ity to the advantage of switching between regions of interest and concluding that patterns may be broken or held over time as the result of noise. We demonstrate that a controller can effectively manipulate the probability of a large fluctuation, thereby modifying escape times exponentially; this demonstrates the potential of optimal control strategies that work in combination with the endemic stochastic environment. To demonstrate this, stochastic simulations and numerical continuation are employed to tie together experimental findings with predictions.

preprint2014arXiv

Inferring unobserved multistrain epidemic sub-populations using synchronization dynamics

A new method is proposed to infer unobserved epidemic sub-populations by exploiting the synchronization properties of multistrain epidemic models. A model for dengue fever is driven by simulated data from secondary infective populations. Primary infective populations in the driven system synchronize to the correct values from the driver system. Most hospital cases of dengue are secondary infections, so this method provides a way to deduce unobserved primary infection levels. We derive center manifold equations that relate the driven system to the driver system and thus motivate the use of synchronization to predict unobserved primary infectives. Synchronization stability between primary and secondary infections is demonstrated through numerical measurements of conditional Lyapunov exponents and through time series simulations.

preprint2014arXiv

Noise Induced Switching and Extinction in Systems with Delay

We consider the rates of noise-induced switching between the stable states of dissipative dynamical systems with delay and also the rates of noise-induced extinction, where such systems model population dynamics. We study a class of systems where the evolution depends on the dynamical variables at a preceding time with a fixed time delay, which we call hard delay. For weak noise, the rates of inter-attractor switching and extinction are exponentially small. Finding these rates to logarithmic accuracy is reduced to variational problems. The solutions of the variational problems give the most probable paths followed in switching or extinction. We show that the equations for the most probable paths are acausal and formulate the appropriate boundary conditions. Explicit general results are obtained for small delay compared to the relaxation rate. We also develop a direct variational method to find the rates. We find that the analytical results agree well with the numerical simulations for both switching and extinction rates.

preprint2014arXiv

Rare Event Extinction on Stochastic Networks

We consider the problem of extinction processes on random networks with a given structure. For sufficiently large well-mixed populations, the process of extinction of one or more state variable components occurs in the tail of the quasi-stationary probability distribution, thereby making it a rare event. Here we show how to extend the theory of large deviations to random networks to predict extinction times. In particular, we use the theory to find the most probable path leading to extinction. We apply the methodology to epidemic models and discover how mean extinction times scale with epidemiological and network parameters in Erdos-Renyi networks. The results are shown to compare quite well with Monte Carlo simulations of the network in predicting both the most probable paths to extinction and mean extinction times.

preprint2013arXiv

Distributed allocation of mobile sensing swarms in gyre flows

We address the synthesis of distributed control policies to enable a swarm of homogeneous mobile sensors to maintain a desired spatial distribution in a geophysical flow environment, or workspace. In this article, we assume the mobile sensors (or robots) have a "map" of the environment denoting the locations of the Lagrangian coherent structures or LCS boundaries. Based on this information, we design agent-level hybrid control policies that leverage the surrounding fluid dynamics and inherent environmental noise to enable the team to maintain a desired distribution in the workspace. We establish the stability properties of the ensemble dynamics of the distributed control policies. Since realistic quasi-geostrophic ocean models predict double-gyre flow solutions, we use a wind-driven multi-gyre flow model to verify the feasibility of the proposed distributed control strategy and compare the proposed control strategy with a baseline deterministic allocation strategy. Lastly, we validate the control strategy using actual flow data obtained by our coherent structure experimental testbed.

preprint2013arXiv

Intervention-Based Stochastic Disease Eradication

Disease control is of paramount importance in public health with infectious disease extinction as the ultimate goal. Although diseases may go extinct due to random loss of effective contacts where the infection is transmitted to new susceptible individuals, the time to extinction in the absence of control may be prohibitively long. Thus intervention controls, such as vaccination of susceptible individuals and/or treatment of infectives, are typically based on a deterministic schedule, such as periodically vaccinating susceptible children based on school calendars. In reality, however, such policies are administered as a random process, while still possessing a mean period. Here, we consider the effect of randomly distributed intervention as disease control on large finite populations. We show explicitly how intervention control, based on mean period and treatment fraction, modulates the average extinction times as a function of population size and rate of infection spread. In particular, our results show an exponential improvement in extinction times even though the controls are implemented using a random Poisson distribution. Finally, we discover those parameter regimes where random treatment yields an exponential improvement in extinction times over the application of strictly periodic intervention. The implication of our results is discussed in light of the availability of limited resources for control.

preprint2013arXiv

Predicting unobserved exposures from seasonal epidemic data

We consider a stochastic Susceptible-Exposed-Infected-Recovered (SEIR) epidemiological model with a contact rate that fluctuates seasonally. Through the use of a nonlinear, stochastic projection, we are able to analytically determine the lower dimensional manifold on which the deterministic and stochastic dynamics correctly interact. Our method produces a low dimensional stochastic model that captures the same timing of disease outbreak and the same amplitude and phase of recurrent behavior seen in the high dimensional model. Given seasonal epidemic data consisting of the number of infectious individuals, our method enables a data-based model prediction of the number of unobserved exposed individuals over very long times.

preprint2013arXiv

Stochastic switching in slow-fast systems: a large fluctuation approach

In this paper we develop a perturbation method to predict the rate of occurrence of rare events for singularly perturbed stochastic systems using a probability density function approach. In contrast to a stochastic normal form approach, we model rare event occurrences due to large fluctuations probabilistically and employ a WKB ansatz to approximate their rate of occurrence. This results in the generation of a two-point boundary value problem that models the interaction of the state variables and the most likely noise force required to induce a rare event. The resulting equations of motion of describing the phenomenon are shown to be singularly perturbed. Vastly different time scales among the variables are leveraged to reduce the dimension and predict the dynamics on the slow manifold in a deterministic setting. The resulting constrained equations of motion may be used to directly compute an exponent that determines the probability of rare events. To verify the theory, a stochastic damped Duffing oscillator with three equilibrium points (two sinks separated by a saddle) is analyzed. The predicted switching time between states is computed using the optimal path that resides in an expanded phase space. We show that the exponential scaling of the switching rate as a function of system parameters agrees well with numerical simulations. Moreover, the dynamics of the original system and the reduced system via center manifolds are shown to agree in an exponentially scaling sense.

preprint2012arXiv

An iterative action minimizing method for computing optimal paths in stochastic dynamical systems

We present a numerical method for computing optimal transition pathways and transition rates in systems of stochastic differential equations (SDEs). In particular, we compute the most probable transition path of stochastic equations by minimizing the effective action in a corresponding deterministic Hamiltonian system. The numerical method presented here involves using an iterative scheme for solving a two-point boundary value problem for the Hamiltonian system. We validate our method by applying it to both continuous stochastic systems, such as nonlinear oscillators governed by the Duffing equation, and finite discrete systems, such as epidemic problems, which are governed by a set of master equations. Furthermore, we demonstrate that this method is capable of dealing with stochastic systems of delay differential equations.

preprint2012arXiv

Coherent Pattern Prediction in Swarms of Delay-Coupled Agents

We consider a general swarm model of self-propelling agents interacting through a pairwise potential in the presence of noise and communication time delay. Previous work [Phys. Rev. E 77, 035203(R) (2008)] has shown that a communication time delay in the swarm induces a pattern bifurcation that depends on the size of the coupling amplitude. We extend these results by completely unfolding the bifurcation structure of the mean field approximation. Our analysis reveals a direct correspondence between the different dynamical behaviors found in different regions of the coupling-time delay plane with the different classes of simulated coherent swarm patterns. We derive the spatio-temporal scales of the swarm structures, and also demonstrate how the complicated interplay of coupling strength, time delay, noise intensity, and choice of initial conditions can affect the swarm. In particular, our studies show that for sufficiently large values of the coupling strength and/or the time delay, there is a noise intensity threshold that forces a transition of the swarm from a misaligned state into an aligned state. We show that this alignment transition exhibits hysteresis when the noise intensity is taken to be time dependent.

preprint2012arXiv

Disease Persistence in Epidemiological Models: The Interplay between Vaccination and Migration

We consider the interplay of vaccination and migration rates on disease persistence in epidemiological systems. We show that short-term and long-term migration can inhibit disease persistence. As a result, we show how migration changes how vaccination rates should be chosen to maintain herd immunity. In a system of coupled SIR models, we analyze how disease eradication depends explicitly on vaccine distribution and migration connectivity. The analysis suggests potentially novel vaccination policies that underscore the importance of optimal placement of finite resources.

preprint2012arXiv

Noise Induced Pattern Switching in Randomly Distributed Delayed Swarm Patterns

We study the effects of noise on the dynamics of a system of coupled self-propelling particles in the case where the coupling is time-delayed, and the delays are discrete and randomly generated. Previous work has demonstrated that the stability of a class of emerging patterns depends upon all moments of the time delay distribution, and predicts their bifurcation parameter ranges. Near the bifurcations of these patterns, noise may induce a transition from one type of pattern to another. We study the onset of these noise-induced swarm re-organizations by numerically simulating the system over a range of noise intensities and for various distributions of the delays. Interestingly, there is a critical noise threshold above which the system is forced to transition from a less organized state to a more organized one. We explore this phenomenon by quantifying this critical noise threshold, and note that transition time between states varies as a function of both the noise intensity and delay distribution.

preprint2012arXiv

Noise, Bifurcations, and Modeling of Interacting Particle Systems

We consider the stochastic patterns of a system of communicating, or coupled, self-propelled particles in the presence of noise and communication time delay. For sufficiently large environmental noise, there exists a transition between a translating state and a rotating state with stationary center of mass. Time delayed communication creates a bifurcation pattern dependent on the coupling amplitude between particles. Using a mean field model in the large number limit, we show how the complete bifurcation unfolds in the presence of communication delay and coupling amplitude. Relative to the center of mass, the patterns can then be described as transitions between translation, rotation about a stationary point, or a rotating swarm, where the center of mass undergoes a Hopf bifurcation from steady state to a limit cycle. Examples of some of the stochastic patterns will be given for large numbers of particles.

preprint2012arXiv

Randomly Distributed Delayed Communication and Coherent Swarm Patterns

Previously we showed how delay communication between globally coupled self-propelled agents causes new spatio-temporal patterns to arise when the delay coupling is fixed among all agents \cite{Forgoston08}. In this paper, we show how discrete, randomly distributed delays affect the dynamical patterns. In particular, we investigate how the standard deviation of the time delay distribution affects the stability of the different patterns as well as the switching probability between coherent states.

preprint2012arXiv

Recruitment dynamics in adaptive social networks

We model recruitment in adaptive social networks in the presence of birth and death processes. Recruitment is characterized by nodes changing their status to that of the recruiting class as a result of contact with recruiting nodes. Only a susceptible subset of nodes can be recruited. The recruiting individuals may adapt their connections in order to improve recruitment capabilities, thus changing the network structure adaptively. We derive a mean field theory to predict the dependence of the growth threshold of the recruiting class on the adaptation parameter. Furthermore, we investigate the effect of adaptation on the recruitment level, as well as on network topology. The theoretical predictions are compared with direct simulations of the full system. We identify two parameter regimes with qualitatively different bifurcation diagrams depending on whether nodes become susceptible frequently (multiple times in their lifetime) or rarely (much less than once per lifetime).

preprint2012arXiv

Robotic Manifold Tracking of Coherent Structures in Flows

Tracking Lagrangian coherent structures in dynamical systems is important for many applications such as oceanography and weather prediction. In this paper, we present a collaborative robotic control strategy designed to track stable and unstable manifolds. The technique does not require global information about the fluid dynamics, and is based on local sensing, prediction, and correction. The collaborative control strategy is implemented on a team of three robots to track coherent structures and manifolds on static flows as well as a noisy time-dependent model of a wind-driven double-gyre often seen in the ocean. We present simulation and experimental results and discuss theoretical guarantees of the collaborative tracking strategy.

preprint2012arXiv

Statistical multi-moment bifurcations in random delay coupled swarms

We study the effects of discrete, randomly distributed time delays on the dynamics of a coupled system of self-propelling particles. Bifurcation analysis on a mean field approximation of the system reveals that the system possesses patterns with certain universal characteristics that depend on distinguished moments of the time delay distribution. Specifically, we show both theoretically and numerically that although bifurcations of simple patterns, such as translations, change stability only as a function of the first moment of the time delay distribution, more complex patterns arising from Hopf bifurcations depend on all of the moments.

preprint2011arXiv

Converging towards the optimal path to extinction

Extinction appears ubiquitously in many fields, including chemical reactions, population biology, evolution, and epidemiology. Even though extinction as a random process is a rare event, its occurrence is observed in large finite populations. Extinction occurs when fluctuations due to random transitions act as an effective force which drives one or more components or species to vanish. Although there are many random paths to an extinct state, there is an optimal path that maximizes the probability to extinction. In this article, we show that the optimal path is associated with the dynamical systems idea of having maximum sensitive dependence to initial conditions. Using the equivalence between the sensitive dependence and the path to extinction, we show that the dynamical systems picture of extinction evolves naturally toward the optimal path in several stochastic models of epidemics.

preprint2011arXiv

Set-based corral control in stochastic dynamical systems: Making almost invariant sets more invariant

We consider the problem of stochastic prediction and control in a time-dependent stochastic environment, such as the ocean, where escape from an almost invariant region occurs due to random fluctuations. We determine high-probability control-actuation sets by computing regions of uncertainty, almost invariant sets, and Lagrangian Coherent Structures. The combination of geometric and probabilistic methods allows us to design regions of control that provide an increase in loitering time while minimizing the amount of control actuation. We show how the loitering time in almost invariant sets scales exponentially with respect to the control actuation, causing an exponential increase in loitering times with only small changes in actuation force. The result is that the control actuation makes almost invariant sets more invariant.

preprint2010arXiv

Maximal Sensitive Dependence and the Optimal Path to Epidemic Extinction

Extinction of an epidemic or a species is a rare event that occurs due to a large, rare stochastic fluctuation. Although the extinction process is dynamically unstable, it follows an optimal path that maximizes the probability of extinction. We show that the optimal path is also directly related to the finite-time Lyapunov exponents of the underlying dynamical system in that the optimal path displays maximum sensitivity to initial conditions. We consider several stochastic epidemic models, and examine the extinction process in a dynamical systems framework. Using the dynamics of the finite-time Lyapunov exponents as a constructive tool, we demonstrate that the dynamical systems viewpoint of extinction evolves naturally toward the optimal path.

preprint2010arXiv

Switching barrier scaling near bifurcation points for non-Gaussian noise

We study noise-induced switching of a system close to bifurcation parameter values where the number of stable states changes. For non-Gaussian noise, the switching exponent, which gives the logarithm of the switching rate, displays a non-power-law dependence on the distance to the bifurcation point. This dependence is found for Poisson noise. Even weak additional Gaussian noise dominates switching sufficiently close to the bifurcation point, leading to a crossover in the behavior of the switching exponent.

preprint2009arXiv

Accurate Noise Projection for Reduced Stochastic Epidemic Models

We consider a stochastic Susceptible-Exposed-Infected-Recovered (SEIR) epidemiological model. Through the use of a normal form coordinate transform, we are able to analytically derive the stochastic center manifold along with the associated, reduced set of stochastic evolution equations. The transformation correctly projects both the dynamics and the noise onto the center manifold. Therefore, the solution of this reduced stochastic dynamical system yields excellent agreement, both in amplitude and phase, with the solution of the original stochastic system for a temporal scale that is orders of magnitude longer than the typical relaxation time. This new method allows for improved time series prediction of the number of infectious cases when modeling the spread of disease in a population. Numerical solutions of the fluctuations of the SEIR model are considered in the infinite population limit using a Langevin equation approach, as well as in a finite population simulated as a Markov process.

preprint2009arXiv

Enhanced vaccine control of epidemics in adaptive networks

We study vaccine control for disease spread on an adaptive network modeling disease avoidance behavior. Control is implemented by adding Poisson distributed vaccination of susceptibles. We show that vaccine control is much more effective in adaptive networks than in static networks due to an interaction between the adaptive network rewiring and the vaccine application. Disease extinction rates using vaccination are computed, and orders of magnitude less vaccine application is needed to drive the disease to extinction in an adaptive network than in a static one.

preprint2009arXiv

Epidemics with Multistrain Interactions: The Interplay Between Cross Immunity and Antibody-Dependent Enhancement

This paper examines the interplay of the effect of cross immunity and antibody-dependent enhancement (ADE) in mutistrain diseases. Motivated by dengue fever, we study a model for the spreading of epidemics in a population with multistrain interactions mediated by both partial temporary cross immunity and ADE. Although ADE models have previously been observed to cause chaotic outbreaks, we show analytically that weak cross immunity has a stabilizing effect on the system. That is, the onset of disease fluctuations requires a larger value of ADE with small cross immunity than without. However, strong cross immunity is shown numerically to cause oscillations and chaotic outbreaks even for low values of ADE.

preprint2009arXiv

Escape Rates in a Stochastic Environment with Multiple Scales

We consider a stochastic environment with two time scales and outline a general theory that compares two methods to reduce the dimension of the original system. The first method involves the computation of the underlying deterministic center manifold followed by a naive replacement of the stochastic term. The second method allows one to more accurately describe the stochastic effects and involves the derivation of a normal form coordinate transform that is used to find the stochastic center manifold. The results of both methods are used along with the path integral formalism of large fluctuation theory to predict the escape rate from one basin of attraction to another. The general theory is applied to the example of a surface flow described by a generic, singularly perturbed, damped, nonlinear oscillator with additive, Gaussian noise. We show how both nonlinear reduction methods compare in escape rate scaling. Additionally, the center manifolds are shown to predict high pre-history probability regions of escape. The theoretical results are confirmed using numerical computation of the mean escape time and escape prehistory, and we briefly discuss the extension of the theory to stochastic control.

preprint2006arXiv

Synchronized dynamics of cortical neurons with time-delay feedback

The dynamics of three mutually coupled cortical neurons with time delays in the coupling are explored numerically and analytically. The neurons are coupled in a line, with the middle neuron sending a somewhat stronger projection to the outer neurons than the feedback it receives, to model for instance the relay of a signal from primary to higher cortical areas. For a given coupling architecture, the delays introduce correlations in the time series at the time-scale of the delay. It was found that the middle neuron leads the outer ones by the delay time, while the outer neurons are synchronized with zero lag times. Synchronization is found to be highly dependent on the synaptic time constant, with faster synapses increasing both the degree of synchronization and the firing rate. Analysis shows that presynaptic input during the interspike interval stabilizes the synchronous state, even for arbitrarily weak coupling, and independent of the initial phase. The finding may be of significance to synchronization of large groups of cells in the cortex that are spatially distanced from each other.