Source author record

Serafim Rodrigues

Serafim Rodrigues 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

12works
12topics
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

12 published item(s)

preprint2021arXiv

Spike-adding and reset-induced canard cycles in adaptive integrate and fire models

We study a class of planar integrate and fire (IF) models called adaptive integrate and fire (AIF) models, which possesses an adaptation variable on top of membrane potential, and whose subthreshold dynamics is piecewise linear (PWL). These AIF models therefore have two reset conditions, which enable bursting dynamics to emerge for suitable parameter values. Such models can be thought of as hybrid dynamical systems. We consider a particular slow dynamics within AIF models and prove the existence of bursting cycles with $N$ resets, for any integer $N$. Furthermore, we study the transition between $N$- and $(N+1)$-reset cycles upon vanishingly small parameter variations and prove (for $N=2$) that such transitions are organised by canard cycles. Finally, using numerical continuation we compute branches of bursting cycles, including canard-explosive branches, in these AIF models, by suitably recasting the periodic problem as a two-point boundary-value problem.

preprint2020arXiv

Towards a new classification of bursting patterns: review & extensions

The mathematical classification of complex bursting oscillations in multiscale excitable systems, seen for example in physics and neuroscience, has been the subject of active enquiry since the early 1980s. This classification problem is fundamental as it also establishes analytical and numerical foundations for studying complex temporal behaviours in multiple timescale models. This manuscript begins by reviewing the seminal works of Rinzel and Izhikevich in classifying bursting patterns of excitable cell models. Moreover, we recall an alternative, yet complementary, mathematical classification approach by Golubitsky, which together with the Rinzel-Izhikevich proposals provide the state-of-the-art foundations to the classification problem. Unexpectedly, while keeping within the Rinzel-Izhikevich framework, we find novel cases of bursting mechanisms not considered before. Moving beyond the state-of-the-art, we identify novel bursting mechanisms that fall outside the current classifications. This leads us towards a new classification, which requires the analysis of both the fast and the slow subsystems of an underlying slow-fast model. This new classification allows the dynamical dissection of a larger class of bursters. To substantiate this, we add a new class of bursters with at least two slow variables, which we denote folded-node bursters, to convey the idea that the bursts are initiated or annihilated via a folded-node singularity. In fact, there are two main families of folded-node bursters, depending upon the phase of the bursting cycle during which folded-node dynamics occurs. If it occurs during the silent phase, we obtain the classical folded-node bursting (or simply folded-node bursting). If it occurs during the active phase, we have cyclic folded-node bursting. We classify both families and give examples of minimal systems displaying these novel types of bursting behaviour.

preprint2020arXiv

Why we should use Topological Data Analysis in Ageing: towards defining the "Topological shape of ageing"

Living systems are subject to the arrow of time; from birth, they undergo complex transformations (self-organization) in a constant battle for survival, but inevitably ageing and disease trap them to death. Can ageing be understood and eventually reversed? What tools can be employed to further our understanding of ageing? The present article is an invitation for biologists and clinicians to consider key conceptual ideas and computational tools (known to mathematicians and physicists), which potentially may help dissect some of the underlying processes of ageing and disease. Specifically, we first discuss how to classify and analyze complex systems, as well as highlight critical theoretical difficulties that make complex systems hard to study. Subsequently, we introduce Topological Data Analysis - a novel Big Data tool - which may help in the study of complex systems since it extracts knowledge from data in a holistic approach via topological considerations. These conceptual ideas and tools are discussed in a rather informal way to pave future discussions and collaborations between mathematicians and biologists studying ageing.

preprint2016arXiv

A modular architecture for transparent computation in Recurrent Neural Networks

Computation is classically studied in terms of automata, formal languages and algorithms; yet, the relation between neural dynamics and symbolic representations and operations is still unclear in traditional eliminative connectionism. Therefore, we suggest a unique perspective on this central issue, to which we would like to refer as to transparent connectionism, by proposing accounts of how symbolic computation can be implemented in neural substrates. In this study we first introduce a new model of dynamics on a symbolic space, the versatile shift, showing that it supports the real-time simulation of a range of automata. We then show that the Goedelization of versatile shifts defines nonlinear dynamical automata, dynamical systems evolving on a vectorial space. Finally, we present a mapping between nonlinear dynamical automata and recurrent artificial neural networks. The mapping defines an architecture characterized by its granular modularity, where data, symbolic operations and their control are not only distinguishable in activation space, but also spatially localizable in the network itself, while maintaining a distributed encoding of symbolic representations. The resulting networks simulate automata in real-time and are programmed directly, in absence of network training. To discuss the unique characteristics of the architecture and their consequences, we present two examples: i) the design of a Central Pattern Generator from a finite-state locomotive controller, and ii) the creation of a network simulating a system of interactive automata that supports the parsing of garden-path sentences as investigated in psycholinguistics experiments.

preprint2015arXiv

Turing Computation with Recurrent Artificial Neural Networks

We improve the results by Siegelmann & Sontag (1995) by providing a novel and parsimonious constructive mapping between Turing Machines and Recurrent Artificial Neural Networks, based on recent developments of Nonlinear Dynamical Automata. The architecture of the resulting R-ANNs is simple and elegant, stemming from its transparent relation with the underlying NDAs. These characteristics yield promise for developments in machine learning methods and symbolic computation with continuous time dynamical systems. A framework is provided to directly program the R-ANNs from Turing Machine descriptions, in absence of network training. At the same time, the network can potentially be trained to perform algorithmic tasks, with exciting possibilities in the integration of approaches akin to Google DeepMind's Neural Turing Machines.

preprint2013arXiv

A method for the reconstruction of unknown non-monotonic growth functions in the chemostat

We propose an adaptive control law that allows one to identify unstable steady states of the open-loop system in the single-species chemostat model without the knowledge of the growth function. We then show how one can use this control law to trace out (reconstruct) the whole graph of the growth function. The process of tracing out the graph can be performed either continuously or step-wise. We present and compare both approaches. Even in the case of two species in competition, which is not directly accessible with our approach due to lack of controllability, feedback control improves identifiability of the non-dominant growth rate.

preprint2013arXiv

Extremum seeking via continuation techniques for optimizing biogas production in the chemostat

We consider the chemostat model with the substrate concentration as the single measurement. We propose a control strategy that drives the system at a steady state maximizing the gas production without the knowledge of the specific growth rate. Our approach separates the extremum seeking problem from the feedback control problem such that each of the two subproblems can be solved with relatively simple algorithms. We are then free to choose any numerical optimization algorithm. We give a demonstration for two choices: one is based on slow-fast dynamics and numerical continuation, the other is a combination of golden-section and Newton iteration. The method copes with non-monotonic growth functions.

preprint2013arXiv

On the electrodynamics of neural networks

We present a microscopic approach for the coupling of cortical activity, as resulting from proper dipole currents of pyramidal neurons, to the electromagnetic field in extracellular fluid in presence of diffusion and Ohmic conduction. Starting from a full-fledged three-compartment model of a single pyramidal neuron, including shunting and dendritic propagation, we derive an observation model for dendritic dipole currents in extracellular space and thereby for the dendritic field potential that contributes to the local field potential of a neural population. Under reasonable simplifications, we then derive a leaky integrate-and-fire model for the dynamics of a neural network, which facilitates comparison with existing neural network and observation models. In particular, we compare our results with a related model by means of numerical simulations. Performing a continuum limit, neural activity becomes represented by a neural field equation, while an observation model for electric field potentials is obtained from the interaction of cortical dipole currents with charge density in non-resistive extracellular space as described by the Nernst-Planck equation. Our work consistently satisfies the widespread dipole assumption discussed in the neuroscientific literature.

preprint2013arXiv

Short-term synaptic plasticity in the deterministic Tsodyks-Markram model leads to unpredictable network dynamics

Short-Term Synaptic Plasticity (STSP) strongly affects the neural dynamics of cortical networks. The Tsodyks and Markram (TM) model for STSP accurately accounts for a wide range of physiological responses at different types of cortical synapses. Here, we report for the first time a route to chaotic behavior via a Shilnikov homoclinic bifurcation that dynamically organises some of the responses in the TM model. In particular, the presence of such a homoclinic bifurcation strongly affects the shape of the trajectories in the phase space and induces highly irregular transient dynamics; indeed, in the vicinity of the Shilnikov homoclinic bifurcation, the number of population spikes and their precise timing are unpredictable and highly sensitive to the initial conditions. Such an irregular deterministic dynamics has its counterpart in stochastic/network versions of the TM model: the existence of the Shilnikov homoclinic bifurcation generates complex and irregular spiking patterns and --acting as a sort of springboard-- facilitates transitions between the down-state and unstable periodic orbits. The interplay between the (deterministic) homoclinic bifurcation and stochastic effects may give rise to some of the complex dynamics observed in neural systems.

preprint2012arXiv

A biophysical observation model for field potentials of networks of leaky integrate-and-fire neurons

We present a biophysical approach for the coupling of neural network activity as resulting from proper dipole currents of cortical pyramidal neurons to the electric field in extracellular fluid. Starting from a reduced threecompartment model of a single pyramidal neuron, we derive an observation model for dendritic dipole currents in extracellular space and thereby for the dendritic field potential that contributes to the local field potential of a neural population. This work aligns and satisfies the widespread dipole assumption that is motivated by the "open-field" configuration of the dendritic field potential around cortical pyramidal cells. Our reduced three-compartment scheme allows to derive networks of leaky integrate-and-fire models, which facilitates comparison with existing neural network and observation models. In particular, by means of numerical simulations we compare our approach with an ad hoc model by Mazzoni et al. [Mazzoni, A., S. Panzeri, N. K. Logothetis, and N. Brunel (2008). Encoding of naturalistic stimuli by local field potential spectra in networks of excitatory and inhibitory neurons. PLoS Computational Biology 4 (12), e1000239], and conclude that our biophysically motivated approach yields substantial improvement.

preprint2012arXiv

Coupling continuous neural networks to the electromagnetic field in nervous tissue

We present a microscopic approach for the coupling of cortical activity, as resulting from proper dipole currents of pyramidal neurons, to the electromagnetic field in extracellular fluid in presence of diffusion and Ohmic conduction. As a result, neural activity becomes represented by a continuous neural field equation, while an observation model for electric field potentials is obtained from the interaction of cortical dipole currents with charge density in non-resistive extracellular space as described by the Nernst-Planck equation.

preprint2010arXiv

Evidence for frequency-dependent extracellular impedance from the transfer function between extracellular and intracellular potentials

We examine the properties of the transfer function F_T = V_m / V_{LFP} between the intracellular membrane potential (V_m) and the local field potential (V_{LFP}) in cerebral cortex. We first show theoretically that, in the subthreshold regime, the frequency dependence of the extracellular medium and that of the membrane potential have a clear incidence on F_T. The calculation of F_T from experiments and the matching with theoretical expressions is possible for desynchronized states where individual current sources can be considered as independent. Using a mean-field approximation, we obtain a method to estimate the impedance of the extracellular medium without injecting currents. We examine the transfer function for bipolar (differential) LFPs and compare to simultaneous recordings of V_m and V_{LFP} during desynchronized states in rat barrel cortex in vivo. The experimentally derived F_T matches the one derived theoretically, only if one assumes that the impedance of the extracellular medium is frequency-dependent, and varies as 1/sqrt(omega) (Warburg impedance) for frequencies between 3 and 500 Hz. This constitutes indirect evidence that the extracellular medium is non-resistive, which has many possible consequences for modeling LFPs.