Margaret Beck

Snakes, Ladders, and Isolas of Localized Patterns

Stable localized roll structures have been observed in many physical problems and model equations, notably in the one-dimensional (1D) Swift–Hohenberg equation. Reflection-symmetric localized rolls are often found to lie on two “snaking” solution branches so that the spatial width of the localized rolls increases when moving along each branch. Recent numerical results by Burke and Knobloch indicate that the two branches are connected by infinitely many “ladder” branches of asymmetric localized rolls. In this paper, these phenomena are investigated analytically. It is shown that both snaking of symmetric pulses and the ladder structure of asymmetric states can be predicted completely from the bifurcation structure of fronts that connect the trivial state to rolls. It is also shown that isolas of asymmetric states may exist, and it is argued that the results presented here apply to 2D stationary states that are localized in one spatial direction.

Using Global Invariant Manifolds to Understand Metastability in the Burgers Equation With Small Viscosity

The large-time behavior of solutions to the Burgers equation with small viscosity is described using invariant manifolds. In particular, a geometric explanation is provided for a phenomenon known as metastability, which in the present context means that solutions spend a very long time near the family of solutions known as diffusive N-waves before finally converging to a stable self-similar diffusion wave. More precisely, it is shown that in terms of similarity, or scaling, variables in an algebraically weighted

Nonlinear Stability of Time-Periodic Viscous Shocks

L^2

Electrical Waves in a One-Dimensional Model of Cardiac Tissue

space, the self-similar diffusion waves correspond to a one-dimensional global center manifold of stationary solutions. Through each of these fixed points there exists a one-dimensional, global, attractive, invariant manifold corresponding to the diffusive N-waves. Thus, metastability corresponds to a fast transient in which solutions approach this “metastable” manifold of diffusive N-waves, followed by a slow decay along this manifold, and, finally, convergence to the self-similar diffusion...

Computing the Maslov index for large systems

In order to understand the nonlinear stability of many types of time-periodic travelling waves on unbounded domains, one must overcome two main difficulties: the presence of embedded neutral eigenvalues and the time-dependence of the associated linear operator. This problem is studied in the context of time-periodic Lax shocks in systems of viscous conservation laws. Using spatial dynamics and a decomposition into separate Floquet eigenmodes, it is shown that the linear evolution for the time-dependent operator can be represented using a contour integral similar to that of the standard time-independent case. By decomposing the resulting Green’s distribution, the leading order behavior associated with the embedded eigenvalues is extracted. Sharp pointwise bounds are then obtained, which are used to prove that time-periodic Lax shocks are linearly and nonlinearly stable under the necessary conditions of spectral stability and minimal multiplicity of the translational eigenvalues. The latter conditions hold, for example, for small-oscillation time-periodic waves that emerge through a supercritical Hopf bifurcation from a family of time-independent Lax shocks of possibly large amplitude.

Using Global Invariant Manifolds to Understand Metastability in the Burgers Equation with Small Viscosity

The electrical dynamics in the heart is modeled by a two-component PDE. Using geometric singular perturbation theory, it is shown that a traveling pulse solution, which corresponds to a single heartbeat, exists. One key aspect of the proof involves tracking the solution near a point on the slow manifold that is not normally hyperbolic. This is achieved by desingularizing the vector eld using a blow-up technique. This feature is relevant because it distinguishes cardiac impulses from, for example, nerve impulses. Stability of the pulse is also shown, by computing the zeros of the Evans function. Although the spectrum of one of the fast components is only marginally stable, due to essential spectrum that accumulates at the origin, it is shown that the spectrum of the full pulse consists of an isolated eigenvalue at zero and essential spectrum that is bounded away from the imaginary axis. Thus, this model provides an example in a biological application reminiscent of a previously observed mathematical phenomenon: that connecting an unstable { in this case marginally stable { front and back can produce a stable pulse. Finally, remarks are made regarding the existence and stability of spatially periodic pulses, corresponding to successive heartbeats, and their relationship with alternans, irregular action potentials that have been linked with arrhythmia.

Superadiabaticity in reaction waves as a mechanism for energy concentration

We address the problem of computing the Maslov index for large linear symplectic systems on the real line. The Maslov index measures the signed intersections (with a given reference plane) of a path of Lagrangian planes. The natural chart parameterization for the Grassmannian of La- grangian planes is the space of real symmetric matrices. Linear system evolu- tion induces a Riccati evolution in the chart. For large order systems this is a practical approach as the computational complexity is quadratic in the order. The Riccati solutions, however, also exhibit singularites (which are traversed by changing charts). Our new results involve characterizing these Riccati sin- gularities and two trace formulae for the Maslov index as follows. First, we show that the number of singular eigenvalues of the symmetric chart represen- tation equals the dimension of intersection with the reference plane. Second, the Cayley map is a diffeomorphism from the space of real symmetric matrices to the manifold of unitary symmetric matrices. We show the logarithm of the Cayley map equals the arctan map (modulo 2i) and its trace measures the angle of the Langrangian plane to the reference plane. Third, the Riccati flow under the Cayley map induces a flow in the manifold of unitary symmetric matrices. Using the natural unitary action on this manifold, we pullback the flow to the unitary Lie algebra and monitor its trace. This avoids singularities, and is a natural robust procedure. We demonstrate the effectiveness of these approaches by applying them to a large eigenvalue problem. We also discuss the extension of the Maslov index to the infinite dimensional case.

A Geometric Construction of Traveling Waves in a Bioremediation Model

The large-time behavior of solutions to the Burgers equation with small viscosity is described using invariant manifolds. In particular, a geometric explanation is provided for a phenomenon known as metastability, which in the present context means that solutions spend a very long time near the family of solutions known as diffusive N-waves before finally converging to a stable self-similar diffusion wave. More precisely, it is shown that in terms of similarity, or scaling, variables in an algebraically weighted

Nonlinear Stability of Semidiscrete Shocks for Two-Sided Schemes

L^2

A geometric theory of chaotic phase synchronization

space, the self-similar diffusion waves correspond to a one-dimensional global center manifold of stationary solutions. Through each of these fixed points there exists a one-dimensional, global, attractive, invariant manifold corresponding to the diffusive N-waves. Thus, metastability corresponds to a fast transient in which solutions approach this “metastable” manifold of diffusive N-waves, followed by a slow decay along this manifold, and, finally, convergence to the self-similar diffusion wave.