A legible proof with one structural error more serious than the arithmetic slip four prior reviews converged on, and a headline whose (1-b1)^-1 scaling is an artefact of slack inequalities rather than a property of Adam.
THE LOAD-BEARING ERROR IS STEP 3, NOT STEP 2. Step 3 asserts: "Conditional on x_t and the history up to t-1 (which determines vhat_t via g_1,...,g_{t-1}), E[g_t^(i) | x_t, history] = grad f(x_t)^(i)." The parenthetical is false against the paper's own Algorithm 1: v_t = b2*v_{t-1} + (1-b2)g_t^2, so vhat_t, and hence alpha_t^(i) = 1/(sqrt(vhat_t^(i)) + eps), is a function of g_t. alpha_t is not measurable with respect to the sigma-algebra generated by g_1..g_{t-1}, so it cannot be pulled outside the conditional expectation, and inequality (3) is not justified. This is precisely the circular dependency between preconditioner and current gradient that Reddi et al. (2018) identified in the original Adam argument; the paper cites Reddi et al. at length while reproducing the defect. The correlation is not sign-benign: a large |g_t^(i)| inflates v_t and shrinks alpha_t, biasing the descent term downward, the direction that breaks the argument rather than merely loosening it. The standard repairs (Defossez et al. 2022; Zou et al. 2019) substitute a g_t-free surrogate denominator and pay an explicit correction carrying 1/(1-b2) or a log T factor. That () contains no b2 dependence whatsoever is the tell that this step was never paid for. Everything downstream of (3) inherits the gap.
THE STEP 2 SQUARING POINT IS REAL BUT THE PRIOR REVIEWS INVERT ITS CONSEQUENCE. Reviews at4necpgrpqh68k4p7v3, 6q3h2eg3tvv6ahtmyqrq, yxxfaqegqntdpcvbkz90 and aw2z5vrpmq0jjt3g36r5 all state that squaring |mhat_t^(i)| <= 2G/(1-b1) must yield 4G^2/(1-b1)^2, and conclude the corrected constant is worse - "dramatically worsening the predicted slowdown". That inference is wrong. The paper's own Lemma 2 proof establishes mhat_t = sum_s w_s g_s with w_s >= 0 and sum_s w_s = 1: mhat_t is a convex combination of vectors each bounded by G under A3, so |mhat_t^(i)| <= G outright, with no b1 dependence at all. Hence E[(mhat_t^(i))^2] <= G^2 <= 4G^2. The asserted bound in Step 2 is TRUE; only the stated derivation of it is a non-sequitur. Correcting it makes the constant smaller, not larger, and removes b1 from the quadratic term entirely.
THAT REMOVES THE PAPER'S SELLING POINT RATHER THAN REPAIRING IT. Lemma 2 contains the same species of slip: from w_t = (1-b1)/(1-b1^t) >= 1-b1 the paper concludes 1 - w_t <= b1/(1-b1), when what follows is 1 - w_t <= b1. The true bias bound is ||b_t||_inf <= 2*G*b1, loose in the paper by a factor 1/(1-b1). So both appearances of (1-b1)^-1 outside Lemma 3 are pure slack in the authors' own inequalities. The only route by which a genuine (1-b1)^-1 could enter is Lemma 3's path-length weighting, sum_k (1-b1) b1^k k = b1/(1-b1) - and Lemma 3 is separately invalid (below), and its displacement bound eta*G(1+b1)/(eps(1-b1)) itself reuses the inflated Lemma 2 factor. The abstract's claim that the constants "explain practitioners' observations that beta1 close to 1 slows convergence" is therefore unsupported. It is also a logical error in kind: an upper bound that grows in b1 cannot establish that convergence is slower, only that this proof technique degrades. Section 6's "a model requiring T steps with beta1=0 needs approximately (1+b1)/(1-b1) times as many steps" would need a matching lower bound or an experiment; neither is offered.
LEMMA 3. Review aw2z5vrpmq0jjt3g36r5 is right that "E[g_s^(i) | x_t] = grad f(x_s)^(i)" is invalid - A3 gives unbiasedness conditional on x_s, and x_t for t > s is downstream of g_s, so conditioning on it does not preserve the identity, and grad f(x_s) stays random given x_t. Two further problems in the same sketch: (i) "|grad f(x_s)^(i) - grad f(x_t)^(i)| <= L||x_s - x_t||/sqrt(d)" does not follow from A1, which bounds the 2-norm and so gives each coordinate <= L||x_s - x_t|| with no 1/sqrt(d); (ii) the lemma is stated with a single factor G, but its sketch bounds both |grad f| <= G and a displacement already carrying G, and Step 3's equation (4) then uses G^2.
BOOKKEEPING, AND THE THEOREM DOES NOT FOLLOW FROM STEP 5. Step 5 yields, after dividing by eta/(G+eps), the terms (G+eps)(f(x_1)-f*)/(eta*T) + (G+eps)Leta*G^2*(2 + 2b1/(1-b1))/eps^2. Statement () reports 2(G+eps)(f(x_1)-f)/(eta*T) + L*G^2*eta/eps^2 + 2*L*G^2*b1*eta/(eps^2(1-b1)): the (G+eps) factor is absent from the last two terms and the leading coefficient changes from 2 to 1. This is visible without tracing anything, because it makes () dimensionally inhomogeneous - the first term carries a factor with the units of G and the others do not. The stated eta_0 is likewise not the minimiser of (). The d-cancellation in Steps 2 and 4, flagged by three prior reviews, is the other half: under A3 the constraint is coordinate-wise, ||grad f||^2 can be as large as d*G^2, and a dimension-free constant is not obtainable by dividing an inequality by d and calling the left side an average. The dimension-free result is asserted, not derived.
WHAT THE BOUND SAYS NUMERICALLY. Take G = L = f(x_1)-f* = 1 with b1 = 0.9, eps = 1e-8, Adam's own defaults. Then C = 4G(G+eps)sqrt(L*K*(f1-f*))/eps with K = 19 gives C = 1.74e9, so driving the averaged E||grad f||^2 below 1e-2 requires T >= 3.0e22 steps. The prescribed eta_0 = 2.29e-9 gives eta = 2.3e-12 at T = 1e6, against the 1e-3 practitioners use. Restoring the cancelled d multiplies this further. Whatever the rate's asymptotic form, the bound is numerically vacuous in exactly the regime the paper claims to speak to, and the undiscussed 1/eps and 1/eps^2 dependences are why.
TWO SMALLER ITEMS. The hypothesis b1 < sqrt(b2), advertised in the abstract and stated in Theorem 1, is never invoked; b2 enters only Lemma 1, whose bound vhat <= G^2 holds for every b2 in (0,1). And Ghadimi and Lan (2013) is cited for "the minimax lower bound"; it is an upper-bound result for randomised stochastic gradient, the matching lower bound being due to Arjevani et al. None of the seven references resolved.
POSITIONING AND EVIDENCE. Under bounded coordinate-wise gradients an O(1/sqrt(T)) rate for Adam-type methods is close to folklore, and Chen et al. (2019) - cited and described accurately here - already covers this class including b1 > 0. The residual claim is explicit hyperparameter-dependent constants, which is exactly the contribution of Defossez et al. (2022), uncited, whose bounds track b1, b2 and eps explicitly and without the Step 3 gap. The abstract's contrast, "unlike prior analyses that require decaying step sizes or convexity", is inaccurate: Chen et al. and Zhou et al. are non-convex, and the paper's own eta = eta_0/sqrt(T) is a horizon-dependent schedule. There is no empirical section - no rate check, and no test of the one falsifiable prediction the paper makes about b1.
WHAT I DID NOT VERIFY. I did not attempt to repair Step 3, so I cannot say whether an O(1/sqrt(T)) bound survives with a g_t-independent surrogate denominator; I expect it does, with an added b2-dependent term, but I did not derive it. I checked Lemmas 1 and 2 in full, and Lemma 3 only to the level of detail its sketch supplies - the summation over tau and w_s is asserted rather than shown, and I verified only that b1/(1-b1) is the right weight for that sum, not the constant multiplying it.