Leaky integrate-and-fire (LIF) network models are commonly used to study how the spiking dynamics of neural networks changes with stimuli, tasks or dynamic network states. However, neurophysiological studies in vivo often rather measure the mass activity of neuronal microcircuits with the local field potential (LFP). Given that LFPs are generated by spatially separated currents across the neuronal membrane, they cannot be computed directly from quantities defined in models of point-like LIF neurons. Here, we explore the best approximation for predicting the LFP based on standard output from point-neuron LIF networks. To search for this best “LFP proxy”, we compared LFP predictions from candidate proxies based on LIF network output (e.g, firing rates, membrane potentials, synaptic currents) with “ground-truth” LFP obtained when the LIF network synaptic input currents were injected into an analogous three-dimensional (3D) network model of multi-compartmental neurons with realistic morphology, spatial distributions of somata and synapses. We found that a specific fixed linear combination of the LIF synaptic currents provided an accurate LFP proxy, accounting for most of the variance of the LFP time course observed in the 3D network for all recording locations. This proxy performed well over a broad set of conditions, including substantial variations of the neuronal morphologies. Our results provide a simple formula for estimating the time course of the LFP from LIF network simulations in cases where a single pyramidal population dominates the LFP generation, and thereby facilitate quantitative comparison between computational models and experimental LFP recordings in vivo.
Leaky integrate-and-fire (LIF) networks are often used to model neural network activity. The spike trains they produce, however, cannot be directly compared to the local field potentials (LFPs) that are measured by low-pass filtering the potential recorded from extracellular electrodes. This is because LFPs are generated by neurons with spatial extensions, while LIF networks typically consist of point neurons. In order to still be able to approximately predict LFPs from LIF network simulations, we here explore simple proxies for computing LFPs based on standard output from LIF network simulations. Predictions from the various LFP proxies were compared with “ground-truth” LFPs computed by means of well-established volume conduction theory where synaptic currents corresponding to the LIF network simulation were injected into populations of multi-compartmental neurons with realistic morphologies. We found that a simple weighted sum of the LIF synaptic currents with a single universally applicable set of weights excellently capture the time course of the LFP signal when the LFP predominantly is generated by a single population of pyramidal cells. Our study therefore provides a simple formula by which the LFP signal can be estimated directly from the LIF network activity, providing a missing quantitative link between simple neural models and LFP measures in vivo.