pith. sign in

arxiv: 2604.17516 · v4 · pith:HVGAHAPXnew · submitted 2026-04-19 · 🌌 astro-ph.IM · astro-ph.GA

Kardashev's Conundrum: Statistical Falsification of the Standard Kardashev Model and the Kardashev--Sagan--Nakamoto Resolution

Sebastian Gurovich This is my paper

Pith reviewed 2026-05-21 08:49 UTC · model grok-4.3

classification 🌌 astro-ph.IM astro-ph.GA
keywords Kardashev scaleenergy productionstatistical modelingBitcoin hashratecivilization typesLandauer limitexponential growthlinear regression
0
0 comments X

The pith

Fitting global energy production to the Kardashev model fails both statistically and physically, requiring normalization by computational hashrate.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper tests whether global primary energy production has grown at the one percent exponential rate conjectured in the standard Kardashev scale. Markov Chain Monte Carlo and linear regression applied to sixty years of data show that a linear model fits better than the exponential and that increments are non-Gaussian. Extrapolating the linear trend to solar luminosity produces a timescale of roughly 1.6 times 10 to the 15 years, far longer than the age of the universe. This leads to the conclusion that energy production alone is dimensionally incomplete for the Kardashev variable. The authors introduce a renormalization dividing energy by Bitcoin hashrate to incorporate information processing capacity.

Core claim

No functional form fitted to P(t) alone can simultaneously satisfy statistical adequacy and physical coherence: the Kardashev variable is dimensionally incomplete. We propose the Kardashev-Sagan-Nakamoto renormalisation B(t) = P(t)/H(t) [J/Hash, the KarNak unit], where H(t) is the annual Bitcoin hashrate. The renormalisation adds no free parameters, is motivated by the Landauer limit, and fulfils Sagan's information-richness requirement.

What carries the argument

The Kardashev-Sagan-Nakamoto renormalisation B(t) = P(t)/H(t), which normalizes energy production by annual Bitcoin hashrate to add the missing information dimension.

If this is right

  • The posterior growth rate from MCMC is 2.01 percent per year with credible interval excluding 1 percent.
  • A linear OLS model is preferred over exponential by WAIC with R squared of 0.987.
  • Year-over-year increments are non-Gaussian and show crisis outliers.
  • Extrapolation to solar luminosity gives a Type II timescale of approximately 1.6E15 years.
  • The renormalized B(t) spans 14 orders of magnitude over 2009-2024.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • This framework may be extended by testing alternative computational proxies such as total transistor production or global data center energy use.
  • Applying the same renormalization to historical data for other energy metrics could reveal whether the 14-order span is robust.
  • The dimensional incompleteness suggests that future models of civilizational progress should include both energy and information processing from the outset.

Load-bearing premise

That the annual Bitcoin hashrate provides a physically motivated proxy for the information-processing capacity needed to normalize energy production, based on the Landauer limit.

What would settle it

Future observations showing consistent energy production growth rates near or above one percent per year, or data demonstrating that energy per hash does not provide a coherent measure across different technological indicators.

Figures

Figures reproduced from arXiv: 2604.17516 by Sebastian Gurovich.

Figure 1
Figure 1. Figure 1: Global energy production 1965–2024 (OWID dataset) in watts, with the standard Kardashev one-percent exponential model [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Zoomed view of the OWID global primary-energy production data (1965–2024) with the linear and exponential best-fit models [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Year-over-year differences ∆Pi = Pi − Pi−1 in global primary-energy production (W) as a function of calendar year 1965– 2024, derived from the OWID direct-method energy dataset. Two stray negative outliers — clearly separated from the general trend of the data — occur at 2008 and 2020, corresponding respectively to the subprime great financial crisis and the COVID-19 pandemic. These are identifiable histor… view at source ↗
read the original abstract

We test the standard Kardashev one-percent exponential conjecture against six decades of global primary-energy production data (1965-2024; Our World in Data). Markov Chain Monte Carlo inference yields a posterior growth rate of r = 2.01 +/- 0.03% per year (95% credible interval [1.94%, 2.08%]), placing the Kardashev 1% value well outside the credible interval. A linear OLS model fits the data with remarkably low dispersion (R^2 = 0.987) and is preferred over the free-rate exponential by the Widely Applicable Information Criterion ({\Delta}WAIC = 5.5). Year-over-year increments are non-Gaussian (Shapiro-Wilk W = 0.925, p = 0.0014; skewness = -0.664) with identifiable crisis outliers (2008, 2020), rejecting the independent-increment multiplicative structure with positive drift required by Kardashev's (1+x)^t geometric series. Extrapolation of the linear model to the solar luminosity yields a Type II civilisational timescale of approximately 1.6E15 years -- approximately 1E5 times both the age of the Universe and the main-sequence lifetime of the Sun -- a physical reductio we term Kardashev's Conundrum. No functional form fitted to P(t) alone can simultaneously satisfy statistical adequacy and physical coherence: the Kardashev variable is dimensionally incomplete. We propose the Kardashev-Sagan-Nakamoto (KSN) renormalisation B(t) = P(t)/H(t) [J/Hash, the KarNak unit], where H(t) is the annual Bitcoin hashrate. The renormalisation adds no free parameters, is motivated by the Landauer limit, and fulfils Sagan's information-richness requirement. Over 2009-2024, B(t) spans 14 orders of magnitude.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The manuscript statistically tests the standard Kardashev 1% exponential growth conjecture against global primary-energy production data (1965-2024). MCMC yields a posterior growth rate of 2.01 ± 0.03 % yr⁻¹, outside the 1% value; an OLS linear model achieves R² = 0.987 and is preferred by WAIC (ΔWAIC = 5.5). Year-over-year increments are non-Gaussian (Shapiro-Wilk p = 0.0014), rejecting the required independent multiplicative structure. Linear extrapolation to solar luminosity produces a Type-II timescale of ~1.6 × 10¹⁵ yr, termed Kardashev’s Conundrum. The authors conclude that no P(t)-only functional form satisfies both statistical adequacy and physical coherence, rendering the Kardashev variable dimensionally incomplete, and propose the Kardashev–Sagan–Nakamoto renormalization B(t) = P(t)/H(t) (J/Hash) using annual Bitcoin hashrate, which spans 14 orders of magnitude over 2009–2024.

Significance. If the statistical preference for linear growth and the physical motivation for the renormalization both hold, the work would supply a data-driven critique of the canonical Kardashev scale and an information-theoretic alternative grounded in the Landauer limit. The explicit use of MCMC posteriors, WAIC model comparison, and reproducible energy data strengthens the falsification claim; however, the resolution’s reliance on a single proxy dataset limits its immediate generality within astrobiology and SETI.

major comments (3)
  1. [extrapolation to solar luminosity] The central claim that the Kardashev variable is dimensionally incomplete rests on the assertion that no functional form fitted to P(t) alone can be both statistically adequate and physically coherent. This conclusion is reached via linear extrapolation of the OLS fit to solar luminosity, but the manuscript does not demonstrate that the linear regime persists over 10¹⁵ yr or rule out saturation, technological transitions, or other physical limits that would invalidate the reductio (see the extrapolation paragraph following the WAIC comparison).
  2. [KSN renormalisation paragraph] The KSN resolution introduces B(t) = P(t)/H(t) motivated by the Landauer limit and Sagan’s information-richness criterion. However, the choice of annual Bitcoin hashrate H(t) as the normalizing proxy is not calibrated against independent computational metrics (e.g., aggregate FLOPS, data-center energy, or global communication volume). Without such cross-validation, the reported 14-order span in B(t) may reflect Bitcoin-specific hardware and adoption dynamics rather than a general information-processing measure, weakening the claim that the renormalization restores physical coherence.
  3. [statistical tests of increments] The rejection of the Kardashev model cites non-Gaussian year-over-year increments (Shapiro-Wilk W = 0.925, p = 0.0014; skewness = −0.664) and identifiable crisis outliers. While this challenges the independent-increment multiplicative structure, the manuscript should explicitly map the statistical test onto the precise assumption in Kardashev’s original geometric series (i.e., that increments are i.i.d. log-normal with positive drift) rather than treating non-Gaussianity as sufficient by itself.
minor comments (2)
  1. [KSN renormalisation paragraph] Define the ‘KarNak unit’ explicitly at first use and clarify whether it is intended as a dimensional replacement or merely a convenient ratio.
  2. [abstract] The abstract states ‘the renormalisation adds no free parameters’; confirm that the Bitcoin hashrate series is treated as an external, fixed dataset with no additional fitting.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive comments that help sharpen the manuscript's claims. We address each major point below and indicate planned revisions.

read point-by-point responses

  1. Referee: [extrapolation to solar luminosity] The central claim that the Kardashev variable is dimensionally incomplete rests on the assertion that no functional form fitted to P(t) alone can be both statistically adequate and physically coherent. This conclusion is reached via linear extrapolation of the OLS fit to solar luminosity, but the manuscript does not demonstrate that the linear regime persists over 10¹⁵ yr or rule out saturation, technological transitions, or other physical limits that would invalidate the reductio (see the extrapolation paragraph following the WAIC comparison).

    Authors: The linear extrapolation serves as a reductio ad absurdum to illustrate that the statistically preferred model for P(t) yields an unphysical timescale, thereby showing that energy production alone is dimensionally insufficient to reach solar luminosity without additional structure. We do not assert that linearity continues indefinitely; rather, the 60-year record shows no saturation, and any future deviation (saturation or transition) would only lengthen the timescale further. We will add a sentence clarifying that the argument is illustrative of the current data's implications rather than a literal long-term forecast. revision: partial

  2. Referee: [KSN renormalisation paragraph] The KSN resolution introduces B(t) = P(t)/H(t) motivated by the Landauer limit and Sagan’s information-richness criterion. However, the choice of annual Bitcoin hashrate H(t) as the normalizing proxy is not calibrated against independent computational metrics (e.g., aggregate FLOPS, data-center energy, or global communication volume). Without such cross-validation, the reported 14-order span in B(t) may reflect Bitcoin-specific hardware and adoption dynamics rather than a general information-processing measure, weakening the claim that the renormalization restores physical coherence.

    Authors: Bitcoin hashrate is adopted as the most transparent, continuously recorded proxy for large-scale computational work, directly tied to energy expenditure and verifiable from public blockchain data since 2009. While independent calibration against aggregate FLOPS or data-center metrics would strengthen generality, such harmonized global datasets do not exist at comparable temporal resolution and precision. The observed 14-order span is consistent with the rapid growth in computational infrastructure and aligns with the Landauer and Sagan motivations. We will expand the text to acknowledge the proxy nature and note the need for future cross-validation as broader metrics become available. revision: partial

  3. Referee: [statistical tests of increments] The rejection of the Kardashev model cites non-Gaussian year-over-year increments (Shapiro-Wilk W = 0.925, p = 0.0014; skewness = −0.664) and identifiable crisis outliers. While this challenges the independent-increment multiplicative structure, the manuscript should explicitly map the statistical test onto the precise assumption in Kardashev’s original geometric series (i.e., that increments are i.i.d. log-normal with positive drift) rather than treating non-Gaussianity as sufficient by itself.

    Authors: We agree that the connection should be stated more explicitly. The Shapiro-Wilk test on year-over-year increments directly assesses the normality of log-increments required by Kardashev's (1 + x)^t geometric series under the assumption of i.i.d. log-normal growth with positive drift. The observed non-Gaussianity and crisis outliers violate this structure. We will revise the relevant paragraph to make this mapping explicit. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation uses external data and independent renormalization

full rationale

The paper fits linear OLS and MCMC exponential models directly to external global primary-energy data (1965-2024, Our World in Data), reports posterior r = 2.01 +/- 0.03% and prefers linear via WAIC and non-Gaussian increments. The 1.6E15-year timescale is explicitly the result of extrapolating the fitted linear model to solar luminosity as a reductio (Kardashev's Conundrum), not an independent first-principles prediction. The KSN resolution defines B(t) = P(t)/H(t) using a separate external Bitcoin hashrate series (2009-2024), adds no free parameters, and invokes the Landauer limit for motivation. No quoted step reduces by construction to a self-definition, fitted parameter renamed as prediction, or load-bearing self-citation; the central claims remain self-contained against the cited datasets and statistical tests.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 1 invented entities

The central claim rests on the assumption that linear extrapolation of fitted energy data remains valid over cosmological timescales and that Bitcoin hashrate is an appropriate stand-in for information-processing capacity; no new free parameters are introduced beyond the standard statistical fits, but the renormalization itself is an invented scaling that lacks independent empirical calibration.

free parameters (2)
  • linear slope of energy production
    Fitted via OLS to 1965-2024 data; directly determines the 1.6E15-year extrapolation.
  • MCMC growth rate posterior
    Posterior mean 2.01% used to contrast with the Kardashev 1% conjecture.
axioms (2)
  • domain assumption Year-over-year energy increments are independent and multiplicative under the Kardashev model
    Invoked when rejecting the geometric series structure on the basis of non-Gaussian increments and crisis outliers.
  • ad hoc to paper Bitcoin hashrate is a valid proxy for computational information processing
    Used to define the KSN renormalization without additional calibration data.
invented entities (1)
  • KSN unit (J/Hash) no independent evidence
    purpose: Renormalized measure of civilizational progress that incorporates information processing
    New scaling introduced to resolve the dimensional incompleteness of raw energy P(t).

pith-pipeline@v0.9.0 · 5905 in / 1891 out tokens · 77694 ms · 2026-05-21T08:49:08.828015+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

  • IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel echoes
    ?
    echoes

    ECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.

    We define the KSN state variable as B(t)=P(t)/H(t) [J Hash⁻¹ ≡ KN], ... motivated by the Landauer limit connecting energy to irreversible computation, and fulfils the information-richness requirement ... The renormalization adds no free parameters

  • IndisputableMonolith/Foundation/AlphaCoordinateFixation.lean J_uniquely_calibrated_via_higher_derivative echoes
    ?
    echoes

    ECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.

    No functional form fitted to P(t) alone can simultaneously satisfy statistical adequacy and physical coherence: the Kardashev variable is dimensionally incomplete.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

13 extracted references · 13 canonical work pages · 1 internal anchor

  1. [1]

    International Journal of Theo- retical Physics 21, 905–940

    The thermodynamics of computation—a review. International Journal of Theo- retical Physics 21, 905–940. doi:10.1007/BF02084158. Blumer, A., Ehrenfeucht, A., Haussler, D., Warmuth, M.K.,

    work page doi:10.1007/bf02084158

  2. [2]

    Information Processing Let- ters 24, 377–380

    Occam’s razor. Information Processing Let- ters 24, 377–380. doi:10.1016/0020-0190(87)90114-1. Box, G.E.P., Cox, D.R.,

    work page doi:10.1016/0020-0190(87)90114-1

  3. [3]

    Journal of the Royal Statistical Society: Se- ries B 26, 211–252

    An analysis of transfor- mations. Journal of the Royal Statistical Society: Se- ries B 26, 211–252. doi:10.1111/j.2517-6161.1964. tb00553.x. Dallaire-Demers, P.L., BTQ Technologies Team,

    work page doi:10.1111/j.2517-6161.1964 1964

  4. [4]

    arXiv preprintarXiv:2603.25519

    Kardashev scale quantum computing for bitcoin min- ing. arXiv preprintarXiv:2603.25519. Freeth, T., Higgon, D., Dacanalis, A., MacDonald, L., Georgakopoulou, M., Wojcik, A.,

    work page arXiv

  5. [5]

    Gray, R.H.,

    URL: https://doi.org/10.1038/s41598-021-84310-w, doi:10.1038/s41598-021-84310-w. Gray, R.H.,

    work page doi:10.1038/s41598-021-84310-w

  6. [6]

    Hastings, W.K.,

    doi:10.3847/1538-3881/ ab792b. Hastings, W.K.,

    work page doi:10.3847/1538-3881/

  7. [7]

    Biometrika , year = 1970, month = apr, volume =

    Monte carlo sampling methods us- ing markov chains and their applications. Biometrika 57, 97–109. doi:10.1093/biomet/57.1.97. Kardashev, N.S.,

    work page doi:10.1093/biomet/57.1.97

  8. [8]

    IBM Journal of Research and Development 5, 183–191

    Irreversibility and heat generation in the computing process. IBM Journal of Research and Development 5, 183–191. doi:10.1147/rd.53.0183. Lares, M., Funes, J.G., Gramajo, L.,

    work page doi:10.1147/rd.53.0183

  9. [9]

    IAU 2015 Resolution B3 on Recommended Nominal Conversion Constants for Selected Solar and Planetary Properties

    IAU 2015 Resolution B3 on Recom- mended Nominal Conversion Constants for Selected Solar and Planetary Properties. arXiv e-prints , arXiv:1510.07674arXiv:1510.07674. Metcalfe, B.,

    work page internal anchor Pith review Pith/arXiv arXiv 2015

  10. [10]

    Computer 46, 26–31

    Metcalfe’s law after 40 years of ether- net. Computer 46, 26–31. doi:10.1109/MC.2013.374. 10 Metropolis, N., Rosenbluth, A.W., Rosenbluth, M.N., Teller, A.H., Teller, E.,

    work page doi:10.1109/mc.2013.374 2013

  11. [11]

    Metropolis, A

    Equation of state calcula- tions by fast computing machines. Journal of Chemical Physics 21, 1087–1092. doi:10.1063/1.1699114. Nakamoto, S.,

    work page doi:10.1063/1.1699114

  12. [12]

    Available: https://doi.org/10.1093/biomet/52.3-4.591

    An analysis of variance test for normality (complete samples). Biometrika 52, 591–611. doi:10.1093/biomet/52.3-4.591. Shklovskii, I.S., Sagan, C.,

    work page doi:10.1093/biomet/52.3-4.591

  13. [13]

    2016.Managing Knowledge Integration Across Boundaries

    Turing, Lovelace, and Babbage, in: The Turing Guide. Oxford University Press. URL:https: //doi.org/10.1093/oso/9780198747826.003.0033, doi:10.1093/oso/9780198747826.003.0033. Watanabe, S.,

    work page doi:10.1093/oso/9780198747826.003.0033