REVIEW 3 major objections 6 minor 26 references
Simulating Rotating Newtonian Universes
T0 review · 3 major / 6 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read A rotating universe, with the maximal shear-free rotation allowed by causality, expands more than 1% faster in the equatorial plane than along its polar axis in Newtonian N-body simulations.
desk verdict First 3D N-body simulations of a rotating universe, but the headline 1% anisotropy is probably an artifact of the ad hoc rescaling the authors themselves call unsatisfactory. 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 central object is a rotating Newtonian Friedmann equation with two scale factors, $a_\parallel$ and $a_\perp$, representing expansion parallel and perpendicular to the rotation axis. The simulations use the StePS N-body code, which compactifies the infinite universe via stereographic projection onto a 4D hypersphere, avoiding the periodic toroidal topology that would otherwise suppress global rotation. To flatten the rotating universe, the authors rescale perpendicular velocities by $s = V_{\perp,\mathrm{nr}}^2 / V_{\perp,\mathrm{r}}^2$, restoring the total kinetic energy of the non-rotating case; this rescaling is the step the paper later calls unsatisfactory.
What would settle it
Run the maximal-rotation simulation without applying the perpendicular-velocity rescaling of Sec. 2.3, leaving the naive rotational velocities; if $H_\perp^2 - H_\parallel^2$ then converges to $\omega_0^2$ at the present time, the reported anisotropy is an artifact of the compensation rather than a property of rotating Newtonian cosmologies.
Extended reading notes
Core claim
The paper's central claim is that a shear-free, rigid rotation of the universe, modeled on the Gödel metric, produces two distinct expansion rates: one along the rotation axis ($a_\parallel$) and one in the equatorial plane ($a_\perp$). In the simulation with the largest permitted rotation, $a_\parallel$ falls below $a_\perp$ by roughly 1% in Einstein–de Sitter and about 2.6% in $\Lambda$CDM at the present time, and the squared Hubble parameters differ by a comparable anisotropy. The authors further find that $H_\perp^2 - H_\parallel^2$ does not converge to the naive expectation $\omega_0^2$, which they interpret as evidence that a curvature-like term is missing from the Newtonian Friedmann equation for a rotating universe, so the exact analytic description is not yet complete.
Load-bearing premise
The load-bearing premise is that rescaling the perpendicular velocities by $s = V_{\perp,\mathrm{nr}}^2 / V_{\perp,\mathrm{r}}^2$ exactly flattens the rotating universe, restoring the non-rotating total kinetic energy; the paper itself concedes this compensation is unsatisfactory and that a curvature-like term is missing, so the 1% anisotropy may be an artifact of that choice.
Editorial extensions
If this is right
- If a global shear-free rotation exists at its maximal allowed rate, the Hubble expansion is anisotropic at the percent level, so precision distance–redshift surveys must account for a direction-dependent scale factor.
- The measured difference $H_\perp^2 - H_\parallel^2$ deviates from $\omega_0^2$, meaning the Newtonian Friedmann equation needs an additional curvature-like term; identifying that term would give a closed analytic model for rotating universes.
- The simulation setup can be extended to perturbed initial conditions, allowing tests of how rotation seeds large-scale structure alongside the expansion anisotropy.
- The predicted equatorial-vs-polar anisotropy is directly testable with future surveys that map cosmic expansion in different directions, for example via supernova distances or the kinetic Sunyaev–Zel'dovich effect.
Reading between the lines
- Because the paper's curvature compensation is acknowledged to be ad hoc, the reported 1% might be an artifact; a rerun without the kinetic-energy rescaling would distinguish a physical anisotropy from a numerical one.
- If the missing curvature term is later identified, the same simulations could be re-analyzed to test whether the anisotropy scales exactly as $\omega_0^2$ or has a different dependence, making a sharp prediction for observations.
- The polar-sector/equatorial-belt measurement scheme could be translated into observable coordinates: a rotating Newtonian universe would produce a dipole-like pattern in the Hubble constant inferred from large-scale structure, which current surveys could search for.
- The authors' approach suggests that Gödel-type rotation, though shear-free, may still be observationally distinguishable through its distinct expansion history, even if it leaves the CMB temperature isotropic.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper presents full 3D Newtonian N-body simulations of a rotating, shear-free universe using the StePS compactified code. The authors set up unperturbed glass initial conditions with a rigid-body angular velocity ω0 about the z-axis, run EdS and ΛCDM cosmologies, and measure scale factors and Hubble parameters in polar sectors (parallel to the rotation axis) and an equatorial belt (perpendicular). For the maximal rotation without closed timelike curves within the horizon, ω0 ≈ 10^-3 Gyr^-1, they report that the perpendicular expansion exceeds the parallel expansion by about 1%, and they interpret the difference through a modified Newtonian Friedmann equation with a curvature-like term. The paper also reports that their initial rescaling prescription used to compensate for this apparent curvature is unsatisfactory, and it defers the full description to future work.
Significance. If the reported anisotropy is robust, the paper is significant: it is the first full 3D numerical realization of a rotating Newtonian universe, it connects to Gödel-type spacetimes and to observational anisotropy claims, and it demonstrates a possible back-reaction of global rotation on the expansion. The use of the StePS code and the stated reproducibility of the simulations are strengths. However, the central quantitative claim, the >1% anisotropy, rests on an ad hoc initial-condition rescaling that the authors themselves call unsatisfactory, and the paper provides no error bars or convergence tests. The qualitative sign of the effect (a∥ < a⊥) appears clear from the figures, but the quantitative precision claimed in the abstract is not yet established.
major comments (3)
- [§2.3, Eq. (2), and Fig. 2] The headline >1% anisotropy is not robust because it depends directly on the perpendicular-velocity rescaling s = V⊥,nr² / V⊥,r² applied at z_init = 63. For the maximal case ω0 = 10^-3 Gyr^-1, Ω_init = ω0(1+z_init)² = 4.1 Gyr^-1, while H(z_init) ≈ 35.4 Gyr^-1 in EdS, giving s ≈ 1/(1+(Ω/H)²) ≈ 0.987 and √s ≈ 0.993. Since a(t) in Eq. (3) is normalized to initial distances, this initial 0.7% suppression of perpendicular velocities is imprinted on the final a⊥/a∥ ratio and is of the same order as the reported ~1% anisotropy. The choice z_init = 63 is not justified; evaluating the same prescription at z = 0 gives s ≈ 0.9998, almost no correction. The paper explicitly concedes in §4 that this method is unsatisfactory and that a curvature-like term is missing from the Newtonian Friedmann equation. The quantitative central claim can therefore not currently be separated from the choice of rescaling prescription. The authors should present results without the rescaling, test sensitivity to z_init, and ideally derive the correction from the equations of motion rather than imposing it ad hoc.
- [§3, Figs. 2 and 3] The claim of 'approximately 1%' anisotropy is made without any error bars or uncertainty quantification. Each cosmology is simulated once, and the spread across opening angles is shown only as a family of curves. The reported deviation is close to the size of the initial-condition rescaling effect discussed above, so the absence of statistical or systematic error estimates is load-bearing for the precision claim. At minimum, the authors should state the variance across the opening-angle ensembles and, if possible, run multiple realizations with different glass configurations to assess particle-shot-noise effects.
- [§4, Eq. (4), and Fig. 3] The relation H⊥² = H∥² + Ω² is asserted without derivation, and Fig. 3 explicitly shows that the simulation does not converge to the expected value ω0² at the present time, which the authors attribute to a missing curvature-like term. Because this relation is used to justify the functional form of the fit c1/a∥² + c2/a∥ and to interpret the back-reaction, the lack of a derivation is a substantive gap. The authors should derive the modified Friedmann equation from the Newtonian equations in the rotating frame, or at least provide the explicit form of the missing term, before using Eq. (4) to interpret the simulation results.
minor comments (6)
- [§2.3, Eq. (1)] The quantities V⊥,nr² and V⊥,r² are described as means of squared velocities, but the notation does not make this explicit. Please define them as ensemble averages or add overbars.
- [§3 and Fig. 3] The text states that H²∥ - H²⊥ is calculated, while Eq. (4) implies H⊥² > H∥² and the figure caption uses an absolute value. Please make the sign convention consistent throughout and state whether Fig. 3 plots |H²∥ - H²⊥| or H²⊥ - H²∥.
- [Fig. 3 caption] The caption text is garbled in places (e.g., '|H2 H2 |', missing superscripts, and axis labels). Please regenerate the figure with clearer mathematical notation and legible axis labels.
- [Table 1] The particle number '15 × 216' is ambiguous. If it means 15 × 2^16, please write it in that form or give the integer value.
- [Abstract and §1] The abstract refers to a 'Gödel-like metric', but the simulations are Newtonian. Please clarify in the introduction whether the Gödel metric is only a motivation for the velocity field or whether a relativistic correspondence is being claimed.
- [§4] The sentence 'H²⊥ - H²∥ does not converge to ω²0 at present' needs a clearer statement of the expected time dependence: since Ω decays as a^-2, the expected present-day value is ω0², so the mismatch shown in Fig. 3 should be stated explicitly as a discrepancy with Eq. (4).
Circularity Check
No significant circularity — the reported >1% anisotropy is a measured N-body output, not a fitted target; the Sec. 2.3 rescaling is a stated initial-condition choice, and the Sec. 4 limitation is an honesty caveat, not a self-referential derivation.
full rationale
The paper's central claim (anisotropy between polar and equatorial expansion rates exceeding 1%) is obtained by integrating N-body trajectories and then measuring scale factors and Hubble parameters from particle positions and velocities. Equation (3) defines a(t) as a mass-weighted ratio of distances; no term in that definition is set equal to the claimed anisotropy, and no output quantity is fed back into the input parameter s. The perpendicular-velocity rescaling in Eq. (2), s = V_perp,nr^2 / V_perp,r^2, is an initial-condition prescription chosen to conserve total kinetic energy; it is computed from initial velocity moments before the simulation runs, not fitted to the final anisotropy, so it does not fall under 'fitted input called prediction.' The c1/a^2 + c2/a curve in Fig. 3 is explicitly presented as a best fit to the data, not as a predicted relation, and the paper openly notes that the square root of H^2_perp - H^2_parallel does not converge to omega0, which is a falsifiable mismatch rather than a construction. The only self-citations are to the open-source StePS code and the compactified-simulation method (Racz et al. 2018, 2019); these are independently published, publicly available tools, not an unverified uniqueness theorem or an ansatz smuggled in by citation. Section 4 concedes that the curvature-compensation method is 'unsatisfactory' and that a curvature-like term is missing from the Newtonian Friedmann equation; while this is a genuine limitation affecting physical robustness and interpretation, it is not circularity, because the simulation output is not defined to equal the rescaling factor or the missing term. The skeptic concern that the initial rescaling already imprints a comparable velocity anisotropy is a legitimate model-robustness criticism, but it does not exhibit an equation in the paper that reduces the prediction to its input by construction. Therefore no circular step is identified, and the appropriate score is 0.
Assumptions & free parameters
free parameters (3)
- ω0 (present-day angular velocity) =
0, 1e-5, 5e-5, 1e-4, 5e-4, 1e-3 Gyr^-1
- s (perpendicular velocity rescaling) =
not reported; computed from initial velocities
- c1, c2 (fit coefficients) =
not reported
assumptions (4)
- domain assumption Newtonian gravity with a rigid, shear-free global rotation adequately models a Gödel-like rotating universe.
- domain assumption The StePS stereographic compactification faithfully represents the infinite universe without boundary artifacts.
- ad hoc to paper The perpendicular velocity rescaling s exactly compensates the apparent curvature from rotation.
- domain assumption The no-CTC condition is equivalent to Ω ≪ H0√a, giving ω0 ≈ 1e-3 Gyr^-1.
Cite this review
Pith. "Pith review of Simulating Rotating Newtonian Universes." pith.science (2026). https://pith.science/paper/O3P64VOD
@misc{pith2026241200594,
author = {Pith},
title = {Pith review of: Simulating Rotating Newtonian Universes},
year = {2026},
howpublished = {\url{https://pith.science/paper/O3P64VOD}},
note = {Machine review of arXiv:2412.00594}
}
abstract
We present the results of a novel type of numerical simulation that realizes a rotating Universe with a shear-free, rigid body rotation inspired by a G\"{o}del-like metric. We run cosmological simulations of unperturbed glasses with various degrees of rotation in the Einstein-de Sitter and the $\Lambda$CDM cosmologies. To achieve this, we use the StePS N-body code capable of simulating the infinite Universe, overcoming the technical obstacles of classical toroidal (periodic) topologies that would otherwise prevent us from running such simulations. Results show a clear anisotropy between the polar and equatorial expansion rates with more than $1\%$ deviation from the isotropic case for maximal rotation without closed timeline curves within the horizon, $\omega_{0} \approx 10^{-3}$ Gyr$^{-1}$; a considerable effect in the era of precision cosmology.
Reference graph
Works this paper leans on
-
[1]
Barrow JD, Juszkiewicz R, Sonoda DH (1985) Universal rotation - How large can it be? Monthly Notices of the Royal Astronomical Society 213:917--943. doi:10.1093/mnras/213.4.917
-
[2]
Visualization of the G\"odel universe
Buser M, Kajari E, Schleich WP (2013) Visualization of the G \"o del universe . New Journal of Physics 15(1):013063. doi:10.1088/1367-2630/15/1/013063, https://arxiv.org/abs/1303.4651 arXiv:1303.4651 [gr-qc]
work page Pith review arXiv 2013
-
[3]
Monthly Notices of the Royal Astronomical Society 162:307
Collins CB, Hawking SW (1973) The rotation and distortion of the Universe . Monthly Notices of the Royal Astronomical Society 162:307. doi:10.1093/mnras/162.4.307
-
[4]
The Astrophysical Journal 886(1):37
Daly RA (2019) Black Hole Spin and Accretion Disk Magnetic Field Strength Estimates for More Than 750 Active Galactic Nuclei and Multiple Galactic Black Holes . The Astrophysical Journal 886(1):37. doi:10.3847/1538-4357/ab35e6, https://arxiv.org/abs/1905.11319 arXiv:1905.11319 [astro-ph.HE]
arXiv 2019
-
[5]
Astrophysical Journal 292:371--394
Davis M, Efstathiou G, Frenk CS, et al (1985) The evolution of large-scale structure in a universe dominated by cold dark matter . Astrophysical Journal 292:371--394. doi:10.1086/163168
doi:10.1086/163168 1985
-
[6]
Annalen der Physik 369(3):253--287
Ewald PP (1921) Die Berechnung optischer und elektrostatischer Gitterpotentiale . Annalen der Physik 369(3):253--287. doi:10.1002/andp.19213690304
-
[7]
Gamow G (1946) Rotating Universe? Nature 158(4016):549. doi:10.1038/158549a0
-
[8]
Reviews of Modern Physics 21(3):447--450
G \"o del K (1949) An Example of a New Type of Cosmological Solutions of Einstein's Field Equations of Gravitation . Reviews of Modern Physics 21(3):447--450. doi:10.1103/RevModPhys.21.447
Show all 26 references
-
[9]
Monthly Notices of the Royal Astronomical Society 142:129
Hawking S (1969) On the rotation of the Universe . Monthly Notices of the Royal Astronomical Society 142:129. doi:10.1093/mnras/142.2.129
1969 doi
-
[10]
Astrophysical Journal Supplement 75:231
Hernquist L, Bouchet FR, Suto Y (1991) Application of the Ewald Method to Cosmological N-Body Simulations . Astrophysical Journal Supplement 75:231. doi:10.1086/191530
1991 doi
-
[11]
U ber eine station \
Lanczos K (1924) \"U ber eine station \"a re Kosmologie im Sinne der Einsteinschen Gravitationstheorie . Zeitschrift fur Physik 21(1):73--110. doi:10.1007/BF01328251
1924 doi
-
[12]
Physical Review D 108(12):123533
Mc Conville R, \'O Colg \'a in E (2023) Anisotropic distance ladder in Pantheon+supernovae . Physical Review D 108(12):123533. doi:10.1103/PhysRevD.108.123533, https://arxiv.org/abs/2304.02718 arXiv:2304.02718 [astro-ph.CO]
2023 arXiv
-
[13]
Astronomy & Astrophysics 636:A15
Migkas K, Schellenberger G, Reiprich TH, et al (2020) Probing cosmic isotropy with a new X-ray galaxy cluster sample through the L _ X -T scaling relation . Astronomy & Astrophysics 636:A15. doi:10.1051/0004-6361/201936602, https://arxiv.org/abs/2004.03305 arXiv:2004.03305 [as...
2020 arXiv
-
[14]
doi:10.1140/epjs/s11734-021-00199-6, https://arxiv.org/abs/2106.03119 arXiv:2106.03119 [astro-ph.CO]
Mohayaee R, Rameez M, Sarkar S (2021) Do supernovae indicate an accelerating universe? European Physical Journal Special Topics 230(9):2067--2076. doi:10.1140/epjs/s11734-021-00199-6, https://arxiv.org/abs/2106.03119 arXiv:2106.03119 [astro-ph.CO]
2021 arXiv
-
[15]
General Relativity and Gravitation 24(2):121--128
Obukhov YN (1992) Rotation in cosmology . General Relativity and Gravitation 24(2):121--128. doi:10.1007/BF00756780
1992 doi
- [16]
-
[17]
Nature 240(5379):298--299
Pathria RK (1972) The Universe as a Black Hole . Nature 240(5379):298--299. doi:10.1038/240298a0
1972 doi
-
[18]
Planck Collaboration , Ade PAR, Aghanim N, et al (2016) Planck 2015 results. XVIII. Background geometry and topology of the Universe . Astronomy & Astrophysics 594:A18. doi:10.1051/0004-6361/201525829, https://arxiv.org/abs/1502.01593 arXiv:1502.01593 [astro-ph.CO]
2016 arXiv
-
[19]
Planck Collaboration , Aghanim N, Akrami Y, et al (2020) Planck 2018 results. VI. Cosmological parameters . Astronomy & Astrophysics 641:A6. doi:10.1051/0004-6361/201833910, https://arxiv.org/abs/1807.06209 arXiv:1807.06209 [astro-ph.CO]
2020 arXiv
-
[20]
Monthly Notices of the Royal Astronomical Society 477(2):1949--1957
R \'a cz G, Szapudi I, Csabai I, et al (2018) Compactified cosmological simulations of the infinite universe . Monthly Notices of the Royal Astronomical Society 477(2):1949--1957. doi:10.1093/mnras/sty695, https://arxiv.org/abs/1711.04959 arXiv:1711.04959 [astro-ph.CO]
2018 arXiv
-
[21]
Astronomy and Computing 28:100303
R \'a cz G, Szapudi I, Dobos L, et al (2019) StePS: A multi-GPU cosmological N-body Code for compactified simulations . Astronomy and Computing 28:100303. doi:10.1016/j.ascom.2019.100303, https://arxiv.org/abs/1811.05903 arXiv:1811.05903 [astro-ph.CO]
2019
-
[22]
doi:10.1103/PhysRevLett.117.131302, https://arxiv.org/abs/1605.07178 arXiv:1605.07178 [astro-ph.CO]
Saadeh D, Feeney SM, Pontzen A, et al (2016) How Isotropic is the Universe? Physical Review Letters 117(13):131302. doi:10.1103/PhysRevLett.117.131302, https://arxiv.org/abs/1605.07178 arXiv:1605.07178 [astro-ph.CO]
2016 arXiv
-
[23]
Astrophysical Journal 143:689
Silk J (1966) Local Irregularities in a G \"O DEL Universe . Astrophysical Journal 143:689. doi:10.1086/148551
1966 doi
-
[24]
Monthly Notices of the Royal Astronomical Society 147:13
Silk J (1970) The instability of a rotating universe . Monthly Notices of the Royal Astronomical Society 147:13. doi:10.1093/mnras/147.1.13
1970 doi
-
[25]
, " * write output.state after.block = add.period write newline
ENTRY address archive author booktitle chapter doi edition editor eid eprint howpublished institution journal key keywords month note number organization pages publisher school series title type url volume year archivePrefix primaryClass adsurl adsnote version label extra.labe...
-
[26]
write newline
" write newline "" before.all 'output.state := FUNCTION add.period duplicate empty 'skip "." * add.blank if FUNCTION if.digit duplicate "0" = swap duplicate "1" = swap duplicate "2" = swap duplicate "3" = swap duplicate "4" = swap duplicate "5" = swap duplicate "6" = swap dupl...
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.