{"id":"c30615c0-e538-4428-92fd-33f5f8235a58","arxiv_id":"2501.12445","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Late-time clusters of classical dipole-conserving fractons map to integrable polygonal billiards for compact interactions, to chaotic motion for non-compact interactions, and generic trajectories exhibit a Janus point with a bidirectional arrow of time.","lead":"Classical dipole-conserving particles, or fractons, tend to cluster at late times. This paper shows the motion inside those clusters is a billiard in a triangle or square, and that trajectories can show a central 'Janus point' with complexity growing in both time directions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 3–1 billiard reduction in §III.B assumes the triangle wall reflects every trajectory, but the alternative transmission branch is never checked; the small-triangle confinement that carries the integrability claim is therefore an unproven assumption, not a derived result.","rationale":"The reader's weakest assumption correctly identifies small-triangle confinement as empirically inferred rather than proven. My concern sharpens and extends this: the missing object is not merely the statistical absence of large-triangle trajectories, but the absence of a branch analysis at the four-fracton walls. The three-fracton problem in Appendix A is treated with a careful discriminant method that sometimes allows transmission and sometimes forces reflection; no such analysis is given in §III.B or §III.C for the four-fracton walls. The reflection rule is simply assumed. Since the claim that the late-time 3–1 and 2–2 clusters are integrable billiards is the paper's primary result, this is the single most load-bearing point. The concrete test I propose would settle it: solve the wall-jump equations and check whether the transmission branch has real solutions for reachable incoming momenta. If the branch check confirms reflection is the only real branch, the small-triangle confinement becomes a theorem rather than a simulation observation, and the central claim is substantially strengthened. If transmission is possible, the paper's strongest claim must be qualified to a subset of initial conditions or to a modified billiard domain, and the verdict would need to become more restrictive than conditional. Because the concern is not yet resolved either way, the reader's conditional verdict remains the appropriate recommendation.","tokens_in":19147,"tokens_out":12219,"duration_ms":140233,"concrete_test":"For the 3–1 wall defined by √12 q3 − √6 q2 = 3, write the full Hamiltonian Eq. (9) with box K, impose the conserved quantities π1 and √2π2 + π3 at the wall together with energy conservation, and solve the jump conditions for both the reflection and transmission branches. Compute the discriminant of the resulting quadratic as a function of incoming (π1, π2, π3) over initial conditions that start in the big-bang Region 3 and form late-time 3–1 states. If the transmission branch is real for any such incoming state, the billiard reduction is incomplete and the integrability claim for 3–1 states fails; if only the reflection branch is real, small-triangle confinement is established analytically. Repeat the same discriminant analysis for the square wall in §III.C.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In §III.B, the 3–1 late-time state is reduced to a billiard in an equilateral triangle by deriving a reflection rule at the wall (e.g., √12 q3 − √6 q2 = 3): π1′ = π1, π2′ = −π2, π3′ = 2√2 π2 + π3. The text explicitly says 'Assuming the trajectory does not escape the triangle,' but it never checks whether the alternative branch of the jump conditions—transmission into the region where the isolated fourth particle becomes coupled—has real solutions. This matters because the three-fracton analysis in Appendix A shows that the analogous boundary jumps have both reflecting and transmitting branches depending on the incoming momenta; the branch is selected by a discriminant, not by assumption. If transmission is possible for any incoming state reachable from the big-bang region, then some 3–1 states leave the small triangle, the billiard domain changes to a larger region with sub-regions where pair terms switch off, and the claimed three integrals of motion and Liouville integrability for 3–1 states fail. The paper's only support for small-triangle confinement is 'This is never observed in simulations,' which is an empirical inference, not a proof. The same gap occurs in the 2–2 square reduction in §III.C, where the reflection rule is similarly stated without checking the transmission branch. Since the central claim of local ergodicity without chaos rests on these two billiard mappings, this missing discriminant check is the most load-bearing unresolved step.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies one-dimensional classical fractons with pair-dependent inertia, continuing earlier work by the same authors. For the four-fracton problem with compact ('box') pair inertia K, it claims that late-time 3–1 and 2–2 cluster states reduce to billiard motion in an equilateral triangle and a square, respectively; because these polygonal billiards are integrable, the cluster dynamics is locally ergodic but not chaotic, while emergent conserved quantities break global ergodicity. The paper also introduces a complexity measure and a Boltzmann entropy, arguing that both increase away from a central 'Janus point', yielding a bidirectional arrow of time. The main technical content is the reduced-coordinate Hamiltonian, reflection rules at cluster boundaries, and numerical simulations of trajectories and Lyapunov exponents.","tokens_in":19458,"tokens_out":9171,"duration_ms":96061,"significance":"If the billiard reduction were rigorously established, this paper would provide a clean Hamiltonian many-body example in which global ergodicity is broken while cluster-level dynamics remains regular (or becomes chaotic for non-compact K). The explicit mapping of interacting fractons to integrable billiards is novel and would be of interest to the statistical-mechanics and fracton communities. The Janus-point/entropy construction also connects fracton dynamics to cosmological arrow-of-time ideas. Strengths include the transparent reduced-coordinate formulation, the concrete reflection rules, the numerical verification of trajectories, and the clear separation of compact and non-compact cases. However, the central reduction relies on an unproven confinement assumption, and the arrow-of-time argument depends on an assumed linear growth law rather than a derivation.","major_comments":[{"comment":"The billiard reduction for the 3–1 state is supported only by the statement 'This is never observed in simulations of 3–1 states' (Section III.B), and the reflection rule is derived after the explicit clause 'Assuming the trajectory does not escape the triangle.' The manuscript never checks whether the boundary jump conditions admit a real transmission branch into the region where the previously isolated fourth fracton becomes coupled. This matters because the analogous three-fracton boundary jumps in Appendix A are resolved by real-solution discriminants (e.g., D1 and D2 in Appendix A), and both reflecting and transmitting branches occur there. The same gap appears in the 2–2 case: the square-wall reflection rule in Section III.C is stated without a discriminant analysis. Since the claimed three integrals of motion and the resulting Liouville integrability in Section IV.A apply only if every trajectory is reflected at the small-triangle or square walls, the confinement and the reflection branch must be proven, not assumed.","section":"III.B and III.C"},{"comment":"The statement 'We now prove this' is not backed by a rigorous proof. The 'three integrals of motion' are not of the same type: the first two (the plane containing the trajectory and the Hamiltonian) are smooth conserved functions, while the third (the finite number of billiard directions) is a discrete invariant coming from the polygon being rational. The cited Liouville integrability theorem for a 2N-dimensional phase space with N integrals does not apply directly to a discrete, non-smooth invariant. The conclusion that equilateral-triangle and square billiards are regular is correct by standard results, but the proof as written is incomplete. Additionally, the phrase 'locally ergodic' is never defined; for these integrable billiards a trajectory is confined to a lower-dimensional invariant set (a fixed direction class and, on unfolding, a torus), so if 'local ergodicity' means exploration of the full cluster phase space, that is false. The authors should state which notion of local ergodicity they intend and prove it for the energy surface restricted to each invariant sector.","section":"IV.A"},{"comment":"The arrow-of-time result rests on the a priori assumption pi = alpha_i + beta_i t (Eq. 22), which is stated as a toy-model input rather than derived from the equations of motion. The complexity P is then shown to grow for t > t_+ (Eq. 24), and the Boltzmann entropy S is defined as a monotone function of P, so the 'second law' in Eq. (31) is a consequence of the choice of macroscopic variable. To support the abstract's claim of a 'generic' bidirectional arrow of time, the paper should either derive Eq. (22) from the dynamics or clearly state that this is a phenomenological model, and it should quantify how generic the observed growth of P is. As it stands, the entropic arrow is tautological once P is observed to increase.","section":"V, Eqs. (22)-(31)"}],"minor_comments":[{"comment":"The abstract and Section IV refer to 'various stadiums', but the paper treats equilateral-triangle and square billiards, not stadium billiards. Please use 'polygonal billiards' or specify the domains.","section":"Abstract, II, IV"},{"comment":"The sentence 'this guarantees only 4 of the 6 exponents to be zero' is imprecise: since Lyapunov exponents in a Hamiltonian system come in pairs (lambda, -lambda), three zero exponents force at least four zeros out of six, leaving room for one nonzero pair. The wording should be corrected.","section":"IV.B"},{"comment":"The arguments based on 'arbitrarily large region' Mathematica RegionPlots are informal; the analytic inequalities in Eqs. (A24)–(A31) are more convincing. Please replace the plot-based claims with the analytic discriminant conditions throughout Appendix A.","section":"Appendix A"},{"comment":"There are recurring typos and notation issues, e.g., 'π′2 = −π2' is written with an unusual prime placement, and the reflection rules in Sections III.B and III.C would be easier to follow if the conserved linear combinations were displayed before the jump rules.","section":"General"},{"comment":"Reference [15] appears incomplete: 'E. Ledberg, Introduction to rational billiards (2018)' lacks a journal or arXiv identifier. Please verify and complete the citation.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The most serious issue is the missing transmission-branch analysis for the four-fracton billiard walls. This is load-bearing for the central claim of local regularity without chaos, but it is likely fixable by performing the discriminant calculation in the style of Appendix A. The entropy/arrow-of-time section is currently more of a phenomenological construction than a theorem. I recommend major revision rather than rejection because the model and the underlying idea are promising, but the manuscript needs a rigorous confinement proof and a clearer statement of the assumptions behind the arrow of time."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The headline: this is a solid, genuinely new contribution to the few-body dynamics of classical dipole-conserving fractons. The billiard reductions for the 3–1 and 2–2 late-time states are the real prize, and they are worked out in enough detail that a referee can check them line by line. The paper also gives a complete classification of three-fracton trajectories that goes beyond the earlier work, and a clean separation between compact and non-compact pair inertia functions.\n\nWhat is actually new: the mapping of 3–1 states to a billiard in an equilateral triangle and 2–2 states to a square, with explicit reflection rules and three conserved quantities; the argument that these polygonal billiards have zero Lyapunov exponents, so local ergodicity coexists with global ergodicity breaking; and the numerical observation of chaotic pairs of Lyapunov exponents for non-compact K. The three-fracton appendix with exact discriminant conditions for reflection vs transmission is careful and reproducible.\n\nWhere I would push back: the billiard reduction for the 3–1 state depends on the claim that late-time trajectories are confined to the 'small triangle' region. That confinement is inferred from simulations ('never observed'), not proved. The stress-test note is right: the paper never checks whether the transmission branch at the triangle wall has real solutions for reachable incoming momenta. The three-fracton analysis in Appendix A shows exactly this kind of branch selection, so the missing discriminant check is a genuine gap, not a stylistic quibble. The same applies to the 2–2 square. This is the most load-bearing soft spot; the integrability conclusion for compact K stands or falls on it.\n\nTwo smaller issues: the Lyapunov exponent numerics have no error bars and the position-momentum vs position-velocity frame subtlety is acknowledged but not fully addressed. And the arrow-of-time section is explicitly a toy model—it assumes linear momentum growth and defines the Boltzmann entropy through the observed complexity P, so the increase of entropy is partly built in. Fine as a heuristic, but it should not be read as a derivation.\n\nThe Liouville integrability phrasing is also a bit loose: polygonal billiards do have zero Lyapunov exponents and non-ergodic motion, but the 'finite number of directions' is a discrete invariant, not a smooth function in involution in the usual sense. That does not change the physics, but a referee should ask for a more precise statement.\n\nBottom line: the paper deserves serious refereeing. The core billiard result is novel, plausible, and checkable, but the missing transmission-branch analysis means I would not accept it as is. The fix is straightforward in principle—do the discriminant calculation for the triangle/square walls the way the three-fracton appendix does.\n\nWho will get value: statistical mechanicians, fracton theorists, and people interested in emergent time arrows. I'd bring it to a reading group, but I wouldn't cite it in my own work until the confinement question is settled.","headline":"A novel billiard reduction for four-fracton late-time states is the real contribution, but the central integrability claim rests on an unproven confinement assumption that a referee should push on.","tokens_in":19971,"tokens_out":3640,"would_cite":false,"duration_ms":33428,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["37J35","37D45","82C05"],"pacs":["05.20.-y","05.45.-a","05.70.Ln"],"model":"deepseek-v4-flash","headline":"Late-time four-fracton clusters are integrable billiards in a triangle or square, so local ergodicity coexists with global non-ergodicity.","keywords":["classical fractons","dipole conservation","Machian dynamics","ergodicity breaking","polygonal billiards","Liouville integrability","Lyapunov chaos","Janus point"],"falsifier":"Run long-time simulations of the four-fracton box-$K$ Hamiltonian from many big-bang-like initial conditions and monitor the reduced coordinates; finding any 3-1 or 2-2 trajectory that enters the large-triangle or outside region (a second particle at rest) or any positive Lyapunov exponent for compact $K$ would refute the integrable-billiard claim.","tokens_in":18954,"feed_emoji":"🎱","tokens_out":10923,"duration_ms":105784,"temperature":0.7,"pith_summary":"The paper studies classical particles that conserve both total momentum and dipole moment, a 'Machian' fracton system in which a particle can move only when another particle is nearby. It claims that the late-time clustered steady states break global ergodicity while keeping a limited form of ergodicity alive inside clusters: for the four-fracton problem with a compact interaction, each two-cluster late-time state maps exactly onto a billiard particle in an equilateral triangle (3-1 state) or a square (2-2 state). Since these are rational polygonal billiards, their motion is Liouville integrable and non-chaotic, so local exploration coexists with conserved quantities that prevent global thermalization. When the pair inertia is non-compact, the same cluster state becomes chaotic, showing that the presence of chaos is controlled by the range of the interaction. The paper also claims that generic non-stationary trajectories have a central 'Janus point' from which a complexity measure and a non-equilibrium Boltzmann entropy grow in both directions of coordinate time, producing a bidirectional arrow of time in a time-reversal-invariant system.","feed_headline":"Four-fracton clusters run as billiards; ergodicity still fails","feed_subtitle":"Late-time 3-1 and 2-2 states map to a triangle and square, staying regular but never thermalizing.","key_machinery":"The machinery is the reduced-coordinate formulation together with the billiard reduction. Coordinates $q_1,q_2,q_3$ and conjugate momenta $\\pi_1,\\pi_2,\\pi_3$ remove the conserved total momentum and dipole moment; the Hamiltonian is a sum of pair-inertia terms, each a momentum-difference squared weighted by a locality function $K(x)$ that is nonzero only when two particles are nearby. In the late-time cluster states only two or three pair terms remain active, so the trajectory is confined to a plane, and the boundaries of the triangular or square allowed region act as mirrors; the paper derives explicit reflection rules under which one reduced-momentum component reverses while another is conserved, making the dynamics exactly a polygonal billiard. Three integrals then follow: the plane containing the trajectory, the conserved kinetic energy, and the finite set of possible velocity directions for a rational polygon, which together imply Liouville integrability. For the Janus-point result, the central object is the complexity $P=\\tfrac12\\sum_{i\\neq j}(p_i-p_j)^2$, whose growth on both sides of the central time is converted into a Boltzmann entropy through a hypersphere-shell phase-space count on a fixed cluster configuration.","core_discovery":"The central claim is that the late-time dynamics of four classical fractons with box (compact) pair inertia $K$ reduces to integrable billiards in reduced coordinates. A 3-1 cluster, with three fractons moving together and one frozen, is confined to an equilateral triangle in the $(q_1,q_2)$ plane; a boundary collision flips one reduced-momentum component while leaving the component parallel to the wall unchanged, which is exactly specular billiard reflection. A 2-2 cluster is likewise confined to a square, with the reflected velocity components swapped in the new coordinates. In both cases the trajectory carries three independent integrals: the conserved plane, the conserved kinetic energy, and the finite set of allowed velocity directions inherited from the rationality of the polygon. Three integrals in a six-dimensional phase space imply Liouville integrability, so the cluster is locally ergodic (it fills the triangle or square) but globally non-ergodic (it cannot explore the full phase space). With non-compact $K$, some 3-1 trajectories reach the 'large triangle' where a second fracton momentarily stops, and these states show a positive Lyapunov exponent. Separately, the paper claims every generic non-stationary trajectory has a Janus point: a time of maximal homogeneity from which complexity $P=\\frac{1}{2}\\sum_{i\\neq j}(p_i-p_j)^2$ and a Boltzmann entropy $S\\simeq k_B\\frac{N-4}{2}\\ln P+\\text{const}$ increase in both temporal directions.","pith_inferences":["Editorial extension: the reflection calculations in the appendix give an algebraic route to proving that large-triangle 3-1 trajectories form a measure-zero set for compact $K$, which would turn the billiard reduction into a theorem rather than an empirical observation.","Editorial extension: for $N>4$, the same reduction should produce billiards in rational polytopes of reduced coordinates; numerically checking whether each cluster type has zero Lyapunov exponents for compact $K$ would test whether integrability is generic for all cluster sizes.","Editorial extension: smoothing the compact box $K$ into a function with a small tail should interpolate between the integrable and chaotic regimes; measuring the first Lyapunov exponent as a function of the tail length would make the compact/non-compact transition quantitative.","Editorial extension: the Janus-point mechanism does not depend on fracton-specific details beyond cluster formation and linear momentum growth, so similar bidirectional arrows of time should appear in other non-ergodic Hamiltonian systems, including dipole-conserving lattice models in their fragmented phases."],"forward_implications":["For box or compact pair inertia, the four-fracton 3-1 and 2-2 late-time states are exactly integrable billiards, so their long-time behavior is determined by the geometry of a triangle or square rather than by generic many-body chaos.","Global ergodicity is broken in a controlled way: emergent conserved quantities confine the dynamics to a low-dimensional subspace even though the cluster motion fills that subspace densely.","Changing the interaction range from compact to non-compact turns the same cluster state from regular to chaotic, so chaos is a property of the pair inertia's support rather than of clustering itself.","Every generic non-stationary trajectory has a central time of maximal homogeneity, with complexity and the non-equilibrium Boltzmann entropy monotone in both time directions, yielding two arrows of time without invoking a past hypothesis.","Clustered steady states are stable to small phase-space perturbations for compact $K$ because nearby trajectories do not separate exponentially, suggesting the clustering found for small $N$ persists for larger $N$."],"supporting_citations":[{"why":"Defines the Machian fracton Hamiltonian and the reduced coordinates that this paper's billiard analysis builds on.","marker":"[1]"},{"why":"Established the late-time attractor and cluster structure that the four-fracton billiard reduction starts from.","marker":"[2]"},{"why":"Cited for the absence of Lyapunov chaos in polygonal billiards, the key regularity property used here.","marker":"[14]"},{"why":"Supplies the finite set of velocity directions in rational billiards, which becomes the third integral of motion.","marker":"[15]"},{"why":"Provides the Liouville-integrability criterion used to conclude that three integrals make the cluster dynamics regular.","marker":"[16]"},{"why":"Introduces the Janus-point and shape-complexity ideas that the paper adapts to fracton trajectories.","marker":"[10]"},{"why":"Supplies the non-equilibrium Boltzmann-entropy construction that the paper follows for its arrow-of-time argument.","marker":"[11]"}],"fun_headline_variants":["Four-fracton billiards: integrable yet non-ergodic","Janus point emerges in classical fracton dynamics","Local chaos, global ergodicity breaking in fractons","Fracton clusters map to triangle and square billiards","Classical fractons: local chaos, broken ergodicity, arrow of time"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The billiard reduction of the 3-1 state rests on the empirical assertion that with box (compact) $K$, no trajectory ever enters the 'large triangle' region of reduced space where a second fracton would be momentarily at rest; if such trajectories exist for some initial conditions, the small-triangle mapping and the no-chaos conclusion would not hold for those states.","fun_headline_variants_meta":{"raw":{"variants":["Four-fracton billiards: integrable yet non-ergodic","Janus point emerges in classical fracton dynamics","Local chaos, global ergodicity breaking in fractons","Fracton clusters map to triangle and square billiards","Classical fractons: local chaos, broken ergodicity, arrow of time"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001009,"raw_usage":{"total_tokens":4295,"prompt_tokens":1006,"completion_tokens":3289,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":622,"completion_tokens_details":{"reasoning_tokens":3200}},"tokens_in":622,"tokens_out":3289,"duration_ms":23532,"temperature":1.0,"reasoning_tokens":3200,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T17:11:40.141100+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run long-time simulations of the four-fracton box-$K$ Hamiltonian from many big-bang-like initial conditions and monitor the reduced coordinates; finding any 3-1 or 2-2 trajectory that enters the large-triangle or outside region (a second particle at rest) or any positive Lyapunov exponent for compact $K$ would refute the integrable-billiard claim.","supporting_citations":[{"cited_title":"This is clearly seen for systems of two fractons, which generically separate out to a fixed distance, and motion comes to a halt, reminiscent of a system with friction","cited_arxiv_id":null,"evidence_quote":"Defines the Machian fracton Hamiltonian and the reduced coordinates that this paper's billiard analysis builds on."},{"cited_title":"However, velocities of individual fractons involve the pair inertia function, and momenta of all other fractons","cited_arxiv_id":null,"evidence_quote":"Established the late-time attractor and cluster structure that the four-fracton billiard reduction starts from."},{"cited_title":"H2−2 = 1 4 ( ˙x2 + ˙y2) = const for 2–2 states","cited_arxiv_id":null,"evidence_quote":"Cited for the absence of Lyapunov chaos in polygonal billiards, the key regularity property used here."},{"cited_title":"galaxies","cited_arxiv_id":null,"evidence_quote":"Supplies the finite set of velocity directions in rational billiards, which becomes the third integral of motion."},{"cited_title":"Generalizing the approach in [1], we derive from the full Hamiltonian Eq","cited_arxiv_id":null,"evidence_quote":"Introduces the Janus-point and shape-complexity ideas that the paper adapts to fracton trajectories."},{"cited_title":"The full Hamiltonian Eq","cited_arxiv_id":null,"evidence_quote":"Supplies the non-equilibrium Boltzmann-entropy construction that the paper follows for its arrow-of-time argument."}],"review_version":1}