Study of invariant manifolds in two different problems : the Hopf-zero singularity and neural synchrony

Author

Castejón i Compnay, Oriol

Director

Baldomá Barraca, Inmaculada

Codirector

Guillamon i Grabolosa, Antoni

Martínez-Seara i Alonso, M. Teresa (Maria Teresa)

Date of defense

2015-07-16

Pages

348 p.



Department/Institute

Universitat Politècnica de Catalunya. Departament de Matemàtica Aplicada I

Abstract

The main object of study of this thesis are invariant manifolds in the field of dynamical systems. We deal with two different and independent topics, namely, the study of exponentially small splitting of invariant manifolds in analytic unfoldings of the Hopf-zero singularity (in Part I) and the applications of dynamical systems in problems inspired by neuroscience (in Part II). In general, this thesis studies both theoretical and applied problems in dynamical systems, using analytical as well as computational tools. In Part I, we consider a certain class of generic unfoldings of the so-called Hopf-zero singularity. One can see that the truncation of the normal form at any finite order of such unfoldings possesses two saddle-focus critical points and, when the parameters lie on a certain curve, they are connected by a one- and a two-dimensional heteroclinic manifolds. However, considering the whole vector field, one expects these heteroclinic connections to be destroyed. This fact can lead to the birth of a homoclinic connection to one of the critical points, producing thus a Shilnikov bifurcation. For the case of $C^\infty$ unfoldings, this was proved by Broer and Vegter during the 80's, but for analytic unfoldings it has remained an open problem. Recently, under some assumptions on the size of the splitting of the heteroclinic connections, Dumortier, Ibáñez, Kokubu and Simó proved the existence of Shilnikov bifurcations in the analytic case. Our study concerns the splitting of the one- and two-dimensional heteroclinic connections. These cannot be detected in the truncation of the normal form at any order, and hence they are exponentially small with respect to one of the perturbation parameters. We give asymptotic formulas of these splittings and, in particular, we prove that under generic conditions the main assumptions made by Dumortier, Ibáñez, Kokubu and Simó hold. In Part II, we deal with tools to provide an accurate prediction of phase variations in an oscillator subject to external stimuli. We construct a method based on the concepts of isochrons, Phase Response Functions (PRF) and Amplitude Response Functions (ARF). In particular, the method can be applied to neurons in a state of repetitive firing. In the special case of a pulse-train periodic stimulus, the application of this theoretical frame leads to a 2D map, one variable controlling phase jumps and the other controlling amplitude jumps. We compare these maps to the classical 1D maps obtained via Phase Response Curves (PRC) and we identify circumstances in which the 2D maps give a more accurate prediction of the synchronization. Moreover, we implement some numerical methods to compute the invariant curves of the 2D maps as well as the dynamics inside these curves. Finally, we compute Arnold tongues of these maps, which allow to determine regions in the parameter space for which the neuron is synchronized to the external input.


El principal objecte d'estudi d'aquesta tesi són les varietats invariants en el camp dels sistemes dinàmics. Considerem dos temes diferents i independents, concretament l'estudi de l'escissió exponencialment petita de varietats invariants en desplegaments analí­tics de la singularitat Hopf-zero (a la Part I), i les aplicacions dels sistemes dinàmics en problemes inspirats per la neurociència (a la Part II). A la Part I, considerem una classe de desplegaments genèrics de l'anomenada singularitat Hopf-zero. Es pot veure que el truncament de la forma normal a qualsevol ordre finit d'aquests desplegaments té dos punts crítics de tipus sella-focus i, quan els paràmetres estan sobre una certa corba, estan connectats per dues varietats heteroclí­niques, una d'unidimensional i una de bidimensional. No obstant, si es considera tot el camp vectorial, s'espera que aquestes connexions heteroclí­niques desapareguin. Això pot causar el naixement d'una òrbita homoclí­nica en un dels dos punts crí­tics, produint així­ el que es coneix com una bifurcació de Shilnikov. En el cas de desplegaments $C^\infty$, això va ser provat per Broer i Vegter durant els anys 80, però el cas de desplegaments analí­tics ha quedat obert. Recentment, sota certes hipòtesis sobre la mida de l'escissió de les connexions heteroclí­niques, Dumortier, Ibáñez, Kokubu i Simó han provat l'existència de bifurcacions de Shilnikov en el cas analític. El nostre estudi consisteix en el càlcul del trencament de les connexions heteroclí­niques. Aquests trencaments no es poden detectar en la forma normal a cap ordre i, per tant, són exponencialment petits en un dels paràmetres de pertorbació. Donem fórmules asimptòtiques d'aquests trencaments i, en particular, provem que sota certes condicions genèriques les principals hipòtesis fetes per Dumortier, Ibàñez, Kokubu i Simó són vàlides. A la Part II, considerem eines per proporcionar una predicció acurada de la variació de fase en un oscil·lador subjecte a estímuls externs. Construïm un mètode basat en els conceptes d'isòcrones, Funcions de Resposta de Fase (PRF, per les seves inicials en anglès) i Funcions de Resposta d'Amplitud (ARF). En particular, el mètode es pot aplicar a neurones en un estat de dispar repetitiu. En el cas especial d'un tren de pulsos periòdic, l'aplicació d'aquest mètode teòric dóna lloc a una aplicació 2D, on una variable controla els canvis de la fase i l'altra els canvis en l'amplitud. Comparem aquestes aplicacions amb les aplicacions 1D clàssiques obtingudes a través de la Corbes de Resposta de Fase (PRC) i identifiquem circumstàncies en què les aplicacions 2D donen una millora substancial de la predicció de sincronització. A més, implementem alguns mètodes numèrics per calcular les corbes invariants de les aplicacions 2D així­ com la dinàmica dins aquestes corbes. Finalment, calculem les llengües d'Arnold corresponents a aquestes aplicacions, que permeten determinar regions en l'espai de paràmetres per a les quals la neurona es sincronitza amb l'estí­mul extern.

Subjects

517 - Analysis

Documents

TOCiC1de1.pdf

4.282Mb

 

Rights

L'accés als continguts d'aquesta tesi queda condicionat a l'acceptació de les condicions d'ús establertes per la següent llicència Creative Commons: http://creativecommons.org/licenses/by-nc-nd/3.0/es/
L'accés als continguts d'aquesta tesi queda condicionat a l'acceptació de les condicions d'ús establertes per la següent llicència Creative Commons: http://creativecommons.org/licenses/by-nc-nd/3.0/es/

This item appears in the following Collection(s)