Source author record

Theodore Kolokolnikov

Theodore Kolokolnikov 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

15works
16topics
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

15 published item(s)

preprint2022arXiv

A ring of spikes

For the Schnakenberg model, we consider a highly symmetric configuration of N spikes whose locations are located at the vertices of a regular N-gon inside either a unit disk or an annulus. We call such configuration a ring of spikes. The ring radius is characterized in terms of the modified Green's function. For a disk, we find that a ring of 9 or more spikes is always unstable with respect to small eigenvalues. Conversely, a ring of 8 or less spikes is stable inside a disk provided that the feed-rate $A$ is sufficiently large. More generally, for sufficiently high feed-rate, a ring of $N$ spikes can be stabilized provided that the annulus is thin enough. As $A$ is decreased, we show that the ring is destabilized due to small eigenvalues first, and then due to large eigenvalues, although both of these thresholds are separated by an asymptotically small amount. For a ring of 8 spikes inside a disk, the instability appears to be supercritical, and deforms the ring into a square-like configuration. For less than 8 spikes, this instability is subcritical and results in spike death.

preprint2022arXiv

Complex oscillatory motion of multiple spikes in a three-component Schnakenberg system

In this paper, we introduce a three-component Schnakenberg model. Its key feature is that it has a solution consisting of N spikes that undergoes a Hopf bifurcation with respect to N distinct modes nearly simultaneously. This results in complex oscillatory dynamics of the spikes, not seen in typical two-component models. For parameter values above the Hopf bifurcations, we derive reduced equations of motion which consist of coupled ordinary differential equations (ODEs) of order 2N for spike positions and their velocities. These ODEs fully describe the slow-time evolution of the spikes near the Hopf bifurcations. We then apply the method of multiple scales to the resulting ODEs to derive long-time dynamics. For a single spike, we find that its long-time motion consists of oscillations near the steady-state whose amplitude can be computed explicitly. For two spikes, the long-time behaviour can be either in-phase or out-of-phase oscillations. Both in and out of phase orbits are stable, coexist for the same parameter values, and the fate of orbit depends solely on the initial conditions. Further away from the Hopf bifurcation point, we offer numerical experiments indicating the existence of highly complex and chaotic oscillations.

preprint2022arXiv

It is better to be semi-regular when you have a low degree

We study the algebraic connectivity for several classes of random semi-regular graphs. For large random semi-regular bipartite graphs, we explicitly compute both their algebraic connectivity and as well as the full spectrum distribution. For an integer $d\in\left[ 3,7\right] $, we find families of random semi-regular graphs that have higher algebraic connectivity than a random $d$-regular graphs with the same number of vertices and edges. On the other hand, we show that regular graphs beat semi-regular graphs when $d\geq8.$ More generally, we study random semi-regular graphs whose average degree is $d$, not necessary an integer. This provides a natural generalization of a $d$-regular graph in the case of a non-integer $d.$ We characterise their algebraic connectivity in terms of a root of a certain 6th-degree polynomial. Finally, we construct a small-world-type network of average degree 2.5 with a relatively high algebraic connectivity. We also propose some related open problems and conjectures.

preprint2021arXiv

Laplace-Beltrami spectrum of ellipsoids that are close to spheres and analytic perturbation theory

We study the spectrum of the Laplace-Beltrami operator on ellipsoids. For ellipsoids that are close to the sphere, we use analytic perturbation theory to estimate the eigenvalues up to two orders. We show that for biaxial ellipsoids sufficiently close to the sphere, the first $N$ eigenvalues have multiplicity at most two, and characterize those that are simple. For the triaxial ellipsoids sufficiently close to the sphere that are not biaxial, we show that at least the first sixteen eigenvalues are all simple. We also give the results of various numerical experiments, including comparisons to our results from the analytic perturbation theory, and approximations for the eigenvalues of ellipsoids that degenerate into infinite cylinders or two-dimensional disks. We propose a conjecture on the exact number of nodal domains of near-sphere ellipsoids.

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.

preprint2016arXiv

Equilibria for an aggregation model with two species

We consider an aggregation model for two interacting species. The coupling between the species is via their velocities, that incorporate self- and cross-interactions. Our main interest is categorizing the possible steady states of the considered model. Notably, we identify their regions of existence and stability in the parameter space. For assessing the stability we use a combination of variational tools (based on the gradient flow formulation of the model and the associated energy), and linear stability analysis (perturbing the boundaries of the species' supports). We rely on numerical investigations for those steady states that are not analytically tractable. Finally we perform a two-scale expansion to characterize the steady state in the limit of asymptotically weak cross-interactions.

preprint2016arXiv

Metastable states for an aggregation model with noise

We study the long-time effect of noise on pattern formation for the aggregation model. We consider aggregation kernels that generate patterns consisting of two delta-concentrations. Without noise, there is a one-parameter family of admissible equilibria that consist of two concentrations whose mass is not necessary equal. We show that when a small amount of noise is added, the heavier concentration "leaks" its mass towards the lighter concentration over a very long time scale, eventually resulting in the equilibration of the two masses. We use exponentially small asymptotics to derive the long-time ODE's that quantify this mass exchange. Our theory is validated using full numerical simulations of the original model -- both of the original stochastic particle system and its PDE limit. Our formal computations show that adding noise destroys the degeneracy in the equilibrium solution and leads to a unique symmetric steady state.

preprint2015arXiv

Weakly Nonlinear Analysis of Vortex Formation in a Dissipative Variant of the Gross-Pitaevskii Equation

For a dissipative variant of the two-dimensional Gross-Pitaevskii equation with a parabolic trap under rotation, we study a symmetry breaking process that leads to the formation of vortices. The first symmetry breaking leads to the formation of many small vortices distributed uniformly near the Thomas-Fermi radius. The instability occurs as a result of a linear instability of a vortex-free steady state as the rotation is increased above a critical threshold. We focus on the second subsequent symmetry breaking, which occurs in the weakly nonlinear regime. At slightly above threshold, we derive a one-dimensional amplitude equation that describes the slow evolution of the envelope of the initial instability. We show that the mechanism responsible for initiating vortex formation is a modulational instability of the amplitude equation. We also illustrate the role of dissipation in the symmetry breaking process. All analyses are confirmed by detailed numerical computations.

preprint2014arXiv

A minimal model of predator-swarm interactions

We propose a minimal model of predator-swarm interactions which captures many of the essential dynamics observed in nature. Different outcomes are observed depending on the predator strength. For a "weak" predator, the swarm is able to escape the predator completely. As the strength is increased, the predator is able to catch up with the swarm as a whole, but the individual prey are able to escape by "confusing" the predator: the prey forms a ring with the predator at the center. For higher predator strength, complex chasing dynamics are observed which can become chaotic. For even higher strength, the predator is able to successfully capture the prey. Our model is simple enough to be amenable to a full mathematical analysis which is used to predict the shape of the swarm as well as the resulting predator-prey dynamics as a function of model parameters. We show that as the predator strength is increased, there is a transition (due to a Hopf bifurcation) from confusion state to chasing dynamics, and we compute the threshold analytically. Our analysis indicates that the swarming behaviour is not helpful in avoiding the predator, suggesting that there are other reasons why the species may swarm. The complex shape of the swarm in our model during the chasing dynamics is similar to the shape of a flock of sheep avoiding a shepherd.

preprint2014arXiv

Collective behaviour of large number of vortices in the plane

We investigate the dynamics of $N$ point vortices in the plane, in the limit of large $N$. We consider {\em relative equilibria}, which are rigidly rotating lattice-like configurations of vortices. These configurations were observed in several recent experiments [Durkin and Fajans, Phys. Fluids (2000) 12, 289-293; Grzybowski {\em et.al} PRE (2001)64, 011603]. We show that these solutions and their stability are fully characterized via a related {\em aggregation model} which was recently investigated in the context of biological swarms [Fetecau {\em et.al.}, Nonlinearity (2011) 2681; Bertozzi {\em et.al.}, M3AS (2011)]. By utilizing this connection, we give explicit analytic formulae for many of the configurations that have been observed experimentally. These include configurations of vortices of equal strength; the $N+1$ configurations of $N$ vortices of equal strength and one vortex of much higher strength; and more generally, $N+K$ configurations. We also give examples of configurations that have not been studied experimentally, including $N+2$ configurations where $N$ vortices aggregate inside an ellipse. Finally, we introduce an artificial ``damping'' to the vortex dynamics, in an attempt to explain the phenomenon of crystalization that is often observed in real experiments. The diffusion breaks the conservative structure of vortex dynamics so that any initial conditions converge to the lattice-like relative equilibrium.

preprint2014arXiv

Drunken robber, tipsy cop: First passage times, mobile traps, and Hopf bifurcations

For a random walk on a confined one-dimensional domain, we consider mean first passage times (MFPT) in the presence of a mobile trap. The question we address is whether a mobile trap can improve capture times over a stationary trap. We consider two scenarios: a randomly moving trap and an oscillating trap. In both cases, we find that a stationary trap actually performs better (in terms of reducing expected capture time) than a very slowly moving trap; however, a trap moving sufficiently fast performs better than a stationary trap. We explicitly compute the thresholds that separate the two regimes. In addition, we find a surprising relation between the oscillating trap problem and a moving-sink problem that describes reduced dynamics of a single spike in a certain regime of the Gray-Scott model. Namely, the above-mentioned threshold corresponds precisely to a Hopf bifurcation that induces oscillatory motion in the location of the spike. We use this correspondence to prove the uniqueness of the Hopf bifurcation.

preprint2014arXiv

Maximizing algebraic connectivity for certain families of graphs

We investigate the bounds on algebraic connectivity of graphs subject to constraints on the number of edges, vertices, and topology. We show that the algebraic connectivity for any tree on $n$ vertices and with maximum degree $d$ is bounded above by $2(d-2) \frac{1}{n}+O(\frac{\ln n}{n^{2}}) .$ We then investigate upper bounds on algebraic connectivity for cubic graphs. We show that algebraic connectivity of a cubic graph of girth $g$ is bounded above by $3-2^{3/2}\cos(π/\lfloor g/2\rfloor) ,$ which is an improvement over the bound found by Nilli [A. Nilli, Electron. J. Combin., 11(9), 2004]. Finally, we propose several conjectures and open questions.

preprint2014arXiv

Mean first passage time for a small rotating trap inside a reflective disk

We compute the mean first passage time (MFPT) for a Brownian particle inside a two-dimensional disk with reflective boundaries and a small interior trap that is rotating at a constant angular velocity. The inherent symmetry of the problem allows for a detailed analytic study of the situation. For a given angular velocity, we determine the optimal radius of rotation that minimizes the average MFPT over the disk. Several distinct regimes are observed, depending on the ratio between the angular velocity $ω$ and the trap size $\varepsilon$, and several intricate transitions are analyzed using the tools of asymptotic analysis and Fourier series. For $ω\sim \mathcal{O}(1)$, we compute a critical value $ω_c>0$ such that the optimal trap location is at the origin whenever $ω< ω_c$, and is off the origin for $ω>ω_c $. In the regime $1 \ll ω\ll \mathcal{O}(\varepsilon^{-1})$ the optimal trap path approaches the boundary of the disk. However as $ω$ is further increased to $\mathcal{O}(\varepsilon^{-1})$, the optimal trap path "jumps" closer to the origin. Finally for $ω\gg \mathcal{O}(\varepsilon^{-1})$ the optimal trap path subdivides the disk into two regions of equal area. This simple geometry provides a good test case for future studies of MFPT with more complex trap motion.

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.

preprint2012arXiv

The Stability of Steady-State Hot-Spot Patterns for a Reaction-Diffusion Model of Urban Crime

The existence and stability of localized patterns of criminal activity are studied for the reaction-diffusion model of urban crime that was introduced by Short et. al. [Math. Models. Meth. Appl. Sci., 18, Suppl. (2008), pp. 1249--1267]. Such patterns, characterized by the concentration of criminal activity in localized spatial regions, are referred to as hot-spot patterns and they occur in a parameter regime far from the Turing point associated with the bifurcation of spatially uniform solutions. Singular perturbation techniques are used to construct steady-state hot-spot patterns in one and two-dimensional spatial domains, and new types of nonlocal eigenvalue problems are derived that determine the stability of these hot-spot patterns to ${\mathcal O}(1)$ time-scale instabilities. From an analysis of these nonlocal eigenvalue problems, a critical threshold $K_c$ is determined such that a pattern consisting of $K$ hot-spots is unstable to a competition instability if $K>K_c$. This instability, due to a positive real eigenvalue, triggers the collapse of some of the hot-spots in the pattern. Furthermore, in contrast to the well-known stability results for spike patterns of the Gierer-Meinhardt reaction-diffusion model, it is shown for the crime model that there is only a relatively narrow parameter range where oscillatory instabilities in the hot-spot amplitudes occur. Such an instability, due to a Hopf bifurcation, is studied explicitly for a single hot-spot in the shadow system limit, for which the diffusivity of criminals is asymptotically large. Finally, the parameter regime where localized hot-spots occur is compared with the parameter regime, studied in previous works, where Turing instabilities from a spatially uniform steady-state occur.