The work we present in this thesis is a series of studies of how heterogeneities in coupling affect the synchronization of coupled neural oscillators. We begin by examining how heterogeneity in coupling strength affects the equilibrium phase difference of a pair of coupled, spiking neurons when compared to the case of identical coupling. This study is performed using pairs of Hodgkin-Huxley and Wang-Buzsaki neurons. We find that heterogeneity in coupling strength breaks the symmetry of the bifurcation diagrams of equilibrium phase difference versus the synaptic rate constant for weakly coupled pairs of neurons. We observe important qualitative changes such as the loss of the ubiquitous in-phase and anti-phase solutions found when the coupling is identical and regions of parameter space where no phase locked solution exists. Another type of heterogeneity can be found by having different types of coupling between oscillators. Synaptic coupling between neurons can either be exciting or inhibiting. We examine the synchronization dynamics when a pair of neurons is coupled with one excitatory and one inhibitory synapse. We also use coupled pairs of Hodgkin-Huxley neurons and Wang-Buzsaki neurons for this work. We then explore the existance of 1:n coupled states for a coupled pair of theta neurons. We do this in order to reproduce an observed effect called quantal slowing. Quantal slowing is the phenomena where jumping between different $1:n$ coupled states is observed instead of gradual changes in period as a parameter in the system is varied. All of these topics fall under the general heading of coupled, non-linear oscillators and specifically weakly coupled, neural oscillators. The audience for this thesis is most likely going to be a mixed crowd as the research reported herein is interdisciplinary. Choosing the content for the introduction proved far more challenging than expected. It might be impossible to write a maximally useful introductory portion of a thesis when it could be read by a physicist, mathematician, engineer or biologist. Undoubtedly readers will find some portion of this introduction elementary. At the risk of boring some or all of my readers we decided it was best to proceed so that enough of the mathematical (biological) background is explained in the introduction so that a biologist (mathematician) is able to appreciate the motivations for the research and the results presented. We begin with a introduction in nonlinear dynamics explaining the mathematical tools we use to characterize the excitability of individual neurons, as well as oscillations and synchrony in neural networks. The next part of the introductory material is an overview of the biology of neurons. We then describe the neuron models used in this work and finally describe the techniques we employ to study coupled neurons.

This dissertation presents analog circuit analysis and design from a nonlinear dynamics perspective. An introduction to fundamental concepts of nonlinear dynamical systems theory is given. The procedure of nondimensionalization is used in order to derive the state space representation of circuits. Geometric tools are used to analyze nonlinear phenomena in circuits, and also to develop intuition about how to evoke certain desired behavior in the circuits. To predict and quantify non-ideal behavior, bifurcation analysis, stability analysis and perturbation methods are applied to the circuits. Experimental results from a reconfigurable analog integrated circuit chip are presented to illustrate the nonlinear dynamical systems theory concepts. Tools from nonlinear dynamical systems theory are used to develop a systematic method for designing a particular class of integrated circuit sinusoidal oscillators. This class of sinusoidal oscillators is power- and area-efficient, as it uses the inherent nonlinearity of circuit components to limit the oscillators' output signal amplitude. The novel design method that is presented is based on nonlinear systems analysis, which results in high-spectral purity oscillators. This design methodology is useful for applications that require integrated sinusoidal oscillators that have oscillation frequencies in the mid- to high- MHz range. A second circuit design example is presented, namely a bandpass filter for front end auditory processing. The bandpass filter mimics the nonlinear gain compression that the healthy cochlea performs on input sounds. The cochlea's gain compression is analyzed from a nonlinear dynamics perspective and the theoretical characteristics of the dynamical system that would yield such behavior are identified. The appropriate circuit for achieving the desired nonlinear characteristics are designed, and it is incorporated into a bandpass filter. The resulting nonlinear bandpass filter performs the gain compression as desired, while minimizing the amount of harmonic distortion. It is a practical component of an advanced auditory processor.

[Truncated abstract] This thesis is split into two parts. In part one; A 50 K dual mode oscillator, the aim of the project was to build a 50 K precision oscillator with frequency stability on the order of 1014 from 1 to 100 seconds. A dual-mode temperature compensation technique was used that relied on a turning point in the frequency-temperature relationship of the difference frequency between two orthogonal whispering gallery modes in a single sapphire crystal. A cylindrical sapphire loaded copper cavity resonator was designed, modelled and built with a turning point in the difference frequency between an E-mode and H-mode pair at approximately 52.5 K . . . The frequencies and Q-factors of whispering spherical modes in the 3-12 GHz range in the fused silica resonator are measured at 6, 77 and 300 K and the Q-factor is used to determine the loss tangent at these temperatures. The frequency and Q-factor temperature dependence of the TM2,1,2 whispering gallery mode at 5.18 GHZ is used to characterise the loss tangent and relative permittivity of the fused silica from 4-300 K. Below 22 K the frequency-temperature dependence of the resonator was found to be consistent with the combined effects of the thermal properties of the dielectric and the influence of an unknown paramagnetic impurity, with a spin resonance frequency at about 138 ± 31 GHz. Below 8 K the loss tangent exhibited a 9th order power law temperature dependence, which may be explained by Raman scattering of Phonons from the paramagnetic impurity ions. A spherical Bragg reflector resonator made from multiple concentric dielectric layers loaded in a spherical cavity that enables confinement of field in the centre of the resonator is described. A set of simultaneous equations is derived that allow the calculation of the required dimensions and resonance frequency for such a resonator and the solution is confirmed using finite element analysis. A spherical Bragg reflector resonator is constructed using Teflon and free-space as the dielectric materials. A Q-factor of 22,000 at 13.87 GHz was measured and found to compare well with the design values.

Made available in DSpace on 2014-06-11T19:22:31Z (GMT). No. of bitstreams: 0 Previous issue date: 2010-08-30Bitstream added on 2014-06-13T19:27:59Z : No. of bitstreams: 1 pereira_mva_me_ilha.pdf: 1675496 bytes, checksum: e8bfb14cdd90155eb3c43096d4c160df (MD5) / Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) / Este trabalho apresenta um Oscilador Controlado por Corrente (CCO) com configuração em anel usando tecnologia CMOS, com melhorias na faixa de operação e ruído de fase. O oscilador proposto tem uma faixa de oscilação de 0,0989 GHz a 1,2 GHz com uma corrente de controle com um intervalo de 0,1 mA a 3 mA com uma potência dissipada de 11,8 mW. A arquitetura apresenta uma melhoria na fase de ruído de -7 dBc / Hz em relação a um oscilador em anel de três estágios (VCO), também apresentado neste trabalho. A estrutura proposta é baseada na mudança da entrada de controle do oscilador e também em modificações nas polarizações dos transistor de carga do estágio de atraso. Estas mudanças, além de aumentar a faixa de operação do oscilador e diminuir o efeito do ruído de fase, também reduzem a variação da amplitude do sinal de saída que acontece a medida que a frequência de operação aumenta ou diminui. Simulações realizadas com ambos os osciladores, confirmam os resultados. / This dissertation presents a Current Controlled Oscillator (CCO-Current-Controlled Oscillator) at ring configuration using CMOS (Complementary Metal-Oxide-Semiconductor) technology, with improvements in operating range and phase noise. The proposed oscillator has an oscillation range of 98.959 MHz to 1.2 GHz with a current control with a range of 0.1 mA to 3 mA with a power dissipation of 11.8 mW. The architecture shows an improvement in phase noise of -7 dBc / Hz when compared with a ring oscillator in three stages (VCO-Voltage- Controlled Oscillator), also presented in this paper. The proposed structure is in the change of input control and also in the polarizations of the load transistor stage of delay. These changes, in modifications increase the operations range of the oscillator, reduce the phase noise and minimize the amplitude variation of the output signal when the frequency operation increase or decrease. Simulations with both oscillators and their comparisons confirm these results.

No description available.

The human sense of hearing is characterized by its exquisite sensitivity, sharp frequency selectivity, and wide dynamic range. These features depend on an active process that in the inner ear boosts vibrations evoked by auditory stimuli. Spontaneous otoacoustic emissions constitute a demonstrative manifestation of this physiologically vulnerable mechanism. In the cochlea, sensory hair bundles transduce sound-induced vibrations into neural signals. Hair bundles can power mechanical movements of their tip, oscillate spontaneously, and operate as tuned nonlinear amplifiers of weak periodic stimuli. Active hair-bundle motility constitutes a promising candidate with respect to the biophysical implementation of the active process underlying human hearing. The responsiveness of isolated hair bundles, however, is seriously hampered by intrinsic fluctuations. In this thesis, we present theoretical and experimental results concerning the noise-imposed limitations of nonlinear amplification by active sensory hair bundles. We analyze the effect of noise within the framework of a stochastic description of hair-bundle dynamics and relate our findings to generic aspects of the stochastic dynamics of oscillatory systems. Hair bundles in vivo are often elastically coupled by overlying gelatinous membranes. In addition to theoretical results concerning the dynamics of elastically coupled hair bundles, we report on an experimental study. We have interfaced dynamic force clamp performed on a hair bundle from the sacculus of the bullfrog with real-time stochastic simulations of hair-bundle dynamics. By means of this setup, we could couple a hair bundle to two virtual neighbors, called cyber clones. Our theoretical and experimental work shows that elastic coupling leads to an effective noise reduction. Coupled hair bundles exhibit an increased coherence of spontaneous oscillations and an enhanced amplification gain. We therefore argue that elastic coupling by overlying membranes constitutes a morphological specialization for reducing the detrimental effect of intrinsic fluctuations.

Understanding the complex dynamics that underpin the response of structures in the occurrence of earthquakes is of paramount importance in ensuring community resilience. The operational continuity of structures is influenced by the performance of nonstructural components, also known as secondary structures. Inherent vulnerability characteristics, nonlinearities and uncertainties in their properties or in the excitation pose challenges that render their response determination as a non-straightforward task. This dissertation settles in the context of mathematical modelling and response quantification of seismically driven secondary systems. The case of bilinear hysteretic, rigid-plastic and free-standing rocking oscillators is first considered, as a representative class of secondary systems of distinct behaviour excited at a single point in the primary structure. The equations governing their full dynamic interaction with linear primary oscillators are derived with the purpose of assessing the appropriateness of simplified analysis methods where the secondary-primary feedback action is not accounted for. Analyses carried out in presence of pulse-type excitation have shown that the cascade approximation can be considered satisfactory for bilinear systems provided the secondary-primary mass ratio is adequately low and the system does not approach resonance. For the case of sliding and rocking systems, much lighter secondary systems need to be considered if the cascade analysis is to be adopted, with the validity of the approximation dictated by the selection of the input parameters. Based on the premise that decoupling is permitted, new analytical solutions are derived for the pulse driven nonlinear oscillators considered, conveniently expressing the seismic response as a function of the input parameters and the relative effects are quantified. An efficient numerical scheme for a general-type of excitation is also presented and is used in conjunction with an existing nonstationary stochastic far-field ground motion model to determine the seismic response spectra for the secondary oscillators at given site and earthquake characteristics. Prompted by the presence of uncertainty in the primary structure, and in line with the classical modal analysis, a novel approach for directly characterising uncertainty in the modal shapes, frequencies and damping ratios of the primary structure is proposed. A procedure is then presented for the identification of the model parameters and demonstrated with an application to linear steel frames with uncertain semi-rigid connections. It is shown that the proposed approach reduces the number of the uncertain input parameters and the size of the dynamic problem, and is thus particularly appealing for the stochastic assessment of existing structural systems, where partial modal information is available e.g. through operational modal analysis testing. Through a numerical example, the relative effect of stochasticity in a bi-directional seismic input is found to have a more prominent role on the nonlinear response of secondary oscillators when compared to the uncertainty in the primary structure. Further extending the analyses to the case of multi-attached linear secondary systems driven by deterministic seismic excitation, a convenient variant of the component-mode synthesis method is presented, whereby the primary-secondary dynamic interaction is accounted for through the modes of vibration of the two components. The problem of selecting the vibrational modes to be retained in analysis is then addressed for the case of secondary structures, which may possess numerous low frequency modes with negligible mass, and a modal correction method is adopted in view of the application for seismic analysis. The influence of various approaches to build the viscous damping matrix of the primary-secondary assembly is also investigated, and a novel technique based on modal damping superposition is proposed. Numerical applications are demonstrated through a piping secondary system multi-connected on a primary frame exhibiting various irregularities in plan and elevation, as well as a multi-connected flexible secondary system. Overall, this PhD thesis delivers new insights into the determination and understanding of the response of seismically driven secondary structures. The research is deemed to be of academic and professional engineering interest spanning several areas including seismic engineering, extreme events, structural health monitoring, risk mitigation and reliability analysis.

The human sense of hearing is characterized by its exquisite sensitivity, sharp frequency selectivity, and wide dynamic range. These features depend on an active process that in the inner ear boosts vibrations evoked by auditory stimuli. Spontaneous otoacoustic emissions constitute a demonstrative manifestation of this physiologically vulnerable mechanism. In the cochlea, sensory hair bundles transduce sound-induced vibrations into neural signals. Hair bundles can power mechanical movements of their tip, oscillate spontaneously, and operate as tuned nonlinear amplifiers of weak periodic stimuli. Active hair-bundle motility constitutes a promising candidate with respect to the biophysical implementation of the active process underlying human hearing. The responsiveness of isolated hair bundles, however, is seriously hampered by intrinsic fluctuations. In this thesis, we present theoretical and experimental results concerning the noise-imposed limitations of nonlinear amplification by active sensory hair bundles. We analyze the effect of noise within the framework of a stochastic description of hair-bundle dynamics and relate our findings to generic aspects of the stochastic dynamics of oscillatory systems. Hair bundles in vivo are often elastically coupled by overlying gelatinous membranes. In addition to theoretical results concerning the dynamics of elastically coupled hair bundles, we report on an experimental study. We have interfaced dynamic force clamp performed on a hair bundle from the sacculus of the bullfrog with real-time stochastic simulations of hair-bundle dynamics. By means of this setup, we could couple a hair bundle to two virtual neighbors, called cyber clones. Our theoretical and experimental work shows that elastic coupling leads to an effective noise reduction. Coupled hair bundles exhibit an increased coherence of spontaneous oscillations and an enhanced amplification gain. We therefore argue that elastic coupling by overlying membranes constitutes a morphological specialization for reducing the detrimental effect of intrinsic fluctuations.

A main thrust of this thesis is to develop and explore linearization-based numeric-analytic integration techniques in the context of stochastically driven nonlinear oscillators of relevance in structural dynamics. Unfortunately, unlike the case of deterministic oscillators, available numerical or numeric-analytic integration schemes for stochastically driven oscillators, often modelled through stochastic differential equations (SDE-s), have significantly poorer numerical accuracy. These schemes are generally derived through stochastic Taylor expansions and the limited accuracy results from difficulties in evaluating the multiple stochastic integrals. We propose a few higher-order methods based on the stochastic version of transversal linearization and another method of linearizing the nonlinear drift field based on a Girsanov change of measures. When these schemes are implemented within a Monte Carlo framework for computing the response statistics, one typically needs repeated simulations over a large ensemble. The statistical error due to the finiteness of the ensemble (of size N, say)is of order 1/√N, which implies a rather slow convergence as N→∞. Given the prohibitively large computational cost as N increases, a variance reduction strategy that enables computing accurate response statistics for small N is considered useful. This leads us to propose a weak variance reduction strategy. Finally, we use the explicit derivative-free linearization techniques for state and parameter estimations for structural systems using the extended Kalman filter (EKF). A two-stage version of the EKF (2-EKF) is also proposed so as to account for errors due to linearization and unmodelled dynamics. In Chapter 2, we develop higher order locally transversal linearization (LTL) techniques for strong and weak solutions of stochastically driven nonlinear oscillators. For developing the higher-order methods, we expand the non-linear drift and multiplicative diffusion fields based on backward Euler and Newmark expansions while simultaneously satisfying the original vector field at the forward time instant where we intend to find the discretized solution. Since the non-linear vector fields are conditioned on the solution we wish to determine, the methods are implicit. We also report explicit versions of such linearization schemes via simple modifications. Local error estimates are provided for weak solutions. Weak linearized solutions enable faster computation vis-à-vis their strong counterparts. In Chapter 3, we propose another weak linearization method for non-linear oscillators under stochastic excitations based on Girsanov transformation of measures. Here, the non-linear drift vector is appropriately linearized such that the resulting SDE is analytically solvable. In order to account for the error in replacing of non-linear drift terms, the linearized solutions are multiplied by scalar weighting function. The weighting function is the solution of a scalar SDE(i.e.,Radon-Nikodym derivative). Apart from numerically illustrating the method through applications to non-linear oscillators, we also use the Girsanov transformation of measures to correct the truncation errors in lower order discretizations. In order to achieve efficiency in the computation of response statistics via Monte Carlo simulation, we propose in Chapter 4 a weak variance reduction strategy such that the ensemble size is significantly reduced without seriously affecting the accuracy of the predicted expectations of any smooth function of the response vector. The basis of the variance reduction strategy is to appropriately augment the governing system equations and then weakly replace the associated stochastic forcing functions through variance-reduced functions. In the process, the additional computational cost due to system augmentation is generally far less besides the accrued advantages due to a drastically reduced ensemble size. The variance reduction scheme is illustrated through applications to several non-linear oscillators, including a 3-DOF system. Finally, in Chapter 5, we exploit the explicit forms of the LTL techniques for state and parameters estimations of non-linear oscillators of engineering interest using a novel derivative-free EKF and a 2-EKF. In the derivative-free EKF, we use one-term, Euler and Newmark replacements for linearizations of the non-linear drift terms. In the 2-EKF, we use bias terms to account for errors due to lower order linearization and unmodelled dynamics in the mathematical model. Numerical studies establish the relative advantages of EKF-DLL as well as 2-EKF over the conventional forms of EKF. The thesis is concluded in Chapter 6 with an overall summary of the contributions made and suggestions for future research.

Le projet HODISS dans le cadre duquel s'effectue nos travaux adresse la problématique de la synchronisation globale des systèmes complexes sur puce (System-on-Chip ou SOCs, par exemple un multiprocesseur monolithique). Les approches classiques de distribution d'horloges étant devenues de plus en plus obsolètes à cause de l'augmentation de la fréquence d'horloge, l'accroissement des temps de propagation, l'accroissement de la complexité des circuits et les incertitudes de fabrication, les concepteurs s’intéressent (pour contourner ces difficultés) à d'autres techniques basées entre autres sur les oscillateurs distribués. La difficulté majeure de cette dernière approche réside dans la capacité d’assurer le synchronisme global du système. Nous proposons un système d'horlogerie distribuée basé sur un réseau d’oscillateurs couplés en phase. Pour synchroniser ces oscillateurs, chacun d'eux est en fait une boucle à verrouillage de phase qui permet ainsi d'assurer un couplage en phase avec les oscillateurs des zones voisines. Nous analysons la stabilité de l'état synchrone dans des réseaux cartésiens identiques de boucles à verrouillage de phase entièrement numériques (ADPLLs). Sous certaines conditions, on montre que l'ensemble du réseau peut synchroniser à la fois en phase et en fréquence. Un aspect majeur de cette étude réside dans le fait que, en l'absence d'une horloge de référence absolue, le filtre de boucle dans chaque ADPLL est piloté par les fronts montants irréguliers de l'oscillateur local et, par conséquent, n'est pas régi par les mêmes équations d'état selon que l'horloge locale est avancée ou retardée par rapport au signal considéré comme référence. Sous des hypothèses simples, ces réseaux d'ADPLLs dits "auto-échantillonnés" peuvent être décrits comme des systèmes linéaires par morceaux dont la stabilité est notoirement difficile à établir. L'une des principales contributions que nous présentons est la définition de règles de conception simples qui doivent être satisfaites sur les coefficients de chaque filtre de boucle afin d'obtenir une synchronisation dans un réseau cartésien de taille quelconque. Les simulations transitoires indiquent que cette condition nécessaire de synchronisation peut également être suffisante pour une classe particulière d'ADPLLs "auto-échantillonnés". / The HODISS project, context in which this work is achieved, addresses the problem of global synchronization of complex systems-on-chip (SOCs, such as a monolithic multiprocessor). Since the traditional approaches of clock distribution are less used due to the increase of the clock frequency, increased delay, increased circuit complexity and uncertainties of manufacture, designers are interested (to circumvent these difficulties) to other techniques based among others on distributed synchronous clocks. The main difficulty of this latter approach is the ability to ensure the overall system synchronization. We propose a clock distribution system based on a network of phase-coupled oscillators. To synchronize these oscillators, each is in fact a phase-locked loop which allows to ensure a phase coupling with the nearest neighboring oscillators. We analyze the stability of the synchronized state in Cartesian networks of identical all-digital phase-locked loops (ADPLLs). Under certain conditions, we show that the entire network may synchronize both in phase and frequency. A key aspect of this study lies in the fact that, in the absence of an absolute reference clock, the loop-filter in each ADPLL is operated on the irregular rising edges of the local oscillator and consequently, does not use the same operands depending on whether the local clock is leading or lagging with respect to the signal considered as reference. Under simple assumptions, these networks of so-called “self-sampled” all-digital phase-locked-loops (SS-ADPLLs) can be described as piecewise-linear systems, the stability of which is notoriously difficult to establish. One of the main contributions presented here is the definition of simple design rules that must be satisfied by the coefficients of each loop-filter in order to achieve synchronization in a Cartesian network of arbitrary size. Transient simulations indicate that this necessary synchronization condition may also be sufficient for a specific class of SS-ADPLLs.