Search

In the work Cho et al. [Jpn. J. Ind. Appl. Math. 33 (2016): 145-166] the authors conjecture that the quadratic nonlinear Schr\"odinger equation (NLS) $i u_t = u_{xx} + u^2 $ for $ x \in \mathbb{T}$ is globally well-posed for real initial data. We identify initial data whose numerical solution blows up in contradiction of this conjecture. The solution exhibits self-similar blowup and potentially nontrivial self-similar dynamics, however the proper scaling ansatz remains elusive. Furthermore, the set of real initial data which blows up under the NLS dynamics appears to occur on a codimension-1 manifold, and we conjecture that it is precisely the stable manifold of the zero equilibrium for the nonlinear heat equation $u_t = u_{xx} + u^2 $. We apply the parameterization method to study the internal dynamics of this manifold, offering a heuristic argument in support of our conjecture., Comment: 21 pages, 7 figures

In the scalar Swift-Hohenberg equation with quadratic-cubic nonlinearity, it is known that symmetric pulse solutions exist for certain parameter regions. In this paper we develop a method to determine the spectral stability of these solutions by associating a Maslov index to them. This requires extending the method of computing the Maslov index introduced by Robbin and Salamon [Topology 32, no.4 (1993): 827-844] to so-called degenerate crossings. We extend their formulation of the Maslov index to degenerate crossings of general order in the case where the intersection is fully degenerate, meaning that if the dimension of the intersection is k, then each of the k crossings is a degenerate one. We then argue that, in this case, this index coincides with the number of unstable eigenvalues for the linearized evolution equation. Furthermore, we develop a numerical method to compute the Maslov index associated to symmetric pulse solutions. Finally, we consider several solutions to the Swift-Hohenberg equation and use our method to characterize their stability.

A variety of nonlinear models of biological systems generate complex chaotic behaviors that contrast with biological homeostasis, the observation that many biological systems prove remarkably robust in the face of changing external or internal conditions. Motivated by the subtle dynamics of cell activity in a crustacean central pattern generator (CPG), this paper proposes a refinement of the notion of chaos that reconciles homeostasis and chaos in systems with multiple timescales. We show that systems displaying relaxation cycles while going through chaotic attractors generate chaotic dynamics that are regular at macroscopic timescales and are thus consistent with physiological function. We further show that this relative regularity may break down through global bifurcations of chaotic attractors such as crises, beyond which the system may also generate erratic activity at slow timescales. We analyze these phenomena in detail in the chaotic Rulkov map, a classical neuron model known to exhibit a variety of chaotic spike patterns. This leads us to propose that the passage of slow relaxation cycles through a chaotic attractor crisis is a robust, general mechanism for the transition between such dynamics. We validate this numerically in three other models: a simple model of the crustacean CPG neural network, a discrete cubic map, and a continuous flow.

In this paper, we consider the dynamics of solutions to complex-valued evolutionary partial differential equations (PDEs) and show existence of heteroclinic orbits from nontrivial equilibria to zero via computer-assisted proofs. We also show that the existence of unbounded solutions along unstable manifolds at the equilibrium follows from the existence of heteroclinic orbits. Our computer-assisted proof consists of three separate techniques of rigorous numerics: an enclosure of a local unstable manifold at the equilibria, a rigorous integration of PDEs, and a constructive validation of a trapping region around the zero equilibrium., Comment: 17 pages, 7 figures

In this paper, we study the dynamics of a class of nonlinear Schr\"odinger equation $ i u_t = \triangle u + u^p $ for $ x \in \mathbb{T}^d$. We prove that the PDE is integrable on the space of non-negative Fourier coefficients, in particular that each Fourier coefficient of a solution can be explicitly solved by quadrature. Within this subspace we demonstrate a large class of (quasi)periodic solutions all with the same frequency, as well as solutions which blowup in finite time in the $L^2$ norm., Comment: 24 pages, 1 figure

Classical results from Sturm-Liouville theory state that the number of unstable eigenvalues of a scalar, second-order linear operator is equal to the number of associated conjugate points. Recent work has extended these results to a much more general setting, thus allowing for spectral stability of nonlinear waves in a variety of contexts to be determined by counting conjugate points. However, in practice, it is not yet clear whether it is easier to compute conjugate points than to just directly count unstable eigenvalues. We address this issue by developing a framework for the computation of conjugate points using validated numerics. Moreover, we apply our method to a parameter-dependent system of bistable equations and show that there exist both stable and unstable standing fronts. This application can be seen as complimentary to the classical result via Sturm-Louiville theory that in scalar reaction-diffusion equations pulses are unstable whereas fronts are stable, and to the more recent result of "Instability of pulses in gradient reaction-diffusion systems: a symplectic approach," by Beck et. al., that symmetric pulses in reaction-diffusion systems with gradient nonlinearity are also necessarily unstable., Comment: 45 pages; 2 figures; fixed typo

In this paper, we study the global dynamics of a class of nonlinear Schr\"odinger equations using perturbative and non-perturbative methods. We prove the semi-global existence of solutions for initial conditions close to constant. That is, solutions will exist for all positive time or all negative time. The existence of an open set of initial data which limits to zero in both forward and backward time is also demonstrated. This result in turn forces the non-existence of any real-analytic conserved quantities. For the quadratic case, we prove the existence of two (infinite) families of nontrivial unstable equilibria and prove the existence of heteroclinic orbits limiting to the nontrivial equilibria in backward time and to zero in forward time. By a time reversal argument, we also obtain heteroclinic orbits limiting to the nontrivial equilibria in forward time and to zero in backward time. The proofs for the quadratic equation are computer-assisted and rely on three separate ingredients: an enclosure of a local unstable manifold at the equilibria, a rigorous integration of the flow (starting from the unstable manifold) and a proof that the solution enters a validated stable set (hence showing convergence to zero)., Comment: 51 pages, 8 figures

This paper develops validated computational methods for studying infinite dimensional stable manifolds at equilibrium solutions of parabolic PDEs, synthesizing disparate errors resulting from numerical approximation. To construct our approximation, we decompose the stable manifold into three components: a finite dimensional slow component, a fast-but-finite dimensional component, and a strongly contracting infinite dimensional "tail". We employ the parameterization method in a finite dimensional projection to approximate the slow-stable manifold, as well as the attached finite dimensional invariant vector bundles. This approximation provides a change of coordinates which largely removes the nonlinear terms in the slow stable directions. In this adapted coordinate system we apply the Lyapunov-Perron method, resulting in mathematically rigorous bounds on the approximation errors. As a result, we obtain significantly sharper bounds than would be obtained using only the linear approximation given by the eigendirections. As a concrete example we illustrate the technique for a 1D Swift-Hohenberg equation., Comment: 57 pages, 3 figures

In this paper, we introduce a method for computing rigorous local inclusions of solutions of Cauchy problems for nonlinear heat equations for complex time values. Using a solution map operator, we construct a simplified Newton operator and show that it has a unique fixed point. The fixed point together with its rigorous bounds provides the local inclusion of the solution of the Cauchy problem. The local inclusion technique is then applied iteratively to compute solutions over long time intervals. This technique is used to prove the existence of a branching singularity in the nonlinear heat equation. Finally, we introduce an approach based on the Lyapunov-Perron method to calculate part of a center-stable manifold and prove that an open set of solutions of the Cauchy problem converge to zero, hence yielding the global existence of the solutions in the complex plane of time., Comment: 42 pages, 5 figures

We propose that the recently defined persistent homology dimensions are a practical tool for fractal dimension estimation of point samples. We implement an algorithm to estimate the persistent homology dimension, and compare its performance to classical methods to compute the correlation and box-counting dimensions in examples of self-similar fractals, chaotic attractors, and an empirical dataset. The performance of the $0$-dimensional persistent homology dimension is comparable to that of the correlation dimension, and better than box-counting., Comment: 24 pages, 13 figures

In this paper, we prove that Wright's equation $y'(t) = - \alpha y(t-1) \{1 + y(t)\}$ has a unique slowly oscillating periodic solution for parameter values $\alpha \in (\tfrac{\pi}{2}, 1.9]$, up to time translation. This result proves Jones' Conjecture formulated in 1962, that there is a unique slowly oscillating periodic orbit for all $ \alpha > \tfrac{\pi}{2}$. Furthermore, there are no isolas of periodic solutions to Wright's equation; all periodic orbits arise from Hopf bifurcations., Comment: 41 pages, 5 figures, 1 table

In this paper, we prove that Wright's equation $y'(t) = - \alpha y(t-1) \{1 + y(t)\}$ has a unique slowly oscillating periodic solution (SOPS) for all parameter values $\alpha \in [ 1.9,6.0]$, up to time translation. Our proof is based on a same strategy employed earlier by Xie [27]; show that every SOPS is asymptotically stable. We first introduce a branch and bound algorithm to control all SOPS using bounding functions at all parameter values $\alpha \in [ 1.9,6.0]$. Once the bounding functions are constructed, we then control the Floquet multipliers of all possible SOPS by solving rigorously an eigenvalue problem, again using a formulation introduced by Xie. Using these two main steps, we prove that all SOPS of Wright's equation are asymptotically stable for $\alpha \in [ 1.9,6.0]$, and the proof follows. This result is a step toward the proof of the Jones' Conjecture formulated in 1962., Comment: 23 pages, 2 figures

Wright's conjecture states that the origin is the global attractor for the delay differential equation $y'(t) = - \alpha y(t-1) [ 1 + y(t) ] $ for all $\alpha \in (0,\tfrac{\pi}{2}]$. This has been proven to be true for a subset of parameter values $\alpha$. We extend the result to the full parameter range $\alpha \in (0,\tfrac{\pi}{2}]$, and thus prove Wright's conjecture to be true. Our approach relies on a careful investigation of the neighborhood of the Hopf bifurcation occurring at $\alpha =\tfrac{\pi}{2}$. This analysis fills the gap left by complementary work on Wright's conjecture, which covers parameter values further away from the bifurcation point. Furthermore, we show that the branch of (slowly oscillating) periodic orbits originating from this Hopf bifurcation does not have any subsequent bifurcations (and in particular no folds) for $\alpha\in(\tfrac{\pi}{2} , \tfrac{\pi}{2} + 6.830 \times 10^{-3}]$. When combined with other results, this proves that the branch of slowly oscillating solutions that originates from the Hopf bifurcation at $\alpha=\tfrac{\pi}{2}$ is globally parametrized by $\alpha > \tfrac{\pi}{2}$., Comment: 45 pages

Biological and physical systems often exhibit distinct structures at different spatial/temporal scales. Persistent homology is an algebraic tool that provides a mathematical framework for analyzing the multi-scale structures frequently observed in nature. In this paper a theoretical framework for the algorithmic computation of an arbitrarily good approximation of the persistent homology is developed. We study the filtrations generated by sub-level sets of a function $f : X \to \mathbb{R}$, where $X$ is a CW-complex. In the special case $X = [0,1]^N$, $N \in \mathbb{N}$ we discuss implementation of the proposed algorithms. We also investigate a priori and a posteriori bounds of the approximation error introduced by our method., Comment: 26 pages; changed title; added revisions