## Spike Talk: Genesis and Neural Coding Scheme Translations

Subham Sahoo, Senior Member, IEEE

Abstract -Although digitalization of future power grids offer several coordination incentives, the reliability and security of information and communication technologies (ICT) hinders its overall performance. In this paper, we introduce a novel architecture Spike Talk via a unified representation of power and information as a means of data normalization using spikes for coordinated control of microgrids. This grid-edge technology allows each distributed energy resource (DER) to execute decentralized secondary control philosophy independently by interacting among each other using power flow along the tie-lines. Inspired from the field of computational neuroscience, Spike Talk basically builds on a fine-grained parallelism on the information transfer theory in our brains, particularly when neurons (modeled as DERs) transmit information (inferred from power streams measurable at each DER) through synapses (modeled as tie-lines). Not only does Spike Talk simplify and address the current bottlenecks of the cyber-physical architectural operation by dismissing the ICT layer, it provides intrinsic operational and cost-effective opportunities in terms of infrastructure development, computations and modeling. Hence, this paper provides a pedagogic illustration of the key concepts and design theories. Since we focus on coordinated control of microgrids in this paper, the signaling accuracy and system performance is studied for several neural coding schemes responsible for converting the realvalued local measurements into spikes.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00773v2-figure-0272dffa-f54c-437b-b2e8-9f937a9a637f.jpg)

Index Terms -Neuromorphic computing, microgrids, hierarchical control, spiking neural network.

## I. INTRODUCTION

D IGITALIZATION of power electronic dominated grids has brought new data-driven directions in health monitoring and coordination of ever-scaling distributed energy resources (DERs). As we move towards decentralized generation, the control and reliability of DERs demand orchestration from information and communication technologies (ICT) to assist with system coordination, as shown in Fig. 1(c). However, these digitalization efforts not only have increased cyber attack vulnerabilities [1], but have added more unreliability metrics in the form of latency, data dropouts, etc. [2]. This mandates an intrinsic coordination principle for power grids that can effectively unify power &amp; information to enhance system reliability.

Apart from the communication reliability issues, the layered hierarchy of the ICT layer overlaying the power grids fabricates strong dependencies on each other, thereby introducing modeling complexity [3]. With increased penetration of ultra-efficient power electronic converters operating at a switching frequency within 100-500 kHz enabling faster control response, it simultaneously mandates upgradation of the ICT infrastructure to level up the coordination efforts. With

Subham Sahoo is an Associate Professor with the Department of Energy, Aalborg University, Denmark. (e-mail: sssa@energy.aau.dk)

Fig. 1. As opposed to the current cyber-physical architecture in (a), an innovative co-transfer technology (b) Spike Talk emanates out of novel protocol, namely publish-subscribe , that supports the inferential exchange of information from power flows in the network. It is a software-defined methodology that offers cognitive abilities to DERs and their coordination by providing numerous benefits without requiring any infrastructural changes.

modeling complexities of the cyber-physical interactions and infrastructural changes still in place as key research questions to be solved, co-transfer technologies that can modulate data over power streams emerged as a monolithic solution to all the unreliability issues in the cyber layer above [4], [5]. The basic principle behind co-transfer technologies is to exploit each DER 1 with a power electronic interface as a transceiver and the tie-line between them as a communication channel. Finally, control parameter adaption is performed at each DER with physical disturbances set as the triggering mechanism.

Differently from the conventional power line communication [6], Power Talk proposed in [7], allows to send information by modulating the bus voltages and corresponding droop parameter and consequently create power deviations, that can be detected and decoded by other converters. This one-way communication model is then extended in [8] for its constellation design and error analysis in the signaling and detection space. From a physical layer standpoint, this technology demystifies the re-design of the primary controller into a communication paradigm with signal transceiver abilities. Another emerging conversion technology, namely

1 Since power electronic converters are equipped in each DER, we use the term DER to represent DC/DC power electronic converter throughout this paper.

Talkative Power [9], [10] embeds information into the pulse width modulation (PWM) stage instead, with a clear rationale of extracting information as harmonics on eliminating the filtering stages. Not only the communication rate can be improved by modifying the switching characteristics, but it does not require additional hardware. An overview behind various co-transfer technologies that transmit power &amp; data along tie-lines can be found in [11]. Although they guarantee good reliability, they are highly dependent on graph theoretic concepts for ensuring low error rates and convergence of information transfer. Furthermore, they follow the request-respond two-way communication protocol (as shown in Fig. 1), which serves as a common push-pull communication approach and limits co-transfer scalability among multiple DERs in large power electronic systems.

From a cyber-physical system perspective, unified transmission of power &amp; information along the physical layer itself in Fig. 1 is capable of solving many issues altogether with: (i) no infrastructural changes, (ii) dismissal of the cyber layer, its inconsistent performance, vulnerabilities and interactions with the physical layer, (iii) simplified modeling. With the requestresponse protocol limiting the scalability of information transfer in large systems, we firstly explore the publish-subscribe communication protocol [12] for power grids in Fig. 1 that facilitates a fluctuation-driven regime, such that the control parameters adaptation using co-transfer technologies are governed solely based on the proximity of selective DERs to the disturbance's spatial location and impact itself. The protocols in Fig. 1 can be explained clearly using the simple example of a newsletter, wherein although Alice, Bob, and Cindy are subscribed to it, they read it only when it interests/concerns any of them. Following up on our preliminary evidence to locally estimate remote signals in power grids using task-aware semantic metrics [13], [14], we can integrate the publishsubscribe based inferential protocol for the same task, only when there is an event. Although the abovementioned work was primarily designed to defend against cyber anomalies in a cyber-physical system setup in Fig. 1, the root cause behind the problem, i.e., communication traffic, jamming, etc. still persists due to the presence of the cyber layer itself. As we transit towards the future power grid architecture as a singlelayer architecture in Fig. 1 enabling independence from the cyber layer, we need an effective tool at each node that can locally synthesize and adapt a decentralized secondary control policy from the power flow data itself.

In this paper, we draw inspiration from the field of computational neuroscience in understanding how several independent sensors act together in our brain to understand perception. This can be attributed to a comprehensive information mechanism that accounts for spatio-temporal data and its processing in our brain. Using a similar operational philosophy without being circumscribed by graph theory, we propose Spike Talk , which treats each DER as an event-driven sensor to infer global information and adapt their local controllers accordingly only using power flows as the medium between them. It is worth notifying that Spike Talk fundamentally differs from Power Talk and Talkative Power due to its new publish-subscribe protocol described above. As a result, the proposed co-transfer technology does not demand transfer of information from other nodes, but should rather be seen as an on-top software-defined control logic that extracts taskoriented data from power flows and its dynamics to adapt DERs' response in the network. More details on a comprehensive comparative analysis between these co-transfer technologies can be obtained in our recent work [15]-[17]. Unlike our previous work [18] in showcasing its energyefficient computational abilities and online learning, we deep dive into the following aspects in this paper:

- 1) Genesis: Since co-transfer of power and information using Spike Talk is fundamentally driven by the information transfer mechanism between neurons in our brains, we first provide a pedagogic illustration to point out the electrical modeling similarities between neurons and DERs. We then use the component-level modeling preliminaries and extend them to a system-level analogy for both neuronal networks and power grids. Firstly, the electrical model of leaky integrate-and-fire (LIF) neuron and a generalized equivalence with converter topology is explored for structural mapping. Secondly, a selective process of task-oriented sparse data is collected for online training, where we leverage the RC dynamics of the LIF neuron and translate the same data collection philosophy for each DER. Using Spike Talk, the entire power electronic grid can be emulated as a multi-sensory perceptive network with DERs at each node inferring via power flows only.
- 2) Neural Coding Schemes: Using the distributed computation facility established at each DER by Spike Talk, the system is adapted only during events represented by sparse dynamic data. Hence, the most critical stage is encoding real analog-valued measurements into sparse binary spikes. This stage affects the following aspects simultaneously:
- · signalling accuracy and system performance,
- · computation time (sparsity in data),
- · computation overheads (online learning).

Given that the secondary controller response is typically expected in seconds [19], Spike Talk is well equipped for this requirement, since its online learning capability can be achieved within general neuronal reaction time of 100-200 ms.

In this paper, we validate these two aspects in an exemplary DC microgrid. With the comparative evaluation of Spike Talk and the existing co-existing technologies already covered in [18], we particularly emphasize on its modeling part and signalling consequences due to different neural coding schemes in this paper. With the modeling preliminaries as a key foundation, other stages of Spike Talk will be introduced later in future scope of this article.

## II. UNRAVELING THE MODELING ANALOGIES

While there are a range of bio-physically accurate neuron models in the field of computational neuroscience from the most complicated ones being Hodgkin-Huxley and Izhikevich models to simple artificial neurons [20], we select the leaky

Fig. 2. Modeling analogies between neurons and DER: (a) neuron membrane and synpases, (b) the electrical modeling of the membrane as a RC circuit, which is retrofit into the output capacitor and equivalent tie-line resistance seen from a DC/DC converter, (c) Similar to the pre-synpatic and post-synaptic connection in (a), two DC DERs (emulated as neurons) are connected to each other via tie-lines (emulated as synpase) to transmit ingrained information via power, (d) structural mapping of each component in the brain and power grid to construct a duality principle.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00773v2-figure-062681c5-98f8-4563-ad1b-e22091a942fd.jpg)

integrate-and-fire (LIF) model since it provides a good tradeoff between the biological accuracy and complexity [21]. Unlike artificial neural networks (ANNs) process the weighted sum of inputs using an activation function, LIF neurons on the other hand integrate the weighted sum input over time with a leakage, much like a RC circuit. In translation, we exploit the same RC dynamics in output of each DER to collect sparse data and incorporate online learning. Before continuing with the sparsity of data, let us understand the modeling similarities between a LIF neuron and DERs in the next subsection.

## A. Component Level Modeling

Our brain has around billions of neurons. Using multiple learning instances in the past, it can sense, perceive and then take dynamic, complex decisions quickly. For a given stimuli, neurons fire to transmit an information using their network to the brain. As shown in Fig. 2(a), two neurons are connected to each other via a synpase, which is responsible for transferring electro-chemical ions as information via their neuro-transmitters. Since the neurons are placed on either side of the synpase, they are called preand post-synpatic neurons . Furthermore, each neuron is surrounded by a thin bi-layer membrane that insulates the intracellular portion from the extracellular medium. Due to this property where two conductive solutions are separated by an insulator, this membrane can be electrically interpreted as a capacitor C . This membrane also exhibit an impermeable property that limits the passage of ions (Na + in Fig. 2(a)) from the extracellular medium inside the neuron. However, there are specific membrane channels that are triggered to open when a current is injected into the neuron, which mimics the electric property of a resistor R . Considering an arbitrary time-varying current I(t) in Fig. 2(b) injected into a neuron, it can be written as:

$$I ( t ) = \frac { V _ { \text{mem} } ( t ) } { R } + C \frac { d V _ { \text{mem} } } { d t } \quad \quad ( 1 ) \quad \text{where} \quad \text{has} \quad \text{that} \quad \,.$$

with τ = RC being the circuit time constant, V mem being the voltage across the membrane. Hence, the LIF neuron membrane can be modeled as a RC circuit, where the first and second term in (1) corresponds to the leakage part and integrate-and-fire part, respectively. In this way, the passive membrane of neuron can be described by a linear differential equation in (1). Considering the derivative dV mem dt substituted into (1) without considering the limit ∆t → 0, we get:

$$\underset { \text{$n$ the $a$} } { \text{thed} } \quad \tau \frac { V _ { \text{mem} } ( t + \Delta t ) - V _ { \text{mem} } ( t ) } { \Delta t } = - V _ { \text{mem} } ( t ) + I ( t ) R \quad ( 2 )$$

Isolating the membrane potential for the next time-step in LHS and assuming input current I(t) = 0 to single out the resistive leaky dynamics, we simplify (1) into:

$$V _ { \text{mem} } ( t + \Delta t ) = ( 1 - \frac { \Delta t } { \tau } ) V _ { \text{mem} } ( t ) \quad \quad ( 3 )$$

Considering the ratio of the subsequent membrane potential values, V mem (t + ∆t) V / mem (t) to be the decay rate β , which can then be mathematically denoted by (1 -∆t /τ ) . However, it is preferred to calculate β using e -∆t /τ , since ∆t &lt;&lt; τ [22].

Remark I: Similar to the electrical RC circuit model of a LIF neuron in (1), we exploit the same network elements in a DC DER in Fig. 2(b) using its output capacitance C o and the equivalent network resistance R eq to formulate a physics-informed threshold based criteria for sparse data collection. The rationale behind using R eq is that it can be adapted/controlled in alignment with power flow dynamics.

## B. System Level Modeling

Using the component modeling for each DER in Remark I as a pre-requisite groundwork, we will now expand it to a system-level control policy in this subsection. For the current carried by the ions in Fig. 2(a), the most critical dynamic variable is the potential difference between the pre-synaptic k and post-synaptic neuron l . This can be attributed to the dynamics of the potential using Kirchhoff's current law:

$$g _ { k l } ( V _ { k } - V _ { 1 } ) = \underbrace { I _ { k } - C _ { k } \frac { d V _ { k } } { d t } } _ { \text{Neuron } k } + \underbrace { C _ { l } \frac { d V _ { 1 } } { d t } - I _ { l } } _ { \text{Neuron } l } \quad ( 4 )$$

where, the term in LHS represent the potential difference between the neuron membranes k and l , whereas the terms in RHS correspond to the capacitive properties of both neurons as well as their synaptic currents I k and I l . It is worth notifying that the g kl in (2) denote the synaptic conductance. Using (1), (2) can be re-written for each neuron as:

$$\text{:ntial} \quad g _ { k l } ( V _ { k } - V _ { 1 } ) = g _ { k } ( V _ { m e m } ^ { k } - V _ { t h } ) - g _ { l } ( V _ { m e m } ^ { l } - V _ { t h } ) \ \ ( 5 )$$

where, g k and g l denote the membrane conductance of neuron k and l , respectively. Using the primary controller modeling preliminaries [23], we can easily translate the electrical equivalent circuit of each neuron in Fig. 2(a) to that of a DC DER. Comparing the multi-neuron and multi-DER system in Fig. 2(a) and 2(c) for a system-level analogy, a duality principle can be inferred, that has been elaborated in Fig. 2(d).

Remark II: With two neurons having their corresponding membrane potential across their surfaces, this can be equivalenced to the voltages across each DER. Furthermore, the synaptic cleft connects the neurons for transmission of information between them, which can be electrically interpreted as synaptic currents. Using the structural mapping in Fig. 2(d), we hypothesize to transmit ingrained information in the current flowing between two DERs by observing their dynamics and response to any disturbances.

Fig. 3. (a) Primary control of a DC DER with R d and ∆R as the static and adaptive droop values, respectively, (b) Single line diagram of the 4-bus DER based microgrid, (c) Spiking feedforward network with LIF neurons.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00773v2-figure-6d5e52af-f641-4edc-98cb-a33322a6ca77.jpg)

Using the structural mapping formalism in Remark II, the current flowing through the line with a tie-line resistance r kl from DER k to DER l operating at a global voltage reference V ref can be given by:

$$r _ { k l } ^ { - 1 } ( V _ { k } - V _ { 1 } ) = \underbrace { \frac { V _ { k } - V _ { r e f } } { R _ { d } ^ { k } + \Delta R _ { k } } } _ { \text{DER} \, k } - \underbrace { \frac { V _ { 1 } - V _ { r e f } } { R _ { d } ^ { 1 } + \Delta R _ { 1 } } } _ { \text{DER} \, l } \quad ( 6 ) \quad \text{power} \\ \text{such $t$} \\ \text{constit}$$

where, R d is a droop parameter for each DER, that can be calculated by ∆V I max with ∆V and I max denoting the maximum allowed voltage deviation and rated current of the DER, respectively. It is worth notifying that we exploited the primary control model in Fig. 3 to simplify (4). Since R d is a static parameter, our objective is to update ∆R to instill a decentralized secondary control policy for each DER, which is corrected in accordance with the power flow dynamics defined as an event . These dynamics are firstly converted into taskoriented spikes, which are then processed by the LIF based neuronal computational unit at each DER to compute ∆R .

Remark III: Explicit network level computational burden with ANN at each node is significantly descended to a neuron level computation using the LIF neuron model for each DER, that not only decreases the computational overhead of the entire system, but also pans out as an energy-efficient and scalable data-driven control methodology for large power electronic systems.

## C. Sparse Data Collection

The data &amp; computational sparsity in Spike Talk can totally be attributed to the binary nature of the data and physicsinformed data collection mechanism. This is achieved simultaneously by collecting spikes that are encoded from real-valued measurements only when a dynamic threshold is reached. Spikes can be simply defined as electrical pulses, that can have multi-modal representations. For example, a single spike can simultaneously represent the voltage and current information. Based on the neuron firing criteria, a cluster of such spikes can serve to reconstruct dynamic profiles of voltages and currents. However, the said criteria can be set up such that spikes are generated in an asynchronous manner. This criteria can be determined analytically as a solution to the LIF neuron in (1).

Applying a step current input in I(t) that switches on at t = t 0 , the general solution can be denoted as:

$$V _ { \text{mem} } ( t ) = I ( t ) R + [ V _ { 0 } - I ( t ) R ] e ^ { - t / \tau } \quad \ ( 7 )$$

where, V 0 is the steady-state voltage. For V mem = 0 V, we can re-write (7) as:

$$V _ { \text{mem} } ( t ) = I ( t ) R [ 1 - e ^ { - t / \tau } ]$$

Replacing e -∆t /τ with β in (8) as per (3), we can exploit a non-physiological assumption to convert the term (1β ) into a learnable weight matrix w to obtain:

$$\ w V _ { \text{mem} } ( t ) = \mathbf I ( t )$$

where, V mem and I denote the vector representations of membrane potential and currents from each neuron. As a result, the membrane voltage difference scaled by the synaptic conductances in w constitute the input current I , which has been shown in the spiking feedforward network in Fig. 3(c). Given that (9) commands good biological plausibility and accuracy, it is also highly congruent to the networked DC power flow theory in a multi-DER microgrid in Fig. 3(b), such that the potential difference between each DER buses constitute the current flow between them. It is worth notifying that (9) is highly congruent to the DC-OPF equation B θ = P [24], where the difference in bus angles θ leads to power flow between them.

With (9) serving as an instantaneous representation of the synaptic current for a given spike, we can analyze the membrane potential of the neighboring DER (post-synaptic neuron) with synaptic current as the input using:

$$V _ { \text{mem} } ( t + 1 ) = \underbrace { \beta \mathbf V _ { \text{mem} } ( t ) } _ { \text{Decay} } + \underbrace { \mathbf w \mathbf X ( t ) } _ { \text{Input} } - \\ \underbrace { S _ { \text{out} } ( \beta \mathbf V _ { \text{mem} } ( t ) + \mathbf w \mathbf X ( t ) ) } _ { \text{Reset-to-zero} } ( 1 0 ) \\... \. \ \mathbf v \wedge \.. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \........................................................................................................................................................................................................ \..................................................................................................................................................................................................... \. \................................................................................................................................................................................................... \... \................................................................................................................................................................................................... ....................................................................................................................................................................................................... ... ................................................................................................................................................................................................... ... \................................................................................................................................................................................................... \. ..................................................................................................................................................................................................... . ..................................................................................................................................................................................................... \. . ................................................................................................................................................................................................... . . ................................................................................................................................................................................................... \... ................................................................................................................................................................................................... \. \. \................................................................................................................................................................................................. \..... \................................................................................................................................................................................................. \... \. \................................................................................................................................................................................................. \. \... \................................................................................................................................................................................................. \. \. \. \............................................................................................................................................................................................... \....... \...............................................................................$$

where, X (t) is an unweighted current from remote DERs. Considering different input current step inputs in (10), it can be seen in Fig. 4(a) that the membrane potential rises upto a resting potential V th2 = 0.4 V, which is where the reset

Fig. 4. LIF neuron model for generating spikes: (a) membrane potential for a synpatic current of 0.5 A only generates spikes with the lower voltage threshold V th1 = 0.4 V, (b) when the synaptic current is increased to 1 A, regularly interspaced spikes are generated for both the voltage thresholds.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00773v2-figure-51a52dc2-cb61-4984-8c2c-50953c245f63.jpg)

mechanism brings it back to zero potential, thereby generating output spikes. In this paper, we have simply used reset-tozero method as the reset mechanism. More details on reset mechanisms with LIF neuron firing can be found in [20]. However for the same input current of 0.5 A introduced as a step input with V th1 raised to 0.8 V, the membrane potential exponentially relaxes to I(t)R , which is smaller than V th1 leading to no output spikes. At the same time, when the input current is increased to 1 A in Fig. 4(b), spikes are generated at different spacing intervals for the resting potential V th1 = 0.8 V. Hence, as opposed to continuously sampled real membrane potential values in Fig. 4, we convert them into periodic interspaced spikes using a neuron firing criteria. Since it returns periodic spikes even for a constant value of input current, it still doesn't sufficiently enforce an eventdriven sparsity. The sparsity in spikes can be further enhanced by using an asynchronous neuron firing criteria, commonly known as neural coding schemes , that will be described in detail in the next section.

## III. NEURAL CODING SCHEMES AND DEPLOYMENT

In this section, we learn about encoding of real-valued signals into binary spikes, which is a common way of data representation in our brain and its perception. These perceptive signatures are basically a statistical reflection of either the signal magnitude or latency, which needs encoding laws that can process big data into few and qualitative batches of spikes.

Fig. 5. Neural coding schemes to represent dynamics in time-domain signals using: (a) rate coding, (b) latency coding, (c) burst coding. It is important to note that burst is a group of short inter-spaced interval (ISI) spikes.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00773v2-figure-008625e1-a94a-4d0e-969c-5eb19274af37.jpg)

## A. Rate Coding

With extensive studies already carried out determining the rate-coding corresponding to the pixel intensity in an image, we conceptualize a similar analogy to time-domain signals based on the magnitude. Since our focus in this paper is on inferring remote information only during events, we can deduce the generation of rate-coded spikes Ω r using:

$$\Omega _ { r } = \begin{cases} 1 \text{ if } \dot { V } _ { k } ( t ) \geq \kappa _ { 1 } e _ { v } ( t ) \\ 0 \text{ else} \end{cases}$$

where, e v (t) = [ V ref -V (t) + I(t)(R k d +∆R(t)) ] is the voltage controller error in Fig. 3(a) and κ 1 is a tuning parameter. As shown in Fig. 5(a), the spike train profile is different with respect to the different dynamic stages. This coding scheme is quite straight-forward, since (11) doesn't necessarily instill any dependency on model parameters.

## B. Latency Coding

Latency coding is derived from a notion of elapsed time to indirectly denote the spike timing based on the feature intensity in a signal. For example, in Fig. 5(b), although only one spike is generated for each time window, it is the elapsed time from the initial set-point t 0 for each spike, that becomes an indirect representative of the signal magnitude. This is a very widely used signal coding technique to process sparse data for spiking neural networks (SNNs) [25]. For a constant current injection in (8), V mem will exponentially reach a steady-state value of I(t)R . Given that a spike is emitted when (8) reaches a threshold of V th , we can then obtain t using:

$$t = \tau [ \ln ( \frac { I R } { I R - V _ { t h } } ) ].$$

Using (12), it can be inferred that the larger input current, the convergence of V mem (t) into V th will be faster with fleeting spikes.

Remark IV: Since (11) is dependent on a nodal variable voltage that varies for each DER, rate coding can generally depict nodal voltage dynamics, whereas the spike generation using latency coding in (12) is primarily dependent on currents flowing between two DERs, it is more appropriate for network level depiction.

Hence, the spikes generation criteria using latency coding with the steady-state feature IR can simply be given by:

$$t ( I R ) = \begin{cases} \tau [ \ln ( \frac { I R } { I R - V _ { t h } } ) ], \text{ if } I R > V _ { t h } \\ \infty, \text{ else} \end{cases} \quad ( 1 3 )$$

## C. Burst Coding

With several bottlenecks in efficiency behind information transmission, burst coding has been an effective tool for fast and energy-efficient inferences using high frequency burst spikes. As shown in Fig. 5(c), burst spikes are a group of low interspaced interval (ISI) spikes next to each other. Basically, burst coding can simply be executed by converting the constant tuning parameter κ in (11) into a burst functionbased parameter κ b to obtain burst spikes Ω b using:

$$\kappa _ { b } ( t ) = \begin{cases} \Xi \kappa _ { b } ( t - 1 ), \text{ if } \Omega _ { b } ( t - 1 ) = 1 \\ 1, \text{ else} \end{cases} \quad ( 1 4 )$$

where, Ξ is a constant parameter. In this way, the activation threshold is adapted during dynamic conditions after every first spike using (14). Although this coding method can control the precision of the transmitted information, it also denotes that significantly more number of spikes needs to be transmitted for the same information that can be characterized by one spike. Hence, the selection criteria behind different coding methods build substantially on a trade-off between precision and the total number of spikes.

## D. Deployment Specifics

As illustrated in Fig. 6(b), Spike Talk operates and learns by utilizing sparse output spike data to label remote features of neighboring DERs. These spikes are generated based on the neural coding schemes described previously, and can be locally deployed as parallel workloads in the controller of each DER. Although the selectivity of qualitative binary spike-based data is achieved post-coding, it is still derived from fixed-time sampled measurements equipped on each DER. For each event in the network (see Fig. 6(a)), the input spikes Ω i are captured using the following coding strategy:

$$\Omega _ { i } = \begin{cases} 1 \text{ if } \dot { I } _ { k } ( t ) \geq \kappa _ { 2 } e _ { k } ( t ) & \text{ where } \\ 0 \text{ else } & \text{ the av} \\ & \text{ to } V _ { r e } \end{cases}$$

where, e (t) = [I i ref -I k (t)] represents the current controller error as depicted in Fig. 6(a), and κ 2 is a tuning parameter. Given that the time constant of the current controller is lower than any other control stage in Fig. 6(a), the state-dependent threshold in (15) justifies the collection of high-resolution input data. This process is highlighted in Fig. 6(b), which is then further refined using various network-based thresholds and DER output RC dynamics to filter semantic output data. The mechanism for online learning and continuous weight adaptation of each DER will be detailed in the next section.

## IV. CONTINUAL ONLINE LEARNING

As already outlined in [18], the learning process in neuroscience is closely related to synaptic plasticity, which can then be attributed to the update of tie-line resistances as per the analogy in Fig. 2(d).

Remark V: From a computational perspective, the notion of plasticity refers to either an increase or decrease of synaptic weights, however from a power network with a fixed topology, the tie-line resistance can instead be updated by adapting the dynamic droop value ∆R in Fig. 3 for each DER.

In comparison to supervised learning that is primarily steered on minimizing the error using backpropagation, Hebbian-based unsupervised learning operates on the relative firing time regime between the pre- and post-synaptic neurons. Going back to the context of DERs in Remark V, whenever there is an increase in the system demand, it increases the power generation from DERs, which is also denoted as an excitatory response. Otherwise, it is called as an inhibitory response corresponding to the decrease in power generation from each DER. Based on their definition and neuronal activity in computational neuroscience, excitatory and inhibitory responses correspond to the increase and decrease in adaptive droop ∆R , respectively.

## A. Voltage Regulation

Voltage regulation as a secondary control objective is widely studied for microgrids and active distribution systems, which is performed differently for centralized or distributed communication topologies [26]. The basic idea behind this objective is to find a corresponding setpoint that allows equivalence of the mean of the DER voltages across the network close to its global reference V ref . Mathematically, it can also be denoted as ∑ N i=1 V i ≈ V ref , with N indicating the total number of DERs. Using Remark V, we can convert the nodal voltage dynamics into spikes via rate coding using (11). This can be simplified by analyzing the online learning framework in Fig. 6(c) performed at DER 1 in Fig. 6(a).

Assuming I 4 = 0, we can approximate the nodal voltages at DERs 2 &amp; 3 in Fig. 6(a) using the network equations and some approximations based on the global voltage regulation objective via:

$$\begin{cases} V _ { 2 } ^ { \prime } ( t ) = V _ { 1 } ( t ) + \psi _ { 1 2 } r _ { d } I _ { 1 } ( t ) \\ V _ { 3 } ^ { \prime } ( t ) = V _ { 1 } ( t ) + \psi _ { 1 3 } ( r _ { a } + r _ { c } ) I _ { 1 } ( t ) \end{cases} \quad \text{(1.6)}$$

where, ψ 12 and ψ 13 denote the regulation coefficient such that the average of all steady-state voltages of the network sum up to V ref . Using (16) as the threshold values for triggering spikes for each DER, we can firstly compute the voltage regulation spikes S (t) v with the help of the rate coding based triggering criteria in (11). Furthermore, it is worth notifying that the sum of all the regulation coefficients in a network must be equal to 1. As shown in Fig. 6(b), different static thresholds V ′ 2 and V ′ 3 are used to formalize the individual spike profile for DER 2 &amp; 3, which then act as an input for the calculation of adaptive droop ∆R .

## B. Power Management

We will now consider different power management schemes, that are particularly driven by either generation/system cost [27], unreliability (electrical and thermal) [28], etc. Since these conditions are dynamic, the energy management scheme augments a certain proportionality value α k such that the currents from the k th DER is shared either equally or disproportionately among other DERs. Due to the low-pass filter disparity accounted into the measurements used in primary controller, there will be a delay in the rise/fall of currents in response to a given disturbance. Hence using Remark IV, we can use latency coding in (12) to investigate the competing rising time between currents with an activation voltage set up for each DER, given by:

$$\begin{cases} V _ { 2 \text{ref} } ^ { \prime } ( t ) = V _ { \text{ref} } - R _ { \text{d} } ^ { 2 } I _ { 1 } ( t ) \\ V _ { 3 \text{ref} } ^ { \prime } ( t ) = V _ { \text{ref} } - R _ { \text{d} } ^ { 3 } I _ { 1 } ( t ) \end{cases}$$

where, R k d = α k R d . Using latency coding, we finally obtain the current sharing spikes S (t) i , as shown in Fig. 6(c). It can be observed that there is a relative firing time difference ∆t = t S v -t S i between the S (t) v and S (t) i layer, that ultimately determines the synaptic weight between the (presynaptic) DER and its (post-synaptic) neighboring DERs.

Fig. 6. Implementation overview of Spike Talk - (a) event detection: each DER will have its local coding mechanism corresponding to each event. Since there is no cyber layer, the controllers are only equipped with primary control scheme, (b) the input spikes are continuously generated for every dynamics, however, the semantic output spikes are filtered using different network based voltage thresholds. These spikes are finally processed in tandem with different control objectives, (c) Continual online learning for each DER - the nodal currents and voltages are firstly being processed into spikes by rate and latency coding based on Remark IV. Finally, the relative timing based plasticity rules [18] were used to adapt the adaptive droop parameter ∆R = a∆w . It is worthy notifying that the update of ∆R is a system-centered approach, which is administered by the power dynamics seen by each DER.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00773v2-figure-6cdb8f03-f104-4aa1-8e50-22f46f6f6591.jpg)

Following the Hebbian unsupervised learning rule, the change in synaptic weight ∆w for k th DER varies as per:

implementation is governed using the power flow dynamics seen by each DER, which is the crux behind Spike Talk .

$$\Delta w _ { k } = \begin{cases} \sum _ { f = 1 } ^ { F } \sum _ { u = 1 } ^ { U } \, \mathbb { A } _ { + } e ^ { \frac { \Delta t } { \tau _ { k } } }, \text{ if } \Delta t > 0 \\ \sum _ { f = 1 } ^ { F } \sum _ { u = 1 } ^ { U } \, \mathbb { A } _ { - } e ^ { \frac { \Delta t } { \tau _ { k } } }, \text{ else}. \end{cases} \quad ( 1 8 ) \quad C. \ T h e$$

where, the magnitudes A + and A -depend on the current synaptic weight. Moreover, f = { 1, 2, 3,... F } , u = { 1, 2, 3,... L } represent the labels of pre-synaptic and post-synaptic spike times. Furthermore, τ k determine the RC time constant of the k th DER. The initial values of R in τ can be obtained and updated using the optimal solution of (9). In this way, the relative firing time difference ∆t is exploited with respect to each DER time constant for inferring the global information for steering energy-efficient data-driven control of microgrids. For incremental deployment of the Hebbian learning policy into microgrids using (18), we modify the primary controller droop using:

$$V _ { \text{ref} } = \mathbf V ( t ) - \mathbf I ( t ) [ \mathbf R _ { \text{d} } - \mathbf a \Delta \mathbf w ( t ) ] \quad \quad ( 1 9 ) \ \text{ using}$$

where, a is a system matrix of a normalization variable. Using (18)-(19), we steer an on-device grid-edge controller for each DER, that drives using the nodal and network equations in (16)-(17) to update their respective weights using the spike timing dependent plasticity (STDP) [18]. Its system-wide

## C. Theoretical Guarantees

Due to the simple nature of Hebbian learning, only the firing time differences between the preand post-synaptic activity may not provide a substantial justification on how the learning process leads to be effective [29]. However, we exploit the microgrid dynamics and the statistical properties of circuit theory in the considered system in correlation with the neuronal circuits. Since Hebbian plasticity exploits these statistical features and domain knowledge in a certain way, it can be categorized as unsupervised learning.

Let us consider a single rate LIF neuron input model in (9), where the synaptic current I(t) is the output y , with the membrane potential V mem (t) of multiple neurons acting as inputs x (t) 1 , ..., x N (t) . Using (9), we can express the output using a response function f , given by:

$$y = f ( \sum _ { i = 1 } ^ { N } w _ { i } x _ { i } ).$$

Since we ascribe Hebbian plasticity in terms of the evolution of the synaptic weight updates, we obtain its rate of change

using (19):

$$\frac { d w _ { i } } { d t } = \eta x _ { i } y = \eta x _ { i } \sum _ { j = 1 } ^ { N } w _ { j } x _ { j }. \quad \quad ( 2 1 ) \quad \text{variat} \quad \text{cases}$$

Considering the vector form of (20) and training over an average ⟨ ... ⟩ over discrete training set, we now use the analogy chart in Fig. 2(d) and Remark V to obtain ∆ R using the weight update policy:

$$\frac { \mathrm d \mathbf w } { \mathrm d t } = \langle \eta \mathbf V \mathbf V ^ { T } \mathbf Y \rangle = \eta \langle \mathbf V \mathbf V ^ { T } \rangle \mathbf Y = \eta \mathbf C \mathbf Y. \quad \ \ ( 2 2 ) \, \begin{array} { c } \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Psi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi\Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phis \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi$$

where, Y is the system admittance matrix with tie-line resistances r kl . Furthermore, C = ⟨ VV T ⟩ denote the correlation matrix. As a result, a system with N DERs can be simply modeled as N coupled linear differential equations. Apart from the stable network models such as Oja's rule [30] or generalized Hebbian algorithm [31], for t →∞ , the solution for synaptic weights can be given by:

$$\mathbf w ( t ) \approx e ^ { \eta \alpha ^ { * } t } \mathbf c ^ { * } \quad \ \ \ \ ( 2 3 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where, α ∗ and c ∗ denote the largest eigenvalue and eigenvector of C , that corresponds to computing the first principal component of the inputs. It should be noted that η can act as a normalization parameter using the total capacity of the system P t . Finally, basic circuit theory performed on the considered system in (22) provides theoretical guarantees of a stable solution in (23), which can be extended to perform a full principal component analysis (PCA) by adding more postsynaptic neurons. In this way, the statistical aspects of the input spikes and the convergence of the output signals using Spike Talk can be justified as a networked unsupervised plasticity-based learning with a significant energy-efficient and ultra-low latency footprint.

## V. PERFORMANCE EVALUATION

Fig. 7. Case I - Performance validation of Remark IV on our hypothesis behind rate and latency coding for processing of spikes corresponding to voltage regulation and current sharing, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00773v2-figure-5b4ef7c1-25c3-4f13-82c5-19f3459c3419.jpg)

In this section, we analyze and evaluate the performance of Spike Talk in the system shown in Fig. 3(b) comprising of N = 4 DERs. The neuro-edge computing tool designed using (19) is integrated into each DER to analyze performance variations under different neural coding schemes. Four test cases were considered, where the modeling of the system in 3(b), generation of spikes using (11)-(14), Hebbian learning using (16)-(19) and the entire framework was carried out in MATLAB/Simulink environment. We consider global voltage regulation as well as proportionate current sharing as the secondary controller objectives for each DER. With these objectives federating the spike generation as well as weight updates, we execute the online learning framework in Fig. 6. The system and control parameters can be found in Appendix.

In Case I, we investigate the performance of Spike Talk in Fig. 7 to test our hypothesis on end-to-end learning using spikes derived from latency coding, whereas edge learning from rate coding. Exploiting the modeling analogy between neuronal architecture and DERs in Fig. 2, we leverage rate coding scheme in (11) with physics-informed firing and reset thresholds in (16). As shown in Fig. 7, the rate coded spikes achieve global voltage regulation with an average around 315 V. These spikes are commonly referred to as the voltage regulation spikes in Fig. 6(c). Since the dynamic droop is accountable for both the objectives, we then formulate proportionate current sharing using (13) as the time-trigger for each DER currents computed locally. Using (13), we exploit the latency of the second spike from the first spike, that basically corresponds to the calculation of the oscillatory frequency for a given disturbance due to heterogeneous control parameters. Hence, we incorporate a pseudo-threshold for each DER, V ′ 2ref and V ′ 3ref in (17), which assists in determining the relative performance of the neighboring DERs and their response times. Since heterogeneous droop parameters imparts different overshoot and rise time to each DER, we curate latency coding to determine these performance metrics. This argument can be clearly verified with the latency coding scheme in Fig. 6(b) where different peaks resulted into different temporal positions of the second spikes. It is worth notifying that the first spike is always triggered to any disturbance, which then translates as t 0 for each event. Finally using the voltage regulation spikes and current sharing spikes , we compute the weight adaptation ∆w k using (18). After customizing ∆ w into ∆ R using (19), we plot the dynamic droop in Fig. 7 that adapts itself based on the disturbance. By visualizing the synaptic conductance of each DER using (5), there is a steady-state difference between them, which is zoomed in Fig. 7 that not only validates our hypothesis in Remark IV, but also validates the modeling analogy between neurons and DER. The steady-state difference between the synaptic conductance can be accounted to the resistive drop in the tie-lines, which consequently demands slightly more current than the overall demand.

Since most event-driven techniques are primarily designed based on the output disturbances for each DER due to the output RC circuit in Fig. 2(b), we introduce a transient disturbance in the input of DER 1 at t = 2 sec in Fig. 7 to determine a correlation factor µ of the system response to input events. Defining the input and output correlation rates ρ i and ρ o respectively as the average rate of spikes over a

Fig. 8. (a) Sensitivity of spikes generation in the output of a DER due to disturbances from the input comprising of the intermittent renewable energy sources - a higher correlation factor µ is adjudged to be providing the same spike density in the output, (b) Correlation detection to determine the correlation of output spikes between DER 1 and 2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00773v2-figure-a228ab00-9d87-41b8-93c1-213ef47656c6.jpg)

given time interval, the correlation factor between the input and output of k th DER can be split using:

$$\kappa ^ { - 1 } \dot { V } _ { k } [ t _ { h } ] = \mu \rho _ { i } I _ { i n } [ t _ { h } ] - \rho _ { o } I _ { o } [ t _ { h } ]. \quad \quad ( 2 4 ) \quad \ \Sigma _ { \text{Sin} }$$

behind computing these variables using the variability and covariability of membrane potentials and spike trains will be covered in the future scope of this work.

With t h being the firing instant in (24), we leverage the voltage dynamics of each DER, which is closely related to the input as well as output current of a DER to approximate the corelationship of spikes. By mapping different correlation rates against each other, we can demonstrate almost a quasi-linear relationship between both ends in Fig. 8(a). It can be seen that the output correlation rate diverges to the lower end away from the linear plot for ρ i within [0.4, 0.7]. This is attributed to the fact that the spike generation gets compromised by the discontinuous mode of operation, which is not aligned with the assigned firing threshold and vary based on the loading condition.

Although online weight adaptation in each DER does not rely on estimating the currents from remote sources, it is still vital to assess the proportion of currents from remote DERs. This aspect can be simply answered by another neuroscience phenomena, known as coincidence detection [32]. This basically means that the spike train correlations between multiple neurons can be quantified by statistical analysis of the spike trains and their correlation coefficients. For example, I 1 in Fig. 8(b) is composed of different proportions of I 2 , I 3 and I 4 , which is impossible to predict without having sufficient global information. As a result, c in Fig. 8(b) is unknown. Since (16) anyway require the regulation coefficients to ensure global voltage regulation within a specified limit, these regulation co-efficients for DER 1 can be determined using:

$$I _ { 1 } = \mu _ { 1 } + \sigma _ { 1 } \underbrace { [ ( 1 - c ) \Psi _ { k } ( t ) + c \Psi _ { c o } ( t ) ] } _ { \psi _ { 1 k } } \quad \quad ( 2 5 ) \quad \text{low la} \quad \text{schem} \quad \text{to inst}$$

where, µ 1 is the temporal average of current I 1 with the next term representing the variance due to Gaussian fluctuations, that comprises of two weighted factors: Ψ k which is independent for the remote k th DERs, and Ψ co which is common for all the DERs. Furthermore, the input correlation coefficient (0 ≤ c ≤ 1) set the relative weight of the proportion of currents from remote DERs, whereas σ 1 is the normalized variance of the total input current. These parameters have been normalized in such a manner that the sum of regulation coefficients naturally is 1, owing to the common variance term in (25) for current from each DER. A detailed proof and analysis

Since coincidence detection is highly dependent on the rate of spikes, we illustrate a case study on correlation-rate relationship in Fig. 8(b). In the case study, a high probability of spikes is seen with an increase in the geometric mean rate of spikes, given by √ v v 1 k , where v k suggest the average rate of output spikes from k th DER. With the correlation ranging from 0 to 0.3, many perturbation techniques [32] can be used to invoke the rate of spikes from each DER without compromising system performance.

In Case II, we swap the coding scheme settings previously used in Case I to evaluate the considered system's performance in Fig. 9. As a result, we use latency coding for global voltage regulation and rate coding for equal current sharing. As it can be seen in Fig. 9, although the coding schemes in itself perform normally, their decoding into the dynamic droop ∆ R using (18) exhibit an oscillatory component, which potentially makes the system operation unstable. This is quite evident from the voltage and current responses for all the DERs, which affirms our hypothesis of the coding assignment for different secondary control objectives in Remark IV. Upon plotting a pseudospectrum estimate of the synaptic conductance g k in Fig. 9, we deduce the strength of the oscillating frequency emanating from relatively mistimed spike positions. The mistimed notion of spikes can be interpreted as ∆ (t + R τ d ) , where τ d represent the synaptic delay that can range from 50-150 ms [33]. As a result of this delay, the interactions between voltage control loops at a system level increases, that results into an unstable condition. Since it directly affects the ultralow latency capability, the suitability of the assigned coding schemes is affected for both control objectives in Remark IV to instill a satisfactory system performance.

As elaborated in the last section, burst coding relies on high frequency events using (14), which inherits an incremental pattern in the signal magnitude in either positive or negative direction over consecutive time intervals. As we deploy each DER with a neuro-edge computational tool modeled as a decentralized secondary controller, high frequency events in power systems do not cater to a response from a secondary controller, which has a typical settling time response in the order of seconds [19]. As a result, the applicability of burst coding is quite limited despite its biological accuracy and reliability of information transfer. However due to its multi-

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00773v2-figure-78f06162-b9ae-44da-b47c-9868cfe54939.jpg)

Fig. 9. Case II - Swapping the neural coding schemes to evaluate the system performance results into instability. The pseudospectrum estimate of the synaptic conductance g k for each DER depicts the strength of the oscillating frequency, which is a byproduct of mistimed temporal spike patterns.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00773v2-figure-70feb280-bd16-48f9-adcd-43818cd3c55a.jpg)

considered system to evaluate system performance in Fig. 10. As it can be seen in Fig. 10, burst coding instills a satisfactory current sharing performance for an increase as well as decrease in load in Case III. It can be argued using the Hebbian plasticity rule in Remark V for DER 1 that an increase in system demand causes an excitatory response with the error excitation in the positive region, which is evident from Fig. 10. Similarly for decrease in load, an inhibitory response can be observed with the voltage error excitation in the negative region at all times. Furthermore, the voltage error is always bounded by the trajectories of ˙ V 1 using (11).

Fig. 10. Case III - Performance evaluation of burst coding in Spike Talk using (14). Although burst coding has a higher degree of biological plausibility, it always reflects on high frequency events, that does not necessarily cater to all the disturbances. Hence, its applicability is limited to high frequency and short inter-spike interval (ISI) events.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00773v2-figure-25395c24-de60-4b6e-b680-e52efab2912d.jpg)

Fig. 11. Case IV - Plug-and-play capability using Spike Talk , when DER 3 is plugged out of the system at t = 4 sec. Since the network topology is modified with outage of DER 3, no spikes are recorded since the neuro-edge computational unit at DER 3 switches to an inert stage with V 3 and I 3 forced to zero. From a system-level computational perspective, the plugged out DER clears away itself from the processing of the output.

fold advantages with sparse data and ingrained energy-efficient computational features, we programmed burst coding into the

In Case IV, the plug-and-play capability of Spike Talk is investigated in Fig. 11. It can be seen that when DER 3 is plugged out of the considered system at t = 4 sec, the rest of the DERs increase their share of current sharing among each other. For a decrease in load at t = 5 sec, the current sharing is still intact with DER 3 as the inert neuron. Using our previous argument in Remark III, it can then be deduced that the network level computational overhead is significantly decreased due to the modeling analogy in Fig. 3(d). It is shown in Fig. 11 that outage of DER 3 is efficiently handled by the Spike Talk architecture, which automatically excludes V w 3 3 from (9) by initializing I 3 = 0. From a computation perspective on a system level, restricted power flow leading to absence of spikes from DER 3 also implies a dropout of computation to be carried out in DER 1. In this way, we ensure a highly accurate, sparse as well as low computational overhead for coordination of DERs only using digital spikes emanating out of real-valued local analog measurements.

## VI. CONCLUSIONS AND FUTURE SCOPE OF WORK

We unravel a novel data-driven co-transfer technology Spike Talk that exploits modality of power and information for sparse data-driven coordinated control of microgrids. This co-transfer principle is inspired from computational neuroscience, where we decipher a modeling analogy between the elements of a brain and power systems. Using this structural complementary mapping, we design a decentralized secondary controller for each DER, that is federated online by power flow in the network. As a result, this completely dismisses the cyber layer and breaks the traditional hierarchy of cyber-physical microgrids to wipe out all the cyber reliability issues and

vulnerabilities. In simple words, the dynamic droop parameter in the primary controller of each DER is adapted by Hebbian unsupervised learning based plasticity rules, where localized input data is selectively collected and processed using the neural coding schemes. These coding schemes introduce a sparse data analysis stage that converts synchronous realvalued measurements into asynchronous digital spikes. With increased investments in re-conductoring and smart conductors [34], Spike Talk will be an enormous grid-edge asset that unlocks power grid's potential in the field of energy efficiency, digitalization and electrification without any compromise in system reliability and performance. As we establish the selection of neuronal weights and thresholds based on each DER's dynamics and maximum capacity respectively, it becomes fairly simple to design and can be incrementally integrated into the controller of each DER.

As a future scope of work, the remaining operational elements of Spike Talk , such as stability, multiplexing of information along tie-lines will be revealed through extensive validation and theoretical guarantees.

## APPENDIX

It should be noted that the system and control parameters are consistent for each DER, unless stated otherwise.

- V ref = 315 V, C 1 = 450 µ F, C 2 = 500 µ F, C 3 = 480 µ F, C 4 = 520 µ F, r a = 0.5 Ω r , b = 0.25 Ω r , c = 0.6 Ω r , d = 0.8

$$\begin{array} { l l$$

## REFERENCES

- [1] S. Sahoo, T. Dragiˇ cevi´, and F. Blaabjerg, 'Cyber security in control of c grid-tied power electronic converters-challenges and vulnerabilities,' IEEE Journal of Emerging and Selected Topics in Power Electronics , vol. 9, no. 5, pp. 5326-5340, 2021.
- [2] E. Espina, R. C´ ardenas-Dobson, J. W. Simpson-Porco, M. Kazerani, and D. S´ aez, 'A Consensus-Based Distributed Secondary Control Optimization Strategy for Hybrid Microgrids,' IEEE Trans. Smart Grid , vol. 14, no. 6, pp. 4242-4255, Nov. 2023.
- [3] A. Platzer, Logical foundations of cyber-physical systems . Springer, 2018, vol. 662.
- [4] P. Popovski, A. M. Fouladgar, and O. Simeone, 'Interactive joint transfer of energy and information,' IEEE Transactions on Communications , vol. 61, no. 5, pp. 2086-2097, 2013.
- [5] E. Hutzell et al. , 'Down the lines: Us army communications in europe, 1942-45,' Ph.D. dissertation, 2017.
- [6] M. Yigit, V. C. Gungor, G. Tuna, M. Rangoussi, and E. Fadel, 'Power line communication technologies for smart grid applications: A review of advances and challenges,' Computer Networks , vol. 70, pp. 366-383, 2014.
- [7] M. Angjelichinoski, C. Stefanovi´ c, P. Popovski, H. Liu, P. C. Loh, and ˇ F. Blaabjerg, 'Multiuser Communication Through Power Talk in DC Microgrids,' IEEE J. Sel. Areas Commun. , vol. 34, no. 7, pp. 20062021, Jul. 2016.
- [8] M. Angjelichinoski, C. Stefanovi´ c, P. Popovski, and F. Blaabjerg, 'Power ˇ talk in dc micro grids: Constellation design and error probability performance,' in 2015 IEEE International Conference on Smart Grid Communications (SmartGridComm) . IEEE, 2015, pp. 689-694.
- [9] S. Monch, C. Kuring, X. Geng, P. A. Hoeher, M. Liserre, R. Quay, and S. Dieckerhoff, 'Efficient talkative power conversion with quasi-squarewave zero-voltage switching hysteretic current control,' IEEE Trans. Power Electron. , 2023.
- [10] Y. Liu, R. Zhu, Y. Leng, and M. Liserre, 'Transient performances analysis of talkative dual-active bridge converter modulated by pwmfsk,' in IECON 2023-49th Annual Conference of the IEEE Industrial Electronics Society . IEEE, 2023, pp. 1-6.
- [11] M. Liserre, H. Beiranvand, Y. Leng, R. Zhu, and P. A. Hoeher, 'Overview of talkative power conversion technologies,' IEEE Open Journ. Power Electron. , vol. 4, pp. 67-80, 2023.
- [12] P. Popovski, Wireless Connectivity: An Intuitive and Fundamental Guide . John Wiley &amp; Sons, 2020.
- [13] K. Gupta, S. Sahoo, and B. K. Panigrahi, 'Delay-aware semantic sampling in power electronic systems,' IEEE Trans. Smart Grid , 2023, early access.
- [14] K. Gupta, S. Sahoo, and B. K. Panigrahi, 'A monolithic cybersecurity architecture for power electronic systems,' IEEE Trans. Smart Grid , 2024, early access.
- [15] Y. Song and S. Sahoo, 'Spike talka sustainable 'one-solution-fitsall' technology for electronic grids,' IEEE Power Electronics Magazine , vol. 11, no. 4, pp. 31-38, 2024.
- [16] Y. Song and S. Sahoo, 'Spike talk in power electronic grids - leveraging post moore's computing laws,' in 2024 IEEE Design Methodologies Conference (DMC) , 2024, pp. 1-6.
- [17] X. Diao, Y. Song, and S. Sahoo, 'Inferring ingrained remote information in ac power flows using neuromorphic modality regime,' in 2024 IEEE International Conference on Communications, Control, and Computing Technologies for Smart Grids (SmartGridComm) , 2024, pp. 86-91.
- [18] X. Diao, Y. Song, S. Sahoo, and Y. Li, 'Neuromorphic event-driven semantic communication in microgrids,' IEEE Transactions on Smart Grid , pp. 1-1, 2024.
- [19] S. Sahoo and S. Mishra, 'A Distributed Finite-Time Secondary Average Voltage Regulation and Current Sharing Controller for DC Microgrids,' IEEE Trans. Smart Grid , vol. 10, no. 1, pp. 282-292, Jan. 2019.
- [20] J. K. Eshraghian, M. Ward, E. O. Neftci, X. Wang, G. Lenz, G. Dwivedi, M. Bennamoun, D. S. Jeong, and W. D. Lu, 'Training spiking neural networks using lessons from deep learning,' Proc. of the IEEE , 2023.
- [21] W. Gerstner, 'A framework for spiking neuron models: The spike response model,' in Handbook of biological physics . Elsevier, 2001, vol. 4, pp. 469-516.
- [22] P. Dayan and L. F. Abbott, Theoretical neuroscience: computational and mathematical modeling of neural systems . MIT press, 2005.
- [23] J. Zhao and F. D¨ orfler, 'Distributed control and optimization in dc microgrids,' Automatica , vol. 61, pp. 18-26, 2015.
- [24] A. R. Bergen, Power systems analysis . Pearson Education India, 2009.
- [25] B. A. Richards, T. P. Lillicrap, P. Beaudoin, Y. Bengio, R. Bogacz, A. Christensen, C. Clopath, R. P. Costa, A. de Berker, S. Ganguli et al. , 'A deep learning framework for neuroscience,' Nat. Neurosc. , vol. 22, no. 11, pp. 1761-1770, 2019.
- [26] S. Sahoo and S. Mishra, 'A Distributed Finite-Time Secondary Average Voltage Regulation and Current Sharing Controller for DC Microgrids,' IEEE Trans. Smart Grid , vol. 10, no. 1, pp. 282-292, Jan. 2019.
- [27] U. Vuyyuru, S. Maiti, C. Chakraborty, and E. I. Batzelis, 'Universal active power control converter for dc-microgrids with common energy storage,' IEEE Open Journ. Indus. Appl. , vol. 2, pp. 21-35, 2021.
- [28] G. Buticchi, M. Andresen, M. Wutti, and M. Liserre, 'Lifetime-based power routing of a quadruple active bridge dc/dc converter,' IEEE Trans. Power Electron. , vol. 32, no. 11, pp. 8892-8903, 2017.
- [29] E. O. Neftci, C. Augustine, S. Paul, and G. Detorakis, 'Event-driven random back-propagation: Enabling neuromorphic deep learning machines,' Front. Neurosc. , vol. 11, p. 324, 2017.
- [30] M. V. Jankovic and H. Ogawa, 'Modulated hebb-oja learning rulea method for principal subspace analysis,' IEEE Trans. Neural Net. , vol. 17, no. 2, pp. 345-356, 2006.
- [31] J. Kaiser, H. Mostafa, and E. Neftci, 'Synaptic plasticity dynamics for deep continuous local learning (decolle),' Front. Neurosc. , vol. 14, p. 424, 2020.
- [32] J. De La Rocha, B. Doiron, E. Shea-Brown, K. Josi´ c, and A. Reyes, 'Correlation between neural spike trains increases with firing rate,' Nature , vol. 448, no. 7155, pp. 802-806, 2007.
- [33] T. P. Lillicrap, A. Santoro, L. Marris, C. J. Akerman, and G. Hinton, 'Backpropagation and the brain,' Nat. Rev. Neurosc. , vol. 21, no. 6, pp. 335-346, 2020.
- [34] S. Marcacci, 'Reconductoring could help solve america's looming grid crisis,' Forbes, April 2024. [Online]. Available: https://www.forbes.com/sites/energyinnovation/2024/04/09/ reconductoring-could-help-solve-americas-looming-grid-problems/
Subham Sahoo
Spike Talk: Genesis and Neural Coding Scheme Translations
Although digitalization of future power grids offer several coordination
incentives, the reliability and security of information and communication
technologies (ICT) hinders its overall performance. In this paper, we introduce
a novel architecture Spike Talk via a unified representation of power and
information as a means of data normalization using spikes for coordinated
control of microgrids. This grid-edge technology allows each distributed energy
resource (DER) to execute decentralized secondary control philosophy
independently by interacting among each other using power flow along the
tie-lines. Inspired from the field of computational neuroscience, Spike Talk
basically builds on a fine-grained parallelism on the information transfer
theory in our brains, particularly when neurons (modeled as DERs) transmit
information (inferred from power streams measurable at each DER) through
synapses (modeled as tie-lines). Not only does Spike Talk simplify and address
the current bottlenecks of the cyber-physical architectural operation by
dismissing the ICT layer, it provides intrinsic operational and cost-effective
opportunities in terms of infrastructure development, computations and
modeling. Hence, this paper provides a pedagogic illustration of the key
concepts and design theories. Since we focus on coordinated control of
microgrids in this paper, the signaling accuracy and system performance is
studied for several neural coding schemes responsible for converting the
real-valued local measurements into spikes.