Network Oscillations Research Articles

Synaptic plasticity: from chimera states to synchronicity oscillations in multilayer neural networks

Synaptic plasticity: from chimera states to synchronicity oscillations in multilayer neural networks

Propagation of transient explosive synchronization in a mesoscale mouse brain network model of epilepsy.

Generalized epileptic attacks, which exhibit widespread disruption of brain activity, are characterized by recurrent, spontaneous, and synchronized bursts of neural activity that self-initiate and self-terminate through critical transitions. Here we utilize the general framework of explosive synchronization (ES) from complex systems science to study the role of network structure and resource dynamics in the generation and propagation of seizures. We show that a combination of resource constraint and adaptive coupling in a Kuramoto network oscillator model can reliably generate seizure-like synchronization activity across different network topologies, including a biologically derived mesoscale mouse brain network. The model, coupled with a novel algorithm for tracking seizure propagation, provides mechanistic insight into the dynamics of transition to the synchronized state and its dependence on resources; and identifies key brain areas that may be involved in the initiation and spatial propagation of the seizure. The model, though minimal, efficiently recapitulates several experimental and theoretical predictions from more complex models and makes novel experimentally testable predictions.

Bistability and the emergence of oscillation in a multiple-loop traffic network

Bistability and the emergence of oscillation in a multiple-loop traffic network

Investigating the Effects of Transcranial Alternating Current Stimulation on Cortical Oscillations and Network Dynamics.

Transcranial electrical brain stimulation techniques like transcranial direct current (tDCS) and transcranial alternating current (tACS) have emerged as potential tools for treating neurological diseases by modulating cortical excitability. These techniques deliver small electric currents to the brain non-invasively through electrodes on the scalp. tDCS uses constant direct current which weakly alters the membrane voltage of cortical neurons, while tACS utilizes alternating current to target and enhance cortical oscillations, though the underlying mechanisms are not fully understood more specifically. To elucidate how tACS perturbs endogenous network dynamics, we simulated spiking neuron network models. We identified distinct roles of the depolarizing and hyperpolarizing phases in driving network activity towards and away from the strong nonlinearity provided by pyramidal neurons. Exploring resonance effects, we found matching tACS frequency to the network's endogenous resonance frequency creates greater entrainment. Based on this, we developed an algorithm to determine the network's endogenous frequency, phase, and amplitude, then deliver optimized tACS to entrain network oscillations. Together, these computational results provide mechanistic insight into the effects of tACS on network dynamics and could inform future closed-loop tACS systems that dynamically tune stimulation parameters to ongoing brain activity.

Functional architecture of intracellular oscillations in hippocampal dendrites

Fast electrical signaling in dendrites is central to neural computations that support adaptive behaviors. Conventional techniques lack temporal and spatial resolution and the ability to track underlying membrane potential dynamics present across the complex three-dimensional dendritic arbor in vivo. Here, we perform fast two-photon imaging of dendritic and somatic membrane potential dynamics in single pyramidal cells in the CA1 region of the mouse hippocampus during awake behavior. We study the dynamics of subthreshold membrane potential and suprathreshold dendritic events throughout the dendritic arbor in vivo by combining voltage imaging with simultaneous local field potential recording, post hoc morphological reconstruction, and a spatial navigation task. We systematically quantify the modulation of local event rates by locomotion in distinct dendritic regions, report an advancing gradient of dendritic theta phase along the basal-tuft axis, and describe a predominant hyperpolarization of the dendritic arbor during sharp-wave ripples. Finally, we find that spatial tuning of dendritic representations dynamically reorganizes following place field formation. Our data reveal how the organization of electrical signaling in dendrites maps onto the anatomy of the dendritic tree across behavior, oscillatory network, and functional cell states.

Refinement of Global Coronal and Interplanetary Magnetic Field Extrapolations Constrained by Remote-sensing and In Situ Observations at the Solar Minimum

Solar magnetic fields are closely related to various physical phenomena on the Sun, which can be extrapolated with different models from photospheric magnetograms. However, the open flux problem (OFP), the underestimation of the magnetic field derived from the extrapolated model, is still unsolved. To minimize the impact of the OFP, we propose three evaluation parameters to quantitatively evaluate magnetic field models and determine the optimal free parameters in the models by constraining the coronal magnetic fields and the interplanetary magnetic fields (IMFs) with real observations. Although the OFP still exists, we find that magnetic field lines traced from the coronal models effectively capture the intricate topological configurations observed in the corona, including streamers and plumes. The OFP is lessened by using the Helioseismic and Magnetic Imager synoptic map instead of the Global Oscillation Network Group daily synoptic maps, and the potential field source surface + potential field current sheet (PFSS+PFCS) model instead of the current sheet source surface (CSSS) model. For Carrington Rotation 2231 at the solar minimum, we suggest that the optimal parameters for the PFSS+PFCS model are R ss = 2.2–2.5 R ☉ and R scs = 10.5–14.0 R ☉, as well as for the CSSS model, are R cs = 2.0–2.4 R ☉, R ss = 11.0–14.7 R ☉, and a = 1.0 R ☉. Despite the IMFs at 1 au being consistent with the measurements by artificially increasing the polar magnetic fields, the IMFs near the Sun are still underestimated. The OFP might be advanced by improving the accuracy of both the weak magnetic fields and polar magnetic fields, especially considering magnetic activities arising from interplanetary physical processes.

Modulation of large rhythmic depolarizations in human large basket cells by norepinephrine and acetylcholine

Rhythmic brain activity is critical to many brain functions and is sensitive to neuromodulation, but so far very few studies have investigated this activity on the cellular level in vitro in human brain tissue samples. This study reveals and characterizes a novel rhythmic network activity in the human neocortex. Using intracellular patch-clamp recordings of human cortical neurons, we identify large rhythmic depolarizations (LRDs) driven by glutamate release but not by GABA. These LRDs are intricate events made up of multiple depolarizing phases, occurring at ~0.3 Hz, have large amplitudes and long decay times. Unlike human tissue, rat neocortex layers 2/3 exhibit no such activity under identical conditions. LRDs are mainly observed in a subset of L2/3 interneurons that receive substantial excitatory inputs and are likely large basket cells based on their morphology. LRDs are highly sensitive to norepinephrine (NE) and acetylcholine (ACh), two neuromodulators that affect network dynamics. NE increases LRD frequency through β-adrenergic receptor activity while ACh decreases it via M4 muscarinic receptor activation. Multi-electrode array recordings show that NE enhances and synchronizes oscillatory network activity, whereas ACh causes desynchronization. Thus, NE and ACh distinctly modulate LRDs, exerting specific control over human neocortical activity.

Mapping responses to focal injections of bicuculline in the lateral parafacial region identifies core regions for maximal generation of active expiration

The lateral parafacial area (pFL) is a crucial region involved in respiratory control, particularly in generating active expiration through an expiratory oscillatory network. Active expiration involves rhythmic abdominal (ABD) muscle contractions during late-expiration, increasing ventilation during elevated respiratory demands. The precise anatomical location of the expiratory oscillator within the ventral medulla’s rostro-caudal axis is debated. While some studies point to the caudal tip of the facial nucleus (VIIc) as the oscillator’s core, others suggest more rostral areas. Our study employed bicuculline (a γ-aminobutyric acid type A [GABA-A] receptor antagonist) injections at various pFL sites (–0.2 mm to +0.8 mm from VIIc) to investigate the impact of GABAergic disinhibition on respiration. These injections consistently elicited ABD recruitment, but the response strength varied along the rostro-caudal zone. Remarkably, the most robust and enduring changes in tidal volume, minute ventilation, and combined respiratory responses occurred at more rostral pFL locations (+0.6/+0.8 mm from VIIc). Multivariate analysis of the respiratory cycle further differentiated between locations, revealing the core site for active expiration generation with this experimental approach. Our study advances our understanding of neural mechanisms governing active expiration and emphasizes the significance of investigating the rostral pFL region.

Mapping responses to focal injections of bicuculline in the lateral parafacial region identifies core regions for maximal generation of active expiration.

The lateral parafacial area (pFL) is a crucial region involved in respiratory control, particularly in generating active expiration through an expiratory oscillatory network. Active expiration involves rhythmic abdominal (ABD) muscle contractions during late-expiration, increasing ventilation during elevated respiratory demands. The precise anatomical location of the expiratory oscillator within the ventral medulla's rostro-caudal axis is debated. While some studies point to the caudal tip of the facial nucleus (VIIc) as the oscillator's core, others suggest more rostral areas. Our study employed bicuculline (a γ-aminobutyric acid type A [GABA-A] receptor antagonist) injections at various pFL sites (-0.2 mm to +0.8 mm from VIIc) to investigate the impact of GABAergic disinhibition on respiration. These injections consistently elicited ABD recruitment, but the response strength varied along the rostro-caudal zone. Remarkably, the most robust and enduring changes in tidal volume, minute ventilation, and combined respiratory responses occurred at more rostral pFL locations (+0.6/+0.8 mm from VIIc). Multivariate analysis of the respiratory cycle further differentiated between locations, revealing the core site for active expiration generation with this experimental approach. Our study advances our understanding of neural mechanisms governing active expiration and emphasizes the significance of investigating the rostral pFL region.

SQUID oscillations in PbTe nanowire networks

SQUID oscillations in PbTe nanowire networks

Dopamine and deep brain stimulation accelerate the neural dynamics of volitional action in Parkinson's disease.

The ability to initiate volitional action is fundamental to human behaviour. Loss of dopaminergic neurons in Parkinson's disease is associated with impaired action initiation, also termed akinesia. Both dopamine and subthalamic deep brain stimulation (DBS) can alleviate akinesia, but the underlying mechanisms are unknown. An important question is whether dopamine and DBS facilitate de novo build-up of neural dynamics for motor execution or accelerate existing cortical movement initiation signals through shared modulatory circuit effects. Answering these questions can provide the foundation for new closed-loop neurotherapies with adaptive DBS, but the objectification of neural processing delays prior to performance of volitional action remains a significant challenge. To overcome this challenge, we studied readiness potentials and trained brain signal decoders on invasive neurophysiology signals in 25 DBS patients (12 female) with Parkinson's disease during performance of self-initiated movements. Combined sensorimotor cortex electrocorticography and subthalamic local field potential recordings were performed OFF therapy (n = 22), ON dopaminergic medication (n = 18) and on subthalamic deep brain stimulation (n = 8). This allowed us to compare their therapeutic effects on neural latencies between the earliest cortical representation of movement intention as decoded by linear discriminant analysis classifiers and onset of muscle activation recorded with electromyography. In the hypodopaminergic OFF state, we observed long latencies between motor intention and motor execution for readiness potentials and machine learning classifications. Both, dopamine and DBS significantly shortened these latencies, hinting towards a shared therapeutic mechanism for alleviation of akinesia. To investigate this further, we analysed directional cortico-subthalamic oscillatory communication with multivariate granger causality. Strikingly, we found that both therapies independently shifted cortico-subthalamic oscillatory information flow from antikinetic beta (13-35 Hz) to prokinetic theta (4-10 Hz) rhythms, which was correlated with latencies in motor execution. Our study reveals a shared brain network modulation pattern of dopamine and DBS that may underlie the acceleration of neural dynamics for augmentation of movement initiation in Parkinson's disease. Instead of producing or increasing preparatory brain signals, both therapies modulate oscillatory communication. These insights provide a link between the pathophysiology of akinesia and its' therapeutic alleviation with oscillatory network changes in other non-motor and motor domains, e.g. related to hyperkinesia or effort and reward perception. In the future, our study may inspire the development of clinical brain computer interfaces based on brain signal decoders to provide temporally precise support for action initiation in patients with brain disorders.

Multipoint study of the rapid filament evolution during a confined C2 flare on 28 March 2022, leading to eruption

This study focuses on the rapid evolution of the solar filament in active region 12975 during a confined C2 flare on 28 March 2022, which finally led to an eruptive M4 flare 1.5 h later. The event is characterized by the apparent breakup of the filament, the disappearance of its southern half, and the flow of the remaining filament plasma into a new, longer channel with a topology very similar to an extreme ultraviolet (EUV) hot channel observed during the flare. In addition, we outline the emergence of the original filament from a sheared arcade and discuss possible drivers for its rise and eruption. We took advantage of Solar Orbiter's favorable position, 0.33 AU from the Sun, and $83. 5^ west of the Sun-Earth line, to perform a multi-point study using the Spectrometer Telescope for Imaging X-rays (STIX) and the Extreme Ultraviolet Imager (EUI) in combination with the Atmospheric Imaging Assembly (AIA) and the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) and Halpha images from the Earth-based Kanzelhöhe Observatory for Solar and Environmental Research (KSO) and the Global Oscillation Network Group (GONG). While STIX and EUI observed the flare and the filament's rise from close up and at the limb, AIA and HMI observations provided highly complementary on-disk observations from which we derived differential emission measure (DEM) maps and nonlinear force-free (NLFF) magnetic field extrapolations. According to our pre-flare NLFF extrapolation, field lines corresponding to both filament channels existed in close proximity before the flare. We propose a loop-loop reconnection scenario based on field structures associated with the AIA 1600 Å flare ribbons and kernels. It involves field lines surrounding and passing beneath the shorter filament channel, and field lines closely following the southern part of the longer channel. Reconnection occurs in an essentially vertical current sheet at a polarity inversion line (PIL) below the breakup region, which enables the formation of the flare loop arcade and EUV hot channel. This scenario is supported by concentrated currents and free magnetic energy built up by antiparallel flows along the PIL before the flare, and by non-thermal X-ray emission observed from the reconnection region. The reconnection probably propagated to involve the original filament itself, leading to its breakup and new geometry. This reconnection geometry also provides a general mechanism for the formation of the long filament channel and realizes the concept of tether cutting. It was probably active throughout the filament's continuous rise phase, which lasted from at least 30 min before the C2 flare until the filament eruption. The C2 flare represents a period of fast reconnection during this otherwise more steady period, during which most of the original filament was reconnected and joined the longer channel. These results demonstrate how rapid changes in solar filament topology can be driven by loop-loop reconnection with nearby field structures, and how this can be part of a long-lasting tether-cutting reconnection process. They also illustrate how a confined precursor flare due to loop-loop reconnection (Type I) can contribute to the evolution towards a full eruption, and that they can produce a flare loop arcade when the contact region between interacting loop systems has a sheet-like geometry similar to a flare current sheet.

Diffusion-driven instabilities in the BT-GN oscillatory carbonylation reaction network.

This study explores the role of diffusion in creating instabilities in the Bruk Temkin-Gorodsky Novakovic (BT-GN) oscillatory carbonylation reaction network. Stoichiometric network analysis and numerical methods revealed the presence of two destabilizing feedback cycles responsible for these instabilities. Analysis of a spatially uniform system showed that the saddle-node bifurcation can be simulated within the reaction network. The introduction of diffusion results in two types of instabilities: one occurs when a spatially uniform system is already unstable, leading to a reaction-diffusion front; and another involves diffusion-driven instabilities where introducing diffusion destabilizes a stable spatially uniform system. Slower PdI2 diffusion plays a key role in inducing these instabilities. Equations describing conditions for the emergence of the instabilities in both cases were derived.

Scalable synchronization cluster in networked chaotic oscillators.

Cluster synchronization in synthetic networks of coupled chaotic oscillators is investigated. It is found that despite the asymmetric nature of the network structure, a subset of the oscillators can be synchronized as a cluster while the other oscillators remain desynchronized. Interestingly, with the increase in the coupling strength, the cluster is expanding gradually by recruiting the desynchronized oscillators one by one. This new synchronization phenomenon, which is named "scalable synchronization cluster," is explored theoretically by the method of eigenvector-based analysis, and it is revealed that the scalability of the cluster is attributed to the unique feature of the eigenvectors of the network coupling matrix. The transient dynamics of the cluster in response to random perturbations are also studied, and it is shown that in restoring to the synchronization state, oscillators inside the cluster are stabilized in sequence, illustrating again the hierarchy of the oscillators. The findings shed new light on the collective behaviors of networked chaotic oscillators and are helpful for the design of real-world networks where scalable synchronization clusters are concerned.

Dynamic assemblies of parvalbumin interneurons in brain oscillations

Brain oscillations are crucial for perception, memory, and behavior. Parvalbumin-expressing (PV) interneurons are critical for these oscillations, but their population dynamics remain unclear. Using voltage imaging, we simultaneously recorded membrane potentials in up to 26 PV interneurons invivo during hippocampal ripple oscillations in mice. We found that PV cells generate ripple-frequency rhythms by forming highly dynamic cell assemblies. These assemblies exhibit rapid and significant changes from cycle to cycle, varying greatly in both size and membership. Importantly, this variability is not just random spiking failures of individual neurons. Rather, the activities of other PV cells contain significant information about whether a PV cell spikes or not in a given cycle. This coordination persists without network oscillations, and it exists in subthreshold potentials even when the cells are not spiking. Dynamic assemblies of interneurons may provide a new mechanism to modulate postsynaptic dynamics and impact cognitive functions flexibly and rapidly.

Resonancelike emergence of chaos in complex networks of damped-driven nonlinear systems.

Characterizing the emergence of chaotic dynamics of complex networks is an essential task in nonlinear science with potential important applications in many fields such as neural control engineering, microgrid technologies, and ecological networks. Here, we solve a critical outstanding problem in this multidisciplinary research field: the emergence and persistence of spatiotemporal chaos in complex networks of damped-driven nonlinear oscillators in the significant weak-coupling regime, while they exhibit regular behavior when uncoupled. By developing a comprehensive theory with the aid of standard analytical methods, a hierarchy of lower-dimensional effective models, and extensive numerical simulations, we uncover and characterize the basic physical mechanisms concerning both heterogeneity-induced and impulse-induced emergence, enhancement, and suppression of chaos in starlike and scale-free networks of periodically driven, dissipative nonlinear oscillators.

Optimal placement of a dissimilar node for chaos suppression in networks

We demonstrate that the presence of a single dissimilar chaotic system suppresses chaos in networks of chaotic oscillators, in a diverse set of network topologies, for sufficiently strong coupling. The key property is determined to be the sum of the path lengths between the dissimilar node and all the other nodes (or its maximum, if coupled to unconnected networks), and there exists a linear relation between this quantity and the critical coupling strength for the onset of a spatiotemporal fixed point. This holds true for a chain with the dissimilar node at different locations, a ring and complete network with one embedded dissimilar node, as well as star networks with a dissimilar hub or dissimilar peripheral node. Furthermore, we show that networks with high average degree and high clustering coefficient are more resilient to the influence of an external dissimilar system. These findings will potentially aid in the design of optimally placed dissimilar nodes for controlling chaos in complex networks.

Energy-Performance Assessment of Oscillatory Neural Networks Based on VO2 Devices for Future Edge AI Computing.

Oscillatory neural network (ONN) is an emerging neuromorphic architecture composed of oscillators that implement neurons and are coupled by synapses. ONNs exhibit rich dynamics and associative properties, which can be used to solve problems in the analog domain according to the paradigm let physics compute. For example, compact oscillators made of VO2 material are good candidates for building low-power ONN architectures dedicated to AI applications at the edge, like pattern recognition. However, little is known about the ONN scalability and its performance when implemented in hardware. Before deploying ONN, it is necessary to assess its computation time, energy consumption, performance, and accuracy for a given application. Here, we consider a VO2-oscillator as an ONN building block and perform circuit-level simulations to evaluate the ONN performances at the architecture level. Notably, we investigate how the ONN computation time, energy, and memory capacity scale with the number of oscillators. It appears that the ONN energy grows linearly when scaling up the network, making it suitable for large-scale integration at the edge. Furthermore, we investigate the design knobs for minimizing the ONN energy. Assisted by technology computer-aided design (TCAD) simulations, we report on scaling down the dimensions of VO2 devices in crossbar (CB) geometry to decrease the oscillator voltage and energy. We benchmark ONN versus state-of-the-art architectures and observe that the ONN paradigm is a competitive energy-efficient solution for scaled VO2 devices oscillating above 100 MHz. Finally, we present how ONN can efficiently detect edges in images captured on low-power edge devices and compare the results with Sobel and Canny edge detectors.

Neuropeptide modulation of bidirectional internetwork synapses.

Oscillatory networks underlying rhythmic motor behaviors, and sensory and complex neural processing, are flexible, even in their neuronal composition. Neuromodulatory inputs enable neurons to switch participation between networks or participate in multiple networks simultaneously. Neuromodulation of internetwork synapses can both recruit and coordinate a switching neuron in a second network. We previously identified an example in which a neuron is recruited into dual-network activity via peptidergic modulation of intrinsic properties. We now ask whether the same neuropeptide also modulates internetwork synapses for internetwork coordination. The crab (Cancer borealis) stomatogastric nervous system contains two well-defined feeding-related networks (pyloric, food filtering, ∼1 Hz; gastric mill, food chewing, ∼0.1 Hz). The projection neuron MCN5 uses the neuropeptide Gly1-SIFamide to recruit the pyloric-only lateral posterior gastric (LPG) neuron into dual pyloric- plus gastric mill-timed bursting via modulation of LPG's intrinsic properties. Descending input is not required for a coordinated rhythm, thus intranetwork synapses between LPG and its second network must underlie coordination among these neurons. However, synapses between LPG and gastric mill neurons have not been documented. Using two-electrode voltage-clamp recordings, we found that graded synaptic currents between LPG and gastric mill neurons (lateral gastric, inferior cardiac, and dorsal gastric) were primarily negligible in saline, but were enhanced by Gly1-SIFamide. Furthermore, LPG and gastric mill neurons entrain each other during Gly1-SIFamide application, indicating bidirectional, functional connectivity. Thus, a neuropeptide mediates neuronal switching through parallel actions, modulating intrinsic properties for recruitment into a second network and as shown here, also modulating bidirectional internetwork synapses for coordination.NEW & NOTEWORTHY Neuromodulation can enable neurons to simultaneously coordinate with separate networks. Both recruitment into, and coordination with, a second network can occur via modulation of internetwork synapses. Alternatively, recruitment can occur via modulation of intrinsic ionic currents. We find that the same neuropeptide previously determined to modulate intrinsic currents also modulates bidirectional internetwork synapses that are typically ineffective. Thus, complementary modulatory peptide actions enable recruitment and coordination of a neuron into a second network.

Filament eruption by multiple reconnections

Context. Filament eruption is a common phenomenon in solar activity, but the triggering mechanism is not well understood. Aims. We focus our study on a filament eruption located in a complex nest of three active regions close to a coronal hole. Methods. The filament eruption is observed at multiple wavelengths: by the Global Oscillation Network Group (GONG), the Solar Terrestrial Relations Observatory (STEREO), the Solar Upper Transition Region Imager (SUTRI), and the Atmospheric Imaging Assembly (AIA) and Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamic Observatory (SDO). Thanks to high-temporal-resolution observations, we were able to analyze the evolution of the fine structure of the filament in detail. The filament changes direction during the eruption, which is followed by a halo coronal mass ejection detected by the Large Angle Spectrometric Coronagraph (LASCO) on board the Solar and Heliospheric Observatory (SOHO). A Type III radio burst was also registered at the time of the eruption. To investigate the process of the eruption, we analyzed the magnetic topology of the filament region adopting a nonlinear force-free-field (NLFFF) extrapolation method and the polytropic global magnetohydrodynamic (MHD) modeling. We modeled the filament by embedding a twisted flux rope with the regularized Biot-Savart Laws (RBSL) method in the ambient magnetic field. Results. The extrapolation results show that magnetic reconnection occurs in a fan-spine configuration resulting in a circular flare ribbon. The global modeling of the corona demonstrates that there was an interaction between the filament and open field lines, causing a deflection of the filament in the direction of the observed CME eruption and dimming area. Conclusions. The modeling supports the following scenario: magnetic reconnection not only occurs with the filament itself (the flux rope) but also with the background magnetic field lines and open field lines of the coronal hole located to the east of the flux rope. This multiwavelength analysis indicates that the filament undergoes multiple magnetic reconnections on small and large scales with a drifting of the flux rope.