Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

The Power-Law Approach of the Consumption Function to its Perfect-Foresight Asymptote

The synthesis proof — the backward-induction spine, certified by its four rigor engines, illustrated on HAFiscal's estimated calibrations

0. The theorem in plain economics: precautionary saving as a discounted stream of premia, built backward (a reader’s guide)

The problem. The consumer solves the standard buffer-stock problem of BST — CRRA utility, permanent and transitory income shocks, geometric discounting — whose complete statement is the paper’s problem section. This document works with the normalized problem under the assumptions listed in §1 ((A1)–(A6), as itemized in statement.md §1; each element deep-linked there to its anchor in the published BST draft).

The question. A very wealthy consumer facing labor-income risk behaves almost like a perfect-foresight consumer: income risk is trivial relative to wealth, so the precautionary motive nearly vanishes. How fast does it vanish as wealth grows? That rate matters computationally — every solution algorithm must assume something above its grid top — and theoretically: the buffer-stock canon proves the ratio version of the approach (BST’s limiting MPC bounds) but neither the level nor the rate (§3.3).

The answer in one sentence. Precautionary saving — the shortfall of consumption below the perfect-foresight rule — dies off as a POWER of wealth ((25), proved in §§5–6: Theorem 1, Paragraph; constructed transparently in §2) — the same functional family as the Pareto laws economists know from wealth, firm-size, and city-size distributions (Gabaix 1999, 2009, 2016) — with an exponent min(1,q)\min(1, \qup) built from patience, returns, and growth ((2)), and (when q>1\qup > 1) a level BB (Theorem 1, (39)) that is nothing more exotic than a growing-perpetuity (Gordon) value.

The route. §1 states the model and the perfect-foresight benchmark; §2 builds the theorem by backward induction from a perfect-foresight resumption date TT — the presentation’s spine; §3 states the claim formally and places it in the literature; §§4–7 are the four rigor engines that certify the spine’s steps; §8 defines and measures the journey the spine describes; §9 cashes the theorem out computationally; §10 tallies the engines; §11 documents figures and reproduction.

The buffer-stock reading. The household’s buffer stock is its end-of-period assets aa. These households are impatient: absent the precautionary motive they would run their assets down to zero (under BST’s maintained zero-income event >0\pZero > 0, which makes the natural borrowing limit zero; with θmin>0\theta_{\min} > 0 they would borrow, toward the natural limit), so the assets they do hold are a buffer against the income risk. The buffer stock has a target a^=m^c(m^)\hat{a} = \hat{m} - \cFunc(\hat{m}): at a^\hat{a}, market resources are expected to reproduce themselves — Ra^+1=m^\Rcal\,\hat{a} + 1 = \hat{m} at ψ1\psi \equiv 1, the asset-side form of BST’s buffer-stock target (§3’s The target has the ψ-general form). In this document’s region the buffer stock stands far above a^\hat{a} and the household runs it down — normalized PF wealth descends at rate ÞΓ\ThornG per period — while the behavior that would hold the buffer stock at its target (maintenance) is target-zone behavior, priced where the journey ends. What the theorem tracks along the descent is precautionary saving: x(m)\psavFunc(m), the eXtra saving relative to the perfect-foresight rule, which along the path satisfies — to leading order; the exact version, with its O(1/wˉ)O(1/\wbar) remainders sized, is §5’s one-step identity —

x(mt)  =  R1Et[x(mt+1)]  +  R1J(wˉt),J(wˉ)  :=  κ(ρ+1)σ22ÞΓwˉ.\psavFunc(\mNrmNow) \;=\; \Rcal^{-1}\,\E_{\prdt}\bigl[\, \psavFunc(\mNrmNxt) \,\bigr] \;+\; \Rcal^{-1}\, \bufAdjFunc(\wbarNow), \qquad \bufAdjFunc(\wbar) \;:=\; \frac{\kap(\rho+1)\sigma^2}{2\,\ThornG\,\wbar}.

J\bufAdjFunc is the Kimball precautionary premium (KPP), and its dating deserves care. It is generated by next period’s shock θt+1\tranShkNxt (σ2=Varθ\sigma^2 = \operatorname{Var}\theta): prudence — Kimball (1990)'s u/u=(ρ+1)/c-\uFunc'''/\uFunc'' = (\rho+1)/c — is evaluated where that shock will strike, tomorrow’s PF position (wealth ÞΓwˉt\ThornG\wbarNow, consumption κÞΓwˉt\kap\ThornG\wbarNow), and multiplies the consumption risk the MPC lets through, 12κ2σ2\tfrac{1}{2}\kap^2\sigma^2; the product is J(wˉt)\bufAdjFunc(\wbarNow), a function of today’s state only because tomorrow’s PF position is today’s state drifted once. It is then priced into today’s saving from one period away — whence its own R1\Rcal^{-1}, the same discount tomorrow’s expected saving carries. Unrolling the recursion, xx is the discounted stock of premia still to be provided along the journey, rung by rung: the KPP is always positive here (the forcing floor of §5, which is why x0x \geq 0), and while the premium itself need not rise monotonically as wealth falls, the sum of the premia still to be provided grows with every backward step of the journey.

The scope of the asset-pricing resemblance — narrower than it first appears, and the rest of this document should be read accordingly. The resemblance is an asymptotic statement: it holds for LARGE wˉ\wbar, and its value is that it proves the power-law tail — the Gordon-perpetuity mechanics of §2 and the engines of §§4–6 are the payoff, and for that purpose the analogy is kept. It is NOT valid where the solution lives: the ergodic set (§8.1) sits at SMALL wˉ\wbar, and there the approximation fails on every count. The displayed formula for J\bufAdjFunc is only the wˉ\wbar \to \infty limit — at smaller wˉ\wbar, especially near wˉh\wbar \approx h (bank balances b=m1\bNrm = m - 1 near zero), the terms the derivation drops are no longer small and no closed form replaces them; even where the formula applies, the premium’s shape (a flow proportional to 1/wˉ1/\wbar) is nothing like a standard payoff stream; and in the limit 0\pZero \downarrow 0 the consumption function’s curvature c\cFunc'' becomes UNBOUNDED near the point where the liquidity constraint would bind (BST proves this), so no smooth pricing approximation can survive at the bottom. That is what necessitates numerical solution of the model in the ergodic region — the EGM methodology with the theorem’s tails (§9); the supplement why the standard asset-pricing toolbox does not replace it records the tool-by-tool case. The pointer: §2 gives the full and careful exposition, deriving the recursion and the premium’s formula by backward recursion from a final period TT at which risk ends, one newly risky period per step. Last, the premium’s positivity is a high-wealth theorem; the buffer-stock adjustment it finances is signed in general — a household whose buffer stock exceeds what it wants provides less, one short of it provides more. With that scope on the table: for the tail, (1) strongly resembles an asset-pricing recursion, and the Gordon-perpetuity mechanics below are the payoff of taking the resemblance seriously — there, and only there.

(Objects, for this section: x(m):=cˉ(m)c(m)\psavFunc(m) := \bar{\cFunc}(m) - \cFunc(m) — precautionary saving, the eXtra saving induced by precaution, the shortfall of consumption below the perfect-foresight rule. The PF-wealth-coordinate form is the gap g(wˉ)g(\wbar) ((7)); wˉ(m):=b+h\wbar(m) := \bNrm + h, with b\bNrm BST’s bank balances; §1 displays the definitions.)

None of this is posited: §2.3 derives the recursion, discount and adjustment included, from two Euler equations and a second-order expansion ((20)), and Lemma 5.1 ((29), §4) is its exact form (the translation box there spells out the precise discount and per-period flow). The adjustment is κ(ρ+1)σ2/(2ÞΓwˉ)\approx \kap(\rho+1)\sigma^2/(2\ThornG\wbar) with wˉ\wbar PF total wealth — the Kimball precautionary premium (KPP) — the precautionary counterpart (Kimball 1990) of the Arrow–Pratt risk premium: absolute risk over PF wealth, hence tiny for the rich; (30), §3.4; the discounting per step is 1/R1/\Rcal, with the state drifting down to ÞΓwˉ\ThornG \wbar (the wealth recursion (8)). Everything in the theorem is bookkeeping on this recursion, and §2.4 does the bookkeeping in a single geometric sum ((21)): the exponent is about which adjustments dominate the present value (near ones or far ones), the amplitude BB is a perpetuity value, and the trichotomy is Gordon-growth’s convergent / r=gr = g / divergent trichotomy.

What is q\qup? (and why min(1,q)\min(1, \qup))

First the raw material — BST’s patience factors, written with the Old English thorn letter: Þ:=(Rβ)1/ρ\Thorn := (R\beta)^{1/\rho} is the absolute patience factor (the growth factor of perfect-foresight consumption); dividing by what consumption must keep up with gives return patience ÞR:=Þ/R\ThornR := \Thorn/R (RIC: ÞR<1\ThornR < 1) and growth patience ÞΓ:=Þ/Γ\ThornG := \Thorn/\Gamma (GIC: ÞΓ<1\ThornG < 1 — the impatience that makes wealth descend). §1 anchors each to its BST definition. Now two per-period log-rates, both familiar:

Their ratio is the star of the theorem (at ψ1\psi\equiv1; in general the root of (26)):

q  =  lnRþg  =  discounting speeddescent speed.\qup \;=\; \frac{\ln \Rcal}{-\GPRte} \;=\; \frac{\text{discounting speed}}{\text{descent speed}}.

The intuition. The household’s state is mm; perfect-foresight wealth is a function of it, wˉ(m):=b+h=m1+h\wbar(m) := \bNrm + h = m - 1 + h — the same state relabeled into PF-wealth units (§1’s convention). The bar is BST’s upper-bound mark, not a mean: wˉ\wbar is not a constant. A change of units, not of position: no consumer ever “begins at m=wˉm = \wbar”. A wealthy household runs its buffer stock down, so its mm falls period by period; measured as wˉ\wbar, the fall is asymptotically exactly geometric — wˉt+1ÞΓwˉt\wbarNxt \approx \ThornG\,\wbarNow, lnwˉ\ln\wbar dropping þg-\GPRte per period — which is exactly why the engines use these units. The journey ends not at a point but in §8’s return region around BST’s buffer-stock target m^\hat{m}: bounded shocks confine the recurrent range of wealth to a compact interval whose upper edge, in PF-wealth units, §8 computes as ζ=C0/(1ÞΓ)\zeta = C_0/(1-\ThornG) (8.1 How rarely is the tail visited? The dual (Kesten) root quantifies how rarely even that edge is visited from above). The time remaining from the current position is therefore τ(wˉ)=ln(wˉ/ζ)/(þg)\trvTime(\wbar) = \ln(\wbar/\zeta)/(-\GPRte) periods — LOGARITHMIC in the current state (made exact by §5.2’s shells and Paragraph's shell-to-wealth conversion; later sections write τ(wˉ)lnwˉ/(þg)\trvTime(\wbar) \approx \ln\wbar/(-\GPRte), absorbing the edge’s additive constant lnζ\ln\zeta into the amplitudes they track). Anything discounted exponentially in TIME therefore becomes a power law in WEALTH:

Rτ(wˉ)=(wˉ/ζ)lnR/(þg)=(wˉ/ζ)q.\Rcal^{-\trvTime(\wbar)} = (\wbar/\zeta)^{-\ln \Rcal/(-\GPRte)} = (\wbar/\zeta)^{-\qup}.

A power law in wealth is exponential discounting read on a logarithmic clock. (The same log-clock conversion generates the Pareto tails of the random-growth literature — there applied to cross-sections of cities and firms, here to a single household’s policy function.) So q\qup answers: by the time the big precautionary action arrives — near the target — how discounted is it, per unit of log-wealth?

Half-lives, and an impatient billionaire. (2)'s ratio carries an inversion worth making explicit: more impatience means a smaller exponent — the gap fades more SLOWLY up the wealth ladder. There is no paradox. An impatient household converges quickly in time, and that is exactly why its gap dies slowly in wealth: even from wˉ=10\wbar = 10 the return region is only decades away, so the destination still shapes behavior far up the ladder. A nearly-patient household is the mirror image — the return temporally remote even at modest wealth, behavior that looks frictionless almost immediately. The per-quarter rates at the recurring calibrations (§3.1), converted to half-lives (ln2/rate\ln 2/\text{rate}, in years):

calibrationdescent þg-\GPRte (%/qtr)half-lifediscounting lnR\ln\Rcal (%/qtr)half-liferatio
HS-mean1.2514\approx 14 yr0.5432\approx 32 yr0.43
COL-TOP0.5233\approx 33 yr0.5134\approx 34 yr0.98
GIC-CAP0.036479\approx 479 yr0.5134\approx 34 yr14\approx 14

Standing at wˉ=10\wbar = 10: the HS-mean household is 46\approx 46 years from the return region and still assigns it relevance weight Rτ=100.430.37\Rcal^{-\trvTime} = 10^{-0.43} \approx 0.37; the ceiling household is 1,600\approx 1{,}600 years away and assigns it 10-14 — economically invisible until the ψ\psi-tilt of (26) is priced in (q=1.47\qup = 1.47: weight 10-1.5 — twelve orders of magnitude of Jensen). Hence: an impatient billionaire is still, at the margin, a buffer-stock consumer — a couple of decades from the return region, behavior visibly shaped by it — while a patient household of far more modest wealth already behaves as if frictionless. Households that converge fast in time are the ones whose gap decays slowly in wealth, and conversely. (The simplified capstone’s race appendix replays this race in slow motion — including a two-line lognormal reading of (26) that reproduces the exact roots of all four calibrations to a few parts in ten thousand.)

A characteristic root, and a saddle path. Economists solve linear rational-expectations models by finding the characteristic roots and discarding explosive solutions (Blanchard–Kahn). Equation (26) is exactly the characteristic equation of the linearized Euler dynamics in log wealth: the homogeneous solutions of the shortfall recursion are the powers wˉq\wbar^{-q} with qq solving (26), and the exponent-detection theorem of §6 (Paragraph) is the corresponding stability selection — compensate by any wrong exponent and the transformed problem explodes or dies; only min(1,q)\min(1, \qup) is the saddle path.

(For readers who do not live in linear rational-expectations models: the ‘homogeneous solutions’ are the family of alternative decay shapes the recursion tolerates on its own before economics is imposed; almost all of them either explode along the induction or violate a boundary requirement, and ‘the saddle path’ is the one shape that survives both tests. §2.4 performs this selection explicitly, with no linear-systems vocabulary required.)

The estimated calibrations. Three calibrations recur in the figures and numerical checks. They arise in the context of the HAFiscal paper (https://llorracc.github.io/HAFiscal-Latest/), which estimates discount-factor distributions by education group for U.S. households: HS-mean (the mean discount factor of the High-School group), COL-TOP (the top atom of the College group’s discount-factor distribution — the most patient estimated type), and GIC-CAP (College fundamentals with the discount factor at the estimation’s cap on β\beta. The cap is a solution-method safeguard, not a point of theoretical interest in itself: for the solution algorithm to perform reliably, the search region for candidate parameters must stop a small gap short of the knife-edge value at which the Growth Impatience Condition fails; in practice HAFiscal sets the cap at 0.9995 times that GIC-failing value. GIC-CAP is thus as patient as the estimation permits — deliberately just below the knife-edge, never at it), plus an auxiliary resonance calibration (code id AUX-RES: College fundamentals with β\beta root-found so q=1\qup = 1 exactly — the Gordon r=gr = g knife-edge, constructed for the theory’s checks rather than reached by the estimation). The parameters are quarterly. One BST ingredient enters silently: the mortality risk BST introduces late (its Mortality section — a constant survival probability L\mathcal{L}, default L=1\mathcal{L} = 1; script L\mathcal{L}, distinct from §6’s fraktur functional L\Lfun). Under the Blanchard–Yaari no-leakage scheme (estates return to survivors as an actuarially fair excess return) mortality enters the solution phase through two substitutions — the discount β^=βL\hat{\beta} = \beta\,\mathcal{L} (HARK’s DiscFacEff) and the annuitized return R^=R/L\hat{R} = R/\mathcal{L} — whose product is invariant, β^R^=βR\hat{\beta}\,\hat{R} = \beta R: the Euler equation, Þ\Thorn, ÞΓ\ThornG, the GIC, and the descent rate are all unchanged, while the budget side shifts (RR/(LΓ)\Rcal \to R/(\mathcal{L}\Gamma), h1/(1LΓ/R)h \to 1/(1 - \mathcal{L}\Gamma/R), κ1ÞL/R\kap \to 1 - \Thorn\mathcal{L}/R, and q\qup rises). At the default L=1\mathcal{L} = 1 everything reduces to the model as written. Where mortality has real consequences is the simulation phase — the ergodic distribution is mortality-financed, and the two mortality channels pull the effective wealth support in opposite directions (killing thins it; the annuity return fattens it — measured in verify_sim_stability_checks.py) — which is why §8.1 and §9 invoke it and why the simulation-phase analysis lives in the ergodic-coverage companion (ergodic_coverage.md). In the frozen verification code these appear under their original identifiers CAL-HS, CAL-CTOP, CAL-CCAP, which are pinned instruments and keep their names.

Translation table.

in the mathematicsin economics
precautionary saving x(m)=cˉ(m)c(m)\psavFunc(m) = \bar{\cFunc}(m) - \cFunc(m), in the PF-wealth coordinate the gap g(wˉ)=κwˉc(m(wˉ))g(\wbar) = \kap \wbar - \cFunc(m(\wbar))the shortfall of consumption below the perfect-foresight (PF) rule
the forcing J(wˉ)\bufAdjFunc(\wbar)the Kimball precautionary premium (KPP) — the coming period’s risk, priced today; the flow the recursion prices
the one-step identity (29)strongly resembles an asset-pricing recursion: xx today = R1Et[x\Rcal^{-1}\,\E_{\prdt}[x tomorrow]] + R1×\Rcal^{-1}\times the adjustment
the backward unroll from TT ((21))the stock of premia still to be provided, accumulated over the remaining risky horizon — one newly risky period per backward step (§2.4)
the eigen-equation (26), root q\qupthe characteristic root: discounting speed ÷ descent speed
min(1,q)\min(1, \qup)the slower-fading of “today’s adjustment” (1/wˉ\propto 1/\wbar) vs “the discounted destination” (wˉq\propto \wbar^{-\qup}) (Paragraph)
the amplitude BBthe growing-perpetuity (Gordon) value of the adjustment stream (Theorem 1)
resonance q=1\qup = 1the Gordon knife-edge r=gr = g: value accrues linearly ⟹ the lnwˉ\ln \wbar law (Paragraph)
the compactified boundary z=0z = 0“infinitely wealthy” treated like a steady state: change coordinates, then linearize there
shellsconstant-percentage bands of PF wealth ≈ periods of the descent
the dual (Kesten) root ζ\zeta^* (8.1 How rarely is the tail visited? The dual (Kesten) root)the Pareto exponent of the wealth distribution’s random-growth machine

1. Model, conditions, and the perfect-foresight benchmark

The model in full (each element deep-linked to its anchor in the published BST draft; the statement is BST’s, so a reader can check any piece against the source).

The symbol correspondence. This corpus (and its HARK consumers) writes relative risk aversion as ρ\rho and permanent income growth as Γ\Gamma; BST’s published draft writes γ\gamma and G\mathcal{G} for the same objects. Everything else is shared: Þ:=(Rβ)1/ρ\Thorn := (R\beta)^{1/\rho} (the absolute patience factor), R:=R/Γ\Rcal := R/\Gamma, κ\kap and κˉ\bar{\kappa} (the limiting MPC bounds), \pZero (the zero-income probability), ψ\psi and θ\theta (permanent and transitory shocks). One structural note: BST separates the transitory shock into the zero-income atom and the employed-state shock (ξ\xi vs θ\theta there); this corpus writes the all-in transitory shock as the single symbol θ\theta — so the corpus’s θ\theta is BST’s ξ\xi, and σ2:=Varθ\sigma^2 := \operatorname{Var}\theta includes the atom.

Preferences. CRRA period utility u(c)=c1ρ/(1ρ)\uFunc(c) = c^{1-\rho}/(1-\rho) with discount factor β>0\beta > 0 (BST, Setup). BST maintains ρ>1\rho > 1, and the existence routing below inherits that; the decay arguments of §§3–5 use only ρ>0\rho > 0 once existence is in hand (assumption (A1) below).

The income process (Friedman–Muth; BST, Assumption assn-shocks, with the atom displayed at its transitory-shock definition). Permanent shocks ψ\psi: E[ψ]=1\E[\psi] = 1, supported on [ψmin,ψmax](0,)[\psi_{\min}, \psi_{\max}] \subset (0, \infty). Transitory shocks: income is zero with probability >0\pZero > 0 — the zero-income event, part of the assumption itself, not an optional garnish — and otherwise equals the employed-state shock (mean one, supported on a compact interval bounded away from zero, per BST) scaled by 1/(1)1/(1-\pZero); the all-in shock θ\theta then satisfies E[θ]=1\E[\theta] = 1 and 0θθmax<0 \leq \theta \leq \theta_{\max} < \infty. Both shocks are bounded and mean-one in levels — this general statement, BST’s own (Assumption assn-shocks’s bounded-support requirement), is what every theory section below argues from; the equiprobable discretizations of lognormal shocks used by the figures are implementations consistent with it (bounded support excludes unbounded specifications but is consistent with all standard discretizations of unbounded distributions onto finite grids). One clause on what the atom buys: the results here do not all need >0\pZero > 0 — with θmin>0\theta_{\min} > 0 existence routes through MST (2020) instead (below) — but >0\pZero > 0 is BST’s maintained specification, it is what the estimated calibrations use, and it is load-bearing for the form of WRIC, for the natural borrowing constraint (the discussion following Assumption assn-shocks), and for the MPC upper bound κˉ=11/ρÞR\bar{\kappa} = 1 - \pZero^{1/\rho}\ThornR.

Budget transitions and normalization — the objects aa, kk, bb, mm, self-contained. BST’s within-period sequence, in levels: the consumer ends period tt holding assets AtA_t (market resources minus consumption); those assets become next period’s kapital, Kt+1=AtK_{t+1} = A_t — the same resources, renamed at the date they earn their return (BST’s kt+1=atk_{t+1} = a_t; without this definition the word “kapital” for date t+1t+1 resources saved at tt would be mysterious); the return makes bank balances Bt+1=RKt+1B_{t+1} = R\,K_{t+1}; adding the period’s labor income gives market resources Mt+1M_{t+1}, from which consumption Ct+1C_{t+1} is chosen. Labor income is BST’s two-part process: permanent income grows by Pt+1=Γψt+1PtP_{t+1} = \Gamma\,\permShkNxt\,P_t (Γ1\Gamma \geq 1 predictable growth, ψ\psi a mean-one permanent shock), and the period’s receipt is Pt+1θt+1P_{t+1}\,\tranShkNxt (θ\theta the mean-one transitory shock of (A3), zero-income event included). Dividing every level by permanent income — lowercase letters are these normalized ratios throughout, and the Γψt+1\Gamma\permShkNxt divisor below is exactly that renormalization crossing the period boundary — the problem is BST’s normalized Bellman equation: choose ct(0,mt]\cNrmNow \in (0, \mNrmNow] (as in the statement page) with transitions

at=mtct,kt+1=at,bt+1=(R/(Γψt+1))kt+1,mt+1=bt+1+θt+1,\aNrmNow = \mNrmNow - \cNrmNow, \qquad \kNrmNxt = \aNrmNow, \qquad \bNrmNxt = \bigl(R/(\Gamma\permShkNxt)\bigr)\,\kNrmNxt, \qquad \mNrmNxt = \bNrmNxt + \tranShkNxt,

where at\aNrmNow is end-of-period assets, kt+1\kNrmNxt next period’s kapital, bt+1\bNrmNxt next period’s bank balances — the returns realized on the kapital, the market-wealth leg of wˉ\wbar ((6)) — and θt+1\tranShkNxt next period’s transitory income. The first-order condition is BST’s normalized Euler equation, ctρ=βREt[(Γψt+1)ρc(mt+1)ρ]\cNrmNow^{-\rho} = \beta R\,\E_{\prdt}[(\Gamma\permShkNxt)^{-\rho} \cFunc(\mNrmNxt)^{-\rho}]. At ψ1\psi \equiv 1 the growth-adjusted return is the constant R=R/Γ\Rcal = R/\Gamma and bt+1=Rat\bNrmNxt = \Rcal\,\aNrmNow.

Human wealth and the perfect-foresight benchmark. Normalized human wealth hhBST’s hNrm, the PDV of labor income including the current period’s — has infinite-horizon limit h=1+R1+R2+=1/(1R1)=R/(RΓ)h = 1 + \Rcal^{-1} + \Rcal^{-2} + \cdots = 1/(1-\Rcal^{-1}) = R/(R-\Gamma) under FHWC (one line from BST’s finite-horizon definition). BST’s perfect-foresight consumption function is cˉ(m)=(b+h)κ\bar{\cFunc}(m) = (\bNrm + h)\,\kap with limiting MPC κ:=1ÞR\kap := 1 - \ThornR, written in the decomposition its finite-horizon form displays — BST’s own b:=m1\bNrm := m - 1 (bank balances) — so cˉ=κwˉ\bar{\cFunc} = \kap\,\wbar in the wˉ\wbar variable displayed just below ((6)), and precautionary saving is x(m):=cˉ(m)c(m)\psavFunc(m) := \bar{\cFunc}(m) - \cFunc(m) — the object of the theorem. In the PF-wealth coordinate the same quantity is the gap g(wˉ):=κwˉc(m(wˉ))g(\wbar) := \kap\wbar - \cFunc(m(\wbar)), with g(wˉ(m))=x(m)g(\wbar(m)) = \psavFunc(m) — the engines’ object, because the PF drift is exactly multiplicative there (at zero gap (8) reads wˉt+1=ÞΓwˉt\wbarNxt = \ThornG\,\wbarNow); economics-facing prose keeps the household’s state mm. The one bookkeeping trap of the correspondence: hh’s leading term is this period’s income, so the PF rule carries the -1; dropping it shifts xx by κ\kap (checked to machine precision).

The perfect-foresight benchmark, displayed (the three defining displays the rest of the document consumes; the conditions named here carry BST anchors in the table below). Under GIC, RIC, and FHWC the normalized consumption function c(m)\cFunc(m) approaches its perfect-foresight asymptote from below:

cˉ(m)=κ(m1+h),κ=1ÞR (BST’s limiting MPC),h=11R1=RRΓ,\bar{\cFunc}(m) = \kap\,(m - 1 + h), \qquad \kap = 1 - \ThornR \ \text{(BST's limiting MPC)}, \qquad h = \frac{1}{1 - \Rcal^{-1}} = \frac{R}{R - \Gamma},

where hh is BST’s human wealth — the PDV of labor income including the current period’s unit of it, which is why the rule (BST’s own, with κ\kap from its limiting-MPC definition) carries the -1. The theorem is about the rate and form of that approach. Write

wˉ(m):=b+h=m1+h(perfect-foresight total wealth, human and market — a function of the state)\wbar(m) := \bNrm + h = m - 1 + h \qquad \text{(perfect-foresight total wealth, human and market — a function of the state)}

— the consumer’s total (human plus market) perfect-foresight wealth, viewed at the decision moment after this period’s returns have been realized on the kapital saved last period (kt=at1\kNrmNow = a_{t-1}): b\bNrm is BST’s bank balances (b=Rkt\bNrm = \Rcal\,\kNrmNow at ψ1\psi \equiv 1 — market resources net of the unit of current income that hh already counts; BST’s own decomposition). The bar follows the convention BST itself states for cˉ\bar{\cFunc} — “the overbar signifies that cˉ\bar{\cFunc} will be an upper bound as we modify the problem to incorporate constraints and uncertainty; analogously, κ\kap is the MPC’s lower bound” — bar = upper bound, underline = lower bound, matching κ\kap; ww itself is unused in BST, so no clash. One caution: bars so often mark CONSTANTS (means, bounds) that wˉ\wbar tempts a constant reading — it is not one: wˉ(m)=m1+h\wbar(m) = m - 1 + h moves one-for-one with the state; after its definition, wˉ\wbar denotes its value — the engines’ coordinate. In this variable the PF rule is simply cˉ=κwˉ\bar{\cFunc} = \kap\,\wbar, and

g(wˉ):=κwˉc(m(wˉ))    0(the gap),g(\wbar) := \kap\, \wbar - \cFunc(m(\wbar)) \;\geq\; 0 \qquad \text{(the gap)},

— the same quantity as precautionary saving: g(wˉ(m))=x(m)g(\wbar(m)) = \psavFunc(m).

The conditions, each with BST’s name, label, and live anchor:

conditionBST labelstatementrole in one clause
finite value of autarkyFVAC0<βΓ1ρE[ψ1ρ]<10 < \beta\Gamma^{1-\rho}\E[\psi^{1-\rho}] < 1value of always spending your permanent income is finite
return impatienceRICÞR:=Þ/R<1\ThornR := \Thorn/R < 1PF consumption is a positive share of PF wealth (κ>0\kap > 0)
weak return impatienceWRIC1/ρÞR<1\pZero^{1/\rho}\ThornR < 1RIC weakened by the atom — enough for existence
finite human wealthFHWCΓ/R<1\Gamma/R < 1 (i.e. R>1\Rcal > 1)the PDV hh converges
growth impatienceGIC (-Raw)ÞΓ:=Þ/Γ<1\ThornG := \Thorn/\Gamma < 1normalized wealth drifts down when high
strong growth impatienceGIC-ModE[Þ/(Γψ)]<1\E[\Thorn/(\Gamma\psi)] < 1delivers the target m^\hat{m}; Jensen ⇒ strictly stronger than GIC for nondegenerate ψ\psi

What the theorems maintain (Stage A / Stage B). Stage A (ψ1\psi \equiv 1): (A1) CRRA ρ>0\rho > 0; (A2) constant R,Γ,βR, \Gamma, \beta; (A3) the income process above with ψ1\psi \equiv 1 — i.i.d. θ\theta, E[θ]=1\E[\theta] = 1, σ2(0,)\sigma^2 \in (0,\infty), support [θmin,θmax][\theta_{\min}, \theta_{\max}] with 0θmin<1<θmax<0 \leq \theta_{\min} < 1 < \theta_{\max} < \infty; θmin=0\theta_{\min} = 0BST’s maintained >0\pZero > 0 zero-income atom (Assumption assn-shocks) — is permitted throughout and is what the estimated calibrations use; (A4) FHWC; (A5) RIC; (A6) GIC(-Raw). Stage B adds i.i.d. ψθ\psi \perp \theta, E[ψ]=1\E[\psi] = 1, suppψ[ψmin,ψmax](0,)\operatorname{supp}\psi \subseteq [\psi_{\min}, \psi_{\max}] \subset (0,\infty), and FVACψ^\psi. Unbounded-above θ\theta is covered by the L4′ extension (E[θk]<\E[\theta^k] < \infty for some k>max(1,q)+1k > \max(1,\qup)+1) — the one deliberate departure from bounded support (8.1 How rarely is the tail visited? The dual (Kesten) root, item 3).

Existence, the atom, and the bounds. With >0\pZero > 0 the non-degenerate limiting solution exists by BST’s contraction theorems under WRIC + FVAC (the Boyd weighted-norm route — the atom is native to BST’s machinery); with θmin>0\theta_{\min} > 0 one may instead invoke MST (2020) Theorem 2.2 (BST’s own remark maps the =0\pZero = 0 normalized problem into the MST framework as a special case). This routing is load-bearing and inherited by every proof in the program. The solution obeys BST’s MPC bounds κmc(m)κˉm\kap m \leq \cFunc(m) \leq \bar{\kappa} m with κˉ=11/ρÞR\bar{\kappa} = 1 - \pZero^{1/\rho}\ThornR, and c(m)/mκ\cFunc(m)/m \to \kap as mm \to \infty — the ratio limit whose level-and-rate refinement is this paper (§3.3).

The target. Under WRIC + FVAC + GIC-Mod there is a unique target m^>0\hat{m} > 0 with Et[mt+1/m]=1\E_{\prdt}[\mNrmNxt/m] = 1 at m^\hat{m} (BST’s buffer-stock-target theorem), computable from the implicit equation (m^c(m^))RE[ψ1]+1=m^(\hat{m} - \cFunc(\hat{m}))\,\Rcal\,\E[\psi^{-1}] + 1 = \hat{m}. An agent failing GIC-Mod has no target — the precise content behind the GIC-CAP name (§0); at the ceiling calibration even the mean log-drift is positive, so normalized wealth drifts up without bound (§8.1). Bounded support of the shocks (Assumption assn-shocks) is what makes the return region around m^\hat{m} (§8) a compact interval; BST’s own bounded-support discussion makes the same point — it is what permits the contraction-rate arguments characterizing the stable target to be made globally rather than only asymptotically.

The imported ladder (all PROVEN and review-hardened in stage_A_proof.md §§1–5 / stage_B_proof.md §§B1–B3; statements only):

(Lemma 5.1 and Corollary 5.2 — the exact one-step law and its monotone comparison sandwich, the induction invariant that certifies the spine — complete the imported ladder; they are stated in §4, where they are consumed.)

2. The spine: equation (1) by backward induction from TT

The theorem’s construction order is backward. Fix a date TT at which the consumer faces — or believes she faces — no further uncertainty, so that she resumes the perfect-foresight rule from TT on; every earlier period faces the transitory risk. Stepping backward from TT, each step adds one newly risky period, and §0’s discounted stock of premia grows exactly as the pricing recursion (1) says it should. This section constructs that recursion — the discount, the drift, and the adjustment all derived, not posited — and unrolls it to the published constants. The setting is the Stage-A world (ψ1\psi \equiv 1: transitory shocks only, per §1); §2.6 records that the construction is ψ\psi-general. Everything here is exact to leading order in 1/wˉ1/\wbar, with the dropped terms itemized in §2.5 and certified by the rigor engines of §§4–7 (the map: (EQ-1) ↦ Lemma 5.1’s (29) and Cor 5.2’s sandwich; the unrolled sum ↦ §5’s shells; the truncation at the log-clock exit ↦ §8’s return region). The spine’s named battery is verify_eq1_backward_checks.py (§11), which implements the backward iteration exactly and checks every displayed constant.

2.1 The terminal date, and the two identities that do all the work

The normalized transitory-only model of §1: CRRA u\uFunc, u(c)=cρ\uFunc'(c) = c^{-\rho}; period budget/transition mt+1=Rat+θt+1\mNrmNxt = \Rcal\, \aNrmNow + \tranShkNxt with at=mtct(mt)\aNrmNow = \mNrmNow - \cFunc\Now(\mNrmNow), R:=R/Γ>1\Rcal := R/\Gamma > 1 (FHWC), E[θ]=1\E[\theta] = 1, σ2:=Varθ\sigma^2 := \operatorname{Var}\theta, bounded θ\theta. The normalized Euler operator is

u(ct(m))=RβΓρ= ÞΓρ by (I1) below  Et[u(ct+1(R(mct(m))+θt+1))].\uFunc'(\cFunc\Now(m)) = \underbrace{R \beta \Gamma^{-\rho}}_{=\ \ThornG^{\rho}\ \text{by (I1) below}}\; \E_{\prdt}\bigl[\, \uFunc'\bigl( \cFunc\Nxt(\, \Rcal(m - \cFunc\Now(m)) + \tranShkNxt \,) \bigr) \,\bigr].

Perfect-foresight (PF) objects ((5), (6)): with κ=1ÞR\kap = 1 - \ThornR (ÞR=(Rβ)1/ρ/R<1\ThornR = (R\beta)^{1/\rho}/R < 1, RIC), h=R/(R1)h = \Rcal/(\Rcal - 1), and perfect-foresight total wealth wˉ=b+h\wbar = \bNrm + h (b\bNrm BST’s bank balances, §1),

cˉ(m)=κwˉ,and the PF-wealth recursionwˉt+1=ÞΓwˉt,ÞΓ=(Rβ)1/ρ/Γ<1 (GIC).\bar{\cFunc}(m) = \kap\,\wbar, \qquad \text{and the PF-wealth recursion} \qquad \wbarNxt = \ThornG\, \wbarNow, \qquad \ThornG = (R\beta)^{1/\rho}/\Gamma < 1 \ \text{(GIC)}.

(The wˉ\wbar-recursion is (8) at zero gap: wˉt+1=R(1κ)wˉ+[R(1h)+h]\wbarNxt = \Rcal(1-\kap)\wbar + [\Rcal(1-h) + h] and the bracket vanishes identically by h=R/(R1)h = \Rcal/(\Rcal-1); R(1κ)=RÞR=ÞΓ\Rcal(1-\kap) = \Rcal\ThornR = \ThornG.)

Two elementary identities carry the whole derivation:

(I1)  RβΓρ=ÞΓρ,(I2)  Rκ=RÞΓ.\text{(I1)}\ \ R\beta\Gamma^{-\rho} = \ThornG^{\,\rho}, \qquad\qquad \text{(I2)}\ \ \Rcal\,\kap = \Rcal - \ThornG.

(I1) is the PF Euler equation itself: ÞΓρ(cˉt+1/cˉt)ρ=1\ThornG^{\rho}\,(\bar{c}_{\prdNxt}/\bar{c}_{\prdt})^{-\rho} = 1 with PF consumption growth cˉt+1/cˉt=ÞΓ\bar{c}_{\prdNxt}/\bar{c}_{\prdt} = \ThornG. (I2) is RIC bookkeeping: κ=1ÞR\kap = 1 - \ThornR and ÞΓ/R=ÞR\ThornG/\Rcal = \ThornR.

The terminal assumption. At date TT the agent faces — or believes she faces — no further uncertainty, and therefore behaves according to the PF rule from TT on:

cT=cˉ,cT=κ,cT=0  (exactly: the PF rule is linear),xT0.\cFunc_T = \bar{\cFunc}, \qquad \cFunc_T' = \kap, \qquad \cFunc_T'' = 0 \ \ \text{(exactly: the PF rule is linear)}, \qquad \psavFunc_T \equiv 0.

Before TT she faces the transitory risk. We iterate backward.

2.2 One step back (T1T-1): the doubly-differentiated Euler subtraction

(The section title unpacked: write the Euler equation twice — once for the agent who knows next period’s shock will be at its mean, once for the agent facing the risk — Taylor-expand both to second order in the shock, and subtract. Everything zeroth- and first-order cancels; the surviving second-order term is the premium, and the surviving first-order-in-the-gap terms are the recursion.)

Write the two Euler equations at T1T-1 — the agent who KNOWS θT=1\theta_T = 1 (she is the PF agent, cˉ\bar{\cFunc}), and the agent who faces the risk (cT1=cˉxT1\cFunc_{T-1} = \bar{\cFunc} - \psavFunc_{T-1}):

u(cˉ(m))=ÞΓρ  u(cT(R(mcˉ(m))+1))(EE-det)\uFunc'(\bar{\cFunc}(m)) = \ThornG^{\rho}\; \uFunc'\bigl( \cFunc_T(\, \Rcal(m - \bar{\cFunc}(m)) + 1 \,) \bigr) \qquad \text{(EE-det)}
u(cˉ(m)xT1(m))=ÞΓρ  E[u(cT(R(mcˉ(m)+xT1(m))+θT))](EE-sto)\uFunc'(\bar{\cFunc}(m) - \psavFunc_{T-1}(m)) = \ThornG^{\rho}\; \E\bigl[\, \uFunc'\bigl( \cFunc_T(\, \Rcal(m - \bar{\cFunc}(m) + \psavFunc_{T-1}(m)) + \theta_T \,) \bigr) \,\bigr] \qquad \text{(EE-sto)}

Let mˉt+1:=R(mcˉ(m))+1\bar{m}_{\prdNxt} := \Rcal(m - \bar{\cFunc}(m)) + 1 be tomorrow’s perfect-foresight market resources (its PF wealth is wˉT=ÞΓwˉ\wbar_T = \ThornG\wbar), and abbreviate x:=xT1(m)x := \psavFunc_{T-1}(m). The stochastic argument is mˉt+1+Rx+(θT1)\bar{m}_{\prdNxt} + \Rcal x + (\theta_T - 1).

Expand (EE-sto) to second order in (θT1)(\theta_T - 1) around E[θ]=1\E[\theta] = 1 — the double differentiation of the spine’s name. With f(θ):=u(cT(mˉt+1+Rx+(θ1)))f(\theta) := \uFunc'(\cFunc_T(\bar{m}_{\prdNxt} + \Rcal x + (\theta - 1))):

f(θ)=u(cT)cT,f(θ)=u(cT)(cT)2+u(cT)cT,f'(\theta) = \uFunc''(\cFunc_T)\cdot \cFunc_T', \qquad f''(\theta) = \uFunc'''(\cFunc_T)\cdot (\cFunc_T')^2 + \uFunc''(\cFunc_T)\cdot \cFunc_T'',

so, taking expectations, the first-order term dies because E[θT1]=0\E[\theta_T - 1] = 0, and at TT the second term of ff'' dies because cT=0\cFunc_T'' = 0 exactly:

E[f(θT)]=u(cT(mˉt+1+Rx))+12σ2u(cT)κ2+(higher moments; bounded in 2.5 below).\E[f(\theta_T)] = \uFunc'(\cFunc_T(\bar{m}_{\prdNxt} + \Rcal x)) + \tfrac{1}{2}\sigma^2\, \uFunc'''(\cFunc_T)\,\kap^2 + \text{(higher moments; bounded in 2.5 below)}.

Expand the remaining argument in Rx\Rcal x (small; verified ex post) and subtract (EE-det). On the left, u(cˉx)u(cˉ)=(u(cˉ))x+O(x2)\uFunc'(\bar{c} - x) - \uFunc'(\bar{c}) = (-\uFunc''(\bar{c}))\,x + O(x^2). The zeroth-order terms cancel by (EE-det) — this is where “the first-order terms disappear”: what survives is first order in xx and second order in (θT1)(\theta_T - 1):

(u(cˉ(m)))x  =  ÞΓρ[u(cT(mˉt+1))κRx  +  12σ2u(cT(mˉt+1))κ2]+h.o.(-\uFunc''(\bar{\cFunc}(m)))\,x \;=\; \ThornG^{\rho}\,\bigl[\, \uFunc''(\cFunc_T(\bar{m}_{\prdNxt}))\,\kap\,\Rcal x \;+\; \tfrac{1}{2}\sigma^2\, \uFunc'''(\cFunc_T(\bar{m}_{\prdNxt}))\,\kap^2 \,\bigr] + \text{h.o.}

Both sides in CRRA, evaluated on the PF ray (cˉ=κwˉ\bar{c} = \kap\wbar, cT(mˉt+1)=κwˉT=κÞΓwˉ\cFunc_T(\bar{m}_{\prdNxt}) = \kap\wbar_T = \kap\ThornG\wbar): u(c)=ρcρ1-\uFunc''(c) = \rho c^{-\rho-1}, u(c)=ρ(ρ+1)cρ2\uFunc'''(c) = \rho(\rho+1)c^{-\rho-2}. Divide through by ρ(κwˉ)ρ1\rho(\kap\wbar)^{-\rho-1} and use (I1) twice (ÞΓρÞΓρ1=ÞΓ1\ThornG^{\rho}\,\ThornG^{-\rho-1} = \ThornG^{-1}):

x[1+κRÞΓ]  =  ÞΓ112κ(ρ+1)σ2(ÞΓwˉ)1.x\,\Bigl[\, 1 + \frac{\kap\Rcal}{\ThornG} \,\Bigr] \;=\; \ThornG^{-1}\cdot \tfrac{1}{2}\kap(\rho+1)\sigma^2\cdot(\ThornG\,\wbar)^{-1}.

By (I2), the bracket is exactly R/ÞΓ\Rcal/\ThornG — the messy marginal-utility coefficients collapse — and:

gT1(wˉ)  =  κ(ρ+1)σ22RÞΓ  wˉ1.g_{T-1}(\wbar) \;=\; \frac{\kap(\rho+1)\sigma^2}{2\,\Rcal\,\ThornG}\; \wbar^{-1}.

This is the time-(T1)(T-1) version of equation (1): a pure power law in wˉ\wbar with exponent 1 and an explicit amplitude. Its reading: 12κ2(ρ+1)σ2/c\tfrac{1}{2}\kap^2(\rho+1)\sigma^2/c is the Kimball precautionary premium in marginal-utility units (prudence u/u=(ρ+1)/c-\uFunc'''/\uFunc'' = (\rho+1)/c times the risk κ2σ2/2\kap^2\sigma^2/2 the MPC passes through); dividing by the local marginal-utility gradient and discounting turns it into consumption units. Note that (19) equals (1/R)[κ(ρ+1)σ2/(2wˉT)](1/\Rcal)\cdot[\kap(\rho+1)\sigma^2/(2\wbar_T)] — the adjustment priced at tomorrow’s PF wealth wˉT=ÞΓwˉ\wbar_T = \ThornG\wbar, then discounted one period at the interest-minus-growth factor R\Rcal. The bracket’s feedback term κR/ÞΓ\kap\Rcal/\ThornG has its own economics: consuming xx less today banks Rx\Rcal x more for tomorrow, which lowers tomorrow’s marginal utility and so offsets part of today’s precautionary wedge; identity (I2) says this offset is exactly what converts the naive discount into the risk-free one.

2.3 The induction step: equation (1) exactly, with its discount and adjustment derived

Now step to T2T-2 (and generically to any t<Tt < T), assuming the shortfall xt+1\psavNxt exists and is small in the bootstrap sense of §2.5 (it is: it was just constructed, of order wˉ1\wbar^{-1}). Repeat the subtraction. Three — not two — pieces survive, all first order in the small quantities {xt,xt+1}\{\psavNow, \psavNxt\} or second order in (θt+11)(\tranShkNxt - 1):

  1. the feedback (as in §2.2): ÞΓρu(cˉ(mˉt+1))κRxt\ThornG^{\rho}\,\uFunc''(\bar{\cFunc}(\bar{m}_{\prdNxt}))\,\kap\,\Rcal\,\psavNow — today’s shortfall banks assets;

  2. the injection (as in §2.2): ÞΓρ12σ2[u(cˉt+1)κ2+u(cˉt+1)ct+1]\ThornG^{\rho}\,\tfrac{1}{2}\sigma^2\,\bigl[\,\uFunc'''(\bar{c}_{\prdNxt})\,\kap^2 + \uFunc''(\bar{c}_{\prdNxt})\,\cFunc\Nxt''\,\bigr]; the c\cFunc'' piece now reads uxt+1=O(wˉρ1xt+1)-\uFunc''\cdot\psavFunc\Nxt'' = O(\wbar^{-\rho-1}\cdot\psavFunc\Nxt'') — two orders in wˉ\wbar below the u\uFunc''' term (§2.5), dropped;

  3. the propagation — new: tomorrow’s consumer is herself short of HER PF rule, so u(ct+1)=u(cˉt+1xt+1)\uFunc'(\cFunc\Nxt) = \uFunc'(\bar{c}_{\prdNxt} - \psavNxt) contributes ÞΓρ(u(cˉ(mˉt+1)))gt+1(mˉt+1)\ThornG^{\rho}\,(-\uFunc''(\bar{\cFunc}(\bar{m}_{\prdNxt})))\,g\Nxt(\bar{m}_{\prdNxt}), the continuation shortfall priced in marginal utility. (E[gt+1(mˉt+1+(θt+11))]=gt+1(mˉt+1)\E[g\Nxt(\bar{m}_{\prdNxt} + (\tranShkNxt - 1))] = g\Nxt(\bar{m}_{\prdNxt}) to leading order; the Jensen correction is second order, §2.5.)

The same division by ρ(κwˉ)ρ1\rho(\kap\wbar)^{-\rho-1}, the same (I1)–(I2) collapse, and — writing everything in PF-wealth coordinates, xt(m)=gt(wˉ)\psavFunc\Now(m) = g\Now(\wbar) with the argument naming the coordinate, mˉt+1\bar{m}_{\prdNxt}'s PF wealth =ÞΓwˉ= \ThornG\wbar — the recursion is

(EQ-1)gt(wˉ)  =  1Rgt+1(ÞΓwˉ)  +  1RJ(wˉ),J(wˉ)  :=  κ(ρ+1)σ22ÞΓwˉ  (the Kimball precautionary premium, KPP)\text{(EQ-1)}\qquad g\Now(\wbar) \;=\; \frac{1}{\Rcal}\, g\Nxt(\ThornG\, \wbar) \;+\; \frac{1}{\Rcal}\, \bufAdjFunc(\wbar), \qquad \bufAdjFunc(\wbar) \;:=\; \frac{\kap(\rho+1)\sigma^2}{2\,\ThornG\, \wbar} \ \ \text{(the Kimball precautionary premium, KPP)}

shortfall today = (discounted at 1/R1/\Rcal) shortfall tomorrow at the deterministically drifted PF wealth ÞΓwˉ\ThornG\wbar, plus the discounted adjustment. This is §0’s pricing recursion (1) — the expectation Et[x(mt+1)]\E_{\prdt}[\psavFunc(\mNrmNxt)] there evaluated at its certainty equivalent, the drifted PF wealth — with every term now derived: the discount is the interest-minus-growth factor 1/R1/\Rcal (not β\beta, not RβΓρR\beta\Gamma^{-\rho} — identities (I1)–(I2) do the conversion into consumption units along the PF ray); the state drifts by ÞΓ\ThornG per step (the wealth recursion (8)); and the flow is the Kimball precautionary premium of §3.4, with wˉJ(wˉ)=cJ:=κ(ρ+1)σ2/(2ÞΓ)\wbar\,\bufAdjFunc(\wbar) = c_J := \kap(\rho+1)\sigma^2/(2\ThornG) identically in this leading-order form — exactly the forcing limit (30) of the rigorous lemma. The rigorous counterpart is Lemma 5.1’s exact identity (29) (§4), (1+ηL(wˉ))Rg(wˉ)=(1+ηR(wˉ))Et[g(wˉt+1)]+J(wˉ)(1 + \etaL(\wbar))\,\Rcal\,g(\wbar) = (1 + \etaR(\wbar))\,\E_{\prdt}[g(\wbarNxt)] + \bufAdjFunc(\wbar), whose monotone comparison form (31)/(32) is (EQ-1) with the honest ±C0\pm C_0 PF-wealth-drift slack; the remainders ηL\etaL, ηR\etaR collect exactly the terms §2.5 drops.

2.4 Unrolling backward: the Gordon trichotomy with the published constants

Iterate (EQ-1) backward from xT0\psavFunc_T \equiv 0, k=Ttk = T - t steps, using J(ÞΓjwˉ)=cJÞΓj/wˉ\bufAdjFunc(\ThornG^{\,j}\wbar) = c_J\,\ThornG^{-j}/\wbar:

xTk(wˉ)  =  j=1kRjJ(ÞΓj1wˉ)  =  cJÞΓwˉj=1k(RÞΓ)j.\psavFunc_{T-k}(\wbar) \;=\; \sum_{j=1}^{k} \Rcal^{-j}\, \bufAdjFunc(\ThornG^{\,j-1}\, \wbar) \;=\; \frac{c_J\,\ThornG}{\wbar}\, \sum_{j=1}^{k}\, (\Rcal\ThornG)^{-j}.

Every claim of the theorem set is now a statement about this geometric sum — the Gordon growing-perpetuity trichotomy, with per-step ratio 1/(RÞΓ)1/(\Rcal\ThornG): each backward step, the adjustment grows by 1/ÞΓ1/\ThornG (the state has drifted down one rung) while discounting compounds by 1/R1/\Rcal. This is §0’s stock of premia literally accumulating, one newly risky period at a time.

(a) RÞΓ>1    q>1\Rcal\ThornG > 1 \iff \qup > 1: the perpetuity converges. kk \to \infty gives

g(wˉ)    cJÞΓwˉ1RÞΓ1  =  κ(ρ+1)σ22(RÞΓ1)wˉ1  =  Bwˉ1,g(\wbar) \;\to\; \frac{c_J\,\ThornG}{\wbar}\cdot\frac{1}{\Rcal\ThornG - 1} \;=\; \frac{\kap(\rho+1)\sigma^2}{2(\Rcal\ThornG - 1)}\, \wbar^{-1} \;=\; B\, \wbar^{-1},

exactly (27) — the number Theorem I (Theorem 1) certifies as the boundary value (39). The gap is wˉ1\wbar^{-1}: the realized exponent is 1=min(1,q)1 = \min(1, \qup).

(b) RÞΓ=1    q=1\Rcal\ThornG = 1 \iff \qup = 1: r=gr = g, the resonance. Every term contributes equally (cJÞΓ/wˉc_J\ThornG/\wbar each). The sum does not run forever: after τ(wˉ)lnwˉ/(þg)\trvTime(\wbar) \approx \ln \wbar/(-\GPRte) steps (the log-clock, (3); þg=lnÞΓ\GPRte = \ln \ThornG, NEGATIVE under GIC, so þg>0-\GPRte > 0 is the per-period descent rate) the drifted PF wealth ÞΓjwˉ\ThornG^{\,j}\wbar has descended to the buffer-stock body and the high-PF-wealth adjustment formula no longer applies. Arithmetic accumulation over the journey:

g(wˉ)    cJÞΓlnwˉ/(þg)wˉ  =  κ(ρ+1)σ22þglnwˉwˉ,g(\wbar) \;\approx\; c_J\,\ThornG\cdot\frac{\ln \wbar/(-\GPRte)}{\wbar} \;=\; \frac{\kap(\rho+1)\sigma^2}{-2\GPRte}\cdot\frac{\ln \wbar}{\wbar},

exactly the ψ1\psi \equiv 1 resonance constant of (40). The ×lnwˉ\times \ln \wbar knife-edge factor of (25) is nothing but r=gr = g arithmetic accumulation.

(c) RÞΓ<1    q<1\Rcal\ThornG < 1 \iff \qup < 1: the perpetuity diverges — the far adjustments dominate. The sum is dominated by its LAST terms: adjustments provided near the end of the descent, at moderate PF wealth, discounted back over the whole journey. Truncating at the log-clock exit jτ(wˉ)j \approx \trvTime(\wbar):

g(wˉ)    1wˉ(RÞΓ)τ(wˉ)×O(1)  =  1wˉwˉ(þglnR)/(þg)×O(1)  =  O(wˉlnR/(þg))  =  O(wˉq),g(\wbar) \;\sim\; \frac{1}{\wbar}\,(\Rcal\ThornG)^{-\trvTime(\wbar)} \times O(1) \;=\; \frac{1}{\wbar}\, \wbar^{(-\GPRte - \ln \Rcal)/(-\GPRte)} \times O(1) \;=\; O\bigl(\wbar^{-\ln \Rcal/(-\GPRte)}\bigr) \;=\; O(\wbar^{-\qup}),

recovering (26)'s exponent at ψ1\psi \equiv 1: q=lnR/(þg)\qup = \ln \Rcal/(-\GPRte) ((2)). Equivalently and more structurally: the homogeneous solutions of (EQ-1) are the self-similar modes wˉq~\wbar^{-\qtilde} with per-step multiplier R1ÞΓq~\Rcal^{-1}\ThornG^{-\qtilde}; self-consistency R1ÞΓq~=1    ÞΓq~=R\Rcal^{-1}\ThornG^{-\qtilde} = 1 \iff \ThornG^{-\qtilde} = \Rcal — the ψ1\psi \equiv 1 eigen-equation (26) — and modes shallower than q\qup die out under backward iteration (λ<1\lambda < 1), modes steeper explode (λ>1\lambda > 1), so the far tail is the slower of {the particular wˉ1\wbar^{-1} injection mode, the homogeneous wˉq\wbar^{-\qup} mode}: exponent min(1,q)\min(1, \qup), which is (25) — the mode selection §6’s Theorem γ-T (Paragraph) certifies. Two features of the sharp statements fall out for free: the amplitude in this regime is set where the truncation happens — in the body, NOT locally at large wˉ\wbar: the maintenance end of the journey, where the descent lands on the hold-at-target (target-zone) behavior of §0’s reading, is for q<1\qup < 1 literally what sets the amplitude (why this case has no local closed-form BB; cf. t1_amplitude_supplement.md) — and the backward steps land on the geometric lattice {ÞΓjwˉ}\{\ThornG^{\,j}\wbar\}, i.e. an arithmetic lattice of spacing þg-\GPRte in lnwˉ\ln \wbar, which is the origin of the log-periodic factor P(lnwˉ)P(\ln \wbar) in the sharp q<1\qup < 1 form (statement.md’s Theorems A2/B2 lattice sum with period þg-\GPRte; the envelope is Paragraph).

(Concretely: feed a trial power-law shape through one backward step of the recursion with the forcing switched off. It comes back as the same shape multiplied by a constant factor that depends on the trial exponent: shallower-than-critical shapes come back smaller — they die out along the induction — while steeper ones come back larger — they explode. ‘Mode’ means nothing more than a shape the step maps to a multiple of itself.)

The finite-horizon corollary, for free. At any finite distance from TT the sum (21) is finite, so the gap is exactly-order wˉ1\wbar^{-1} at ANY finite horizon regardless of q\qup — the Method of Moderation’s “the slope is exactly 1 at finite horizon; min(1,q)\min(1, \qup) only in the limit” (§9).

2.5 Why the constants come out exact: the 1/wˉ1/\wbar order count

The expansion parameter is NOT σ\sigma: it is the size of the shock relative to PF wealth. Fix the θ\theta-distribution and let wˉ\wbar grow. With q~:=min(1,q)\qtilde := \min(1, \qup) and x=O(wˉq~)x = O(\wbar^{-\qtilde}) the bootstrap class, the dropped terms are:

Table 1:Error bookkeeping: every dropped term is subleading in 1/wˉ1/\wbar

dropped termsize relative to the kept terms
E[(θ1)3]uκ3/6\E[(\theta-1)^3]\cdot \uFunc''''\,\kap^3/6 (third moment)O(1/wˉ)O(1/\wbar)
uct+1\uFunc''\cdot \cFunc\Nxt'' in ff'' (gap curvature)O(wˉ(q~+1))O(\wbar^{-(\qtilde+1)}) — an order below even O(1/wˉ)O(1/\wbar)
u\uFunc' linearization curvature O(x2)O(x^2)O(x)=O(wˉq~)O(x) = O(\wbar^{-\qtilde})
Jensen in propagation (E[x(mˉt+1+θ1)]\E[\psavFunc(\bar{m}_{\prdNxt} + \theta - 1)])O(σ2x/x)=O(1/wˉ2)O(\sigma^2\,x''/x) = O(1/\wbar^2) relative
evaluation-point shifts (aa vs aˉ\bar{a})O(x/wˉ)O(x/\wbar) relative

All are O(1/wˉ)O(1/\wbar) or better relative to the retained terms — for a FIXED specification; the constants are NOT uniform over the 0\pZero \downarrow 0 family, in whose limit c\cFunc'' becomes unbounded near the constraint kink (BST) and the threshold below grows accordingly; nothing in this section is claimed for small wˉ\wbar — uniformly on wˉwˉ0\wbar \geq \wbar_0 over the bootstrap class — the same inventory Lemma 5.1 certifies as (1+ηL,R)(1 + \eta_{L,R}) factors with η=O(1/wˉ)\eta = O(1/\wbar) ((29), §4). So the truncation at second order is not an approximation that happens to be decent: every neglected term is subleading in 1/wˉ1/\wbar, which is why the construction lands on the exact (27), the exact resonance constant, and the exact eigen-root — the same constants the rigorous chain certifies (Theorem 1, Paragraph). “Second order around E[θ]\E[\theta]” is the whole story at the top because prudence (u\uFunc''') is the lowest-order derivative that sees risk at all, and 1/wˉ1/\wbar scaling demotes everything beyond it.

What the backward telling does not prove — and where the rigor engines take over: that the true infinite-horizon cc is the TT \to -\infty limit of this construction with uniform control (value-iteration convergence with the stated uniformity is part of the rigorous chain); that the dropped O(1/wˉ)O(1/\wbar) remainders cannot conspire (Lemma 5.1’s ηL\etaL, ηR\etaR bounds plus the comparison engine (31)/(32) handle exactly this; §4); and the sharp forms at q1\qup \leq 1 (Theorems A2/B2 and the γ-engine statements of §6).

2.6 The spine is ψ\psi-general (pointer)

Nothing above is special to ψ1\psi \equiv 1. The W3 spike (eq1_backward_psi_spike.md; battery verify_eq1_psi_spike_checks.py, ALL PASS, cross-run byte-identical) re-runs the doubly-differentiated subtraction with permanent shocks and finds three exact ψ\psi-cancellations: (K1) the ψρ\psi^{-\rho} Euler weight exactly absorbs the ψ\psi-dependence of next-period consumption’s level along the PF ray; (K2) the ψ\psi-weighted innovation mean is exactly zero, so the first-order term dies ψ\psi-generally and the adjustment keeps its Kimball form with σ2σB2\sigma^2 \to \sigma_B^2 ((42)); (K3) the asset-feedback coefficient is ψ\psi-invariant, so the discount stays exactly 1/R1/\Rcal. The mode multiplier becomes λψ(q~)=E[ψ1+q~]ÞΓq~/R\lambda_\psi(\qtilde) = \E[\psi^{1+\qtilde}]\,\ThornG^{-\qtilde}/\Rcal; self-consistency λψ(q~)=1\lambda_\psi(\qtilde) = 1 is the eigen-equation (26) exactly, and the unrolled amplitude at q>1\qup > 1 is (28) exactly. §7 states the rigorous Stage-B results; the derivation itself is not imported.

3. The formal claim and its place in the literature

3.1 The claim

(A reading guide for the displays that follow: they are one claim specialized three ways. The exponent is always the smaller of 1 and q_∞. When q_∞ exceeds 1 an amplitude LEVEL exists — the Gordon value, with its ψ-general form. Exactly at q_∞ = 1 the level is replaced by the logarithmic resonance form. Below 1 the amplitude is set in the body, and a log-periodic wobble can ride on it.)

§2 built the approach; here is its formal statement. An infinite-horizon consumer with CRRA utility (coefficient ρ\rho) faces the Friedman–Muth income process of §1 — i.i.d. mean-one permanent shocks ψ\psi and transitory shocks θ\theta, whose specification includes the zero-income atom — under the buffer-stock conditions GIC, RIC, and FHWC ((A1)–(A6), §1). The statement lives in the perfect-foresight wealth coordinate (6), with the gap g(wˉ)g(\wbar) of (7). The gap decays as a power law with exponent min(1,q)\min(1, \qup),

g(wˉ)wˉmin(1,q)(×lnwˉ exactly at the knife-edge q=1),g(\wbar) \asymp \wbar^{-\min(1,\, \qup)} \qquad (\times \ln \wbar \ \text{exactly at the knife-edge } \qup = 1),

where q\qup is the unique positive root of the eigen-equation

E[ψ1+q]=RÞΓq(ψ1:  q=lnRþg,þg:=lnÞΓ),\E[\psi^{1+q}] = \Rcal\, \ThornG^{\,q} \qquad\Bigl(\psi \equiv 1:\ \ \qup = \frac{\ln \Rcal}{-\GPRte},\quad \GPRte := \ln \ThornG\Bigr),

with sharp closed-form amplitudes. At q>1\qup > 1,

wˉg(wˉ)B=κ(ρ+1)σ22(RÞΓ1)(transitory-only; σ2=Varθ),\wbar\, g(\wbar) \to B = \frac{\kap(\rho+1)\sigma^2}{2(\Rcal\ThornG - 1)} \qquad \text{(transitory-only; } \sigma^2 = \operatorname{Var}\theta\text{)},
wˉg(wˉ)Bψ=κ(ρ+1)σB22(RÞΓE[ψ2])(permanent shocks; σB2 below),\wbar\, g(\wbar) \to B_\psi = \frac{\kap(\rho+1)\sigma_B^2}{2(\Rcal\ThornG - \E[\psi^2])} \qquad \text{(permanent shocks; } \sigma_B^2 \text{ below)},

and at q=1\qup = 1, (wˉ/lnwˉ)g(wˉ)κ(ρ+1)σB2/(2E[ψ2]L(1))(\wbar/\ln \wbar)\, g(\wbar) \to \kap(\rho+1)\sigma_B^2 / (2\,\E[\psi^2]\,\Lfun'(1)) (at ψ1\psi\equiv 1: κ(ρ+1)σ2/(2þg)\kap(\rho+1)\sigma^2/(-2\GPRte); (40)). At q<1\qup < 1 the sharp form carries a periodic factor PP (Theorems A2/B2 of the original chain; the envelope is Paragraph). Figure 1 shows the object and the three slopes.

Left panel, the consumption function and its perfect-foresight asymptote visibly merging by m equals 5000 on a log scale, with the gap shaded; right panel, the gap versus wealth on log-log axes for the three estimated calibrations, each with its min(1,q_∞) reference slope.

Figure 1:The gap and its power law. Left: c(m)\cFunc(m) and the PF asymptote κwˉ=κ(b+h)\kap\,\wbar = \kap(\bNrm+h) for HS-mean, on m[m^,5000]m \in [\hat{m}, 5000] (log–log, from the target m^\hat{m} up): the two rules visibly merge — at the right edge human wealth (h=184.6h = 184.6 for this calibration) is 3.6%\approx 3.6\% of PF total wealth wˉ\wbar — with the gap xx shaded. Right: xx vs wˉ\wbar on log–log axes for the three estimated calibrations, with dashed reference slopes min(1,q)-\min(1,\qup).

How to read this figure. Left: the true consumption function and the perfect-foresight line, plotted to high enough wealth (m=5000m = 5000, where human wealth is 3.6%\approx 3.6\% of PF total wealth) that the two visibly merge — the shaded wedge between them (precautionary saving xx) is what the theorem is about. Right: log–log paper, where a power law is a straight line — exactly as in the Pareto plots economists draw for wealth or city-size distributions. The slope is the exponent min(1,q)-\min(1, \qup). A shallower line means prudence fades more slowly with PF wealth: for HS-mean (slope -0.38) the precautionary motive remains quantitatively alive far up the wealth ladder, while at the patience ceiling (slope -1) it dies off at the fastest admissible rate.

Why measured slopes bend: the purification rate. At any finite depth the plotted lines are not yet perfectly straight, and the theorem also says how fast they straighten. To leading order the gap is not one decaying thing but a mixture of two dying components: the homogeneous component of the pricing recursion (1), fading at the eigen-root rate q\qup of (26) per ee-fold of wˉ\wbar, and the component continually re-fed by the KPP forcing, fading at rate 1 per ee-fold (the premium itself shrinks like 1/wˉ1/\wbar). The asymptotic power law belongs to whichever component dies more slowly — that is the min(1,q)\min(1,\qup) in the theorem — and the other is a contaminant still visibly present at any realistic depth. Because the ratio of two exponentials decays at the difference of their rates, the mixture purifies at the mode-gap rate 1q|1-\qup|: the contaminant’s share shrinks by the factor e1qe^{-|1-\qup|} per ee-fold of wˉ\wbar, so 1/1q1/|1-\qup| is the number of ee-folds one must climb to shrink the measured slope’s error by a factor of ee. At COL-TOP (q=0.694\qup = 0.694, mode gap 0.31) a unit of purification costs 3.3\approx 3.3 ee-folds — a ×26\times 26 climb in PF wealth — which is why the fitted deep slope there is still 7%7\% shy of q-\qup even at wˉ106\wbar \sim 10^{6} (digits and fit windows in figures/FIGURES.md), while at HS-mean (q=0.381\qup = 0.381, mode gap 0.62) the same unit costs only 1.6\approx 1.6 ee-folds and the fitted slope agrees with q-\qup to 1.5%1.5\%. The same arithmetic is the speed limit on any tail attached at a finite grid top: what a knot-based extrapolation measures at the top node is the mixture, not the asymptote, and the measured slope’s excess over the true exponent dies off no faster than e1qe^{-|1-\qup|} per ee-fold beyond the knot. Near the resonance q=1\qup = 1 the mode gap closes and purification stalls — only the logarithmic factor separates the components — so there the knot’s measured slope stays informative arbitrarily far out.

All computational figures use calibrations estimated in the HAFiscal paper (full provenance in figures/FIGURES.md; discount factors are survival-adjusted, as there): HS-mean, based on the estimated parameters of a representative high-school-educated household (q=0.381\qup = 0.381); COL-TOP, the most patient college-educated type as estimated (the top of the estimated discount-factor distribution, q=0.694\qup = 0.694); and GIC-CAP, the college patience ceiling (q=1.474\qup = 1.474): its β\beta is capped a hair below the boundary at which BST’s raw growth impatience condition (GIC) would fail; at that patience, BST’s Strong Growth Impatience condition (GIC-Mod, §1) — whose failure marks the frontier beyond which an agent is so patient that no target level of wealth exists and wealth grows without bound — has already failed (§8.1). The three span the trichotomy: the two estimated calibrations live in the q<1\qup<1 regime (exponent q\qup itself), the patience ceiling in the q>1\qup>1 regime (exponent 1). Two clearly-labeled auxiliary calibrations appear only where a regime cannot be exhibited at estimated patience (the exact-resonance panel; the resolved periodic-factor inset).

3.2 Situating the result: function-approach, not distribution-tail

A reader who knows the modern heavy-tails literature in economics should be told immediately what this theorem is not. The celebrated power laws of Toda (2014), Benhabib–Bisin–Zhu (2011, 2015), Stachurski–Toda (2019), Beare–Toda (2022), and Gouin-Bonenfant–Toda (2023) concern the cross-sectional distribution of the level of wealth (or income, or firm size): Pareto tails of a stationary distribution, with the tail exponent solving a Kesten-type root equation for the forward dynamics. This paper’s object is different in kind: it is the consumption function itself — a policy function, not a distribution — and the claim is that c(m)\cFunc(m) approaches the perfect-foresight solution according to a power-law decay process in the individual decision-maker’s wealth as that wealth goes to infinity.

The distinction is sharp, not cosmetic. In the transitory-only model (ψ1\psi \equiv 1, bounded θ\theta) under GIC, the ergodic distribution of wealth has compact support — there is no Pareto tail of any kind, a fact whose econ-packaged form is precisely Stachurski–Toda’s (2019) impossibility theorem (risk-free saving + constant β\beta with βR<1\beta R < 1 ⇒ wealth inherits the income tail). Yet the consumption gap is exactly a power law. Formally, the two literatures’ root equations are Mellin duals: the distributional-tail equation E[(ÞΓ/ψ)ζ]=1\E[(\ThornG/\psi)^\zeta] = 1 may fail to have any root (it does fail at ψ1\psi \equiv 1) while the primal equation (26) always has one — the level R>1\Rcal > 1 is what saves it. The full kinship map (which constructions could and could not transport a distributional-tail theorem to the primal gap) is alt_proof_econlit.md §6; the one live bridge candidate (a geometrically-killed running maximum, nearest published statement Beare–Seo–Toda 2022) is recorded there with its four blocking gaps. Section 8.1 How rarely is the tail visited? The dual (Kesten) root develops the dual root explicitly — with the estimated calibrations’ numbers — for readers meeting it for the first time.

(Mellin duality, in one sentence: the Mellin transform is the change of variables to log wealth, which turns multiplying wealth by a factor into adding a step — under it the two literatures’ root equations are the same equation seen from the two sides, which is why one side can have a closed form while the other does not.)

3.3 The near-miss frontier: what the literature proves, what is new

The buffer-stock/income-fluctuation canon comes remarkably close to this theorem without stating it. (Throughout, lemmas of this program are cited by their source labels — L0, L1, …, with primes marking strengthened variants such as L2′; γ-prefixed labels belong to the compactified engine of §5.) What is published:

known resultsourcewhat it giveswhat it lacks
ratio limit c(m)/mκ\cFunc(m)/m \to \kapBST (Limiting MPC Bounds); MST (2020) Prop 2.5; Ma–Toda (2021) Thm 2.3; Ma–Toda (2022) Thm 2.6x(m)=o(m)\psavFunc(m) = o(m)no level, no rate
derivative limit c(m)κc'(m) \downarrow \kapBST (in-text squeeze; not a numbered theorem — flagged)xx eventually monotoneno rate
rich-end growth Et[ct+1/ct]Þ\E_{\prdt}[\cNrmNxt/\cFunc\Now] \to \ThornBST (Asymptotic Consumption Growth Factors)behavior \to PF benchmarka limit, not a rate
sign: approach from belowCarroll–Kimball (1996) concavity; Leland (1968)/Sandmo (1970)x0x \geq 0
level convergence x0x \to 0nowhere (verified sweep, 2026-07-07)proven here (L3; §1’s imported ladder)
the rate min(1,q)\min(1,\qup) + amplitudesnowherethe theorem

Two hygiene notes from the verified survey: Ma–Toda (2022) states explicitly that it does not characterize the rate; and the Euler-iteration existence canon (Li–Stachurski 2014, MST 2020 Assumption 2.3, Ma–Toda 2021) requires E[u(Y)]<\E[\uFunc'(Y)] < \infty, which fails at the zero-income atom (CRRA E[θρ]=\Rightarrow \E[\theta^{-\rho}] = \infty when >0\pZero > 0) — so existence for the model with the atom routes through BST’s WRIC/FVAC machinery, not through that canon (§1).

3.4 The forcing is the Kimball precautionary premium (KPP)

The engine below runs on a one-step recursion whose inhomogeneous term J(wˉ)\bufAdjFunc(\wbar) is, to leading order, κ(ρ+1)σ2/(2ÞΓ)wˉ1\kap(\rho+1)\sigma^2/(2\ThornG) \cdot \wbar^{-1} — economically, (κ\kap times) the Kimball precautionary premium (KPP) — the precautionary counterpart (Kimball 1990) of the Arrow–Pratt risk premium — for the coming period’s normalized income noise, evaluated at PF-wealth scale wˉ\wbar: Kimball’s (1990) form ψKησ2/2\psi_K \approx \eta\,\sigma^2/2 with CRRA relative prudence ηc=ρ+1\eta \cdot c = \rho+1. Its sign in this region (the gap is positive: uncertainty depresses consumption below the PF line) is Leland (1968)/Sandmo (1970); its small-risk quadratic template is Pratt (1964); its consumption interpretation is Drèze–Modigliani (1972).

4. Rigor engine I: the induction invariant

Every step of the spine leaned on a leading-order recursion; the four engines now make it rigorous, and the first is the induction invariant itself: the exact form of (EQ-1). Lemma 5.1 rewrites the Euler equation exactly in gap form — (EQ-1) ↦ (29) — with the spine’s leading-order adjustment replaced by the exact forcing J(wˉ)\bufAdjFunc(\wbar) (two-sided bounds; limit (30) == the cJc_J of §2.3), and with everything §2.5’s table drops carried honestly as the certified remainder weights ηL,ηR=O(1/wˉ)\etaL, \etaR = O(1/\wbar). Corollary 5.2 is the sandwich that makes §2.4’s unrolling honest: its two branches (31)/(32) are (EQ-1) with the stochastic PF-wealth step bracketed by its worst-case displacement ±C0\pm C_0, so the backward iteration proceeds without knowing the realized path. Both are PROVEN and review-hardened in stage_A_proof.md §§1–5 / stage_B_proof.md §§B1–B3 (statements only, completing §1’s imported ladder):

Economic reading of (31)(32): today’s gap = (discounted at 1/R1/\Rcal) tomorrow’s gap at the deterministically-shrunk PF wealth ÞΓwˉ\ThornG \wbar, plus the discounted adjustment R1J(wˉ)\Rcal^{-1}\bufAdjFunc(\wbar), the coming period’s risk priced today — the buffer-stock journey, one rung at a time.

Formally, the sandwich (31)(32) strongly resembles a pricing recursion — shortfall = discount × continuation shortfall + adjustment, with the Kimball adjustment as the flow (§0’s buffer-stock reading). Readers who internalize this one line can predict every result in §§5–8 by asking the corresponding asset-pricing question: does the perpetuity converge? what happens at r=gr = g? which adjustments carry the value?

5. Rigor engine II: the unrolling made uniform — the compactified boundary fixed point

§2.4 stepped backward one period at a time and truncated the sum at the log-clock exit; making that honest for the infinite-horizon policy — uniformly in the number of backward steps, with the stochastic drift and the η\eta-remainders live — is this engine’s job. Its device is a change of coordinates under which the whole backward journey fits on one compact picture: the shells of §5.2 are §2.4’s backward steps seen in the state space (one shell = one rung of the descent = one tick of the log-clock), and the boundary point z=0z = 0 is “TT \to -\infty seen from the state space” — the place where infinitely many backward steps accumulate, compactified so that “take the limit” becomes “evaluate at the boundary”. The scalar equation at that boundary (§5.4) is §2.4’s geometric sum, resummed as a fixed point.

(Proof bodies: alt_proof_compactified.md §§1–3, refuter-reviewed (RF1/RF2, zero broken); labels below are that document’s. One notational bridge: that document writes the trial exponent as ss; here it is q~\qtilde (and precautionary saving is xx, §0).)

5.1 The coordinate, and “assuming the power law in the limit” as a choice of units

Map z:=1/wˉz := 1/\wbar. The half-line [wˉ0,)[\wbar_0, \infty) becomes (0,zˉ](0, \bar{z}], and infinite wealth becomes the finite boundary point z=0z = 0 — adjoin it. (This is the one-dimensional Poincaré-compactification chart of dynamical systems; Dumortier–Llibre–Artés 2006, ch. 5 — decorative ancestry, no theorem of theirs is consumed.) Then choose units that embed the conjectured power law: the compensated gap

(The idea in one line: adjoin ‘wealth equals infinity’ to the state space as an honest boundary point, so that ‘take the limit as wealth grows’ becomes ‘evaluate at the boundary point’ — statements about limits become statements about continuity at a place the analysis can stand on.)

Wq~(z):=wˉq~g(wˉ),q~=the compensation (trial) exponent.W_{\qtilde}(z) := \wbar^{\qtilde}\, g(\wbar), \qquad \qtilde = \text{the compensation (trial) exponent}.

“Assume the power law holds in the limit” is thereby demoted from a hypothesis to a choice of coordinates: nothing is assumed — §6’s Theorem II (6. Rigor engine III: the trichotomy certified — Theorem II, resonance, and the q<1\qup < 1 sharp forms) shows the compactified problem itself detects the right exponent (Wq~W_{\qtilde} has a finite nonzero boundary behavior only at q~=min(1,q)\qtilde = \min(1, \qup)). Figure 2 gives the geometry.

Left, the wealth-space ladder marching down toward the return region around the target; right, the compactified interval with shells accumulating at the boundary point z equals zero, arrows pointing away from the boundary.

Figure 2:The ladder and the compactified shells. Left: in wealth space the recursion marches down toward the return region around the target, providing an adjustment σ2/wˉ\asymp \sigma^2/\wbar per rung. Right: in z=1/wˉz = 1/\wbar the shells accumulate at the boundary point z=0z=0 (“PF total wealth == \infty”); the referred point zz/ÞΓz \mapsto z/\ThornG points away from the boundary, so the recursion at shell n+1n{+}1 always reads off the already-controlled shell nn.

How to read this figure. Left: economics. A household richer than its buffer-stock target saves less than income-plus-interest, so its wealth-to-income ratio drifts down about þg-\GPRte percent per period (the GIC drift); each rung of the descent calls for an adjustment proportional to 1/(PF wealth)1/\text{(PF wealth)}. Right: the same journey after the change of variables z=1/wˉz = 1/\wbar, which puts “infinitely wealthy” ON the page as the point z=0z = 0 — the way a phase diagram puts a steady state on the page. The arrows show why the argument is not circular: the recursion always maps a point toward the region where the solution is already controlled, so information flows from known territory to the boundary, never from an assumption about the boundary.

5.2 The shell equation

Partition (0,zˉ](0, \bar{z}] into geometric shells SnS_n (log-width þg-\GPRte), accumulating at the boundary; enlarge each by the constant ζ:=C0/(1ÞΓ)\zeta := C_0/(1-\ThornG) — the exact fixed point of the one-step displacement map (the unique ζ\zeta with ÞΓζ+C0=ζ\ThornG\zeta + C_0 = \zeta) — so the referred point of every wˉ\wbar in shell n+1n{+}1 lands exactly in shell nn (Lemma γ0; no band slack, no self-reference: the map zz/ÞΓz \mapsto z/\ThornG points away from the boundary, so the recursion at shell n+1n{+}1 always reads off the already-controlled shell nn). In these units the imported sandwich (31)(32) becomes the shell equation (Lemma γ1).

(‘Referred point’ means: the wealth level the one-step map sends the current level to. The enlargement constant is chosen precisely so that this image always lands a full shell down — never straddling a shell boundary — and that non-straddling is what lets the induction walk down one shell at a time without circular reasoning.)

For compensation q~\qtilde, its upper and lower branches read

Wq~(wˉ)(1+D^wˉ)[λ(q~)Wq~(ÞΓwˉC0)+Fq~(wˉ)]W_{\qtilde}(\wbar) \leq \Bigl(1 + \tfrac{\hat{D}}{\wbar}\Bigr)\bigl[\, \lambda(\qtilde)\, W_{\qtilde}(\ThornG \wbar - C_0) + F_{\qtilde}(\wbar) \,\bigr]
Wq~(wˉ)(1D^wˉ)[λ(q~)Wq~(ÞΓwˉ+C0)+Fq~(wˉ)],W_{\qtilde}(\wbar) \geq \Bigl(1 - \tfrac{\hat{D}}{\wbar}\Bigr)\bigl[\, \lambda(\qtilde)\, W_{\qtilde}(\ThornG \wbar + C_0) + F_{\qtilde}(\wbar) \,\bigr],

with boundary-continuous coefficients: the multiplier λ(q~)=ÞΓqq~\lambda(\qtilde) = \ThornG^{\,\qup - \qtilde} (at ψ1\psi \equiv 1 — an algebraic identity via R=ÞΓq\Rcal = \ThornG^{-\qup}) and the forcing Fq~(wˉ)=wˉq~J(wˉ)/R(1+O(1/wˉ))F_{\qtilde}(\wbar) = \wbar^{\qtilde} \bufAdjFunc(\wbar)/\Rcal \cdot (1 + O(1/\wbar)), whose boundary limit at q~=1\qtilde = 1 is F(0):=cJ/RF(0) := c_J/\Rcal by (30), consumed rate-free. There is also the exact identity

(1+ηL)RW1(wˉ)=(1+ηR)Et[wˉwˉt+1W1(wˉt+1)]+wˉJ(wˉ),(1 + \etaL)\,\Rcal\, W_1(\wbar) = (1 + \etaR)\,\E_{\prdt}\Bigl[\tfrac{\wbar}{\wbarNxt}\, W_1(\wbarNxt)\Bigr] + \wbar \bufAdjFunc(\wbar),

which is never used to produce a limit — only to evaluate one after it exists.

5.3 The two core lemmas (the entire engine — 76 lines in the source)

Lemma γ2 (boundedness). W1W_1 is bounded on the compactified domain. Idea: iterate the sup of (35) inward shell by shell; with λ(1)=1/(RÞΓ)<1\lambda(1) = 1/(\Rcal\ThornG) < 1 and bounded forcing, the sup recursion Mn+1λ(1+δn)Mn+Fˉ(1+δn)M_{n+1} \leq \lambda(1+\delta_n)M_n + \bar{F}(1+\delta_n) with δn<\sum\delta_n < \infty is trapped by max(anchor,Fˉ/(1λ))\max(\text{anchor}, \approx\bar{F}/(1-\lambda)). Self-contained — Theorem A1 is not imported.

(Intuition for boundedness: a contraction factor strictly below one means each shell’s worst case is pulled toward the injection’s scale instead of compounding on itself; the summability of the perturbations says their total lifetime effect is finite — too weak to undo the contraction.)

Lemma γ3 (boundary stability — the engine). For two-sided shell recursions Wn+1λnWn+FnW_{n+1} \lessgtr \lambda_n W_n + F_n with λnλ<1\lambda_n \to \lambda < 1 and FnFF_n \to F_\infty: both lim sup\limsup and lim inf\liminf of WnW_n equal F/(1λ)F_\infty/(1-\lambda). Proof (the four lines): let M:=lim supMnM := \limsup M_n; passing to the limsup in the sup-recursion gives MλM+FM \leq \lambda M + F_\infty, so MF/(1λ)M \leq F_\infty/(1-\lambda); the mirrored inf recursion gives mF/(1λ)m \geq F_\infty/(1-\lambda); mMm \leq M closes it. \blacksquare

(The proof in words: the recursion says each shell’s worst case is at most a contraction factor strictly below one times the previous shell’s worst case, plus a bounded injection. Worst cases therefore cannot ratchet upward: in the limit the largest recurring value M satisfies M ≤ λM + F̄, which pins M ≤ F̄/(1−λ). Running the same argument from below on the best cases gives a matching floor, and the two one-sided bounds meet at the Gordon value.)

5.4 Theorem I (the owner’s one number): the boundary value at q>1\qup > 1

Proof (assembled). γ2 gives boundedness; γ0/γ1 put the shell sequences in γ3’s hypotheses (λnλ<1\lambda_n \to \lambda < 1 by boundary continuity of the coefficients, FnF(0)F_n \to F(0) by (30), the (1±D^/wˉ)(1 \pm \hat{D}/\wbar) factors are the δn\delta_n with δn<\sum \delta_n < \infty because 1/wˉn\sum 1/\wbar_n is geometric); γ3 forces lim sup=lim inf=F(0)/(1λ)\limsup = \liminf = F(0)/(1-\lambda), i.e. the limit exists. Only then is (37) passed to the limit (dominated convergence; ηL,ηR0\etaL, \etaR \to 0, wˉ/wˉt+1ÞΓ1\wbar/\wbarNxt \to \ThornG^{-1} uniformly, W1(wˉt+1)W(0)W_1(\wbarNxt) \to W(0) a.s.-and-bounded), which yields (38) — the order matters, and is what makes the argument non-circular (refuter-verified explicitly, RF1 §audit-1). \blacksquare

(Unpacking the chain for readers who have not memorized the labels: γ2 says the compensated gap never explodes on any shell — its worst cases stay bounded; γ0 and γ1 are housekeeping — the injections converge and the weights settle; γ3 then converts ‘bounded, with settling inputs’ into an actual limit; and the limit is the Gordon value. Each lemma is stated and proved just above.)

This is the promised punchline: at the boundary point “PF wealth == \infty” the entire Euler equation degenerates to one linear scalar equation, and the finite number representing the exact limit is its unique solution. The unrolled form B=k0λkF(0)B = \sum_{k \geq 0} \lambda^k F(0) is the geometric resolvent — the same number every other proof in the program computes by its own route (§10). (Rigor record: the source suite measures the deep-tail match to the closed form at 0.9% with the predicted approach rate wˉ(q1)\wbar^{-(\qup-1)} fitted to 0.597 vs 0.6.) Figure 3 shows both faces: panel (a) the GIC-CAP calibration (q=1.474\qup = 1.474, Stage B), where the boundary value Bψ=357B_\psi = 357 exists and the boundary equation pins it — but the approach outruns reachable PF wealth (W/Bψ=0.42W/B_\psi = 0.42 at the top of the estimation grid, 0.91 even at 130×130\times the grid top; the knife-edge window of §9); panel (b) a labeled illustration calibration where the same convergence completes in-window (match 1.3%; inset slope -0.598 vs -0.6).

The Gordon-growth reading of BB is §2.4’s unrolling, resummed: adjustment growth 1/ÞΓ1/\ThornG per rung of the descent, discounting 1/R1/\Rcal per period, convergence exactly when RÞΓ>1\Rcal\ThornG > 1 — the Gordon condition “discount rate exceeds growth rate”, which is q>1\qup > 1 (case (a) of §2.4, (21)) — and B=cJÞΓ/(RÞΓ1)B = c_J\ThornG/(\Rcal\ThornG - 1): the exact closed form of Theorem I, reached there by backward accumulation and here as a boundary fixed point. The entire theorem can be read as: the precautionary shortfall of the wealthy is the Gordon value of their future adjustments.

Two panels of the compensated gap versus the compactified coordinate; left the GIC-CAP calibration en route to a boundary value far above the computational grid; right an illustration calibration completing its convergence to B.

Figure 3:The compactified boundary value. (a) The GIC-CAP calibration: the boundary equation pins W(0)=Bψ=357W(0) = B_\psi = 357, but at þg=3.6104-\GPRte = 3.6\cdot 10^{-4} the approach outruns reachable PF wealth. (b) An illustration calibration (ÞΓ=0.9\ThornG = 0.9, q=1.6\qup = 1.6; not HAFiscal): the same convergence, completed in-window.

How to read this figure. The vertical axis is PF wealth × shortfall — call it the amplitude. If the shortfall truly behaves like B/wˉB/\wbar, this product levels off at the constant BB; moving leftward along the horizontal axis (z=1/wˉ0z = 1/\wbar \to 0) means growing richer. Panel (b), the labeled illustration economy: the amplitude does level off, at exactly the closed-form perpetuity value — the marked dot at the boundary. Panel (a), the patience-ceiling calibration GIC-CAP: the same leveling is provably underway, but its destination (Bψ=357B_\psi = 357) lies far above anything visible on the computational grid (shaded region) — at the grid top the amplitude has covered only 42% of the way. Practical moral: on any feasible grid, a fitted constant would badly understate BψB_\psi; the top-of-grid boundary condition should use the closed form, not a fit.

6. Rigor engine III: the trichotomy certified — Theorem II, resonance, and the q<1\qup < 1 sharp forms

§2.4 classified the unrolled sum by its per-step ratio 1/(RÞΓ)1/(\Rcal\ThornG); Theorem II certifies that classification at the boundary, including the two members the spine could only sketch — the resonance constant, and the q<1\qup < 1 sharp forms in which §2.4(c)'s backward lattice {ÞΓjwˉ}\{\ThornG^{\,j}\wbar\} becomes rigorous shell/renewal bookkeeping.

Test a trial decay exponent q~\qtilde against the recursion: a xwˉq~x \propto \wbar^{-\qtilde} mode is multiplied, per backward step, by the discount 1/R1/\Rcal and by the drift-rescaling ÞΓq~\ThornG^{-\qtilde} (the state contracts to ÞΓwˉ\ThornG\wbar, so a steeper trial exponent harvests a larger rescaling). The two effects exactly offset when q~=q\qtilde = \qup — that balance is what defines q\qup — and for any other trial exponent the mismatch compounds at the rate λ(q~)=ÞΓqq~\lambda(\qtilde) = \ThornG^{\,\qup-\qtilde} per step: modes shallower than q\qup die out under backward iteration (λ<1\lambda < 1), modes steeper explode (λ>1\lambda > 1), which is why the fixed point selects q\qup (general ψ\psi: λB(q~)=eL(q~)\lambda_B(\qtilde) = e^{\Lfun(\qtilde)} — so “boundary multiplier =1= 1” is literally Lemma A5’s eigen-equation (26)). Everything in the theorem is the elementary classification of the shell recursion by (λ,limFn)(\lambda, \lim F_n):

regimemultiplierforcingboundary behaviorresult
q>1\qup > 1, q~=1\qtilde=1λ<1\lambda < 1FncJ/R>0F_n \to c_J/\Rcal > 0contracting fixed point — one numberTheorem I: W(0)=BW(0) = B
q=1\qup = 1, q~=1\qtilde=1λ=1\lambda = 1FncJ/R>0F_n \to c_J/\Rcal > 0linear escape — a slope, not a valueTheorem γ-R
q<1\qup < 1, q~=q\qtilde=\qupλ=1\lambda = 1Fn0F_n \to 0 geometricallymultiplier exactly 1, fading forcing — a limit setγ-C

In Gordon-growth language the three rows are §2.4’s trichotomy, certified: the convergent perpetuity (a), the r=gr = g knife-edge whose linear accrual is the lnwˉ\ln \wbar law (b), and the divergent case (c), in which the “value” is dominated by the latest, largest adjustments — the journey’s end near the target, whose imprint is the periodic fine structure.

Resonance (Theorem γ-R; PROVEN). At q=1\qup = 1 (RÞΓ=1\Rcal\ThornG = 1), the shell recursion is Wn+1=Wn(1+O(δn))+FnW_{n+1} = W_n(1 + O(\delta_n)) + F_n with FncJ/RF_n \to c_J/\Rcal: by the perturbed Stolz–Cesàro lemma (γ-R1), Wn/ncJ/RW_n/n \to c_J/\Rcal; converting shells to PF wealth (n(wˉ)=ln(wˉ/wˉa)/(þg)n(\wbar) = \ln(\wbar/\wbar_a)/(-\GPRte)),

(Stolz–Cesàro is the discrete l’Hôpital rule: when the per-step increments of a sequence settle to a constant, the running average settles to the same constant. At resonance the discounting and the descent cancel exactly, so every backward step contributes the SAME undiscounted amount — after n steps the sum is about n times that amount, and since the journey takes about log-wealth-many steps, a factor of log wealth survives in the law. That is the entire origin of the knife-edge logarithm.)

wˉlnwˉg(wˉ)    cJR(þg)=κ(ρ+1)σ22þg\frac{\wbar}{\ln \wbar}\, g(\wbar) \;\to\; \frac{c_J}{\Rcal\,(-\GPRte)} = \frac{\kap(\rho+1)\sigma^2}{-2\GPRte}

— the sharp resonance constant of Theorem B-res’s ψ1\psi\equiv1 corollary, re-proven in a few lines (measured to 0.07% in the source suite). Figure 4(b) exhibits it at the resonance patience level of the College fundamentals (β\beta root-found so q=1.000000\qup = 1.000000 exactly, sitting between the most patient estimated College type (COL-TOP) and the patience ceiling (GIC-CAP)): WnW_n grows linearly in the tilted shell index with per-step increment cJB/R\to c_J^B/\Rcal (deepest measured increment 1.101 vs theory 1.113). This is a knife-edge, not a neighborhood: the theorem is a fixed-q\qup statement, and the crossover to it is non-uniform on any window with q1lnwˉ=O(1)|\qup - 1| \ln \wbar = O(1) (Remarks 3.4/4.1 of the source; the two deliberately-registered FAILs of the source’s falsifier suite live exactly there, adjudicated honest by two refuter panels).

q<1\qup < 1: the boundary is a circle, not a point (γ-C; envelope PROVEN, sharp form imported). At q~=q\qtilde = \qup the map zz/ÞΓz \mapsto z/\ThornG shifts log PF wealth by exactly þg-\GPRte, so it leaves unchanged the position within the multiplicative cycle: the fractional part φ={lnwˉ/(þg)}\varphi = \{\ln \wbar/(-\GPRte)\}. In log PF wealth the dynamics repeat with period þg-\GPRte, so what matters asymptotically is the position modulo þg-\GPRte — a circle of circumference þg-\GPRte rather than a single boundary point. PF-wealth levels sharing the same position in the cycle (equal log PF wealth modulo þg-\GPRte) form a family that the shell recursion maps to itself; on each such family the recursion neither contracts nor expands (its multiplier is exactly 1) while its forcing fades geometrically. Proven self-contained here: the two-sided envelope converges to a bounded positive limit set [m,M][m_\infty, M_\infty], 0<mM<0 < m_\infty \leq M_\infty < \infty (γ-C1), and the period-by-period smoothing of the position-in-the-cycle distribution contributed by the continuous part of the shock is geometrically summable (γ-C2). Upgrading to convergence separately at each position of the cycle needs a uniform bound on how much the limit can vary under small shifts of the position (a modulus of continuity in the cycle variable) — GAP-γ-equicont, honestly open for this engine (it is L8/L9 doubling territory); the sharp statement — wˉqg(wˉ)=P({lnwˉ/(þg)})+o(1)\wbar^{\qup} g(\wbar) = P(\{\ln \wbar/(-\GPRte)\}) + o(1), with PP the log-periodic factor, whose non-constancy is now directly measured (first resolved tone: amplitude 7.81047.8\cdot10^{-4} at þg=0.693-\GPRte = 0.693) — is Theorem A2 of the original chain, imported. Figure 4(c) makes the practical point sharply: at HS-mean’s estimated patience (þg=0.0125-\GPRte = 0.0125) the periodic component provably has no effect at any attainable precision: its total amplitude — the oscillation osc(P):=supPinfP\operatorname{osc}(P) := \sup P - \inf P of the log-periodic factor — satisfies osc(P)/Pˉeα/(þg)\operatorname{osc}(P)/\bar{P} \sim e^{-\alpha/(-\GPRte)}, far below any numerical floor. That is, the boundary circle collapses numerically to a point — which is exactly why a constant-coefficient power-law extrapolant works in the HAFiscal/HARK code; the inset shows the one calibration in the program with a numerically resolved tone (þg=0.69-\GPRte = 0.69, amplitude 8104\approx 8\cdot10^{-4}, correlation 0.9999 between the measured wiggle and its predicted position in the cycle) — the honest evidence that the circle is real.

(Where the circle comes from: below the resonance the journey’s endgame depends on where within one descent-step the log-clock lands — its fractional part — and positions that differ by a whole step are equivalent. ‘Position modulo one step’ is a point on a circle, and the log-periodic wobble in the law is the shadow that circle casts.)

(Plainly, what is missing: a proof that the wobble of the compensated gap within a single shell cannot speed up without bound as the shells accumulate — ‘equicontinuity’ is the analysts’ name for that uniform-wobble bound, and the cited doubling estimates are the standard tool that delivers it. No economic content is in doubt; only this uniformity step is open for this engine.)

Three panels: the GIC-CAP calibration climbing toward its boundary value; the exact-resonance calibration growing linearly with a Cesaro inset; the HS-mean cycle-position plot, flat, with a resolved-tone inset.

Figure 4:The trichotomy at the boundary. (a) q>1\qup > 1 patience ceiling (GIC-CAP): WnBψW_n \nearrow B_\psi, en route. (b) q=1\qup = 1 exactly (AUX-RES): linear growth, Wn/ncJB/RW_n/n \to c_J^B/\Rcal. (c) q<1\qup < 1 (HS-mean): the periodic component is undetectably small at the estimated calibration’s þg-\GPRte; inset — the resolved tone at þg=0.69-\GPRte = 0.69.

How to read this figure. The three Gordon cases in motion. (a) Convergent case (q>1\qup > 1): the amplitude climbs toward its perpetuity value — a plateau (here honestly en route; see Figure 3a). (b) The knife-edge r=gr = g case (q=1\qup = 1 exactly): no plateau exists; instead value accrues linearly with the rungs of the descent — just as a growing annuity with discount rate equal to growth rate is worth (number of periods) × (per-period flow) — and the slope is the resonance constant (inset: the running average converging to it). This linear accrual is precisely where the extra lnwˉ\ln \wbar in the theorem comes from. (c) The q<1\qup < 1 case, plotted against the fractional part of log PF wealth: theory says the limit is a periodic function on this axis (a circle, §6), but at HS-mean’s estimated patience the predicted oscillation is eα/(þg)\sim e^{-\alpha/(-\GPRte)} — astronomically below any solver’s resolution — so the limit is flat for every practical purpose. The inset shows the one deliberately coarse calibration where the wiggle becomes visible: the circle is real, just irrelevant at realistic patience.

Wrong-exponent detection (Theorem γ-T; PROVEN). Compensate with any other q~\qtilde and the boundary behavior degenerates at the exact geometric per-shell rates ÞΓmin(1,q)q~\ThornG^{\,\min(1,\qup)-\qtilde}: q~<min(1,q)Wq~0\qtilde < \min(1,\qup) \Rightarrow W_{\qtilde} \to 0; q~>\qtilde > the true exponent Wq~\Rightarrow W_{\qtilde} \to \infty. Two dials govern min(1,)\min(1,\cdot): the multiplier λ(q~)\lambda(\qtilde) crosses 1 at q~=q\qtilde = \qup, and the forcing wˉq~J(wˉ)\wbar^{\qtilde}\bufAdjFunc(\wbar) (wˉq~1\asymp \wbar^{\qtilde-1}) stops fading at q~=1\qtilde = 1 — so the only compensation with finite nonzero boundary behavior is q~=min(1,q)\qtilde = \min(1, \qup). The compactified problem finds the exponent by itself (Figure 5, on HS-mean: trial exponents q0.15\qup - 0.15 / q+0.15\qup + 0.15 drift at measured log-slopes -0.144 / +0.156 against theory -0.15 / +0.15, while q~=q\qtilde = \qup stays flat at +0.006; detection rates verified to ±0.012(þg)\pm 0.012\cdot(-\GPRte) per shell in the source suite, and reconfirmed on fresh calibrations by the RF2 panel). This doubles as the numerical diagnostic of record: the grid-depth migration of the fitted local exponent toward min(1,q)\min(1,\qup) is a theorem-backed signature, not a stylized fact.

The compensated gap for three trial exponents on log-log axes; only the true exponent stays flat, the low trial decays, the high trial grows.

Figure 5:The compactified problem detects the exponent (HS-mean): only q~=min(1,q)\qtilde = \min(1,\qup) gives a finite nonzero boundary value.

How to read this figure. The closest econometric analogy is choosing the right detrending. Multiply the shortfall by a trial power wˉq~\wbar^{\qtilde} and watch what happens: too small a q~\qtilde and the product still dies away (under-compensated); too large and it explodes (over-compensated); exactly right and it levels off. Only q~=min(1,q)\qtilde = \min(1, \qup) levels off — the solved model picks its own exponent, with no fitting involved. This is also §0’s saddle-path selection made visual: wrong exponents put you on explosive or degenerate solution branches.

7. Rigor engine IV: permanent shocks — Theorem III and the human-wealth-revaluation channel

(Proof body: alt_proof_compactified.md §7 (γ-B1/γ-B2/γ-B3, Theorem γ-B); Stage-B foundations imported from stage_B_proof.md §§B1–B3.)

With permanent shocks the exact recursion becomes the master identity

ψt+1wˉt+1=ÞΓwˉ+Wt+1+Rg(wˉ),W:=ψ(θ+h1)h,E[W]=0,\permShkNxt\, \wbarNxt = \ThornG\, \wbar + W_{\prdNxt} + \Rcal\, g(\wbar), \qquad W := \psi(\theta + h - 1) - h, \quad \E[W] = 0,
σB2:=Var(W)=E[ψ2]σθ2+h2σψ2.\sigma_B^2 := \operatorname{Var}(W) = \E[\psi^2]\,\sigma_\theta^2 + h^2\,\sigma_\psi^2.

(Watch the typographical neighbors: lowercase wˉ\wbar is PF total wealth, capital Wt+1W_{\prdNxt} the mean-zero income-revaluation shock — and neither is the compensated gap Wq~()W_{\qtilde}(\cdot) of §5.) The h2σψ2h^2\sigma_\psi^2 term is the human-wealth revaluation channel: a permanent shock reprices the entire future income stream — all hh units of it — so the precautionary forcing is larger than in the transitory-only model. The ψ\psi-weighted one-step relation (B5.2) makes the compactified multiplier λB=E[ψ2]/(RÞΓ)=eL(1)<1    q>1\lambda_B = \E[\psi^2]/(\Rcal\ThornG) = e^{\Lfun(1)} < 1 \iff \qup > 1, and the identical γ2/γ3 pattern (boundedness via the stopped expectation-unroll Lemma γ-B2, which covers the general bounded-ψ\psi case including the empirically relevant ψminÞΓ\psi_{\min} \leq \ThornG; key identity E[ψ2(ÞΓ/ψ)]=ÞΓ\E[\psi^2 (\ThornG/\psi)] = \ThornG) delivers:

(In words: unroll the expectation of the recursion only up to the random time the state first exits the good region — ‘stopping’ there is legitimate because the displayed identity says the drift is exactly right on average — and boundedness follows without ever tracking an individual path.)

The spine reaches the same constants. §2.6’s ψ\psi-general one-step subtraction (the W3 spike, eq1_backward_psi_spike.md; battery verify_eq1_psi_spike_checks.py) derives the mode condition E[ψ1+q]=RÞΓq\E[\psi^{1+q}] = \Rcal\,\ThornG^{\,q} — the eigen-equation (26) — and the amplitude (28) directly from the backward telling, with the discount still exactly 1/R1/\Rcal by the (K3) feedback invariance; Theorem III is the rigorous certificate of that construction at q>1\qup > 1.

Stage-B resonance and q<1\qup < 1 are OPEN for this engine (the shell index becomes a random walk whose occupation control is exactly the original chain’s tilt toolkit); the sharp results are Theorems B-res, B2, B2-arith of stage_B_proof.md, imported.

(The obstruction, plainly: with permanent shocks the shell index no longer marches down deterministically — each period it moves by a random amount, a random walk with downward drift. Controlling how much time such a walk spends at each level — its ‘occupation’ — is exactly what the exponential-tilting toolkit does for the original Kesten chain; porting that control is the outstanding engineering, not any new economics.)

8. The journey, measured: the economics of the proof

§2.4’s unrolling stopped at the log-clock exit — the backward steps end where the descent reaches the buffer-stock body. Read forward (route α, alt_proof_econlit.md, which prices the same numbers as a path sum), the gap at PF wealth wˉ\wbar is the expected present discounted value of the Kimball precautionary premium (KPP)s the household will still provide along its buffer-stock journey from wˉ\wbar down toward the return region around the target m^\hat{m} — the precise object behind the phrase “toward the target”, and the destination this section defines before measuring where along the journey the gap is earned. In BST’s language the destination is the individual target m^\hat{m} (unique under WRIC + FVAC + GIC-Mod per BST’s buffer-stock-target theorem, and computable from BST’s implicit equation (m^c(m^))RE[ψ1]+1=m^(\hat{m} - \cFunc(\hat{m}))\,\Rcal\,\E[\psi^{-1}] + 1 = \hat{m}) or the pseudo-target mˇ\check{m} (GIC-Raw). What the unrolled argument itself consumes is only this: a compact wealth interval containing the target, on which the descent’s bookkeeping can stop. Concretely (Stage A), take the PF-wealth levels wˉζ\wbar \leq \zeta with ζ:=C0/(1ÞΓ)\zeta := C_0/(1-\ThornG) — §5.2’s enlargement constant: above ζ\zeta the recursion (8) moves strictly down almost surely (wˉt+1ÞΓwˉ+C0<wˉ\wbarNxt \leq \ThornG\wbar + C_0 < \wbar), from inside the interval PF wealth cannot exit upward (ÞΓζ+C0=ζ\ThornG\zeta + C_0 = \zeta), the descent from any wˉ\wbar enters it in ln(wˉ/ζ)/(þg)\approx \ln(\wbar/\zeta)/(-\GPRte) periods, and the target must lie inside (above ζ\zeta PF wealth falls almost surely, contradicting Et[mt+1/m]=1\E_{\prdt}[\mNrmNxt/m] = 1 at m^\hat{m}). The region is compact because the transitory shock has bounded support — part of the maintained income-process assumption (BST, Assumption assn-shocks; §1): bounded θ\theta, with the sandwich xxˉx \leq \bar{x}, is what makes the one-step displacement bound C0C_0 — hence the upper edge ζ\zeta — finite and computable from the primitives (C0=max(1,θmax1+Rxˉ)C_0 = \max(1,\, \theta_{\max} - 1 + \Rcal\bar{x}) works, read off (8)). With unbounded transitory shocks the region is not compact — a single large draw passes any point, no finite C0C_0 exists — and the renewal argument as given needs the bound; the decay theorem itself still holds via the L4′ moment extension, which replaces the compact-interval bookkeeping with moment bounds (8.1 How rarely is the tail visited? The dual (Kesten) root, item 3). With permanent shocks the interval is entered and re-entered rather than absorbing (8.1 How rarely is the tail visited? The dual (Kesten) root, item 2 quantifies the excursions above it). The journey picture is exactly Carroll’s (1997) narrative made quantitative. Where along the journey is the gap earned? §2.4 answered in principle: the head of the journey at q>1\qup > 1, every rung equally at resonance, the journey’s end at q<1\qup < 1. Figure 6 shows it measured. The moral for this project: the descent spends ~89% of its lifetime inside any sensibly-sized grid — the far field is visited too rarely to deserve gridpoints, and what happens there is exactly what the theory prices. An in-grid representation plus the theorem’s tail is the whole solution; no ad-hoc external criterion (a grid top, an extrapolation form, a tolerance) does any work.

Measured on the estimated calibrations (the bookkeeping is the ψ2\psi^2-tilted Wald chain — step L(1)\Lfun'(1), weight λB\lambda_B per rung — which the figure build verified is the only unroll-consistent accounting at estimated patience): the journey’s-end/boundary share is 89% for HS-mean (q=0.38\qup = 0.38, end-dominated), 50% for COL-TOP (q=0.69\qup = 0.69, the transition), and 2% for the ceiling GIC-CAP (q=1.47\qup = 1.47, head-dominated) — with all three per-rung weights λB\lambda_B within 1% of 1: estimated patience sits close to the resonance knife-edge, which is the deep reason the historical exponential top-of-grid heuristic was locally harmless even though it is asymptotically wrong (the gap cannot fade faster than the 1/wˉ1/\wbar forcing floor — stage-A Lemma 6.1 / stage-B Lemma B-6.1).

Per-rung contribution shares along the journey for the three estimated calibrations, showing end-domination, near-flatness, and head-domination.

Figure 6:Where the gap is earned along the journey (the three estimated calibrations; ψ²-tilted Wald bookkeeping). Each panel unrolls the pricing recursion rung by rung (one rung = one model period) and plots each rung’s share of the summed rung values; the terminal boundary term is excluded from the normalization and its share is reported separately inside each panel. The moral for this project: the descent spends ~89% of its lifetime inside any sensibly-sized grid — the far field is visited too rarely to deserve gridpoints, and what happens there is exactly what the theory prices. An in-grid representation plus the theorem’s tail is the whole solution; no ad-hoc external criterion (a grid top, an extrapolation form, a tolerance) does any work.

How to read this figure (the “when is the shortfall earned?” decomposition). Recall the buffer-stock reading (§0): the shortfall a wealthy household exhibits today is the present value of the stream of premia it will still provide along its descent toward the target. This figure decomposes that present value by when the premia are provided. The horizontal axis is progress along the descent, rung by rung, rescaled to [0,1][0,1]: 0 = the first rung — today, at PF wealth wˉ=105\wbar = 10^5; 1 = the last rung providing an adjustment, just above the unroll’s cutoff at wˉ=600\wbar = 600 — well inside the approach to the target but still above the return region around m^\hat{m} (8.1 How rarely is the tail visited? The dual (Kesten) root quantifies how rarely the lowest rungs are ever visited). The height is that rung’s share of the summed rung values — the plotted shares sum to one — while the terminal boundary term (the continuation value parked at the cutoff) is excluded from the normalization and reported separately in each panel. Three readings:

The punchline for HAFiscal: all three estimated-calibration profiles are close to flat, because the estimated economies sit near the r=gr = g knife-edge (above, the λB\lambda_B-within-1%-of-1 measurement). That single fact is the deepest reason top-of-grid behavior is delicate at realistic patience: the value of the shortfall is spread almost evenly over the whole journey, so no part of the wealth domain is safely negligible.

8.1 How rarely is the tail visited? The dual (Kesten) root

The theorem describes c(m)\cFunc(m) at wealth levels the agent rarely — in one important case never — holds. This subsection quantifies “rarely”, because the correct quantifier is itself an object worth knowing: the tail exponent of a Kesten process. (Evidence for every number below: verify_dualroot_checks.py, ALL PASS, same 7-atom ψ\psi discretization as the figures.)

Economists have met this machine under another name: random growth. Multiply a variable by an i.i.d. factor each period (plus a small additive perturbation); if the factor shrinks on average but sometimes exceeds one, the stationary distribution develops a Pareto upper tail whose exponent solves E[factorζ]=1\E[\text{factor}^\zeta] = 1. This is the mechanism behind Zipf’s law for cities — Gabaix (1999) shows that Gibrat-style proportional random growth of city populations, with a small friction keeping cities from vanishing, produces a power-law size distribution, with the Zipf exponent 1 emerging in the small-friction limit — and the standard account of Pareto tails in wealth and firm size. Gabaix (2009) surveys both the empirical catalog (income and wealth, cities, firms, stock returns, trading volume, executive pay) and the generating theory, including precisely the exponent equation E[Aζ]=1\E[A^\zeta] = 1 below; Gabaix (2016) is the introduction written for the general economist, arguing these are true quantitative laws — empirically established and theoretically understood. Everything below is that familiar machine, applied to a single household’s normalized perfect-foresight wealth rather than to a cross-section of cities or firms.

Rearranging the master identity (41), normalized perfect-foresight wealth follows

wˉt+1=At+1wˉ+Bt+1,At+1:=ÞΓψt+1,Bt+1:=Wt+1+Rg(wˉ)ψt+1  (bounded: xxˉ),\wbarNxt = A_{\prdNxt}\, \wbar + B_{\prdNxt}, \qquad A_{\prdNxt} := \frac{\ThornG}{\permShkNxt}, \qquad B_{\prdNxt} := \frac{W_{\prdNxt} + \Rcal\, g(\wbar)}{\permShkNxt} \ \ (\text{bounded: } x \leq \bar{x}),

a stochastic recursion with i.i.d. multiplier — the workhorse of the wealth-distribution literature (Kesten 1973; Goldie 1991; econ packagings Benhabib–Bisin–Zhu 2011, Toda 2014; the Gouin-Bonenfant–Toda 2023 toolkit’s getZeta computes exactly the object below). The dual (Kesten) root is the unique ζ>0\zeta^* > 0 solving

E[Aζ]  =  E[(ÞΓ/ψ)ζ]  =  1,\E[A^{\zeta^*}] \;=\; \E\Bigl[\bigl(\ThornG/\psi\bigr)^{\zeta^*}\Bigr] \;=\; 1,

which exists precisely when the multiplier (i) contracts on average, E[lnA]<0\E[\ln A] < 0, and (ii) occasionally expands, Pr(A>1)>0\Pr(A > 1) > 0. When it exists, the stationary distribution of wˉ\wbar has a power tail with exponent ζ\zeta^*: Pr(wˉ>t)\Pr(\wbar > t) decays essentially like tζt^{-\zeta^*} (the upper bound is an elementary Chernoff/supermartingale estimate; the matching lower bound is the Kesten–Goldie theorem, with the usual non-arithmeticity caveat — a discretized ψ\psi technically sits in the lattice case, where Kevei’s 2017 version supplies the log-periodic analogue). Duality made precise: ζ\zeta^* and the theorem’s q\qup solve the two members of one Mellin family — the dual equation (45) governs the tail of the wealth distribution; the primal equation (26) governs the approach of the consumption function. §3.2’s obstruction is the case where the dual root fails to exist while the primal always does.

(The caveat, plainly: the clean constant-times-power tail needs the random log-steps of the descent NOT to be all multiples of one common unit. A discretized permanent shock violates that — the ‘lattice’ case — and then the constant is replaced by a bounded periodic wobble in log wealth: the same log-periodic phenomenon as the resonance circle, and equally harmless for the exponent.)

The reachability taxonomy (answering: with random shocks, what can “unreachable” even mean?):

  1. ψ1\psi \equiv 1, bounded θ\theta — literally unreachable. AÞΓ<1A \equiv \ThornG < 1: no expansion, ever; wˉt+1ÞΓwˉ+C0\wbarNxt \leq \ThornG \wbar + C_0 gives an absorbing interval (the return region of §8), so the ergodic support is compact — zero probability beyond a finite point, not merely small. This is the strongest form of §3.2’s contrast: no wealth tail of any kind, yet the gap is exactly a power law.

  2. Permanent shocks, even bounded (the HAFiscal case) — reachable but power-law rare. The 7-atom estimated ψ\psi has ψmin=0.9217<ÞΓ\psi_{\min} = 0.9217 < \ThornG at all three calibrations, so Pr(A>1)>0\Pr(A > 1) > 0: runs of small permanent-income realizations (ψ<ÞΓ\psi < \ThornG, so A>1A > 1) shrink the denominator of the normalized ratio and so inflate normalized perfect-foresight wealth without bound — excursions to high normalized PF wealth are produced by bad income luck, not good (the household is income-poor and ratio-rich) — the support is unbounded even though every shock is bounded. “Behaviorally reachable” is then a measure-and-time statement: tail mass tζ\sim t^{-\zeta^*}, excursions mean-reverting at rate E[lnA]|\E[\ln A]| per period.

    calibrationE[lnA]\E[\ln A]dual root ζ\zeta^*Pr(wˉ>10×typical)\Pr(\wbar > 10\times\text{typical}) (Chernoff scale)
    HS-mean-0.01139.1961010\sim 6\cdot10^{-10}
    COL-TOP-0.00393.177104\sim 7\cdot10^{-4}
    GIC-CAP (the patience ceiling)+0.0009\mathbf{+0.0009}noneno stationary wˉ\wbar from impatience alone

    Patience fattens the tail (9.23.29.2 \to 3.2 \to none) — the familiar patience-to-Pareto mechanism of the wealth-distribution literature, here seen from inside a single agent’s problem.

  3. Unbounded θ\theta — any level is possible. One large transitory draw passes any point, so there is no finite maximal wealth at all; the ergodic tail then inherits the income tail (Stachurski–Toda 2019). The theorem is indifferent: the L4′ extension covers unbounded θ\theta.

The patience-ceiling exception. At the GIC-CAP ceiling the mean log-drift is positive: under the reduced i.i.d. dynamics that calibration’s normalized perfect-foresight wealth drifts up forever — impatience alone never pulls it back, and no dual root exists. (As in item 2, the upward drift is a drift of the NORMALIZED ratio, fed by small permanent-income realizations that shrink its denominator — bad income luck, not good.) The cap’s construction pins raw growth impatience to hold — barely (ÞΓ=0.99964<1\ThornG = 0.99964 < 1) — while the ψ\psi-modified conditions fail (E[ln(ÞΓ/ψ)]=+0.0009>0\E[\ln(\ThornG/\psi)] = +0.0009 > 0 and E[ÞΓ/ψ]1.002>1\E[\ThornG/\psi] \approx 1.002 > 1): the raw, log-drift, and expectation-of-ratio conditions are separated by Jensen gaps, and the ceiling calibration sits between raw and modified — an estimated-parameters instance of BST’s condition-lattice distinctions, the case in which GIC-Mod fails while GIC holds, so only the pseudo-target survives. In (2)'s race language the cap’s side effect on the exponent is stark: pinning ÞΓ\ThornG a hair under 1 kills the descent rate — þg=3.6×104-\GPRte = 3.6\times10^{-4} per quarter, a wealth-ratio half-life of roughly 479 years against a discounting half-life of 34 (§0’s half-life table) — and since the exponent is the ratio discounting/descent, a dead descent rate is precisely what hands the race to discounting: the deterministic ratio is 14\approx 14, and only the ψ\psi-tilt of (26) drags the realized root down to 1.47 — supercritical still. The cap and the exponent condition thus pull on the same dial in opposite directions: the cap keeps the solution method safely inside raw growth impatience (§3.1), and does that job; but proximity to the raw-GIC knife-edge is exactly what pushes the exponent past its own boundary — which is why the simplified capstone maintains RÞΓ<1\Rcal\ThornG < 1 (its standing assumption, GIC-policy) instead of proximity to the GIC as the discipline on patience. What truncates that calibration’s ergodic normalized PF wealth in the full model (the distribution whose extreme upper quantile fixes the top of the estimation grid) is mortality with replacement — the reset mechanism of the wealth-tail literature (Toda 2014) — and within an expected working life its upward drift moves normalized PF wealth only 15%\approx 15\%.

Practical conclusion, sharpened. The asymptotic region — the range wˉwˉ0\wbar \geq \wbar_0 on which the gap law’s leading-order form governs and the remainders ηL\etaL, ηR\etaR are O(1/wˉ)O(1/\wbar)-small — is visited with stationary probability ranging from t9t^{-9} (HS-mean) down to “only via the mortality-reset tail” (GIC-CAP) — and in the ψ1\psi \equiv 1 benchmark, never. The theorem’s practical role is therefore exactly as §9 frames it: a boundary condition and extrapolation form for the function, not a description of states to simulate.

9. The computational payoff

An in-grid representation plus the theorem’s tail is the whole solution; no ad-hoc external criterion (a grid top, an extrapolation form, a tolerance) does any work.

  1. The boundary value is a boundary condition. Solving on a compactified grid z[0,zˉ]z \in [0, \bar{z}] (equivalently, a top-of-grid asymptotic condition) is standard numerical practice; Theorem I upgrades the practice to principle: the exact boundary data at z=0z = 0 is W(0)=BW(0) = B (q>1\qup > 1), i.e. c(m)κ(b+h)B/m\cFunc(m) \approx \kap(\bNrm+h) - B/m at the top of the grid — with BB known in closed form from the calibration.

  2. The extrapolation form of record. The gap extrapolant xCwˉq=C(b+h)qx \approx C\,\wbar^{-q} = C\,(\bNrm+h)^{-q}, q=min(1,q)q = \min(1, \qup), is the asymptotically correct form — the theorem behind the decay_extrap='powerlaw' flip in the HAFiscal/HARK code (RECONCILED-002; HARK PR fix-aggshock-pf-decay-extrap), replacing the historical exponential heuristic, which decays faster than the forcing floor (stage-A Lemma 6.1 / stage-B Lemma B-6.1) permits as an asymptotic form.

  3. The knife-edge window, quantified on HAFiscal’s own numbers. At the estimated calibrations the growth patience factor is close to 1 (þg0-\GPRte \to 0), the rescaling is near-additive, and the asymptotic region begins only at astronomically large PF wealth. GIC-CAP makes this concrete: þg=3.6104-\GPRte = 3.6\cdot10^{-4}, q=1.474\qup = 1.474, Bψ=357B_\psi = 357 — and the compensated gap has covered only 42% of the way to its boundary value at the top of the estimation grid, 91% even at 130×130\times the grid top (Figure 3a). That is why the exponential heuristic was locally harmless on historical grids (0.22%\leq 0.22\% policy effect at the audited window) even though it is asymptotically wrong. The plateau-onset scale lnwˉlnwˉc+O(1)/(q1)\ln \wbar \gtrsim \ln \wbar_c + O(1)/(\qup-1) quantifies this.

  4. A built-in diagnostic. Theorem γ-T’s detection property is a practical grid test: compensate the measured gap by trial exponents; only min(1,q)\min(1, \qup) stabilizes. The grid-depth migration of fitted exponents toward min(1,q)\min(1, \qup) is the theorem-backed convergence signature for solver validation. The migration’s pace is the purification rate of Figure 1: the fitted exponent’s residual dies like e1qe^{-|1-\qup|} per ee-fold of grid depth — slowest near resonance.

  5. The computer-science face. In algorithmic terms the sandwich (31)(32) is a perturbed divide-and-conquer recurrence, and (26) is its Akra–Bazzi critical-exponent equation iaibip=1\sum_i a_i b_i^{p} = 1 at p=qp = -\qup; the trichotomy including the resonance log is the master theorem’s three cases (Akra–Bazzi 1998; Leighton 1996, whose perturbation tolerance — shifts of the recurrence argument up to size t/log1+εtt/\log^{1+\varepsilon} t at argument tt — our bounded shifts satisfy with room). That proof is carried out in full in alt_proof_econcs.md, with hypothesis-match tables against the machine-formalized statements (Isabelle AFP Akra_Bazzi, Eberl 2015/2017; Lean 4 mathlib4 Computability/AkraBazzi, with its honest mismatches recorded) and executed interval-arithmetic certificates for q\qup, the denominator signs, and BB/BψB_\psi — the certification practice a formalization effort would inherit.

  6. The finite-horizon slope, for free. §2.4’s construction makes the truncated sum (21) finite at any finite distance from TT, so the gap is exactly-order wˉ1\wbar^{-1} at ANY finite horizon regardless of q\qup — the theorem behind the Method of Moderation’s finite-horizon practice: the correct top-of-grid slope is exactly 1 for a finite-horizon solver, and min(1,q)\min(1, \qup) is an infinite-horizon limit statement, not a finite-TT one.

(The bridge, for economists who have not met divide-and-conquer recurrences: an algorithm that splits a problem into a pieces of a fraction b the size obeys exactly our kind of one-step recursion, level by level; its ‘critical exponent’ equation asks at what power the per-level contributions neither die out nor compound — the same question our root equation asks about backward steps of the journey.)

10. One theorem, four engines — and one telling

enginedocumentfacedistinctive reachreview
backward induction from TT (the spine, §2)this document, from eq1_backward_derivation.mdeconomics / constructionthe organizing telling: derives (EQ-1) and unrolls it to every published constantbattery verify_eq1_backward_checks.py: ALL PASS
renewal / implicit renewal + tiltstage_A_proof.md, stage_B_proof.mdclassical probabilityeverything, incl. sharp PP at q<1\qup<1 (A2/B2/B2-arith), general-ψ\psi resonance (B-res), unbounded θ\theta (L4′)R1–4, RB1–4, RC1–4: zero broken
journey path-sum (PV of the adjustment stream)alt_proof_econlit.mdeconomicsthe literature-anchored narrative; general-ψ\psi BψB_\psi with no margin conditionRD1/RD2: zero broken
perturbed master theoremalt_proof_econcs.mdcomputer scienceorder + constants by induction; formalization anchors + certificates; Stage-B under ψmin>ÞΓ\psi_{\min} > \ThornGRE1/RE2: zero broken
compactified boundary fixed pointalt_proof_compactified.mdgeometry / computation (this core)shortest engine; the scalar boundary equation; exponent detection; Cesàro resonance; general-ψ\psi BψB_\psiRF1/RF2: zero broken

The intellectual division of labor: the spine constructs the theorem backward from the perfect-foresight resumption date and makes every constant visible in one geometric sum; the compactified core proves it the simplest way; the renewal chain proves the most; the path-sum proof tells the same story forward in the economics literature’s language; the master-theorem proof ports it to the community whose tools certify and formalize it.

11. Figures and reproduction

All figures are produced by make_final_proof_figures.py (single deterministic run, ~140 s; 12/12 pre-registered checks; log committed as make_final_proof_figures_out.txt); per-figure provenance and headline numbers are in figures/FIGURES.md, including the calibration table with source-line citations into the HAFiscal estimation code (HS-mean q=0.3813\qup = 0.3813; COL-TOP q=0.6942\qup = 0.6942; GIC-CAP q=1.4735\qup = 1.4735, Bψ=356.6B_\psi = 356.6; resonance auxiliary q=1.000000\qup = 1.000000) and the honest i.i.d.-reduction caveat (HAFiscal’s persistent unemployment spells are mapped to an i.i.d. atom at income 0.7 with the education-specific ergodic unemployment rate). The underlying solvers are the review-pack EGM libraries (review/R4_egm_lib.py, review/RB4_egm_lib_B.py), run under the program’s guards (longdouble; tail points quoted only where x/c > 1e-10; large-(þg)(-\GPRte) calibrations for amplitude exhibits). The theorem-level verification of record is the pre-registered falsifier suites of the four source documents (all committed with outputs; 31/31, 22/22, 13/15+2-honest-FAILs, and the six refuter packs’ independent batteries).

The spine’s named instrument is verify_eq1_backward_checks.py (pre-registered gates E1–E6 frozen in its header; committed verify_eq1_backward_checks_out.txt), which implements §2’s backward iteration exactly — pointwise Euler solves against the PF continuation, assembled in difference form — and checks the one-step law (19), (EQ-1) as an identity along the iteration, the unrolled amplitude against (27), the q<1\qup < 1 exponent, the resonance slope, and determinism; its ψ\psi-general companion is verify_eq1_psi_spike_checks.py (§2.6, committed verify_eq1_psi_spike_checks_out.txt).

References

Buffer-stock / income-fluctuation foundations. Carroll (1997), QJE 112(1):1–55. Carroll, Theoretical Foundations of Buffer Stock Saving — cite as NBER WP 10867 (2004) and forthcoming, Quantitative Economics (no volume/pages exist yet; maintained draft at llorracc.github.io/BufferStockTheory). Carroll–Kimball (1996), Econometrica 64(4):981–992. Carroll–Holm–Kimball (2021), JET 195:105276. Schechtman–Escudero (1977), JET 16(2):151–166. Deaton (1991), Econometrica 59(5):1221–1248. Chamberlain–Wilson (2000), RED 3(3):365–395. Li–Stachurski (2014), JEDC 45:353–365. Ma–Stachurski–Toda (2020), JET 187:105003. Ma–Toda (2021), JET 192:105193; Ma–Toda (2022), J. Math. Econ. 98:102562. Szeidl (unpublished ms., 2008/2013) — flagged, not load-bearing. Harmenberg (2021), JEDC 129:104185. Zeldes-style zero-income event as adopted in BST’s income process.

Precautionary-saving canon. Kimball (1990), Econometrica 58(1):53–73; Kimball (1993), Econometrica 61(3):589–611 (decorative). Pratt (1964), Econometrica 32(1/2):122–136. Arrow (1965/1971), Aspects/Essays in the Theory of Risk-Bearing (lecture volumes). Leland (1968), QJE 82(3):465–473. Sandmo (1970), REStud 37(3):353–360. Drèze–Modigliani (1972), JET 5(3):308–335.

Heavy tails in economics (kinship — the wealth-LEVEL literature; §3.2). Benhabib–Bisin–Zhu (2011), Econometrica 79(1):123–157; (2015), JET 159:489–515. Toda (2014), JET 154:310–348; (2019), JME 104:101–113. Stachurski–Toda (2019), JET 182:1–24 (+ 2020 corrigendum) — the econ-packaged compact-support obstruction. Beare–Toda (2022), Econometrica 90(4):1811–1833. Beare–Seo–Toda (2022), Econometric Theory 38(5):986–1013 (+ corrigendum ET 41(5), 2025). Gouin-Bonenfant–Toda (2023), Quantitative Economics 14(1):201–233 (+ the Pareto-extrapolation toolkit).

Power laws in the random-growth canon (§0 and 8.1 How rarely is the tail visited? The dual (Kesten) root). Gabaix (1999), “Zipf’s Law for Cities: An Explanation,” QJE 114(3):739–767 — shows that proportional random growth (Gibrat’s law) with a small friction preventing collapse to zero generates a power-law city-size distribution, and that the Zipf exponent of 1 emerges in the small-friction limit; the template for “random growth ⟹ Pareto tail.” Gabaix (2009), “Power Laws in Economics and Finance,” Annu. Rev. Econ. 1:255–294 — the survey of the empirical power laws (income and wealth, cities, firms, stock returns, trading volume, executive pay) and of the random-growth/Kesten theory that generates them, including the exponent equation E[Aζ]=1\E[A^\zeta]=1 used in 8.1 How rarely is the tail visited? The dual (Kesten) root. Gabaix (2016), “Power Laws in Economics: An Introduction,” J. Econ. Perspect. 30(1):185–206 — the economist-facing introduction arguing that these are true quantitative laws, both empirically established and theoretically understood.

Probability engines (used by the original chain). Kesten (1973), Acta Math. 131:207–248. Goldie (1991), Ann. Appl. Probab. 1(1):126–166. Athreya–McDonald–Ney (1978), Ann. Probab. 6(5):788–797. Grincevičius (1975); Grey (1994), AAP 4(1):169–183; Vervaat (1979); Kevei (2017), J. Appl. Probab. 54(3):732–749 (+ 2016 ECP note). Buraczewski–Damek–Mikosch (2016), Springer (Thms 2.4.4/2.4.7). Jelenković–Olvera-Cravioto (2012), AAP 44(2):528–561.

Divide-and-conquer / formalization (the CS face; §9.5). Akra–Bazzi (1998), Comput. Optim. Appl. 10(2):195–210. Leighton (1996), MIT ms. Eberl (2017), J. Automated Reasoning 58(4):483–508

Compactification ancestry. Dumortier–Llibre–Artés (2006), Qualitative Theory of Planar Differential Systems, Springer, ch. 5.