{"id":"6b7bb6a7-0240-48e2-a356-9a247bca728b","arxiv_id":"2412.00594","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"First 3D N-body simulations of a rotating, shear-free Newtonian universe show more than 1% anisotropy between polar and equatorial expansion rates at maximal allowed rotation.","lead":"This paper runs the first 3D cosmological simulations of a universe that rotates like a rigid body, using a code that avoids periodic boundaries. It finds that the expansion rate differs by more than 1% between directions parallel and perpendicular to the rotation axis, a potentially observable effect in precision cosmology.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The >1% anisotropy may be an artifact of the Sec. 2.3 curvature-compensation rescaling: at z_init=63 the rescaling already suppresses the perpendicular Hubble flow by ~0.7%, comparable to the reported effect, and the paper concedes the method is unsatisfactory.","rationale":"The reader's conditional verdict already identifies the curvature-compensation rescaling as the weakest point, and my analysis agrees: the concern is not the existence of rotation in the simulation but the fact that the quantitative prediction is controlled by an initial-condition rescaling whose physical basis is conceded to be unsatisfactory. The specific quantitative addition here is that the rescaling is not a small correction: it changes the perpendicular Hubble flow by ~0.66% at z_init=63 for the maximal rotation case, i.e., the same order as the reported 1% effect. Since the simulation uses an unperturbed glass and the scale factors are ratios to the initial configuration, this kinematic imprint alone could produce much of the claimed anisotropy. The proposed control run, which evaluates the same rescaling at z=0 instead of z_init, would settle whether the effect is robust. Because the paper is transparent, provides a reproducible code, and presents the result as a first step, a conditional verdict remains appropriate; no change to the reader's verdict is needed.","tokens_in":7680,"tokens_out":19842,"duration_ms":212376,"concrete_test":"Rerun the omega0=10^-3 Gyr^-1 EdS and Lambda-CDM cases with the same scaling prescription (Eq. 2) but evaluate s at z=0 instead of z_init=63: s(z=0) ~ 0.9998, whereas s(z_init) ~ 0.987. If the final polar/equatorial scale-factor anisotropy drops from the reported ~1-2.6% to well below 1%, the headline is an artifact of applying the rescaling at the initial epoch; if it remains >1%, the concern is refuted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim depends directly on the ad hoc curvature compensation in Sec. 2.3. After adding rigid rotation, all perpendicular velocity components are multiplied by sqrt(s), with s = V_perp,nr^2 / V_perp,r^2. For the maximal case (omega0=10^-3 Gyr^-1, z_init=63), omega_init = omega0(1+z)^2 = 4.1 Gyr^-1 and H(z_init) ~ 35.4 Gyr^-1, so s ~ 0.987 and the perpendicular Hubble flow is reduced by ~0.66% relative to the parallel Hubble flow. Because a(t) in Eq. (3) is measured against initial distances, this initial velocity anisotropy is imprinted on the final a_perp/a_parallel ratio and is already of the same order as the headline >1% anisotropy. Moreover, the choice of epoch z_init=63 is arbitrary: evaluating the same formula at z=0 gives s ~ 0.9998, i.e. almost no correction. The paper explicitly states (Sec. 4) that this method is unsatisfactory and that a curvature-like term is missing from the Newtonian Friedmann equation. The simulation therefore cannot currently separate a physical rotational back-reaction from the choice of the rescaling prescription.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":7945,"tokens_out":3646,"duration_ms":135855,"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":[{"comment":"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.","section":"§2.3, Eq. (2), and Fig. 2"},{"comment":"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.","section":"§3, Figs. 2 and 3"},{"comment":"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.","section":"§4, Eq. (4), and Fig. 3"}],"minor_comments":[{"comment":"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.","section":"§2.3, Eq. (1)"},{"comment":"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²∥.","section":"§3 and Fig. 3"},{"comment":"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.","section":"Fig. 3 caption"},{"comment":"The particle number '15 × 216' is ambiguous. If it means 15 × 2^16, please write it in that form or give the integer value.","section":"Table 1"},{"comment":"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.","section":"Abstract and §1"},{"comment":"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).","section":"§4"}],"recommendation":"major_revision","confidential_remarks":"The qualitative direction of the anisotropy is likely real, but the quantitative headline claim is not yet supported because the initial-condition rescaling in Sec. 2.3 introduces a systematic effect of the same order as the reported signal and is explicitly called unsatisfactory by the authors. This is fixable with additional simulations and a derived correction, so I recommend major revision rather than rejection, but the revision must contain new analysis, not just new text."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Kati — quick take. This is the first paper I know of that runs actual 3D N-body simulations of a shear-free rotating universe, using the StePS compactified infinite-volume scheme. That is a real step beyond the earlier Gödel-universe visualization work by Buser et al., and the method alone earns a read. The qualitative result — polar scale factor ends up smaller than equatorial — is internally consistent and clearly presented.\n\nThe quantitative headline is another matter. The >1% anisotropy is not robust. The Sec. 2.3 rescaling that adjusts perpendicular velocities to match the non-rotating kinetic energy is ad hoc, and the authors say so themselves in Sec. 4: they call it unsatisfactory and admit a curvature-like term is missing from the Newtonian Friedmann equation. The stress-test arithmetic checks out: at z_init=63, omega0=1e-3 Gyr^-1 gives omega_init ~ 4.1 Gyr^-1 and H(z_init) ~ 35.4 Gyr^-1, so s ~ 0.987 and the perpendicular Hubble flow is suppressed by ~0.65% at the start — already comparable to the reported final anisotropy. The choice of z_init=63 looks arbitrary; at z=0 the same formula gives s ~ 0.9998. So the headline number may be mostly a rescaling artifact rather than a physical rotational back-reaction. There are also no error bars or convergence tests, and the asserted relation H_perp^2 = H_parallel^2 + Omega^2 is not derived; the authors themselves show it does not hold and then fit a c1/a^2 + c2/a curve, adding two free parameters post hoc.\n\nNone of this is a takedown. The paper is honest about its gaps, the method is new, and the code is open-source with parameters given, so the result is in principle reproducible. The StePS self-citations are legitimate — it is the code they actually use. I would not cite the quantitative claim as evidence, but I would send this to a competent referee rather than desk-reject it. It is a legitimate first step, and the referee can ask for error estimates, convergence tests, and a proper treatment of the curvature-like term. Bring it to reading group if you want a case study in how an acknowledged ad hoc step can contaminate a headline number.","headline":"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.","tokens_in":8518,"tokens_out":3095,"would_cite":false,"duration_ms":100820,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["rotating universe","shear-free rotation","cosmological anisotropy","N-body simulation","Gödel metric","expansion rate","compactified simulations","Newtonian cosmology"],"falsifier":"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.","tokens_in":7446,"feed_emoji":"🌀","tokens_out":7663,"duration_ms":68312,"temperature":0.7,"pith_summary":"This paper reports N-body simulations of a universe that rotates as a rigid, shear-free whole, in the spirit of Gödel's rotating spacetime. The central result is that at the maximum angular velocity that avoids closed timelike curves within the horizon, $\\omega_0 \\approx 10^{-3}$ Gyr$^{-1}$, the cosmic expansion becomes direction-dependent: the equatorial expansion rate exceeds the polar one by more than 1% in both the Einstein–de Sitter and $\\Lambda$CDM cosmologies. If correct, this means a global rotation of the kind still allowed by CMB isotropy constraints would leave a measurable imprint on the expansion history, relevant for precision cosmology. The authors use a compactified infinite-universe simulation technique that avoids the toroidal topology of standard N-body boxes, and they measure scale factors in polar sectors versus an equatorial belt to isolate the anisotropy.","feed_headline":"Rotating universes expand 1% faster at the equator","feed_subtitle":"Shear-free rotation at the largest allowed rate makes equatorial expansion outpace the polar direction.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the compactified infinite-universe simulation method (stereographic projection onto a 4D hypersphere) that makes global rotation simulable without toroidal topology.","marker":"Rácz et al. (2018)"},{"why":"Introduces the StePS N-body code used for the simulations.","marker":"Rácz et al. (2019)"},{"why":"The rotating, shear-free spacetime that inspires the rotation prescription used here.","marker":"Gödel (1949)"},{"why":"Constrains anisotropic Bianchi cosmologies, leaving shear-free Gödel-like rotation as the permitted class of rotating models.","marker":"Planck Collaboration et al. (2016)"},{"why":"Argues that shear-free global rotation remains compatible with observations, providing the theoretical motivation for the model.","marker":"Obukhov (1992)"}],"fun_headline_variants":["Rotating universe: poles lag equator by 1%","Gödel-inspired spin skews cosmic expansion","Rotation makes universe expand unevenly","Simulated rotating cosmos shows 2.6% anisotropy","Shear-free spin alters Hubble expansion rates"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Rotating universe: poles lag equator by 1%","Gödel-inspired spin skews cosmic expansion","Rotation makes universe expand unevenly","Simulated rotating cosmos shows 2.6% anisotropy","Shear-free spin alters Hubble expansion rates"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000114,"raw_usage":{"total_tokens":1021,"prompt_tokens":854,"completion_tokens":167,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":470,"completion_tokens_details":{"reasoning_tokens":96}},"tokens_in":470,"tokens_out":167,"duration_ms":2832,"temperature":1.0,"reasoning_tokens":96,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T05:11:52.332983+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[],"review_version":1}