Determination of modular multilevel converters admittances and their impacts on HVDC power system stability

This paper deals with the determination of the complex admittances and trans-admittances of power converters in general and in particular of two modular multilevel converters (MMCs) connected by a high-voltage direct current (HVDC) link to two distinct AC systems. With the term trans-admittances we mean the cross admittance “seen” between the points of common coupling of the two MMCs in the two AC grids. These quantities show how the HVDC link couples the two seemingly separated AC systems. The accurate computation of admittances and trans-admittances requires careful consideration of several frequency harmonics and frequency up-conversion/down-conversion phenomena through the HVDC link. To determine them, we exploit the periodic small-signal (PAC) analysis that we briefly introduce from a theoretical perspective. Results show how AC and DC instabilities can be induced by MMCs and how they propagate through the HVDC link

Shooting method based modular multilevel converter initialization for electromagnetic transient analysis

This paper presents an accurate and efficient initialization strategy for modular multilevel converters (MMCs) based on the shooting method, a numerical technique aimed at deriving the periodic steady-state operating condition of any circuit. This technique is compatible with MMC models of different levels of detail and whose control scheme may include modulation strategies and capacitor voltage balancing algorithms. Electromagnetic transient simulations of the NORDIC32 power system modified by adding a high-voltage direct current link with 128-level MMCs prove that the proposed initialization strategy allows starting simulations close to steady-state, thereby significantly limiting initialization transients and their corresponding extra CPU time

Electric power-system's global-inertia estimation

This work presents a technique to estimate on-line the global inertia of an electric power system by exploiting the footprint of the principal frequency system dynamics. This method can estimate the inertia provided as a whole by synchronous machines, as well as by converter-interfaced generators controlled to emulate the behavior of the former through virtual inertia. Probing tones are injected by a grid-forming converter-interfaced generator and its virtual rotor speed is used to extrapolate its footprint. As a result, the method requires neither measuring the active power exchange of each synchronous generator nor extrapolating their rotor speeds. Since the proposed technique is entirely data driven, it does not require any model of the power system generators/prime-movers/controllers and of the interconnecting grid. The method is comprehensively tested on a modified version of the IEEE 39-BUS system and a dynamic version of the IEEE 118-BUS system, both containing grid-forming converter-interfaced generators

On the Firing Rate Dependency of the Phase Response Curve of Rat Purkinje Neurons In Vitro

Synchronous spiking during cerebellar tasks has been observed across Purkinje cells: however, little is known about the intrinsic cellular mechanisms responsible for its initiation, cessation and stability. The Phase Response Curve (PRC), a simple input-output characterization of single cells, can provide insights into individual and collective properties of neurons and networks, by quantifying the impact of an infinitesimal depolarizing current pulse on the time of occurrence of subsequent action potentials, while a neuron is firing tonically. Recently, the PRC theory applied to cerebellar Purkinje cells revealed that these behave as phase-independent integrators at low firing rates, and switch to a phase-dependent mode at high rates. Given the implications for computation and information processing in the cerebellum and the possible role of synchrony in the communication with its post-synaptic targets, we further explored the firing rate dependency of the PRC in Purkinje cells. We isolated key factors for the experimental estimation of the PRC and developed a closed-loop approach to reliably compute the PRC across diverse firing rates in the same cell. Our results show unambiguously that the PRC of individual Purkinje cells is firing rate dependent and that it smoothly transitions from phase independent integrator to a phase dependent mode. Using computational models we show that neither channel noise nor a realistic cell morphology are responsible for the rate dependent shift in the phase response curve

The what and where of adding channel noise to the Hodgkin-Huxley equations

One of the most celebrated successes in computational biology is the Hodgkin-Huxley framework for modeling electrically active cells. This framework, expressed through a set of differential equations, synthesizes the impact of ionic currents on a cell's voltage -- and the highly nonlinear impact of that voltage back on the currents themselves -- into the rapid push and pull of the action potential. Latter studies confirmed that these cellular dynamics are orchestrated by individual ion channels, whose conformational changes regulate the conductance of each ionic current. Thus, kinetic equations familiar from physical chemistry are the natural setting for describing conductances; for small-to-moderate numbers of channels, these will predict fluctuations in conductances and stochasticity in the resulting action potentials. At first glance, the kinetic equations provide a far more complex (and higher-dimensional) description than the original Hodgkin-Huxley equations. This has prompted more than a decade of efforts to capture channel fluctuations with noise terms added to the Hodgkin-Huxley equations. Many of these approaches, while intuitively appealing, produce quantitative errors when compared to kinetic equations; others, as only very recently demonstrated, are both accurate and relatively simple. We review what works, what doesn't, and why, seeking to build a bridge to well-established results for the deterministic Hodgkin-Huxley equations. As such, we hope that this review will speed emerging studies of how channel noise modulates electrophysiological dynamics and function. We supply user-friendly Matlab simulation code of these stochastic versions of the Hodgkin-Huxley equations on the ModelDB website (accession number 138950) and http://www.amath.washington.edu/~etsb/tutorials.html.Comment: 14 pages, 3 figures, review articl

Simple, Fast and Accurate Implementation of the Diffusion Approximation Algorithm for Stochastic Ion Channels with Multiple States

The phenomena that emerge from the interaction of the stochastic opening and closing of ion channels (channel noise) with the non-linear neural dynamics are essential to our understanding of the operation of the nervous system. The effects that channel noise can have on neural dynamics are generally studied using numerical simulations of stochastic models. Algorithms based on discrete Markov Chains (MC) seem to be the most reliable and trustworthy, but even optimized algorithms come with a non-negligible computational cost. Diffusion Approximation (DA) methods use Stochastic Differential Equations (SDE) to approximate the behavior of a number of MCs, considerably speeding up simulation times. However, model comparisons have suggested that DA methods did not lead to the same results as in MC modeling in terms of channel noise statistics and effects on excitability. Recently, it was shown that the difference arose because MCs were modeled with coupled activation subunits, while the DA was modeled using uncoupled activation subunits. Implementations of DA with coupled subunits, in the context of a specific kinetic scheme, yielded similar results to MC. However, it remained unclear how to generalize these implementations to different kinetic schemes, or whether they were faster than MC algorithms. Additionally, a steady state approximation was used for the stochastic terms, which, as we show here, can introduce significant inaccuracies. We derived the SDE explicitly for any given ion channel kinetic scheme. The resulting generic equations were surprisingly simple and interpretable - allowing an easy and efficient DA implementation. The algorithm was tested in a voltage clamp simulation and in two different current clamp simulations, yielding the same results as MC modeling. Also, the simulation efficiency of this DA method demonstrated considerable superiority over MC methods.Comment: 32 text pages, 10 figures, 1 supplementary text + figur

Human Pluripotent Stem-Cell-Derived Cortical Neurons Integrate Functionally into the Lesioned Adult Murine Visual Cortex in an Area-Specific Way

The transplantation of pluripotent stem-cell-derived neurons constitutes a promising avenue for the treatment of several brain diseases. However, their potential for the repair of the cerebral cortex remains unclear, given its complexity and neuronal diversity. Here, we show that human visual cortical cells differentiated from embryonic stem cells can be transplanted and can integrate successfully into the lesioned mouse adult visual cortex. The transplanted human neurons expressed the appropriate repertoire of markers of six cortical layers, projected axons to specific visual cortical targets, and were synaptically active within the adult brain. Moreover, transplant maturation and integration were much less efficient following transplantation into the lesioned motor cortex, as previously observed for transplanted mouse cortical neurons. These data constitute an important milestone for the potential use of human PSC-derived cortical cells for the reassembly of cortical circuits and emphasize the importance of cortical areal identity for successful transplantation. Espuny-Camacho et al. show that transplanted ESC-derived human cortical neurons integrate functionally into the lesioned adult mouse brain. Transplanted neurons display visual cortical identity and show specific restoration of damaged cortical pathways following transplantation into the visual but not the motor cortex, suggesting the importance of areal-identity match for successful cortical repair

The location of the axon initial segment affects the bandwidth of spike initiation dynamics

The dynamics and the sharp onset of action potential (AP) generation have recently been the subject of intense experimental and theoretical investigations. According to the resistive coupling theory, an electrotonic interplay between the site of AP initiation in the axon and the somato-dendritic load determines the AP waveform. This phenomenon not only alters the shape of AP recorded at the soma, but also determines the dynamics of excitability across a variety of time scales. Supporting this statement, here we generalize a previous numerical study and extend it to the quantification of the input-output gain of the neuronal dynamical response. We consider three classes of multicompartmental mathematical models, ranging from ball-and-stick simplified descriptions of neuronal excitability to 3D-reconstructed biophysical models of excitatory neurons of rodent and human cortical tissue. For each model, we demonstrate that increasing the distance between the axonal site of AP initiation and the soma markedly increases the bandwidth of neuronal response properties. We finally consider the Liquid State Machine paradigm, exploring the impact of altering the site of AP initiation at the level of a neuronal population, and demonstrate that an optimal distance exists to boost the computational performance of the network in a simple classification task. Copyright

High Bandwidth Synaptic Communication and Frequency Tracking in Human Neocortex

Neuronal firing, synaptic transmission, and its plasticity form the building blocks for processing and storage of information in the brain. It is unknown whether adult human synapses are more efficient in transferring information between neurons than rodent synapses. To test this, we recorded from connected pairs of pyramidal neurons in acute brain slices of adult human and mouse temporal cortex and probed the dynamical properties of use-dependent plasticity. We found that human synaptic connections were purely depressing and that they recovered three to four times more swiftly from depression than synapses in rodent neocortex. Thereby, during realistic spike trains, the temporal resolution of synaptic information exchange in human synapses substantially surpasses that in mice. Using information theory, we calculate that information transfer between human pyramidal neurons exceeds that of mouse pyramidal neurons by four to nine times, well into the beta and gamma frequency range. In addition, we found that human principal cells tracked fine temporal features, conveyed in received synaptic inputs, at a wider bandwidth than for rodents. Action potential firing probability was reliably phase-locked to input transients up to 1,000 cycles/s because of a steep onset of action potentials in human pyramidal neurons during spike trains, unlike in rodent neurons. Our data show that, in contrast to the widely held views of limited information transfer in rodent depressing synapses, fast recovering synapses of human neurons can actually transfer substantial amounts of information during spike trains. In addition, human pyramidal neurons are equipped to encode high synaptic information content. Thus, adult human cortical microcircuits relay information at a wider bandwidth than rodent microcircuits

Accurate and Fast Simulation of Channel Noise in Conductance-Based Model Neurons by Diffusion Approximation

Stochastic channel gating is the major source of intrinsic neuronal noise whose functional consequences at the microcircuit- and network-levels have been only partly explored. A systematic study of this channel noise in large ensembles of biophysically detailed model neurons calls for the availability of fast numerical methods. In fact, exact techniques employ the microscopic simulation of the random opening and closing of individual ion channels, usually based on Markov models, whose computational loads are prohibitive for next generation massive computer models of the brain. In this work, we operatively define a procedure for translating any Markov model describing voltage- or ligand-gated membrane ion-conductances into an effective stochastic version, whose computer simulation is efficient, without compromising accuracy. Our approximation is based on an improved Langevin-like approach, which employs stochastic differential equations and no Montecarlo methods. As opposed to an earlier proposal recently debated in the literature, our approximation reproduces accurately the statistical properties of the exact microscopic simulations, under a variety of conditions, from spontaneous to evoked response features. In addition, our method is not restricted to the Hodgkin-Huxley sodium and potassium currents and is general for a variety of voltage- and ligand-gated ion currents. As a by-product, the analysis of the properties emerging in exact Markov schemes by standard probability calculus enables us for the first time to analytically identify the sources of inaccuracy of the previous proposal, while providing solid ground for its modification and improvement we present here