Source author record

Michael J. Ward

Michael J. Ward 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

14works
10topics
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

14 published item(s)

preprint2026arXiv

The Effective Reactivity for Capturing Brownian Motion by Partially Reactive Patches on a Spherical Surface

We analyze the trapping of diffusing ligands, modeled as Brownian particles, by a sphere that has $N$ partially reactive boundary patches, each of small area and arbitrary shape, on an otherwise reflecting boundary. For such a structured target, the partial reactivity of each boundary patch is characterized by a Robin boundary condition, with a local boundary reactivity $κ_i$ for $i=1,\ldots,N$. For any spatial arrangement of well-separated patches on the surface of the sphere, the method of matched asymptotic expansions is used to derive explicit results for the capacitance $C_{\rm T}$ of the structured target, which is valid for any $κ_i>0$. This target capacitance $C_{\rm T}$ is defined in terms of a Green's matrix, which depends on the spatial configuration of patches, the local reactive capacitance $C_i(κ_i)$ of each patch and another coefficient that depends on the local geometry near a patch. The analytical dependence of $C_{i}(κ_i)$ on $κ_i$ is uncovered via a spectral expansion over Steklov eigenfunctions. For circular patches, the latter are readily computed numerically and provide an accurate fully explicit sigmoidal approximation for $C_{i}(κ_i)$. In the homogenization limit of $N\gg 1$ identical uniformly-spaced patches with $κ_i=κ$, we derive an explicit scaling law for the effective capacitance and the effective reactivity of the structured target that is valid in the limit of small patch area fraction. From a comparison with numerical simulations, we show that this scaling law provides a highly accurate approximation over the full range $κ>0$, even when there is only a moderately large number of reactive patches.

preprint2023arXiv

Symmetry-Breaking Bifurcations for Compartmental Reaction Kinetics Coupled by Two Bulk Diffusing Species with Comparable Diffusivities in 2-D

For a 2-D coupled PDE-ODE bulk-cell model, we investigate symmetry-breaking bifurcations that can emerge when two bulk diffusing species are coupled to two-component nonlinear intracellular reactions that are restricted to occur only within a disjoint collection of small circular compartments, or "cells", of a common small radius that are confined in a bounded 2-D domain. Outside of the union of these cells, the two bulk species with comparable diffusivities and bulk degradation rates diffuse and globally couple the spatially segregated intracellular reactions through Robin boundary conditions across the cell boundaries, which depend on certain membrane reaction rates. In the singular limit of a small common cell radius, we construct steady-state solutions for the bulk-cell model and formulate a nonlinear matrix eigenvalue problem that determines the linear stability properties of the steady-states. For a certain spatial arrangement of cells for which the steady-state and linear stability analysis become highly tractable, we construct a symmetric steady-state solution where the steady-states of the intracellular species are the same for each cell. As regulated by the ratio of the membrane reaction rates on the cell boundaries, we show for various specific prototypical intracellular reactions, and for a specific two-cell arrangement, that our 2-D coupled PDE-ODE model admits symmetry-breaking bifurcations from this symmetric steady-state, leading to linearly stable asymmetric patterns, even when the bulk diffusing species have comparable or possibly equal diffusivities. Overall, our analysis shows that symmetry-breaking bifurcations can occur without the large diffusivity ratio requirement for the bulk diffusing species as is well-known from a Turing stability analysis applied to a spatially uniform steady-state for typical two-component activator-inhibitor systems.

preprint2023arXiv

The Emergence of Spatial Patterns for Compartmental Reaction Kinetics Coupled by Two Bulk Diffusing Species with Comparable Diffusivities

Originating from the pioneering study of Alan Turing, the bifurcation analysis predicting spatial pattern formation from a spatially uniform state for diffusing morphogens or chemical species that interact through nonlinear reactions is a central problem in many chemical and biological systems. From a mathematical viewpoint, one key challenge with this theory for two component systems is that stable spatial patterns can typically only occur from a spatially uniform state when a slowly diffusing "activator" species reacts with a much faster diffusing "inhibitor" species. However, from a modeling perspective, this large diffusivity ratio requirement for pattern formation is often unrealistic in biological settings since different molecules tend to diffuse with similar rates in extracellular spaces. As a result, one key long-standing question is how to robustly obtain pattern formation in the biologically realistic case where the time scales for diffusion of the interacting species are comparable. For a coupled 1-D bulk-compartment theoretical model, we investigate the emergence of spatial patterns for the scenario where two bulk diffusing species with comparable diffusivities are coupled to nonlinear reactions that occur only in localized "compartments", such as on the boundaries of a 1-D domain. The exchange between the bulk medium and the spatially localized compartments is modeled by a Robin boundary condition with certain binding rates. As regulated by these binding rates, we show for various specific nonlinearities that our 1-D coupled PDE-ODE model admits symmetry-breaking bifurcations, leading to linearly stable asymmetric steady-state patterns, even when the bulk diffusing species have equal diffusivities. Depending on the form of the nonlinear kinetics, oscillatory instabilities can also be triggered. Moreover, the analysis is extended to treat a periodic chain of compartments.

preprint2020arXiv

Localized Signaling Compartments in 2-D Coupled by a Bulk Diffusion Field: Quorum Sensing and Synchronous Oscillations in the Well-Mixed Limit

We analyze oscillatory instabilities for a coupled PDE-ODE system modeling the communication between localized spatially segregated dynamically active signaling compartments that are coupled through a passive extracellular bulk diffusion field in a bounded 2-D domain. Each signaling compartment is assumed to secrete a chemical into the extracellular medium (bulk region) and it can also sense the concentration of this chemical in the region around its boundary. This feedback from the bulk region, resulting from the entire collection of cells, in turn modifies the intracellular dynamics within each cell. In the limit where the signaling compartments are circular disks with a small common radius $ε\ll 1$ and where the bulk diffusivity is asymptotically large, a matched asymptotic analysis is used to reduce the dimensionless PDE-ODE system into a nonlinear ODE system with global coupling. For Sel'kov reaction kinetics, this ODE system for the intracellular dynamics and the spatial average of the bulk diffusion field is then used to investigate oscillatory instabilities in the dynamics of the cells that are triggered due to the global coupling. In particular, numerical bifurcation software on the ODEs is used to study the overall effect of coupling defective cells (cells that behave differently from the remaining cells) to a group of identical cells. Moreover, when the number of cells is large, the Kuramoto order parameter is computed to predict the degree of phase synchronization of the intracellular dynamics. Quorum sensing behavior, characterized by the onset of collective behavior in the intracellular dynamics as the number of cells increases above a threshold, is also studied. Our analysis shows that the cell population density plays a dual role of triggering and then quenching synchronous oscillations in the intracellular dynamics.

preprint2020arXiv

Pattern Formation and Oscillatory Dynamics in a Two-Dimensional Coupled Bulk-Surface Reaction-Diffusion System

On a two-dimensional circular domain, we analyze the formation of spatio-temporal patterns for a class of coupled bulk-surface reaction-diffusion models for which a passive diffusion process occurring in the interior bulk domain is linearly coupled to a nonlinear reaction-diffusion process on the domain boundary. For this coupled PDE system we construct a radially symmetric steady state solution and from a linearized stability analysis formulate criteria for which this base state can undergo either a Hopf bifurcation, a symmetry-breaking pitchfork (or Turing) bifurcation, or a codimension-two pitchfork-Hopf bifurcation. For each of these three types of bifurcations, a multiple time-scale asymptotic analysis is used to derive normal form amplitude equations characterizing the local branching behavior of spatio-temporal patterns in the weakly nonlinear regime. Among the novel aspects of this weakly nonlinear analysis are the two-dimensionality of the bulk domain, the systematic treatment of arbitrary reaction kinetics restricted to the boundary, the bifurcation parameters which arise in the boundary conditions, and the underlying spectral problem where both the differential operator and the boundary conditions involve the eigenvalue parameter. The normal form theory is illustrated for both Schnakenberg and Brusselator reaction kinetics, and the weakly nonlinear results are favorably compared with numerical bifurcation results and results from time-dependent PDE simulations of the coupled bulk-surface system. Overall, the results show the existence of either subcritical or supercritical Hopf and symmetry-breaking bifurcations, and mixed-mode oscillations characteristic of codimension-two bifurcations. Finally, the formation of global structures such as large amplitude rotating waves is briefly explored through PDE numerical simulations.

preprint2020arXiv

Spot Patterns in the 2-D Schnakenberg Model with Localized Heterogeneities

A hybrid asymptotic-numerical theory is developed to analyze the effect of different types of localized heterogeneities on the existence, linear stability, and slow dynamics of localized spot patterns for the two-component Schnakenberg reaction-diffusion model in a 2-D domain. Two distinct types of localized heterogeneities are considered: a strong localized perturbation of a spatially uniform feed rate and the effect of removing a small hole in the domain, through which the chemical species can leak out. Our hybrid theory reveals a wide range of novel phenomena such as, saddle-node bifurcations for quasi-equilibrium spot patterns that otherwise would not occur for a homogeneous medium, a new type of spot solution pinned at the concentration point of the feed rate, spot self-replication behavior leading to the creation of more than two new spots, and the existence of a creation-annihilation attractor with at most three spots. Depending on the type of localized heterogeneity introduced, localized spots are either repelled or attracted towards the localized defect on asymptotically long time scales. Results for slow spot dynamics and detailed predictions of various instabilities of quasi-equilibrium spot patterns, all based on our hybrid asymptotic-numerical theory, are illustrated and confirmed through extensive full PDE numerical simulations.

preprint2020arXiv

Stable Asymmetric Spike Equilibria for the Gierer-Meinhardt Model with a Precursor Field

Precursor gradients in a reaction-diffusion system are spatially varying coefficients in the reaction-kinetics. Such gradients have been used in various applications, such as the head formation in the Hydra, to model the effect of pre-patterns and to localize patterns in various spatial regions. For the 1-D Gierer-Meinhardt (GM) model we show that a simple precursor gradient in the decay rate of the activator can lead to the existence of stable, asymmetric, two-spike patterns, corresponding to localized peaks in the activator of different heights. This is a qualitatively new phenomena for the GM model, in that asymmetric spike patterns are all unstable in the absence of the precursor field. Through a determination of the global bifurcation diagram of two-spike steady-state patterns, we show that asymmetric patterns emerge from a supercritical symmetry-breaking bifurcation along the symmetric two-spike branch as a parameter in the precursor field is varied. Through a combined analytical-numerical approach we analyze the spectrum of the linearization of the GM model around the two-spike steady-state to establish that portions of the asymmetric solution branches are linearly stable. In this linear stability analysis a new class of vector-valued nonlocal eigenvalue problem (NLEP) is derived and analyzed.

preprint2020arXiv

Synchrony and Oscillatory Dynamics for a 2-D PDE-ODE Model of Diffusion-Sensing with Small Signaling Compartments

We analyze a class of cell-bulk coupled PDE-ODE models, motivated by quorum and diffusion sensing phenomena in microbial systems, that characterize communication between localized spatially segregated dynamically active signaling compartments that have a permeable boundary. Each cell secretes a signaling chemical into the bulk region at a constant rate and receives a feedback of the bulk chemical from the entire collection of cells. This global feedback, which activates signaling pathways within the cells, modifies the intracellular dynamics according to the external environment. The cell secretion and global feedback are regulated by permeability parameters across the cell membrane. For arbitrary reaction-kinetics within each cell, the method of matched asymptotic expansions is used in the limit of small cell radius to construct steady-state solutions of the PDE-ODE model, and to derive a globally coupled nonlinear matrix eigenvalue problem (GCEP) that characterizes the linear stability properties of the steady-states. In the limit of large bulk diffusivity an asymptotic analysis of the PDE-ODE model leads to a limiting ODE system for the spatial average of the concentration in the bulk region that is coupled to the intracellular dynamics within the cells. Results from the linear stability theory and ODE dynamics are illustrated for Sel'kov reaction-kinetics, where the kinetic parameters are chosen so that each cell is quiescent when uncoupled from the bulk medium. For various specific spatial configurations of cells, the linear stability theory is used to construct phase diagrams in parameter space characterizing where a switch-like emergence of intracellular oscillations can occur through a Hopf bifurcation.

preprint2020arXiv

Weakly Nonlinear Analysis of Peanut-Shaped Deformations for Localized Spots of Singularly Perturbed Reaction-Diffusion Systems

Spatially localized 2-D spot patterns occur for a wide variety of two component reaction-diffusion systems in the singular limit of a large diffusivity ratio. Such localized, far-from-equilibrium, patterns are known to exhibit a wide range of different instabilities such as breathing oscillations, spot annihilation, and spot self-replication behavior. Prior numerical simulations of the Schnakenberg and Brusselator systems have suggested that a localized peanut-shaped linear instability of a localized spot is the mechanism initiating a fully nonlinear spot self-replication event. From a development and implementation of a weakly nonlinear theory for shape deformations of a localized spot, it is shown through a normal form amplitude equation that a peanut-shaped linear instability of a steady-state spot solution is always subcritical for both the Schnakenberg and Brusselator reaction-diffusion systems. The weakly nonlinear theory is validated by using the global bifurcation software {\em pde2path} [H.~Uecker et al., Numerical Mathematics: Theory, Methods and Applications, {\bf 7}(1), (2014)] to numerically compute an unstable, non-radially symmetric, steady-state spot solution branch that originates from a symmetry-breaking bifurcation point.

preprint2020arXiv

Weakly Nonlinear Theory for Oscillatory Dynamics in a One-Dimensional PDE-ODE Model of Membrane Dynamics Coupled by a Bulk Diffusion Field

We study the dynamics of systems consisting of two spatially segregated ODE compartments coupled through a one-dimensional bulk diffusion field. For this coupled PDE-ODE system, we first employ a multi-scale asymptotic expansion to derive amplitude equations near codimension-one Hopf bifurcation points for both in-phase and anti-phase synchronization modes. The resulting normal form equations pertain to any vector nonlinearity restricted to the ODE compartments. In our first example, we apply our weakly nonlinear theory to a coupled PDE-ODE system with Sel'kov membrane kinetics, and show that the symmetric steady state undergoes supercritical Hopf bifurcations as the coupling strength and the diffusivity vary. We then consider the PDE diffusive coupling of two Lorenz oscillators. It is shown that this coupling mechanism can have a stabilizing effect, characterized by a significant increase in the Rayleigh number required for a Hopf bifurcation. Within the chaotic regime, we can distinguish between synchronous chaos, where both the left and right oscillators are in-phase, and chaotic states characterized by the absence of synchrony. Finally, we compute the largest Lyapunov exponent associated with a linearization around the synchronous manifold that only considers odd perturbations. This allows us to predict the transition to synchronous chaos as the coupling strength and the diffusivity increase.

preprint2015arXiv

The dynamics of localized spot patterns for reaction-diffusion systems on the sphere

In the singularly perturbed limit corresponding to a large diffusivity ratio between two components in a reaction-diffusion (RD) system, quasi-equilibrium spot patterns are often admitted, producing a solution that concentrates at a discrete set of points in the domain. In this paper, we derive and study the differential algebraic equation (DAE) that characterizes the slow dynamics for such spot patterns for the Brusselator RD model on the surface of a sphere. Asymptotic and numerical solutions are presented for the system governing the spot strengths, and we describe the complex bifurcation structure and demonstrate the occurrence of imperfection sensitivity due to higher order effects. Localized spot patterns can undergo a fast time instability and we derive the conditions for this phenomena, which depend on the spatial configuration of the spots and the parameters in the system. In the absence of these instabilities, our numerical solutions of the DAE system for $N = 2$ to $N = 8$ spots suggest a large basin of attraction to a small set of possible steady-state configurations. We discuss the connections between our results and the study of point vortices on the sphere, as well as the problem of determining a set of elliptic Fekete points, which correspond to globally minimizing the discrete logarithmic energy for $N$ points on the sphere.

preprint2014arXiv

Slowly varying control parameters, delayed bifurcations and the stability of spikes in reaction-diffusion systems

We present three examples of delayed bifurcations for spike solutions of reaction-diffusion systems. The delay effect results as the system passes slowly from a stable to an unstable regime, and was previously analysed in the context of ODE's in [P.Mandel, T.Erneux, J.Stat.Phys, 1987]. It was found that the instability would not be fully realized until the system had entered well into the unstable regime. The bifurcation is said to have been "delayed" relative to the threshold value computed directly from a linear stability analysis. In contrast, we analyze the delay effect in systems of PDE's. In particular, for spike solutions of singularly perturbed generalized Gierer-Meinhardt (GM) and Gray-Scott (GS) models, we analyze three examples of delay resulting from slow passage into regimes of oscillatory and competition instability. In the first example, for the GM model on the infinite real line, we analyze the delay resulting from slowly tuning a control parameter through a Hopf bifurcation. In the second example, we consider a Hopf bifurcation on a finite one-dimensional domain. In this scenario, as opposed to the extrinsic tuning of a system parameter through a bifurcation value, we analyze the delay of a bifurcation triggered by slow intrinsic dynamics of the PDE system. In the third example, we consider competition instabilities of the GS model triggered by the extrinsic tuning of a feed rate parameter. In all cases, we find that the system must pass well into the unstable regime before the onset of instability is fully observed, indicating delay. We also find that delay has an important effect on the eventual dynamics of the system in the unstable regime. We give analytic predictions for the magnitude of the delays as obtained through analysis of certain explicitly solvable nonlocal eigenvalue problems. The theory is confirmed by numerical solutions of the full PDE systems.

preprint2013arXiv

Logarithmic Expansions and the Stability of Periodic Patterns of Localized Spots for Reaction-Diffusion Systems in $\R^2

The linear stability of steady-state periodic patterns of localized spots in $\R^2$ for the two-component Gierer-Meinhardt (GM) and Schnakenburg reaction-diffusion models is analyzed in the semi-strong interaction limit corresponding to an asymptotically small diffusion coefficient $\eps^2$ of the activator concentration. In the limit $\eps\to 0$, localized spots in the activator are centered at the lattice points of a Bravais lattice with constant area $|Ω|$. To leading order in $ν={-1/\log\eps}$, the linearization of the steady-state periodic spot pattern has a zero eigenvalue when the inhibitor diffusivity satisfies $D={D_0/ν}$, for some $D_0$ independent of the lattice and the Bloch wavevector $\kb$. From a combination of the method of matched asymptotic expansions, Floquet-Bloch theory, and the rigorous study of certain nonlocal eigenvalue problems, an explicit analytical formula for the continuous band of spectrum that lies within an ${\mathcal O}(ν)$ neighborhood of the origin in the spectral plane is derived when $D={D_0/ν} + D_1$, where $D_1={\mathcal O}(1)$ is a de-tuning parameter. The periodic pattern is linearly stable when $D_1$ is chosen small enough so that this continuous band is in the stable left-half plane $\mbox{Re}(λ)<0$ for all $\kb$. Moreover, for both the Schnakenburg and GM models, our analysis identifies a model-dependent objective function, involving the regular part of the Bloch Green's function, that must be maximized in order to determine the specific periodic arrangement of localized spots that constitutes a linearly stable steady-state pattern for the largest value of $D$. From a numerical computation, it is shown within the class of oblique Bravais lattices that a regular hexagonal lattice arrangement of spots is optimal for maximizing the stability threshold in $D$.

preprint2010arXiv

The Stability and Dynamics of Localized Spot Patterns in the Two-Dimensional Gray-Scott Model

The dynamics and stability of multi-spot patterns to the Gray-Scott (GS) reaction-diffusion model in a two-dimensional domain is studied in the singularly perturbed limit of small diffusivity $ε$ of one of the two solution components. A hybrid asymptotic-numerical approach based on combining the method of matched asymptotic expansions with the detailed numerical study of certain eigenvalue problems is used to predict the dynamical behavior and instability mechanisms of multi-spot quasi-equilibrium patterns for the GS model in the limit $ε\to 0$. A differential algebraic ODE system for the collective coordinates $S_j$ and ${\mathbf x}_j$ for $j=1,...,k$ is derived, which characterizes the slow dynamics of a spot pattern. Instabilities of the multi-spot pattern due to the three distinct mechanisms of spot self-replication, spot oscillation, and spot annihilation, are studied by first deriving certain associated eigenvalue problems by using singular perturbation techniques. From a numerical computation of the spectrum of these eigenvalue problems, phase diagrams representing in the GS parameter space corresponding to the onset of spot instabilities are obtained for various simple spatial configurations of multi-spot patterns. In addition, it is shown that there is a wide parameter range where a spot instability can be triggered only as a result of the intrinsic slow motion of the collection of spots. The hybrid asymptotic-numerical results for spot dynamics and spot instabilities are validated from full numerical results computed from the GS model for various spatial configurations of spots.