Main limit first: reduction, b1 = 0.9, b2 = 0.999, x = 1e-4, five to ten starts per lr. At lr 0.275 and 0.3 the groups merge, so a "cycle" is no longer a period and I do not read those two lr. Several lines below were found after seeing the data and are marked post hoc.
Your numbers
The slope of ln(1/rho) on F is -0.0295 (lr 0.275) to -0.0578 (lr 0.15), your "-0.030 to -0.058". Between lr, d ln M / df = +0.0489 between f = 0.35 and 0.45. |S_F| is 0.0578, 0.0577, 0.0565 at lr 0.15, 0.175, 0.2, then 0.0486, 0.0411. Your sign argument holds: with ln(1/rho) = h(F) = -1/2 ln[(A_e F + a_e (1 - F)) / (A_l F + a_l (1 - F))], the slope inside a record is h'(F) and across lr M is about h(f), so d ln M / df = h'(f), the same sign. We observe opposite signs.
One correction to how S_F reads. h'(F) = -1/2 [(A_e - a_e) / v_e - (A_l - a_l) / v_l] with v_i = A_i F + a_i (1 - F). It is the difference of two relative burst excesses, each over the current mix, not the contrast between the amplitudes. S_F < 0 means the relative excess of e3 is larger than that of l3.
Does rho move inside a calm
No. Summed |change of ln rho|, calm against bursts, on regenerated records: lr 0.1: 0.00016 against 0.08138 (ratio 0.0020); lr 0.2: 0.00016 against 0.36341 (0.0004); lr 0.3: 0.00008 against 0.76665 (0.0001). The standard deviation of the change over a calm is 0.00000. So d_g is the increment of burst g.
d_g against rho_g
lr slope beta r R2 R2 of quintile means cycle (steps) slope of the model -(1 - b2^cycle)
0.1 -0.0707 -0.188 0.035 0.919 426 -0.347
0.15 -0.0599 -0.171 0.030 0.915 364 -0.306
0.2 -0.0497 -0.157 0.025 0.946 296 -0.256
0.225 -0.0351 -0.132 0.017 0.920 256 -0.226
0.25 -0.0313 -0.125 0.016 0.912 224 -0.201
0.275 -0.0672 -0.183 0.033 0.702 236 (merging)
0.3 -0.1975 -0.314 0.099 0.879 500 (merging)
The mean of d_g given rho_g lies on a line (quintile means at 0.9 or more for seven of nine lr). The cloud is not a clean line: |r| = 0.13 to 0.19 for lr <= 0.275. The simple slope falls from 0.0707 to 0.0313 over lr 0.1 to 0.25. Between lr 0.2 and 0.225 it drops 29 percent (local exponent of |beta| -2.97 against -0.70, -0.06, -1.00, -0.25 and -1.06 at the other intervals up to 0.25). The starts spread by about 15 percent (lr 0.2: -0.0583 to -0.0421; lr 0.225: -0.0408 to -0.0270), so that drop is about two sigma. A two-segment fit of ln|beta| over the nine lr puts the kink at 0.26, pulled up by the two merging points. So the slope does not clearly move at 0.2.
Adding the burst length, ln sqrt(v_l3) and the calm before to the regression takes R2 from 0.025 to 0.984 at lr 0.2 (0.035 to 0.983 at lr 0.1; 0.016 to 0.964 at lr 0.25). The map is nearly exact once these are in. The partial coefficient of ln rho is -0.3522, -0.3261, -0.3029, -0.2749, -0.2501, -0.2191, -0.1962 for lr 0.1 to 0.25, against the model -(1 - 0.999^cycle) = -0.3469, -0.3264, -0.3055, -0.2802, -0.2560, -0.2261, -0.2005: ratio 0.97 to 1.015. It is smooth in lr and follows the cycle length. That follows from the identity below; it is not an independent confirmation.
The exact map
Between two instants just before successive bursts (i0 = a_g - 1, i2 = a_(g+1) - 1, T steps), v_i(i2) = theta v_i(i0) + I_i, with theta = b2^T and I_i the net injection (cycle, calm included). Then rho_(g+1)^2 = omega_g rho_g^2 + (1 - omega_g) rho*_g^2, omega_g = theta v_l(i0) / v_l(i2), rho*^2 = I_e / I_l. This holds with no linearity assumption. I computed it on five regenerated records per lr. It holds to 7e-16 on every cycle, which is algebra since I defined I from the v's; the content is that both injections are positive in 100 percent of the cycles at every lr.
- omega: mean 0.6542, 0.6735, 0.6944, 0.7199, 0.7433, 0.7741, 0.7988 for lr 0.1 to 0.25, against b2^cycle 0.6541, 0.6736, 0.6943, 0.7198, 0.7432, 0.7739, 0.7986. That is stationarity of v_l (the ratio v_l(i0) / v_l(i2) averages to 1), not a test. Its sd is 0.014 to 0.021 (0.088 at lr 0.275, 0.212 at lr 0.3). Its local exponent rises from 0.240 to 0.344 as lr passes 0.2, because the cycle shortens faster (the number of groups grows with exponent 0.77, then 1.2).
- Target: the weighted mean of ln rho* equals the weighted mean of ln rho (stationarity again). Over lr 0.15 to 0.25, post hoc since the two merging points leave: exponent of 1/rho 0.02265 against 0.01367 (0.2 to 0.25), step 0.0090; of 1/rho* 0.02248 against 0.01314, step 0.0093. The step is in the target; omega does not enter. On the window I had written down (0.2 to 0.3, which contains the merging lr), the unweighted target step is 0.0041 against 0.0089 for M, so my prediction fails as written.
What moves with lr (post hoc)
The slope gamma of ln rho*_g on the length of burst g, per step:
lr 0.1 0.125 0.15 0.175 0.2 0.225 0.25
gamma -7.8e-6 +1.9e-5 +7.5e-6 +1.16e-4 +2.96e-4 +4.73e-4 +7.38e-4
corr -0.017 +0.029 +0.008 +0.099 +0.220 +0.337 +0.567
Up to lr 0.15 the target does not depend on the length of its burst (gamma within +-2e-5 across starts). From lr 0.175 it does, and increasingly. This is an amplitude that moves with lr, as you require, and it turns on just before the step at 0.2 to 0.225. It does not read as a contrast shrinking. Both injections depend on the length: the slope of ln I_e and ln I_l on length is -3.7e-3 and -3.9e-3 at lr 0.175, -8.9e-3 and -9.5e-3 at 0.2, -1.31e-2 and -1.41e-2 at 0.225, -2.04e-2 and -2.19e-2 at 0.25, against slopes within +-1e-3 up to lr 0.15. The longer the burst (by my cutoff), the less energy it injects into each of the two coordinates, and the two differ by 6 percent. The target keeps only that difference. So the thing the formula has to reproduce is the shared dependence, with the e3 / l3 difference second.
What I cannot claim: that this term explains the step. I tried to split the exponent of 1/rho* into a level at a fixed reference length and the length term. The result depends on the reference length: at 100 steps the exponents at the six intervals up to 0.25 are 0.019, 0.021, 0.036, 0.057, 0.062, 0.108; at 120 steps 0.017, 0.022, 0.022, 0.030, 0.031, 0.058; at 140 steps 0.014, 0.024, 0.008, 0.003, 0.001, 0.007. The sign of the contribution changes within the range of lengths I have (107 to 131 steps). I drop that split.
Where I missed
Held: rho constant in a calm (under 0.2 percent of the change); the multiple regression (R2 +0.15 or more); the injections positive; the target carries the step on lr <= 0.25. Missed: |r| of the simple line at or above 0.3 (it is 0.13 to 0.19); the simple slope within 30 percent of -(1 - b2^cycle) (it is 0.15 to 0.20 of it; the partial coefficient is within 3 percent, which I did not predict); a kink of the slope in [0.19, 0.22] (the fit puts it at 0.26); the unweighted target step in 0.005 to 0.013 on the window 0.2 to 0.3 (0.0041). The identity check (7e-16) and the omega agreement with b2^cycle are bookkeeping.
What I have not tested
x other than 1e-4, b1 other than 0.9, b2 other than 0.999 for this map. Whether the dependence on length is the same thing as the dependence on F found at round 90. Why the injections fall with length, and why from lr 0.175. The merged regime. A fit of A_e, A_l, a_e, a_l with the injections as inputs.
One question back
Since rho does not move inside a calm and the exact map has omega = b2^cycle smooth in lr, with the step in the burst's target; and since both injections of a burst start to fall with its length at lr 0.175 (slope -3.7e-3 against under 1e-3 up to 0.15) while the e3 / l3 difference is 6 percent of that, is the next measurement the energy a burst injects per step, its peak amplitude and its length taken separately (the slopes above only give the injection against length), to see why a longer burst injects less and why from lr 0.175, before any formula with fixed amplitudes: yes or no?
Scripts: verifier_tour91_carte.py (modes calme, carte, injection, injection_synthese).
Notebook: section "VRAIE CRITIQUE DE DIPANKARSARKAR, 04/11/2026 (tour 91, PAS simulée)".