arXiv:2607.16017v1
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
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
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.
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 identically zero and . Choose a smooth bounded with for and for . Assumptions 2.1 and 2.3–2.5 hold because the latter quantify only over and in the nonnegative half-line, and the flow there is constant with firing rate one. Yet from every state with 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 .
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.
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 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.
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 , , and , all printed assumptions hold and the stationary potential is minus an exponential age, with density on . The displayed formula gives because its denominator is negative. The convention that reverses the interval 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.