Subthreshold firing in Mott nanodevices.

Resistive switching, a phenomenon in which the resistance of a device can be modified by applying an electric field1-5, is at the core of emerging technologies such as neuromorphic computing and resistive memories6-9. Among the different types of resistive switching, threshold firing10-14 is one of the most promising, as it may enable the implementation of artificial spiking neurons7,13,14. Threshold firing is observed in Mott insulators featuring an insulator-to-metal transition15,16, which can be triggered by applying an external voltage: the material becomes conducting ('fires') if a threshold voltage is exceeded7,10-12. The dynamics of this induced transition have been thoroughly studied, and its underlying mechanism and characteristic time are well documented10,12,17,18. By contrast, there is little knowledge regarding the opposite transition: the process by which the system returns to the insulating state after the voltage is removed. Here we show that Mott nanodevices retain a memory of previous resistive switching events long after the insulating resistance has recovered. We demonstrate that, although the device returns to its insulating state within 50 to 150 nanoseconds, it is possible to re-trigger the insulator-to-metal transition by using subthreshold voltages for a much longer time (up to several milliseconds). We find that the intrinsic metastability of first-order phase transitions is the origin of this phenomenon, and so it is potentially present in all Mott systems. This effect constitutes a new type of volatile memory in Mott-based devices, with potential applications in resistive memories, solid-state frequency discriminators and neuromorphic circuits.

A Mean Field to Capture Asynchronous Irregular Dynamics of Conductance-Based Networks of Adaptive Quadratic Integrate-and-Fire Neuron Models.

Mean-field models are a class of models used in computational neuroscience to study the behavior of large populations of neurons. These models are based on the idea of representing the activity of a large number of neurons as the average behavior of mean-field variables. This abstraction allows the study of large-scale neural dynamics in a computationally efficient and mathematically tractable manner. One of these methods, based on a semianalytical approach, has previously been applied to different types of single-neuron models, but never to models based on a quadratic form. In this work, we adapted this method to quadratic integrate-and-fire neuron models with adaptation and conductance-based synaptic interactions. We validated the mean-field model by comparing it to the spiking network model. This mean-field model should be useful to model large-scale activity based on quadratic neurons interacting with conductance-based synapses.

Operando characterization of conductive filaments during resistive switching in Mott VO2.

Vanadium dioxide (VO2) has attracted much attention owing to its metal-insulator transition near room temperature and the ability to induce volatile resistive switching, a key feature for developing novel hardware for neuromorphic computing. Despite this interest, the mechanisms for nonvolatile switching functioning as synapse in this oxide remain not understood. In this work, we use in situ transmission electron microscopy, electrical transport measurements, and numerical simulations on Au/VO2/Ge vertical devices to study the electroforming process. We have observed the formation of V5O9 conductive filaments with a pronounced metal-insulator transition and that vacancy diffusion can erase the filament, allowing for the system to "forget." Thus, both volatile and nonvolatile switching can be achieved in VO2, useful to emulate neuronal and synaptic behaviors, respectively. Our systematic operando study of the filament provides a more comprehensive understanding of resistive switching, key in the development of resistive switching-based neuromorphic computing.

EEG-fMRI in awake rat and whole-brain simulations show decreased brain responsiveness to sensory stimulations during absence seizures.

In patients suffering absence epilepsy, recurring seizures can significantly decrease their quality of life and lead to yet untreatable comorbidities. Absence seizures are characterized by spike-and-wave discharges on the electroencephalogram associated with a transient alteration of consciousness. However, it is still unknown how the brain responds to external stimuli during and outside of seizures. This study aimed to investigate responsiveness to visual and somatosensory stimulation in Genetic Absence Epilepsy Rats from Strasbourg (GAERS), a well-established rat model for absence epilepsy. Animals were imaged under non-curarized awake state using a quiet, zero echo time, functional magnetic resonance imaging (fMRI) sequence. Sensory stimulations were applied during interictal and ictal periods. Whole-brain hemodynamic responses were compared between these two states. Additionally, a mean-field simulation model was used to explain the changes of neural responsiveness to visual stimulation between states. During a seizure, whole-brain responses to both sensory stimulations were suppressed and spatially hindered. In the cortex, hemodynamic responses were negatively polarized during seizures, despite the application of a stimulus. The mean-field simulation revealed restricted propagation of activity due to stimulation and agreed well with fMRI findings. Results suggest that sensory processing is hindered or even suppressed by the occurrence of an absence seizure, potentially contributing to decreased responsiveness during this absence epileptic process.

Modeling the relationship between neuronal activity and the BOLD signal: contributions from astrocyte calcium dynamics.

Functional magnetic resonance imaging relies on the coupling between neuronal and vascular activity, but the mechanisms behind this coupling are still under discussion. Recent experimental evidence suggests that calcium signaling may play a significant role in neurovascular coupling. However, it is still controversial where this calcium signal is located (in neurons or elsewhere), how it operates and how relevant is its role. In this paper we introduce a biologically plausible model of the neurovascular coupling and we show that calcium signaling in astrocytes can explain main aspects of the dynamics of the coupling. We find that calcium signaling can explain so-far unrelated features such as the linear and non-linear regimes, the negative vascular response (undershoot) and the emergence of a (calcium-driven) Hemodynamic Response Function. These features are reproduced here for the first time by a single model of the detailed neuronal-astrocyte-vascular pathway. Furthermore, we analyze how information is coded and transmitted from the neuronal to the vascular system and we predict that frequency modulation of astrocytic calcium dynamics plays a key role in this process. Finally, our work provides a framework to link neuronal activity to the BOLD signal, and vice-versa, where neuronal activity can be inferred from the BOLD signal. This opens new ways to link known alterations of astrocytic calcium signaling in neurodegenerative diseases (e.g. Alzheimer's and Parkinson's diseases) with detectable changes in the neurovascular coupling.

Mean-field based framework for forward modeling of LFP and MEG signals.

The use of mean-field models to describe the activity of large neuronal populations has become a very powerful tool for large-scale or whole brain simulations. However, the calculation of brain signals from mean-field models, such as the electric and magnetic fields, is still under development. Thus, the emergence of new methods for an accurate and efficient calculation of such brain signals is currently of great relevance. In this paper we propose a novel method to calculate the local field potentials (LFP) and magnetic fields from mean-field models. The calculation of LFP is done via a kernel method based on unitary LFP's (the LFP generated by a single axon) that was recently introduced for spiking-networks simulations and that we adapt here for mean-field models. The calculation of the magnetic field is based on current-dipole and volume-conductor models, where the secondary currents (due to the conducting extracellular medium) are estimated using the LFP calculated via the kernel method and the effects of medium-inhomogeneities are incorporated. We provide an example of the application of our method for the calculation of LFP and MEG under slow-waves of neuronal activity generated by a mean-field model of a network of Adaptive-Exponential Integrate-and-Fire (AdEx) neurons. We validate our method via comparison with results obtained from the corresponding spiking neuronal networks. Finally we provide an example of our method for whole brain simulations performed with The Virtual Brain (TVB), a recently developed tool for large scale simulations of the brain. Our method provides an efficient way of calculating electric and magnetic fields from mean-field models. This method exhibits a great potential for its application in large-scale or whole-brain simulations, where calculations via detailed biological models are not feasible.