REVIEW 3 major objections 5 minor 2 cited by
Thermodynamic description of worldwide distribution of energy and carbon emission
T0 review · 3 major / 5 minor · reviewed 2026-07-09 · glm-5.2
Pith's one-line read Rayleigh-Jeans condensation fits 40 years of energy and CO2 data
desk verdict RJ thermalization applied to country-level energy/CO2 distributions: visually close fits but no quantitative benchmarks against alternatives read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The Rayleigh-Jeans distribution rho_m = T / (E_m - mu) with two implicit constraint equations (sum of rho_m = 1 and sum of E_m * rho_m = E). The RJE model uses energy levels E_k = (e^{ak/N} - 1) / (e^a - 1) with parameter a controlling the density of states. The rescaled energy parameter epsilon = E/B (ratio of average to maximum energy) is matched to the Gini coefficient, and a is fit to minimize the geometric distance between theoretical and empirical Lorenz curves. Lorenz curves are constructed by computing cumulative probability h(m) and cumulative wealth fraction w(m), and Pareto curves are the complementary cumulative distribution C(w_m) = 1 - h(m).
What would settle it
If the rescaled Lorenz and Pareto curves for energy consumption or CO2 emission were to shift significantly over time (e.g., the Gini coefficient moving outside the 0.87-0.89 range for a sustained period), or if the RJE model could not reproduce the curves for a future decade with the same structural relationship between epsilon and G, the claim that a stable thermodynamic equilibrium governs the distribution would be undermined.
Extended reading notes
Core claim
The Rayleigh-Jeans distribution rho_m = T / (E_m - mu), derived from two conservation constraints (total energy and probability norm), reproduces the Lorenz and Pareto curves for country-level energy consumption, electricity generation, and CO2 emission over a 40-50 year period with Gini coefficients near 0.88. The inequality in these distributions is attributed to RJ condensation at low-energy modes, which naturally creates a large phase of energy-poor countries and a tiny phase of energy-rich countries. The rescaled curves remain stable over decades even as absolute values change, which the authors interpret as evidence of an underlying thermodynamic equilibrium.
Load-bearing premise
The paper assumes that countries behave like nonlinear coupled oscillators in a closed system where total energy and total probability are conserved, leading to dynamical thermalization. In reality, global energy is continuously injected from external sources (fossil fuels, solar, nuclear) and dissipated through consumption, so the conservation laws required for the Rayleigh-Jeans distribution to emerge from thermalization may not hold in the way the physical analogy requires
Editorial extensions
If this is right
- If the RJ condensation mechanism genuinely governs energy distribution, then policy interventions that redistribute energy between countries without changing the ratio epsilon (average to maximum energy) would leave the rescaled distribution unchanged, even if absolute emissions drop significantly.
- The stability of the rescaled curves over 40-50 years suggests that the parameter epsilon is a structural invariant of the international system, potentially more informative than absolute emission targets for understanding distributional dynamics.
- The same framework applied to wealth (in companion work) and now to energy and CO2 suggests a universal condensation mechanism across different economic quantities, raising the question of whether any exchange-based distribution with two conserved quantities inevitably produces RJ-type inequality.
- The RJE model's ability to recover effective per-country energy values from the fitted Lorenz curve (with most countries matching within small relative error, except top consumers like China and US at roughly 20% deviation) could be used to identify countries whose energy consumption deviates from the thermodynamic prediction, flagging structural anomalies in the global distribution.
Reading between the lines
- The framework treats the international system as effectively closed with respect to energy and probability norm, but real energy systems are open: energy is continuously injected from fossil fuels and solar sources and dissipated through use. If the system is not genuinely conservative, the RJ distribution may be fitting a steady-state pattern without the mechanism that generates it being thermody
- The two-parameter RJE model (epsilon and a) fitting Lorenz curves for ~200 countries may be a flexible enough functional form that good fits are not strongly discriminating. A testable prediction would be whether the fitted parameter a shows meaningful trends or correlations with independent structural features of the global economy, rather than being a free fit parameter each year.
- If the condensation mechanism is real, perturbing the system (e.g., a major country drastically changing its energy share) should produce relaxation back to the same rescaled distribution over some characteristic timescale, which could in principle be measured from historical data after large shocks.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript applies the Rayleigh-Jeans (RJ) thermalization framework—previously used by the authors for wealth distributions—to worldwide distributions of energy consumption, electricity generation, and CO2 emission across countries over a 40–50 year period. The central claim is that the RJ extended (RJE) model, with parameters ε (matched to the Gini coefficient) and a (fit by minimizing geometric distance to the Lorenz curve), provides an excellent description of the observed Lorenz and Pareto curves. The paper presents visual comparisons for several years and quantities, along with world maps of effective RJE values and their deviations from data.
Significance. The paper addresses a legitimate empirical question: whether the stable, high-inequality distributions of energy consumption and CO2 emission across countries can be captured by a thermodynamic functional form. The dataset spans a substantial time period and the rescaled stability of the Lorenz curves is a genuine empirical observation worth reporting. The RJE model produces visually close fits to the data. However, the significance of this finding is substantially weakened by the fitting protocol: one of two free parameters is set directly from the Gini coefficient of the target data, and no quantitative goodness-of-fit metrics or comparisons with alternative two-parameter distributions are provided. The authors acknowledge this concern in Section 4 but do not address it quantitatively.
major comments (3)
- Section 2, paragraph following Eq. (1): The fitting procedure uses two free parameters—ε is set to match the Gini coefficient G of the target data, and a is optimized to minimize geometric distance to the same Lorenz curve. Since G is a scalar summary of the Lorenz curve being fitted, one parameter reproduces one summary statistic and the other fits the remaining shape of a smooth monotonic curve from (0,0) to (1,1). With this structure, many two-parameter families (lognormal, gamma, beta, etc.) would likely produce visually similar fits. The claim that the RJE model 'provides an excellent description' (Section 4) is not established without (i) quantitative goodness-of-fit metrics (RMSE, R², KS statistic, etc.) and (ii) comparison against alternative two-parameter distributions fitted with the same protocol. This is load-bearing for the central claim and must be addressed.
- Section 3, Figs. 2 and 4: The 'nearly perfect agreement' and 'very good agreement' assessments are based solely on visual inspection of overlaid curves. No quantitative fit metrics are reported anywhere in the paper. Without such metrics, the strength of the agreement—and whether the RJE form specifically outperforms simpler alternatives—cannot be assessed. At minimum, the authors should report a quantitative distance metric for each fit and confirm that it is small relative to the curve-to-curve variability across years.
- Section 1 and Section 2 (paragraph containing Eq. (1)): The structural premise that countries can be modeled as nonlinear coupled oscillators undergoing dynamical thermalization with conserved total energy and probability norm is not demonstrated. Energy is continuously injected (fossil fuels, solar, etc.) and dissipated, not conserved in a closed system. The probability norm conservation (Σρ_m = 1) is imposed by construction. Without a microscopic dynamics argument or empirical evidence that international energy/CO2 dynamics satisfy these conservation laws in the relevant sense, the RJ distribution functions as a phenomenological ansatz rather than a derived consequence. The authors should either provide such an argument or reframe the claim accordingly. As stated, the physical mechanism claim is unsupported.
minor comments (5)
- Section 3.1: The Gini coefficient range is stated as '0.894 ≤ G ≤ 0.972' but the individual values listed (0.894, 0.876, 0.872, 0.875, 0.879) have a maximum of 0.894, not 0.972. This appears to be a typo.
- Section 2: The RJE spectrum E_k = (e^{ak/N} − 1)/(e^a − 1) is introduced without much motivation beyond 'the density of states is decreasing at high energies.' A brief physical or empirical justification for this specific functional form would help readers assess its plausibility.
- Fig. 3, right panel: The logarithmic difference Δ between RJE and data values is shown, but no summary statistic (e.g., median or RMS of |Δ|) is reported. A quantitative summary would strengthen the claim that 'for most countries these differences are rather small.'
- Section 4: The authors acknowledge the circularity concern ('a critical mind can argue...') but do not respond to it beyond asserting universality. A direct response with the quantitative metrics requested above would be more effective.
- References [10, 11] are central to the methodology but appear to be by the same authors and very recent (2025–2026). The dependence on these works for procedural details is heavy; ensuring their availability at production stage would aid reproducibility.
Circularity Check
Two-parameter curve fit with one parameter matched to Gini; not strictly circular but weak evidence
full rationale
The paper's RJE model has two free parameters: ε is set so that the model's Gini coefficient equals the data's Gini, and a is optimized to minimize the geometric distance between the model and data Lorenz curves (Section 2). The Gini coefficient is a single scalar summary of the Lorenz curve, so matching it does not by construction reproduce the full curve—many distinct Lorenz curves share the same Gini. The second parameter a then selects among the remaining family of curves. This is a standard two-parameter curve fit, not a case where the output equals the input by construction. The Pareto curve agreement is not independent evidence since it is a reparameterization of the same cumulative distribution, but this is a weakness of evidence rather than circularity. The methodology (RJ distribution, RJE spectrum, fitting protocol) is imported from self-citations [10, 11] by the same authors, which is load-bearing for the approach. However, the method is applied to genuinely new datasets (energy consumption, electricity, CO2 emission) not analyzed in those prior works. The paper does not invoke a uniqueness theorem, does not rename a known result, and does not claim first-principles derivation of the RJ distribution itself. The authors explicitly acknowledge the limitation: 'a critical mind can argue that the fact that the ENTH theory fits well the real Lorenz curves may be useful but not a decisive argument.' No specific step reduces to its inputs by construction. The circularity is at most minor: the self-citation chain provides the method, and one fit parameter is a scalar summary of the target curve, but the central descriptive claim retains independent content because the full Lorenz curve shape is not determined by the Gini alone. Score 2 reflects the load-bearing self-citation chain without strict circularity in the fitting procedure. No specific circular steps are identified because none meet the threshold of being reducible to inputs by construction. The absence of quantitative goodness-of-fit metrics and comparison to alternative two-parameter distributions is a correctness risk, not a circularity finding.
Assumptions & free parameters
free parameters (3)
- ε (rescaled energy = E/B = average energy / bandwidth) =
0.0105-0.0166 for energy; 0.00672-0.0128 for CO2
- a (RJE spectrum parameter controlling density of states) =
3.44-4.82 across datasets
- N (number of thermalized agents) =
≥10000
assumptions (4)
- ad hoc to paper Countries can be represented as nonlinear coupled oscillators undergoing dynamical thermalization
- domain assumption Total system energy and probability norm are conserved during interactions between countries
- ad hoc to paper The density of energy states follows the RJE form E_k = (e^{ak/N} - 1)/(e^a - 1)
- domain assumption Lorenz curves should be constructed for total energy consumption (not per capita) over countries
invented entities (1)
-
Thermalized agents representing countries (oscillator modes)
Cite this review
Pith. "Pith review of Thermodynamic description of worldwide distribution of energy and carbon emission." pith.science (2026). https://pith.science/paper/WSOINA3O
@misc{pith2026260707315,
author = {Pith},
title = {Pith review of: Thermodynamic description of worldwide distribution of energy and carbon emission},
year = {2026},
howpublished = {\url{https://pith.science/paper/WSOINA3O}},
note = {Machine review of arXiv:2607.07315}
}
read the original abstract
Based on public data, we analyze the distributions of energy and carbon emission over world countries on a scale of the last 40-50 years using their presentation via Lorenz and Pareto curves. These curves in rescaled format remain remarkably stable on this time period being characterized by high values of the Gini coefficient indicating a strong inequality of energy distribution. To explain these distributions, we introduce the ENergy Thermalization Hypothesis (ENTH) according to which these distributions result from the Rayleigh-Jeans (RJ) thermalization and condensation of agents representing different countries. We show that this hypothesis provides an excellent description of Lorenz and Pareto curves obtained from data on the above time period. It also gives natural grounds for inequality relating it to the RJ condensation at low energy states. We additionally trace parallels with the wealth inequality in the world.
Figures
Forward citations
Cited by 2 Pith papers
-
Thermodynamic statistics of given names in USA and France
Name popularity in the US and France is extremely unequal and has stayed that way for over a century; the authors fit this pattern with a two-parameter Rayleigh-Jeans condensation model borrowed from wealth physics.
-
Thermodynamic theory of voting and EU elections
Vote shares in EU and French elections are shown to follow a Rayleigh-Jeans thermal distribution, with the Gini coefficient set by a fitted energy parameter.
Reference graph
Works this paper leans on
-
[1]
International Energy Agency. Global Energy Review 2025.https://www.iea.org/reports/ global-energy-review-2025(Accessed June 25, 2026)
work page 2025
-
[2]
EMBER. Global Electricity Review 2025.https://ember-energy.org/countries-and-regions/ world/(Accessed June 25, 2026)
work page 2025
-
[3]
EDGAR - Emissions Database for Global Atmospheric Research. GHG emissions of all world countries 2025 report.https://edgar.jrc.ec.europa.eu/report_2025(Accessed June 25, 2026)
work page 2025
-
[4]
U.S. Energy Information Administration (2026); Energy Institute - Statistical Review of World Energy (2025)
work page 2026
-
[5]
Ember (2026), Energy Institute - Statistical Review of World Energy (2025) – with major processing by Our World In Data [6]https://ourworldindata.org/(Accessed June 25, 2026)
work page 2026
-
[6]
Towards future infrastructures for sustainable multi-energy systems: A review
Guelpa E., Bischi A., Verda V ., Chertkov M., Lund H. Towards future infrastructures for sustainable multi-energy systems: A review. Energy 2019; 184: 2. 10
work page 2019
-
[7]
Statistical physics, Elsvier, Amsterdam (1980)
Landau L.D., and Lifshitz E.M.. Statistical physics, Elsvier, Amsterdam (1980)
work page 1980
-
[8]
Kolmogorov spectra of turbulence I,Springer , Berlin, (1992)
Zakharov V .E., L’vov V .S., and Falkovich G. Kolmogorov spectra of turbulence I,Springer , Berlin, (1992)
work page 1992
Show all 34 references
-
[9]
M., and Shepelyansky D
Frahm K. M., and Shepelyansky D. L. Wealth thermalization hypothesis and social networks. arXiv:2506.17720 [cond-mat.stat-mech] (2025); J. Stat. Phys. 2026; 193: 64
2025 arXiv
-
[10]
M., Ermann L
Frahm K. M., Ermann L. and Shepelyansky D. L. Thermodynamic description of wealth inequality in the world. arXiv:2606.17965[cond-mat.stat-mech] (2026)
2026 arXiv
-
[11]
Capital in the Twenty-First Century
Piketty T. Capital in the Twenty-First Century. Belknap Press of Harvard University Press, Cambridge, MA (2014(
2014
-
[12]
World Inequality Report 2022
Chancel L., Piketty T., Saez E., and Zucman G.et al.. World Inequality Report 2022. World Inequality Lab https://wir2022.wid.world(Accessed on June 25, 2026)
2022
-
[13]
Social and Economic Inequality
Roach B., Rajkanikar P.J., Goodwin N., and Harris J. Social and Economic Inequality. An ECI Teaching Mod- ule on Social and Environmental Issues, Economics in Context Initiative, Global Development Policy Center, Boston University (2003),https://www.bu.edu/eci/files/2023/05/In...
2003
-
[14]
Observation of geometric parametric instability induced by the periodic spatial self-imaging of multimode waves
Krupa K., Tonello A., Barthelemy A., Shalaby B.M., Bendahmane A., Millot G., and Wabnitz S. Observation of geometric parametric instability induced by the periodic spatial self-imaging of multimode waves. Phys. Rev. Lett. 2016; 116: 183901
2016
-
[15]
Condensation of classical nonlinear waves
Connaughton C., Josserand C., Picozzi A., Pomeau Y ., and Rica S. Condensation of classical nonlinear waves. Phys. Rev. Lett. 2005; 95: 263901
2005
-
[16]
Optical wave turbulence: Towards a unified nonequilibrium thermodynamic formulation of statistical nonlinear optics
Picozzi A., Garnier J., Hansson T., Suret P., Randoux S., Millot G., and Christodoulides D.N. Optical wave turbulence: Towards a unified nonequilibrium thermodynamic formulation of statistical nonlinear optics. Phys. Rep. 2014; 542: 1
2014
-
[17]
Classical Rayleigh-Jeans conden- sation of light waves: observation and thermodynamic characterization
Baudin K., Fusaro A., Krupa K., Garnier J., Rica S., Millot G., and Picozzi A. Classical Rayleigh-Jeans conden- sation of light waves: observation and thermodynamic characterization. Phys. Rev. Lett. 2020; 125: 244101
2020
-
[18]
Thermalization of orbital angular momentum beams in multimode optical fibers
Podivilov E.V ., Mangini F., Sidelnikov O.S., Ferraro M., Gervaziev H., Kharenko D.S., Zitelli M., Fedoruk M.P., Babin S.A., and Wabnitz S. Thermalization of orbital angular momentum beams in multimode optical fibers. Phys. Rev. Lett. 2022; 128: 243901
2022
-
[19]
Dynamical thermalization, Rayleigh-Jeans condensate, vortexes and wave collapse in quantum chaos fibers and fluid of light
Ermann L., Chepelianskii A.D., and Shepelyansky D.L. Dynamical thermalization, Rayleigh-Jeans condensate, vortexes and wave collapse in quantum chaos fibers and fluid of light. J. Phys. A: Math. Theor. 2026; 59: 055702
2026
-
[20]
Condensation transition in polydisperse hard rods
Evans M.R., Majumdar S.N., Pagonabarraga I., and Trizac E. Condensation transition in polydisperse hard rods. J. Chem. Phys. 2010; 132: 014102
2010
-
[21]
Constraint driven condensation in large fluctuations of linear statistics
Szavits-Nossan J., Evans M.R., and Majumdar S.N. Constraint driven condensation in large fluctuations of linear statistics. Phys. Rev. Lett. 2014; 112: 020602
2014
-
[22]
On the concentration of large deviations for fat tailed distributions, with application to financial data
Filiasi M., Livan G., Marsili M., Peressi M., Vesselli E., and Zarinelli E. On the concentration of large deviations for fat tailed distributions, with application to financial data. J. Stat. Mech. 2014; P09030
2014
-
[23]
Methods of measuring the concentration of wealth
Lorenz M.O. Methods of measuring the concentration of wealth. Quarterly Publications of the American Sta- tistical Association, 1905; 9: 209 (New Series, No. 70);https://doi.org/10.1080/15225437.1905. 10503443(Accessed June 25, 2026)
1905 doi
-
[24]
Sulla misura della concentrazione e della variabilit‘a dei caratteri
Gini C. Sulla misura della concentrazione e della variabilit‘a dei caratteri. Atti del Reale Istituto Veneto di Scienze, Lettere ed Arti 1914; 73: 1203; English translation in Metron - Int. J. Statisrics 2005; 63: 3;https: //www.dss.uniroma1.it/RePec/mtn/articoli/2005-1-1.pdf(...
1914
-
[25]
1906 Manuale di Economia Politica 1906: EGEA – Università Bocconi Editore: Milan, 2006, ISBN 88-8350-084-9
Pareto V . 1906 Manuale di Economia Politica 1906: EGEA – Università Bocconi Editore: Milan, 2006, ISBN 88-8350-084-9
1906
-
[26]
and Yakovenko V
Dragulescu A. and Yakovenko V . Statistical mechanics of money. Eur. Phys. J. B. 2000; 17: 723
2000
-
[27]
and Rosser J.B
Yakovenko V .M. and Rosser J.B. Colloquium: Statistical mechanics of money, wealth, and income. Rev. Mod. Phys. 2009; 81: 1703
2009
-
[28]
Global Inequality in Energy Consumption from 1980 to 2010
Lawrence S., Liu Q, and Yakovenko V .M. Global Inequality in Energy Consumption from 1980 to 2010. Entropy 21013; 15: 5565
1980
-
[29]
and Yakovenko V .M
Semieniuk G. and Yakovenko V .M. Historical evolution of global inequality in carbon emissions and footprints versus redistributive scenarios. J. Cleaner Production 2020; 264: 121420
2020
-
[30]
and Wang H
Boghosian B.M., Devitt-Lee A., Johnson M., Li J., Marcq J.A. and Wang H. Oligarchy as a phase transition: the effect of wealth-attained advantage in a Fokker–Planck description of asset exchange. Physica A 2017; 476: 15
2017
-
[31]
The inescapable casino
Boghosian B.M. The inescapable casino. Sci. American, 2019; November: 71
2019
-
[32]
and Shepelyansky D.L
Frahm K.M. and Shepelyansky D.L. Nonlinear perturbation of Random Matrix Theory. Phys. Rev. Lett. 2023; 131: 077201
2023
-
[33]
and Shepelyansky D.L
Frahm K.M. and Shepelyansky D.L. Dynamical thermalization and turbulence in social stratification models. 2026; arXiv:2603.24190 [cond-mat.stat-mech]
2026
-
[34]
List of countries by energy consumption and production.https://en.wikipedia
Wikipedia contributors. List of countries by energy consumption and production.https://en.wikipedia. org/wiki/List_of_countries_by_energy_consumption_and_production(Accessed June 25, 2026). 12
2026
Reviewed July 9, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.