arXiv:2607.16017v1

On large networks of integrate-and-fire neurons with short-term synaptic plasticity

Quentin Cormier, Eva Löcherbach, Valentin Schmutz

math.PR60H1060K3545D0537A3060G55

Abstract

This work studies the mean-field limit of large networks of interacting stochastic leaky integrate-and-fire (LIF) neurons subject to short-term synaptic depression (STD). The macroscopic dynamics of this system is governed by a two-dimensional, non-linear McKean-Vlasov equation that couples the evolution of the neurons' membrane potentials with a synaptic depression variable. We investigate the long-time behavior of this limit system. To this end, we introduce an auxiliary linearized Markov process by freezing the interaction non-linearity to a constant. By exploiting the regeneration of the membrane potential at spike times, we are able to explicitly compute the conditional expectation of the synaptic depression variable, conditionally on the potential value, under the invariant measure of this two-dimensional linear process. This is a crucial ingredient to study time-dependent local perturbations thereof. As a consequence we are able to identify an analytic criterion guaranteeing the local stability of any invariant probability measure of the fully non-linear system. This stability criterion is formulated in terms of the zeros of the Laplace transform of a specific linear response function. Finally, we provide numerical examples demonstrating that the two-dimensional framework induces a richer spectrum of long-time dynamics than purely one-dimensional models. For example, synaptic depression can lead to low-frequency oscillations around a unique, unstable invariant measure where the oscillations are much slower than the neurons' firing rates.

AI-generated audit

Audit summary

Audited against arXiv v1

Not a correctness certificate. A “Correct” result may include yellow typos or minor formal corrections that do not affect substantive soundness. It means this audit found no unresolved substantive error under the stated criteria; it does not replace expert scrutiny or formal verification.

Current report

Detailed mathematical audit

Generated August 18, 2026
01Statements2 reported findingsContains wrong statements

The auxiliary recurrence theorem and the local-stability theorem are false under the stated assumptions, which constrain the flow only for nonnegative initial potentials. Smooth examples with absorbing negative states satisfy every printed assumption but have many invariant measures and no local attraction.

Theorem 2.6Incorrect

Uniqueness and positive Harris recurrence fail for allowed negative states

Pages 7–8 · Assumptions 2.1, 2.3–2.5 and Theorem 2.6 · arXiv:2607.16017v1

Take bb identically zero and α=0\alpha=0. Choose a smooth bounded ff with f(v)=0f(v)=0 for v1v\le-1 and f(v)=1f(v)=1 for v0v\ge0. Assumptions 2.1 and 2.3–2.5 hold because the latter quantify only over vv and v~\widetilde v in the nonnegative half-line, and the flow there is constant with firing rate one. Yet from every state (v,1)(v,1) with v1v\le-1 the process never jumps and stays fixed. Thus there are infinitely many invariant Dirac measures, contradicting positive Harris recurrence and uniqueness on the stated state space R×[0,1]\mathbb{R}\times[0,1].

Theorem 2.9Not able to verify

The local-stability theorem depends on the false global recurrence input

Pages 9–10 and 23–25 · Theorem 2.9 and its proof · arXiv:2607.16017v1

The proof uses the global exponential contraction and unique invariant measure asserted in Theorem 2.6 to define the decay range for the response kernel and to control every nearby initial law. The counterexample to Theorem 2.6 shows that Assumptions 2.3–2.5 do not supply those inputs on the stated state space. Because Theorem 2.9 additionally refers to a decay constant inherited from that failed input, the same example is not recorded as a separate formal counterexample to its full premise; rather, the local-stability conclusion remains unsupported unless an accessibility hypothesis is added and the argument is rechecked.

02Proofs2 reported findingsContains incorrect or incomplete proofs

The coupling proof establishes estimates only after access to the nonnegative region, while the theorem ranges over all real potentials. It then invokes ‘classical arguments’ for positive Harris recurrence without proving accessibility or tightness.

Proposition 3.12Incorrect as written

The coupling estimate does not cover the stated state space

Pages 15–17 · Proposition 3.12 and Proposition 3.13 · arXiv:2607.16017v1

The uniform coupling-time bound in Proposition 3.13 is taken only over nonnegative initial potentials. Proposition 3.12 immediately applies it to all v,v~Rv,\widetilde v\in\mathbb{R} and concludes positive Harris recurrence. The counterexample in the statement audit shows that negative absorbing states need never reach the controlled region, so this extension is invalid and cannot be repaired without adding a substantive accessibility hypothesis.

Proposition 2.7Incorrect as written

The invariant-density formula has the wrong sign for negative drift

Pages 8 and 18–20 · Proposition 2.7 and its proof · arXiv:2607.16017v1

For b=0b=0, f=1f=1, and α=1\alpha=-1, all printed assumptions hold and the stationary potential is minus an exponential age, with density eve^v on v0v\le0. The displayed formula gives ev-e^v because its denominator b(v)+αb(v)+\alpha is negative. The convention that reverses the interval [0,σ)[0,\sigma) does not turn a negative Lebesgue density into a probability density. An absolute-value/orientation treatment and a corresponding revision of the proof are required.

03Novelty0 reported findingsNo non-novelty findings

No non-novelty findings.

Detailed audit reportFull reasoning, manuscript locations, and references.
Open report PDF ↗

Author response

Challenge an audit finding

Local workflow preview

A listed author may submit formal evidence that an audit is inaccurate. The response would be considered in a fresh AI re-evaluation; it would not edit the audit automatically.

Paper
arXiv:2607.16017v1
Authors listed
Quentin Cormier, Eva Löcherbach, Valentin Schmutz
Audit date
August 18, 2026
  1. 01Establish identityMatch an authenticated scholarly identity to this paper.
  2. 02Submit evidenceIdentify the finding and give a formal mathematical response.
  3. 03Re-evaluateA separate agent checks the response and records a disposition.
Recommended production method

Authenticate with ORCID, then require an exact arXiv match

MathAudit should accept the identity only when ORCID OAuth authenticates the claimant's iD and this exact arXiv paper appears in arXiv's public authority feed for that iD. A matching name alone is not sufficient.

ORCID OAuth and arXiv authority-record lookup are not connected in this local prototype.

Email fallback for papers without a linked ORCID

A production fallback could send a one-time link only when the submitted address matches an independently maintained author-contact allowlist for this paper. MathAudit must return the same message for every address so the form cannot reveal which contacts are on that list.

This demonstration does not send, store, or compare the address.

Structured response preview

This form remains unavailable until production identity verification succeeds. Nothing entered here is submitted.

This panel never establishes authorship in the local prototype. A production result should be described narrowly as an authenticated ORCID match or control of a separately allowlisted author-contact mailbox.