Pool Size and Synaptic Scaling

Paired controls separate fixed summed I→E coupling from fixed expected synaptic strength as the inhibitory pool grows.

Abstract

We tested inhibitory pool size under two controls: fixed summed inhibition and fixed expected per-synapse strength. Across 3 seeds, increasing I size from 16 to 256 changed mean E firing from 6.01 to 6 Hz at fixed summed coupling, versus 16.62 to 6 Hz at fixed synaptic strength. Pool-size effects depended on scaling. These untrained-network rates do not by themselves establish gamma oscillations.

Results

Pool-size scaling response

At nominal GI→E=G_{I \to E} = 2 μS, E rates changed from 6.01 to 6 Hz across 16–256 I neurons. At wˉI→E≈\bar{w}_{I \to E} \approx 7.81 nS, E rates fell from 16.62 to 6 Hz and I rates from 120.41 to 34.41 Hz. All tested fixed-mean levels showed this decrease, consistent with stronger summed inhibition and reduced E–I feedback (Fig. 1).

Four panels of E and I firing rates for inhibitory pools of 16, 64 and 256 neurons. Fixed summed coupling gives nearly flat rates; fixed expected synaptic strength gives falling rates.
Figure 1.

Reanalysed population firing rates. (A) Excitatory (E) rate with fixed expected summed coupling GI→EG_{I \to E}; (B) E rate with fixed expected synaptic strength wˉI→E\bar{w}_{I \to E}; (C) inhibitory (I) rate with fixed summed coupling; (D) I rate with fixed synaptic strength. Markers are means ±1 sample standard deviation across 3 seeds, with 8 trials per seed. Shared conditions reused the same simulations.

Methods

We compared two scaling controls using recorded outputs from simulations of untrained, dense recurrent excitatory–inhibitory networks, without additional simulation.

  1. Initialize fan-in-normalized weights. For an I→E matrix with 1024 excitatory columns and NIN_{I} inhibitory rows, each weight was WkjIE=GdrawNIW_{k j}^{I E} = \frac{G_{\text{draw}}}{N_{I}}, with Gdraw=max⁡(0,X)G_{\text{draw}} = \max ( 0 , X ) and X∼N(μinit,σinit2)X \sim \mathcal{N} ( \mu_{\text{init}} , \sigma_{\text{init}}^{2} ) a Gaussian draw of mean μinit\mu_{\text{init}} and standard deviation σinit=0.1μinit\sigma_{\text{init}} = 0.1 \mu_{\text{init}}, in μS. Defining expected summed coupling GI→E=E[Gdraw]G_{I \to E} = \mathcal{E} [ G_{\text{draw}} ] gives

    E[WkjIE]=wˉI→E=GI→ENI,E[∑k=1NIWkjIE]=GI→E.(1)\mathcal{E} [ W_{k j}^{I E} ] = \bar{w}_{I \to E} = \frac{G_{I \to E}}{N_{I}} , \quad \mathcal{E} [ \sum_{k = 1}^{N_{I}} W_{k j}^{I E} ] = G_{I \to E} . \tag{1}

    Here E\mathcal{E} averages over weight initialization, kk indexes inhibitory neurons and jj excitatory neurons; wˉI→E\bar{w}_{I \to E} is the expected conductance of one synapse, not an identical realised weight. An inhibitory volley gives

    Δginh,post=E,j=∑k∈IactiveWkjIE,(2)\Delta g_{\text{inh,post=E} , j} = \sum_{k \in \mathcal{I}_{\text{active}}} W_{k j}^{I E} , \tag{2}

    where Iactive\mathcal{I}_{\text{active}} is the active inhibitory set and Δginh,post=E,j\Delta g_{\text{inh,post=E} , j} the inhibitory conductance increment at E neuron jj, in μS; both weights and participation therefore enter the increment.

  2. Apply paired pool-size controls. We swept NI∈{N_{I} \in \{16, 64, 256}\}. Fixed-summed controls used parent means μinit∈{\mu_{\text{init}} \in \{1, 2, 4}\} μS, so wˉI→E∝1NI\bar{w}_{I \to E} \propto \frac{1}{N_{I}}; fixed-mean-synapse controls scaled μinit\mu_{\text{init}} with NIN_{I}, giving nominal wˉI→E∈{\bar{w}_{I \to E} \in \{3.90625, 7.8125, 15.625}\} nS. Nominal values approximate the post-clamp expectations; the arms coincide at 256 I neurons, with one additional shared condition at 64 I neurons and nominal 1 μS summed coupling.

  3. Drive and measure the networks. Nominal E→I summed coupling stayed at 1 μS; input weights used parent mean 1.2 μS, 95% initial zeros, compensation for zeroing and fan-in normalization. Each network received 784 independent 25 Hz input channels, implemented as Bernoulli spikes per 0.1 ms timestep, for 500 ms and 8 trials, with seeds 40, 41, 42. We included all conditions and averaged spike counts over the full duration, trials and neurons within each population; overlapping controls reused measurements, giving 3 seeds at each of 18 conditions from 42 distinct simulations. The reanalysis used these rates, not raw spike trains; population rates alone do not establish gamma oscillations [1].

  4. Display pool-size comparisons. We displayed the recorded population-rate responses for fixed total and fixed per-synapse coupling with their shared controls and seed aggregation.

Network diagrams

PING with variable inhibitory pool; illustrative model schematic.
Figure A1. PING with variable inhibitory pool. Illustrative schematic of the documented model; it is not a recovered historical execution bundle or a measurement. XX denotes input spike channels, EE/II excitatory/inhibitory populations and RR the class readout. Every weighted pathway is labelled WXYW_{XY}, with target XX first and source YY second. Arrowheads show excitatory or feedforward drive; bars show inhibition. Recording and score reductions are omitted.

Parameters

Model values are read from original native snapshots through the selected presentation’s compute ancestry. Paired values are E / I. Initializer parent means and SDs describe the original weight distributions, in µS for conductance projections and dimensionless for the current-based readout; they are not measured trained-weight statistics. Missing snapshots or fields are unavailable; current recipes and historical copied tables are not substituted. The recorded configuration identifies the input, numerical and analysis choices of that computation; it does not establish a separate upstream training protocol.

Parameter16 I cells256 I cells
Network
E / I neuronsUnavailable / UnavailableUnavailable / Unavailable
Output neuronsUnavailableUnavailable
Neurons
Capacitance, E / I (nF)Unavailable / UnavailableUnavailable / Unavailable
Leak conductance, E / I (µS)Unavailable / UnavailableUnavailable / Unavailable
Resting voltage, E / I (mV)Unavailable / UnavailableUnavailable / Unavailable
Threshold, E / I (mV)Unavailable / UnavailableUnavailable / Unavailable
Reset, E / I (mV)Unavailable / UnavailableUnavailable / Unavailable
Refractory period, E / IUnavailable / UnavailableUnavailable / Unavailable
Synapses and weights
Input→E synapseUnavailable / UnavailableUnavailable / Unavailable
E→I synapseUnavailable / UnavailableUnavailable / Unavailable
I→E synapseUnavailable / UnavailableUnavailable / Unavailable
E→Output synapseUnavailable / UnavailableUnavailable / Unavailable
Input→E initializer parent mean / SDUnavailable / UnavailableUnavailable / Unavailable
E→I initializer parent mean / SDUnavailable / UnavailableUnavailable / Unavailable
I→E initializer parent mean / SDUnavailable / UnavailableUnavailable / Unavailable
E→Output initializer parent mean / SDUnavailable / UnavailableUnavailable / Unavailable
Native parameter stageUnavailableUnavailable
Input and simulation
Timestep (ms)UnavailableUnavailable
Duration (ms)500500
Evaluation samplesUnavailableUnavailable
Seeds40 / 41 / 4240 / 41 / 42
Input rate (Hz)2525
Post-burn exclusion (ms)UnavailableUnavailable
Readout and training
Output neuron modelUnavailableUnavailable
EpochsUnavailableUnavailable
Checkpoint roleUnavailableUnavailable
Checkpoint purposeUnavailableUnavailable
Analysis
Frequency search (Hz)UnavailableUnavailable
Count-bin width (ms)UnavailableUnavailable
Autocorrelation maximum lag (ms)UnavailableUnavailable

Dataset

Upstream: —

Downstream: —

Current runset

RunDateDurationSizeEvidenceStorageRan onVerification
exp047-r003-present<1s140.3 KBCurrentR2LocalLocal: absent; R2: catalog-verified
exp047-r002-analyse<1s17.3 KBCurrentR2LocalLocal: absent; R2: catalog-verified
exp047-r001-compute4m 49s36.1 KBCurrentR2UnknownLocal: absent; R2: catalog-verified

References

  1. G. Buzsáki and X.-J. Wang. “Mechanisms of Gamma Oscillations.” Annual Review of Neuroscience 35, 203–225 (2012). doi:10.1146/annurev-neuro-062111-150444