Petri Net Description of Biological Neural Circuits for Fast Hardware Prototyping
Authors: Carlo daCunha, Rodrigo Pena, Marcos Turqueti
Organizations: Helen and John C. Hartmann Department of Electrical Engineering, New Jersey Institute of Technology, 323 Martin Luther King Jr. Blvd., Newark, 07102, NJ, USA. · School of Applied Engineering and Technology, New Jersey Institute of Technology, 323 Martin Luther King Jr. Blvd., Newark, 07102, NJ, USA. · Department of Biological Sciences, Florida Atlantic University, 5353 Parkside Drive, Jupiter, 33458, FL, USA. · Electronic Systems Group, Lawrence Berkeley National Laboratory, 1 Cyclotron Road., Berkeley, 94720, CA, USA.
Current approaches to simulating biological neural circuits, whether on general-purpose hardware or dedicated neuromorphic platforms, remain constrained by fixed-timestep numerical integration, hardware-imposed precision limits, and an inability to guarantee timing correctness for event-driven spiking dynamics under real-time constraints. Here, we propose a Petri net description of biological neural circuits that overcomes these limitations by modeling neurons, synapses, and spike events as a T-timed Petri net with formally verifiable timing semantics, enabling deadline-guaranteed real-time execution and analytically tractable correspondence to continuous-time leak-integrate-and-fire dynamics, independent of the underlying integration timestep. To test the model, we present the results of three simulated microcircuits: feedback inhibition, lateral inhibition, and hierarchical feature detector. The Petri neuron reproduces the expected dynamical signatures of each circuit while providing formally bounded timing guarantees throughout, with worst-case response times matching analytical predictions across all three cases.
Figures & tables
Figure 1: Phases of a biological action potential. Depolarization (rapid voltage rise past threshold), repolarization (voltage returns toward rest), hyperpolarization (transient overshoot below resting potential), and the subsequent return to rest via ionic pump activity. The refractory period spans repolarization and hyperpolarization, during which the neuron is unresponsive (absolute) or requires elevated excitation to fire (relative). This cycle motivates the input-gating, threshold-firing, and recovery transitions of the Petri neuron in Fig. 2 .
Figure 2: Petri net representation of the Petri neuron, mapping accumulation, threshold firing, axonal propagation, and refractory recovery onto places and transitions.
Tinput
Tspike
Tprop
Trec
Tleak
pacc
1/0/+1
0/θ/−θ
0/0/0
0/0/0
0/1/−1
prdy
0/0/0
0/1/−1
0/0/0
1/0/+1
0/0/0
ppre
0/0/0
1/0/+1
0/1/−1
0/0/0
0/0/0
pout
0/0/0
0/0/0
1/0/+1
0/0/0
0/0/0
prec
0/0/0
1/0/+1
0/0/0
0/1/−1
0/0/0
Table 1: Incidence matrix with entries C+/C−/C . θ denotes the current token count of pacc at firing time (flush arc).
Figure 3: Firing rate as a function of input excitation I/Ith for the LIF model (Eq. 3 ) and the Petri neuron (Eq. 5 ), parameterized via the Padé mapping. Agreement is close for spontaneous and moderate firing ( f<50 Hz), with relative error shrinking further as I increases into the high-frequency regime.
Machine
Nom. (ms)
Mean
Std
p99
ϵ (max)
AMD
1
9.7
2.3
16.9
95.9
5
171.4
221.6
717.2
901.9
10
207.4
227.7
728.2
917.7
50
198.6
234.7
728.6
885.0
100
252.0
240.8
730.0
905.3
RP2040
1
0.03
0.17
1.00
1.00
Table 2: AMD Linux POSIX timer and RP2040 time jitter. All values in μ s. ϵ denotes the worst-case absolute error ( N=10,000 ).
Figure 4: Feedback inhibition network. Chevron indicates an excitatory synapse, while the small black disks indicate inhibitory synapses.
Parameter
E
I
R [M Ω ]
30–60
337
C [ ×10−10 F]
1.3–4
0.44–0.74
τ [ms]
8–12
15–25
L [mm]
1–3
4–8
Table 3: Neuronal parameters of excitatory (E) and inhibitory (I) neurons in a typical feedback-inhibition network.
Figure 5: Inter-event interval histogram for the firing of E and I neurons in a feedback inhibition network simulated using Petri neuron.
Figure 6: Lateral inhibition microcircuit motif with excitatory neurons EA,B and inhibitory interneurons IA,B .
Figure 7: IEIH of excitatory neurons and their firing rates as a function of time.
Figure 8: The hierarchical feature detector consisting of a 2×2 grid, simple neurons S that capture lines, and complex neurons C that indicate the presence of horizontal and vertical lines.
Spiking neural networks process information through discrete, time-dependent neural events, aiming for efficient computation and biologically inspired learning.
We study the expressivity of spiking neural networks, which provide a natural framework for asynchronous, event-driven computation complementary to conventional feedforward neural networks. We consider the time-to-first-spike model in a setting for which the input-output map is continuous and piecewise linear, with affine pieces governed by causal feasibility constraints that determine which presynaptic spikes occur before a neuron fires. We first show that each neuron's firing time admits a maxout-like representation with exponentially many, highly constrained affine pieces. We then formalize causal regions as polyhedral regions with fixed causal sets and derive upper and lower bounds on the maximal number of causal regions in both shallow and multilayer feedforward spiking networks. Our theoretical and experimental results show that spiking networks can generate richer partitions of the input space than conventional feedforward ReLU networks.
Manjot Singh, Guido Montúfar, Gitta Kutyniok
Department of Mathematics, Ludwig-Maximilians-Universität München · Munich Center for Machine Learning · Departments of Mathematics and Statistics & Data Science, University of California, Los Angeles +2
Spiking neural networks (SNNs) promise low-power event-driven computation for temporally rich tasks, but commonly used neuron models often trade off gradient-based trainability, dynamical richness, and high activity sparsity. These limitations are acute in regression, where approximation error, noise and spike discretization can severely degrade continuous-valued outputs. Indeed, many state-of-the-art (SOTA) SNNs rely on simple phenomenological dynamics trained with surrogate gradients and offer limited control over spiking diversity and sparsity. To overcome such limitations, we introduce multi-timescale conductance spiking networks, a gradient-trainable framework in which neural dynamics emerge from shaping the current-voltage (I-V) curve by tuning fast, slow and ultra-slow conductances. This parametrization allows systematic control over excitability, can be implemented efficiently in analog circuits, and yields rich firing regimes including tonic, phasic and bursting responses within a single model. We derive a discrete-time formulation of these differentiable dynamics, enabling direct backpropagation through time without surrogate-gradient approximations. To probe both trainability and accuracy, we evaluate feedforward networks of these neurons at the predictability limit of Mackey-Glass time-series regression and compare them to baseline LIF and SOTA AdLIF networks. Our model outperforms LIF and AdLIF networks, while exhibiting substantially sparser activity from both communication and computational perspectives. These results highlight multi-timescale conductance spiking neurons as a promising building block for energy-aware temporal processing and neuromorphic implementation.
Alex Fulleda-Garcia, Saray Soldado-Magraner, Josep Maria Margarit-Taulé
Instituto de Microelectrónica de Barcelona (IMB-CNM), Consejo Superior de Investigaciones Científicas (CSIC) Cerdanyola del Vallès, Spain · Department of Neurobiology, University of California Los Angeles (UCLA) Los Angeles, USA · Joint Senior authors
Short-term plasticity (STP) is fundamental to temporal information processing in biological neural systems but remains difficult to realize efficiently in neuromorphic hardware. Memristive electrochemical random-access memory (ECRAM) devices naturally exhibit non-equilibrium ionic dynamics that produce transient conductance modulation; however, these behaviors are typically treated as undesirable variability or tolerated as side effects in memory-centric computing paradigms. In this work, we instead transform these volatile dynamics from a tolerated device artifact into a computational resource through a cross-layer device-circuit-system co-design framework. We introduce a delay-feedback leaky integrate-and-fire (LIF) neuron architecture co-designed with ECRAM synapses that exploits activity-dependent conductance modulation with negligible additional circuit overhead. The architecture integrates ECRAM-based synapses with a tunable delay-feedback spike-generation path, enabling transient device dynamics to directly modulate neuron excitability and synaptic efficacy. We used experimentally characterized ECRAM devices exhibiting transient conductance modulation (1.5 KOhms per spike) to develop a compact behavioral model suitable for circuit-level simulation. Circuit simulations demonstrate two key STP behaviors -- synaptic facilitation and intrinsic excitability modulation -- while consuming 2 pJ per spike, and the same device-driven mechanisms extend across multiple neuron topologies. Network-level analysis further demonstrates frequency-selective spike processing, allowing individual synapses to act as tunable temporal filters within spiking neural networks. This work demonstrates that non-equilibrium ECRAM dynamics can serve as a native hardware substrate for STP and temporal computation in neuromorphic circuits.
Alex Currie, Sean Borkholder, Nithil Harris Manimaran +4
RAMLab, Department of Electrical and Microelectronic Engineering, Rochester Institute of Technology, Rochester NY, USA · Iontronic and Nanoelectronics Lab, School of Physics and Astronomy, Rochester Institute of Technology, Rochester NY, USA · Brain Lab, Department of Computer Engineering, Rochester Institute of Technology, Rochester NY, USA