{"id":"1dcb4b90-89ea-47b0-a17e-0a4812f90799","arxiv_id":"2607.27244","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"A shell with N internal particles has a piecewise coordinate distribution set by particle masses and a momentum distribution set by degrees of freedom, with universal RMS displacement Σm_i² R² / 3(M+m_tot)².","lead":"Exact coordinate and momentum distributions are derived for a movable shell containing a finite number of particles, under energy and momentum conservation. Comparing the two distributions can reveal how many internal particles and degrees of freedom the system has, and the shell wanders chaotically even in vacuum.","discovery_kind":"extension","skeptic_critique":{"model":"moonshotai/kimi-k3","headline":"The projected density in Eq. 6 has the wrong exponent sign: the inverse Laplace step Eq. 5→6 yields (…)^{−1/2}, not √(…). Final formulas are nonetheless correct — the error is in the written derivation, not the conclusions.","rationale":"The reader flagged ergodicity as the weakest assumption and called the derivations \"standard microcanonical projections carried through carefully.\" The ergodicity point is real but openly acknowledged by the authors, numerically probed at N=3, 5, 15, and the N=2 failure is demonstrated in the paper itself (Fig. 3) — it is a stated limitation, not a hidden one, and the headline claims are explicitly restricted to N>2 with MD support. The more load-bearing issue sits inside the analytic machinery the reader endorsed: the inverse Laplace step Eq. 5→6 flips the exponent sign, making Eq. 6 (and Eq. 14) incorrect as printed and internally inconsistent with Eq. 17's exponent and with Eq. 20's N=2 divergent behavior. This is the single most load-bearing concern because Eq. 6 is the fountainhead of every subsequent formula — yet it does not overturn the conclusions, because (a) the correct −1/2 exponent is independently derivable by elementary delta-function elimination or density-of-states scaling, (b) the final exponents, normalizations, the universal RMS formula (verified here in the M→∞ and N=1 limits), and the 2^N−1 segment count all check out, and (c) the paper's own simulations confirm the final distributions. That is precisely a CONDITIONAL situation: the paper should be accepted contingent on correcting the §III derivation (replace the formal momentum-Laplace step with a Fourier treatment or direct elimination, fix Eq. 6 to the −1/2 power, and adjust Eq. 14's integrand notation), since as printed the chain from Eq. 1 to Eq. 19 does not close. I disagree with the reader's identification of the weakest point: the acknowledged ergodicity assumption is less load-bearing than an unflagged sign error in the central derivation, and the reader's \"correctness_risk: low / derivations carried through carefully\" assessment is not accurate as the text stands — though the error is repairable without changing any result, which is why the adjustment is to CONDITIONAL rather than REJECT.","tokens_in":33986,"tokens_out":1486,"duration_ms":474534,"concrete_test":"Compute the projected density for the minimal case (shell + N=2 particles, D=1; three momenta total) by direct delta-function elimination over the two particle momenta: ρ(P) = Σ_roots 1/|∂(g_E,g_P)/∂(p_1,p_2)| ∝ [2(E_tot−P²/2M)(m_1+m_2) − P²]^{−1/2}, diverging at the allowed-momentum edges. Compare with Eq. 6 applied to the same system (√, vanishing at edges) and with Eq. 20 at D=1, N=2 (exponent −1/2, divergent, as the paper itself remarks). If direct elimination gives −1/2, Eq. 6 and Eq. 14 must be corrected; then verify Eqs. 17–19 follow only from the corrected exponent.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The whole paper rests on the projected microcanonical density Eq. 6, obtained by Laplace-transforming Eq. 1 in (E_tot, P_tot) and integrating out two momenta. Two problems. (1) The Laplace transform in momentum is only formal: integrating e^{−s2·p} against the Gaussian energy factor produces e^{+a s2²/2s1}, which has no inverse Laplace transform; the legitimate route is a Fourier transform in momentum, which does reproduce Eq. 5's structure. (2) Decisively, the step Eq. 5 → Eq. 6 is wrong as printed: L^{-1}[s1^{−1/2} e^{−B s1}] = (E−B)^{−1/2}/√π, an inverse square root, whereas Eq. 6 writes +√(2(E_tot−E1)(m_i+m_j) − (P_tot−P1)²). The correct exponent is −1/2, as direct delta-function elimination confirms: integrating δ(E−H)δ(P_tot−Σp) over p_i, p_j, each quadratic root contributes 1/|f′|, giving [2(E−E1)(m_i+m_j) − (P_tot−P1)²]^{−1/2}.\n\nThe error propagates into Eq. 14 (integrand written as √) but is silently corrected by Eq. 17: the iteration rule of Eq. 16 adds 1/2 to the exponent per momentum integrated, so starting from Eq. 14's +1/2 one would arrive at exponent (N−1)/2 after the N−2 x-integrations, not the N/2−3/2 written in Eq. 17. The chain only closes with the correct −1/2 start: −1/2 + (DN−2)/2 = (DN−3)/2, matching Eqs. 19–20. Crispest internal contradiction: for D=1, N=2, Eq. 6 IS the shell momentum marginal (only P remains after elimination), and it vanishes at the boundary (√), while the paper's own Eq. 20 at N=2 diverges at the boundary (exponent −1/2), as the text itself notes. Same object, opposite exponents.\n\nThe headline results survive: the (DN−3)/2 exponent, the normalization constants, and Eq. 22 all check out via an independent density-of-states/coarea derivation, the N=3, D=1 constant-distribution prediction, and the MD comparisons. Section IV is unaffected: it needs only that the momentum integral of the projected density is a coordinate-independent constant, true for either sign. So this is a derivation-level flaw requiring repair, not a resu","agreement_with_reader":"disagree"},"referee_report":{"model":"moonshotai/kimi-k3","summary":"The paper studies an isolated composite object — a movable 1D shell of mass M containing N internal particles — in the center-of-mass frame with fixed total energy and zero total momentum. From the microcanonical measure projected onto the constraint shell (Eqs. 1–6), the authors derive: (i) the shell coordinate distribution ρ(X) as the volume of admissible internal configurations (Eqs. 7–8), shown to consist of 2^N − 1 polynomial branches with an explicit closed form for N = 3 (Eq. 9) and a simple universal second moment ⟨X²⟩ = (Σmᵢ²)R²/[3(M+m_tot)²] (Eq. 12); (ii) the shell momentum distribution ρ(P) ∝ (E_tot − (M+m_tot)P²/2Mm_tot)^{(DN−3)/2} (Eq. 20), depending only on M, m_tot, and the number of internal degrees of freedom, implying a non-equipartition mean shell energy (Eq. 22). Molecular-dynamics simulations for colliding and non-colliding particles agree with the theory for N > 2, while N = 2 fails reproducibly, signaling an extra integral of motion (Fig. 3). The coordinate/momentum distributions together are proposed as a diagnostic of internal particle number versus dimensionality.","tokens_in":18107,"tokens_out":4847,"duration_ms":105973,"significance":"If the results hold, this is a clean, parameter-free contribution to the statistical mechanics of small systems: exact closed-form marginals under simultaneous energy, momentum, and center-of-mass constraints; a universal RMS formula (Eq. 12) that is independent of individual masses and total energy and is verified against independent random mass draws and time averages (Figs. 5–7); and a concrete, falsifiable diagnostic distinguishing internal particle number from dimensionality. The non-equipartition result (Eq. 22) extends the known non-Maxwellian microcanonical phenomenology (Ray & Graben, Ref. [22]) to a structured composite object of direct relevance to nanopeapods and rotaxanes. The numerical cross-checks are extensive and the N = 2 failure is reported honestly. These strengths make the paper worth publishing once the derivation is repaired.","major_comments":[{"comment":"Eq. 6 (§III) has the wrong exponent sign. The inverse Laplace transform of Eq. 5 gives L⁻¹[s₁^(−1/2) e^(−B s₁)] = (E−B)^(−1/2)/√π — an inverse square root — whereas Eq. 6 writes the positive square root. Direct delta-function elimination of pᵢ, pⱼ from δ(E−H)δ(P_tot−Σp) confirms the −1/2 exponent: each quadratic root contributes 1/|f′|, yielding [2(E−E₁)(mᵢ+mⱼ)−(P_tot−P₁)²]^(−1/2). The error is not cosmetic: for D=1, N=2, Eq. 6 IS the shell momentum marginal and as printed vanishes at the support boundary, while the paper's own Eq. 20 at D=1, N=2 has exponent −1/2 and diverges there — a direct internal contradiction. The chain Eq. 14 → 16 → 17 corroborates this: the iteration rule (Eq. 16) adds 1/2 per x-integration, so starting from Eq. 14's +1/2 one would reach (N−1)/2 in Eq. 17, not the written N/2−3/2; only the correct −1/2 start closes the chain. The final results (Eqs. 17–22) are c","section":"§III, Eq. 6; §V, Eqs. 14–17"},{"comment":"Eqs. 2–5: the Laplace transform in the momentum variable is only formal. Integrating e^(−s₂·p) against the Gaussian energy factor produces e^(+a s₂²), which has no inverse Laplace transform (even as a tempered distribution the step is unjustified). The legitimate route is a Fourier transform in the momentum delta function — which does reproduce the structure of Eq. 5 — or a direct delta-function calculation. The Θ(P(p,q)) factor and the claim that nonnegativity of total momentum holds 'without loss of generality' should also be re-examined in this light. Since the paper advertises exact distributions, the derivation should be made rigorous; the conclusions are unaffected.","section":"§III, Eqs. 2–5"},{"comment":"All 'exact' distributions are conditional on the assumption (stated in §IV) that a typical trajectory uniformly fills the constraint shell. For colliding particles this is well supported, but for non-colliding 1D motion with N>2 it is only numerically supported, and the paper itself demonstrates the assumption failing at N=2 within the same model family (Fig. 3). The abstract and conclusions should state explicitly that exactness is conditional on ergodicity on the shell, and the numerical evidence for the absence of extra integrals at N>2 should be strengthened or systematized (e.g., convergence in simulation time, incommensurate vs. near-commensurate mass ratios, or a direct diagnostic such as momentum-set discreteness as used for N=2).","section":"§III–IV; Figs. 3–4"}],"minor_comments":[{"comment":"Reference [26] (Ramshaw 2025) appears in the bibliography but is never cited in the text (the list jumps [25] → [27,28]). Either cite or remove it.","section":"References"},{"comment":"Several equation cross-references are malformed: 'The distribution 7 is non-zero' (missing 'Eq.'); 'Eq.(12) )'; 'Eq20'; 'Eq.9' appears without space; and there is a doubled period after '2³ − 1 segments. .' in §IV.","section":"Throughout"},{"comment":"Fig. 6 caption: 'for cases N = 5(a)' is misprinted as 'or cases'; Fig. 8 caption repeats 'as well as the values of the integrals of motion' twice.","section":"Figs. 6, 8"},{"comment":"Eq. 19 is the D=2 special case of Eq. 20 but this is not stated where Eq. 19 appears; one sentence clarifying that the N y-integrations have already been performed would help the reader follow the exponent bookkeeping.","section":"§V, Eq. 19"},{"comment":"Typesetting of the max/min expressions in Eq. 8 and the x_Max bounds in §IV ('M ax', 'x{1}_Max') is garbled and should be normalized.","section":"§IV, Eq. 8"},{"comment":"Consider adding a second PACS number (e.g., 05.20.-y, classical statistical mechanics) alongside 05.45.-a.","section":"Front matter"}],"recommendation":"major_revision","confidential_remarks":"The final formulas (Eqs. 9, 12, 20–22) appear correct and are well supported by the simulations; the required changes are concentrated in the §III/§V derivation and in qualification of the ergodicity claim, so I expect a revision to be straightforward. I verified the exponent-sign concern independently (both via the Laplace inversion and via direct delta-function elimination) and it is genuine — the written Eq. 6 contradicts the paper's own Eq. 20 at D=1, N=2. Self-citations [36,37] supply only the projection technique, which seems appropriate. Novelty relative to the authors' prior energy/angular-momentum work looks adequate, since the coordinate-distribution construction (2^N−1 branches, Eq. 12) is new."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The useful core here is concrete: under E and P conservation they get the shell coordinate law as the volume of admissible internal configs (2^N−1 polynomial branches, mass-ratio dependent, plateaus when a heavy particle pins the geometry), the momentum law ρ(P)∝(E_tot−(M+m_tot)P²/(2M m_tot))^{(DN−3)/2}, and the universal ⟨X²⟩=(Σ m_i²)R²/[3(M+m_tot)²]. Comparing ρ(X) to ρ(P) really does separate particle number/masses from dimensionality. That package is new relative to their earlier energy-distribution papers and is the reason to read it.\n\nWhat they do well is carry the microcanonical projection through to closed forms (especially the full N=3 coordinate expression) and check it with MD for colliding and non-colliding cases when N>2. The RMS formula is simple, parameter-free, and matches random-mass and equal-mass runs. Non-equipartition of shell energy is stated cleanly. Section IV on the “internal Brownian” motion without a bath is the right physical framing for rotaxanes/nanopeapods-type models.\n\nSoft spots, in proportion. The printed step from the Laplace form to Eq. 6 has the wrong root sign (should be inverse square root); the iteration that produces Eqs. 19–20 only closes if you start from −1/2. Final exponents, normalizations, the N=3 flat momentum distribution, and the MD all match the correct density-of-states answer, so this is a derivation write-up bug, not a wrong headline result—but a referee will catch it and it needs fixing. Ergodicity for non-colliding particles is assumed; they themselves show N=2 has an extra integral and discrete momenta, and for N>2 they only have numerical support. No code/data shipped. The “infer internal structure from ρ(X) vs ρ(P)” claim is plausible and not operationally demonstrated on a realistic observable.\n\nThis is for people who do small-system statistical mechanics, confined ideal gases, or need analytic benchmarks for shell/cluster MD. Not field-changing, but solid and usable. I would send it to peer review; after the Eq. 6 repair and a clearer ergodicity caveat it is an accept-level note.","headline":"Exact shell distributions and a clean universal RMS formula for a movable container with N internal particles; derivation write-up has a fixable sign error, but the final results and numerics hold.","tokens_in":18978,"tokens_out":602,"would_cite":false,"duration_ms":22857,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["05.45.-a"],"model":"grok-4.5","headline":"A movable shell enclosing finitely many particles has exact, mass-dependent coordinate and momentum distributions that reveal the internal structure and produce medium-free chaotic motion.","keywords":["small systems","microcanonical ensemble","coordinate distribution","momentum distribution","movable shell","energy equipartition violation","finite degrees of freedom","chaotic shell motion"],"falsifier":"For N>2 non-colliding particles, measure the shell’s long-time position histogram and check whether it collapses onto the predicted multi-branch volume formula (and whether ⟨X²⟩ equals the mass-squared expression) independent of initial conditions.","tokens_in":18526,"feed_emoji":"⚛️","tokens_out":798,"duration_ms":15003,"temperature":0.7,"pith_summary":"Small isolated systems with a movable outer shell and a finite number of particles inside cannot be treated with ordinary thermodynamics. This paper derives the exact probability distributions for the shell’s position and its momentum when total energy and total momentum are conserved. The position distribution splits into a finite number of polynomial branches set by the particle masses and can contain flat plateaus; its width is given by a simple universal formula involving only the sum of the squared masses. The momentum distribution is a single power-law function whose exponent counts the internal degrees of freedom and automatically violates energy equipartition. Because one distribution depends on particle number while the other depends on dimensionality, comparing the two tells an observer how many particles are inside and how they move. The same impacts that generate these distributions also drive persistent, Brownian-like wandering of the shell even in perfect vacuum.","feed_headline":"Shell position and momentum reveal hidden particle count","feed_subtitle":"Exact distributions and a universal rms formula diagnose the interior even in vacuum","key_machinery":"The projected microcanonical density obtained by eliminating two momenta via Laplace transform (or Liouville), which reduces ρ(X) to a pure configuration-volume integral with mass-dependent limits and reduces ρ(P) to an iterated momentum integral that yields the explicit power-law form.","core_discovery":"Under conservation of total energy and zero total momentum, the shell coordinate distribution equals the (N−1)-dimensional volume of admissible internal configurations and therefore consists of 2^N−1 polynomial branches whose shapes change with mass ratios, while the shell momentum distribution is exactly proportional to (E_tot − (M+m_tot)P²/(2M m_tot)) raised to the power (DN−3)/2; the resulting mean-square displacement of the shell is universally ⟨X²⟩ = (Σ m_i²) R² / [3(M+m_tot)²].","pith_inferences":["The same volume-construction technique should apply to nested shells or to shells free to move in two or three dimensions, yielding higher-dimensional analogues of the branch structure.","If the shell is weakly coupled to an external bath, the vacuum rms formula supplies the natural noise floor that must be exceeded before ordinary Brownian motion dominates.","Discrete momentum spectra seen for N=2 suggest that few-body integrable shells could be used as calibrated momentum filters in nanoscale devices."],"forward_implications":["Shell position statistics alone fix the number and mass ratios of internal particles, independent of motion dimensionality.","Shell momentum statistics alone fix the total number of internal degrees of freedom and the reduced mass.","Comparing the two distributions therefore diagnoses both particle number and dimensionality of an unseen interior.","The shell executes persistent chaotic wandering whose rms width is set solely by internal masses, even in perfect vacuum.","Average shell energy is systematically lower than equipartition whenever the shell is heavier than a single internal degree of freedom."],"fun_headline_variants":["Shell coordinate branches and momentum spectra count hidden particles","Vacuum shell chaos and rms formula diagnose internal particle number","Exact shell distributions reveal composite system interior structure","Position plateaus plus momentum power law expose particle count","Shell coords stay dimension-free while momenta track degrees of freedom"],"cache_read_input_tokens":128,"weakest_assumption_plain":"A typical trajectory is assumed to fill the entire phase-space region allowed by the conserved quantities uniformly; without that filling the analytic histograms would not match the long-time motion.","fun_headline_variants_meta":{"raw":{"variants":["Shell coordinate branches and momentum spectra count hidden particles","Vacuum shell chaos and rms formula diagnose internal particle number","Exact shell distributions reveal composite system interior structure","Position plateaus plus momentum power law expose particle count","Shell coords stay dimension-free while momenta track degrees of freedom"]},"model":"grok-4.5","effort":"low","cost_usd":0.004718,"raw_usage":{"total_tokens":1336,"prompt_tokens":770,"num_sources_used":0,"completion_tokens":77,"cost_in_usd_ticks":47184000,"prompt_tokens_details":{"text_tokens":770,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":489,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":770,"tokens_out":77,"duration_ms":10158,"temperature":1.0,"reasoning_tokens":489,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-31T22:13:09.572318+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"For N>2 non-colliding particles, measure the shell’s long-time position histogram and check whether it collapses onto the predicted multi-branch volume formula (and whether ⟨X²⟩ equals the mass-squared expression) independent of initial conditions.","supporting_citations":[],"review_version":1}