Fix some $$p\in [0,1]$$
p
∈
[
0
,
1
]
and a positive integer n. The discrete Bak–Sneppen model is a Markov chain on the space of zero-one sequences of length n with periodic boundary conditions. At each moment of time a minimum element (typically, zero) is chosen with equal probability, and it is then replaced alongside both its neighbours by independent Bernoulli(p) random variables. Let $$\nu ^{(n)}(p)$$
ν
(
n
)
(
p
)
be the probability that an element of this sequence equals one under the stationary distribution of this Markov chain. It was shown in Barbay and Kenyon (in Proceedings of the Twelfth Annual ACM-SIAM Symposium on Discrete Algorithms (Washington, DC, 2001), pp. 928–933, SIAM, Philadelphia, PA, 2001) that $$\nu ^{(n)}(p)\rightarrow 1$$
ν
(
n
)
(
p
)
→
1
as $$n\rightarrow \infty $$
n
→
∞
when $$p>0.54\dots $$
p
>
0.54
⋯
; the proof there is, alas, not rigorous. The complimentary fact that $$\displaystyle \limsup _{n\rightarrow \infty } \nu ^{(n)}(p)< 1$$
lim sup
n
→
∞
ν
(
n
)
(
p
)
<
1
for $$p\in (0,p')$$
p
∈
(
0
,
p
′
)
for some $$p'>0$$
p
′
>
0
is much harder; this was eventually shown in Meester and Znamenski (J Stat Phys 109:987–1004, 2002). The purpose of this note is to provide a rigorous proof of the result from Barbay and Kenyon (in Proceedings of the Twelfth Annual ACM-SIAM Symposium on Discrete Algorithms (Washington, DC, 2001), pp. 928–933, SIAM, Philadelphia, PA, 2001), as well as to improve it, by showing that $$\nu ^{(n)}(p)\rightarrow 1$$
ν
(
n
)
(
p
)
→
1
when $$p>0.45$$
p
>
0.45
. (Our method, in fact, shows that with some finer tuning the same is true for $$p>0.419533$$
p
>
0.419533
.)