Pith. sign in

REVIEW 4 major objections 5 minor 2 cited by

A table-top assembly of vibrating granular disks spontaneously organizes into a rotating triangular lattice — a classical continuous spacetime crystal — and melts in three stages as the packing fraction is lowered.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-02 18:27 UTC pith:J6GFGEJK

load-bearing objection Careful experiment, real melting phenomenology, but the 'spontaneous continuous time crystal' label overreaches the current evidence. the 4 major comments →

arxiv 2603.09649 v1 pith:J6GFGEJK submitted 2026-03-10 cond-mat.soft cond-mat.stat-mech

Three-stage melting of a macroscopic continuous spacetime crystal

Guoqing Liu , Jimin Bai , Matteo Baggioli , Jie Zhang This is my paper
classification cond-mat.soft cond-mat.stat-mech
keywords continuous spacetime crystaltime-translation symmetry breakingtwo-dimensional meltingactive granular matterhexatic phasetopological defectsspontaneous symmetry breakingdirectional persistence
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper claims that a dense, circularly confined monolayer of vibrated granular disks spontaneously forms a continuous spacetime crystal: a triangular lattice that rotates rigidly with a period of roughly 4.7–5.5 hours at high packing fraction. It then reports the first experimental melting sequence for such a phase, driven by lowering the packing fraction. Spatial and temporal order melt at different packing fractions and through different mechanisms: temporal order is destroyed by loss of directional persistence as many-body interactions weaken, while spatial order is destroyed by proliferation of topological defects, passing through a hexatic phase. If correct, this is direct evidence that breaking spatial and temporal translation symmetry can be decoupled in a driven classical many-body system, and it extends spacetime-crystal phenomena into a macroscopic, table-top regime.

Core claim

The central discovery is that active granular disks confined in a circle self-organize at high packing fraction into a state in which particles sit on a triangular lattice while the whole lattice rotates coherently for nearly a day. The rotation period is about five hours, six orders of magnitude longer than the vertical drive; the motion shows a gapless phase-fluctuation mode, resists strong injected noise, and appears with random onset times and phases across nominally identical small systems. Lowering the packing fraction melts the time-crystalline order first (around φ ≈ 0.709), leaving a spatially hexatic phase in a coexistence regime, then melts the remaining spatial order (around φ ≈

What carries the argument

The carrying object is the spacetime-crystalline phase itself: a rigid-body-rotating two-dimensional triangular lattice characterized by two measured order parameters — the spatial crystalline fraction (Bragg-peak weight in the static structure factor) and the time-crystalline fraction (spectral weight of the emergent oscillation peak). Temporal melting is tracked by directional persistence and non-affine displacement; spatial melting is tracked by hexatic correlations and Voronoi-based topological defects such as dislocations, disclinations, and defect clusters. The rotation is explained as emergent flocking under circular confinement, where accumulated collisions create effective mutual at

Load-bearing premise

The load-bearing premise is that the observed rotation is a spontaneous breaking of continuous time-translation symmetry rather than an externally selected limit cycle: the paper's evidence is random onset times and phases across seven 15-cm replicas, yet only five of seven rotated within the 15-hour window, their periods spread from 1.68 to 4.95 hours, and at the highest packing fraction all rotated counterclockwise, which the supplementary text attributes to minor experimen

What would settle it

Measure residual tilt and local anisotropy of the plate while running the seven small replicas with their positions randomly shuffled; if the rotation direction and phase of each replica track the plate's measured asymmetry rather than varying randomly across shuffled runs, the 'spontaneous' label is falsified.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • A driven, dissipative classical system can sustain a spacetime crystal for macroscopic times, with a collective rotation period near 10^4 seconds set by many-body interactions rather than by the 10^-2 second drive.
  • Spatial and temporal order behave as independent axes: at intermediate packing fractions the system is time-disordered but spatially hexatic, showing that temporal order can melt before spatial order.
  • The melting route is three-stage — spacetime crystal, time-coexistence (hexatic), space-coexistence, fluid — with critical packing fractions 0.734, 0.709, and 0.687.
  • Temporal rigidity has a many-body origin: the time-crystalline phase remains mostly intact under maximum injected noise, consistent with a spontaneously broken symmetry rather than a fragile driven oscillation.
  • A gapless Goldstone mode accompanies the broken time symmetry, with phase-fluctuation structure factor scaling roughly as 1/q^2 and a linearly dispersing, weakly damped mode, giving the phase the rigidity expected of a symmetry-broken state.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • An extension the authors do not pursue: reshaping the confining boundary from circular to flower-like should suppress the rotating crystal and replace the three-stage sequence with ordinary two-dimensional melting, which would directly test the flocking mechanism.
  • If the residual-bias explanation is right, better leveling and randomized manufacturing asymmetries should increase the fraction of replicas that begin rotating within a fixed window and randomize the high-density rotation direction; the distribution of onset waiting times could be measured and compared with nucleation-like statistics.
  • The separate order parameters introduced here could be exported to other driven many-body systems, including simulations of active Brownian disks at low activity-to-diffusion ratio, to look for the same three-stage decoupling of spatial and temporal melting.
  • The three-stage route suggests a two-parameter phase diagram in which packing fraction and noise/activity are varied independently, potentially revealing a window where temporal order is lost while a true spatial crystal, rather than a hexatic, survives.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. This manuscript reports experiments on vertically vibrated granular disks with ratchet legs confined in a circular boundary. At high packing fraction φ=0.835, the disks form a triangular lattice that undergoes coherent rigid-body rotation with a period of several hours (~5 h), which the authors interpret as a classical continuous spacetime crystal spontaneously breaking both spatial and continuous temporal translational symmetry. Decreasing φ, they identify four phases: a spacetime crystal (φ>0.734), a 'T-coexistence' region (0.709–0.734), an 'S-coexistence' region (0.687–0.709), and a fluid. They claim that spatial and temporal order melt separately and through distinct mechanisms: spatial order via proliferation of topological defects, and temporal order via decay of directional persistence caused by progressive weakening of many-body interactions.

Significance. If substantiated, this would be a striking table-top demonstration of spontaneous breaking of continuous time-translation symmetry in a classical dissipative system and the first experimental account of spacetime-crystal melting, with implications for out-of-equilibrium phase transitions. Strengths include macroscopic visualization, persistence for almost a day, robustness to injected acoustic noise, multiple nominally identical replicas, and a phase diagram built from independent spatial and temporal probes. The manuscript also contains honest statements about the role of experimental imperfections in fixing chiral order. However, the core identification of a spontaneous time crystal and the claimed distinct melting mechanisms rest on evidence that is not yet fully secured; the concerns below are load-bearing for the central claims.

major comments (4)
  1. [Main text, 'Experimental observation'; SM1 and SM3] The spontaneous breaking of continuous time-translation symmetry is not established because the rotation chirality is attributed to an external bias. At φ=0.835, RS, MRS, and CS particles all rotate counterclockwise, and this is explicitly attributed to 'minor imperfections in the experimental setup that are amplified by many-body interactions' (SM3). Under such a bias, random onset times and phases among only 5 of 7 replicas (SM1, periods ranging 1.68–4.95 h) can arise from stochastic nucleation in a biased potential rather than from Mexican-hat phase selection. No quantitative test (e.g., Rayleigh test for uniform phase distribution) is reported, and two replicas never rotate within the window. Please provide a statistical test using all seven replicas with appropriate censoring, measure or bound the bias, or perform experiments with reversed or removed bias to show that the phase is s
  2. [SM11, Eqs. (S13)–(S14)] The claimed Goldstone mode is not a clean discriminator. φ_i(t) is defined as the residual after subtracting each particle's best-fit linear rotation, so S_φ(q) ~ 1/q^2 is expected for ordinary 2D displacement fluctuations and does not uniquely indicate a temporal Goldstone mode. The dynamic structure factor S_φ(q,ω) is constructed from the equal-time spatial covariance of these residuals; a linearly dispersing mode with γ→0 at q→0 could be an artifact of detrending and the normal-mode procedure. Please compare against a control (e.g., a non-rotating but spatially ordered configuration, or a synthetic model with prescribed phase noise) and show that the gapless mode is genuinely a phase mode rather than a phonon artifact.
  3. [Figure 3D; Abstract; Figure 5] The claim that spatial and temporal order melt through 'distinct mechanisms' is not directly evidenced. Temporal melting is characterized by the decay of directional persistence, which is the very quantity plotted in Fig. 3D; the statement that this is 'caused by the progressive weakening of many-body interactions' is an interpretation, as interaction strength is never measured. Please provide a direct measurement of interaction strength/collision rate as a function of φ, or a mechanistic model with a measurable interaction parameter that yields the observed decay of directional persistence. As written, the mechanism for temporal melting restates the order parameter.
  4. [Methods, Eqs. (19)–(24) and thresholds] The phase boundaries φ1, φ2, and φ3 depend on several free parameters: the angular window l_W=10° in Eq. (19), the hexatic threshold |ψ6|>0.64, the MSD criterion MSD/D^2=2.5, and the non-affine threshold log10(Dmin^2/D^2)<0. Please include a sensitivity analysis showing that the three-stage melting scenario and the reported critical packing fractions are robust under reasonable variations of these parameters. Without it, the separation of spatial and temporal melting may be an artifact of the operational definitions rather than a physical decoupling.
minor comments (5)
  1. [Section 'Melting of 2D spatial order'] Typo: 'we turn into the the melting' should be 'we turn to the melting.' Also, there are spacing issues in the typeset text with 'V oronoi' in the Methods and SKM sections.
  2. [Figure 1, panels E/F] The main text reports T≈4.70 h and a range 4.7–5.5 h, while Fig. 1F gives f=5.5×10^-5 Hz, corresponding to a period of 5.05 h. Please ensure consistency and report the uncertainty of the period estimate.
  3. [Methods, Eq. (17)] The conversion of mode frequencies from inverse length to physical frequency via V0 is dimensionally motivated but not justified. Please explain why V0 is the appropriate velocity scale for the phonons after subtracting the global rotation, and how the result depends on the choice of velocity measure (e.g., mean speed vs. RMS speed).
  4. [Main text, Eq. (44) and definition of G(t)] The definition G(t)=⟨˜y(t)˜y(0)⟩−⟨˜y(t)⟩⟨˜y(0)⟩ is unusual for a signal already normalized to [−1,1] by min-max scaling; if the mean is zero by construction, the subtracted term may be negligible. Please clarify the purpose of the subtraction and whether the same normalization is used in the envelope fitting.
  5. [References] A number of references are to arXiv preprints or in-press articles with 2026 dates (e.g., refs. [34], [51]). Please verify final publication status and update if possible.

Circularity Check

0 steps flagged

No significant circularity: the melting phase diagram is built from independent operational order parameters, and the self-citations are not load-bearing.

full rationale

The central phase diagram rests on two independent operational observables: the spatial crystalline fraction (normalized Bragg-peak angular weight, Methods Eqs. 18-21) and the time-crystalline fraction (normalized spectral power around the dominant frequency peak, Methods Eqs. 22-24). These are computed from different data (instantaneous positions vs. particle-trajectory Fourier spectra) and are not algebraically linked, so the observation that they drop at different packing fractions is an empirical result rather than a construction. The directional-persistence and non-affine-motion analyses are additional independent probes; the statement that "temporal order is lost through the decay of directional persistence" is a phenomenological description of the same loss of coherent tangential motion, but no fitted parameter is renamed as a prediction and no equation reduces one order parameter to another. The SM11 Goldstone-mode analysis defines the phase residual via Eq. S13, φ_i(t)=θ_total,i(t)-θ_fit,i(t), and the observed 1/q² phase spectrum is therefore partly a consistency check on the same phase variable that defines the time-crystalline order; however, the paper does not rely on the Goldstone mode as the sole evidence for spontaneous time-translation symmetry breaking, and the random-phase and noise-robustness tests are separate, not derived from that analysis. The self-citations (refs. 44, 50, 51, 60) support ancillary mechanistic interpretations such as the flocking origin of rotation and the kinetic-temperature scaling, which are externally published experiments/simulations rather than unverified uniqueness claims; they are not load-bearing for the three-stage melting scenario. Concerns about external bias (consistent counterclockwise chirality, only 5/7 small replicas rotating within 15 h, period scatter, and no quantitative phase-uniformity test) are validity threats to the spontaneous-symmetry-breaking interpretation, not circularity in the paper's derivation chain. Overall, no step reduces by construction to its own inputs.

Axiom & Free-Parameter Ledger

4 free parameters · 6 axioms · 0 invented entities

No new particles, forces, or mediators are introduced. The central claim rests on operational definitions of order and on interpreting the rotation as spontaneous; both involve arbitrary thresholds and interpretive leaps. The main free parameters are analysis thresholds whose values are not justified by a convergence study.

free parameters (4)
  • angular window l_W for spatial crystalline fraction = 10 degrees
    Width of the angular window around each Bragg peak used to define the spatial crystalline fraction (Methods). The value is chosen by hand and no convergence study is reported.
  • hexatic threshold |ψ6| > 0.64 = 0.64
    Threshold for classifying a particle as structurally ordered in the structural percolation analysis (Methods). Arbitrary choice affecting the percolation curves.
  • MSD criterion MSD/D^2 = 2.5 = 2.5
    Defines the characteristic lag time t* used in the directional persistence measure (Methods). Arbitrary choice affecting the persistence values.
  • non-affine threshold log10(Dmin^2/D^2) < 0 = 0
    Threshold for identifying dynamically ordered particles in the dynamical percolation analysis (Methods). Arbitrary choice affecting cluster statistics.
axioms (6)
  • domain assumption The vibrated granular system reaches a nonequilibrium steady state with energy injection balanced by dissipation and no slow aging over the ~20 h observation window.
    All phase assignments assume each packing fraction represents a stationary phase; no hysteresis or equilibration tests are reported.
  • ad hoc to paper The 100 Hz vertical vibration does not impose in-plane periodic forcing, so the ~5 h rotation is a spontaneous low-frequency symmetry breaking rather than a response to the drive.
    Asserted in the main text; the drive frequency is six orders of magnitude faster, but the shaker still provides a lab-frame reference and energy input.
  • ad hoc to paper Random onset time and phase across seven replicas establish spontaneous breaking of continuous time-translation symmetry.
    Evidence is 5 of 7 systems rotating within 15 h with periods 1.68-4.95 h; the randomness could also reflect extrinsic variability rather than a true Goldstone phase.
  • domain assumption Standard 2D melting phenomenology (KTHNY) applies qualitatively to this out-of-equilibrium active system.
    Used throughout to interpret hexatic order and defect proliferation; the authors note the system is beyond equilibrium KTHNY but compare to active Brownian disk simulations.
  • domain assumption Normal-mode reconstruction from the displacement/phase covariance matrix, rescaled by a characteristic velocity or rate, yields physical phonon and Goldstone spectra.
    Methods and SM11; the absolute frequency scale depends on the arbitrary rescaling factor V0 or Gamma.
  • domain assumption The time-crystalline fraction, defined as normalized spectral power in the dominant oscillation peak, is a valid order parameter for temporal crystalline order.
    Operational definition in Methods; the phase diagram is built on this metric.

pith-pipeline@v1.3.0-alltime-deepseek · 25908 in / 17564 out tokens · 162486 ms · 2026-08-02T18:27:57.195190+00:00 · methodology

0 comments
read the original abstract

A spacetime crystal is a phase of matter that spontaneously develops periodic order in both space and time. Spacetime crystals have been experimentally observed in microscopic quantum many-body systems and, very recently, in a mesoscopic nematic liquid crystal. However, the melting process of a spacetime crystal and its underlying physical mechanisms have not yet been experimentally reported. Here, we present a direct observation of a classical continuous spacetime crystal melting in a table-top experiment with macroscopic active granular disks in 2+1 spacetime dimensions. The spacetime crystal is characterized by the spontaneous formation of a coherent, rigid-body rotation of a 2D triangular lattice that persists for almost a day and remains remarkably robust to noise. By tuning the disk packing fraction, we observe a complex three-stage melting process involving a spatially hexatic phase and multiple coexistence regions. Importantly, we show that spatial and temporal crystalline orders melt separately through distinct mechanisms: spatial order is destroyed by the proliferation of topological defects, while temporal order is lost through the decay of directional persistence caused by the progressive weakening of many-body interactions. Our results demonstrate that the spontaneous breaking of spatial and temporal translational symmetries can be decoupled, leading to the emergence of exotic out-of-equilibrium classical phases of matter.

Figures

Figures reproduced from arXiv: 2603.09649 by Guoqing Liu, Jie Zhang, Jimin Bai, Matteo Baggioli.

Figure 1
Figure 1. Figure 1: Macroscopic continuous spacetime crystal. A. Experimental setup: Particles are placed on an aluminum-alloy plate and confined within a circular boundary. Vertical vibrations are applied normal to the plate using an electromagnetic shaker. Inset: schematic of the ratcheted particles. B. Single-particle translational motion: Particles initially placed near the center exhibit random, disordered trajectories, … view at source ↗
Figure 2
Figure 2. Figure 2: Three-stage melting of spacetime order. A. Phase diagram showing, from right to left, the spacetime crystal, time-coexistence, space-coexistence, and fluid phases. The temporal and spatial crystalline fractions decrease independently across the coexistence regimes. The critical packing fractions separating the phases are ϕ1 = 0.687, ϕ2 = 0.709, and ϕ3 = 0.734. One representative packing fraction is selecte… view at source ↗
Figure 3
Figure 3. Figure 3: Dynamical melting of time-crystalline order. A. Dynamical coexistence. The dynamical parameter log10 D 2 min/D 2  quantifies local non-affine displacements. In the spacetime crystalline phase, this quantity remains small, indicating highly ordered and coherent dynamics. In the fluid phase, no ordered domains are present. In the time-coexistence (T-coexistence) region, dynamically ordered and disordered do… view at source ↗
Figure 4
Figure 4. Figure 4: Structural correlations and topological defects. A. Spatial distribution of the local hexatic order parameter |ψ6| in the S-coexistence phase at ϕ = 0.662, where fragmented ordered domains coexist with a disordered fluid. Here, |ψ6| represents the modulus of ψ6. B. Distribution of ψ6 in the complex plane for the configuration shown in panel A. The central peak corresponds to disordered fluid particles, whi… view at source ↗
Figure 5
Figure 5. Figure 5: Spatiotemporal phase diagram of spatial and temporal order. A summary of the different phases in terms of their spatial and temporal crystalline order is shown. Remarkably, spatial and temporal symmetries are either both preserved or both broken in the low-packing-fraction fluid phase and in the high-packing-fraction spacetime crystalline state. In contrast, in the two intermediate coexistence regions, spa… view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Phase Time Crystals and Pairing in Binary Active Chiral Systems

    cond-mat.stat-mech 2026-06 unverdicted novelty 7.0

    Binary active chiral particle systems under out-of-phase circular drives form phase time crystals including bound-pair triangular lattices along with stripes, glasses, and phase-separated states.

  2. Emergence of Time Semicrystals in Holographic Driven-Dissipative Systems

    hep-th 2026-04 unverdicted novelty 5.0

    A holographic driven-dissipative model exhibits a time semicrystal phase with persistent discrete subharmonic peaks on a continuous spectrum, critical scaling at the transition from discrete time crystals, and log-per...

Reference graph

Works this paper leans on

43 extracted references · 13 canonical work pages · cited by 2 Pith papers

  1. [2]

    M., Lubensky, T

    Chaikin, P. M., Lubensky, T. C. & Witten, T. A.Principles of condensed matter physics, vol. 10 (Cambridge university press Cambridge, 1995). 3.Frenkel, D. Order through entropy.Nat. materials14, 9–12 (2015)

    1995

  2. [4]

    Wood, W. W. & Jacobson, J. D. Preliminary results from a recalculation of the monte carlo equation of state of hard spheres.The J. Chem. Phys.27, 1207–1208 (1957). 5.Alder, B. J., Wainwright, T. E.et al.Phase transition for a hard sphere system.The J. chemical physics27, 1208 (1957). 6.Strandburg, K. J. Two-dimensional melting.Rev. Mod. Phys.60, 161–207, ...

    work page doi:10.1103/revmodphys.60.161 1957

  3. [7]

    Kosterlitz, J. M. & Thouless, D. J. Ordering, metastability and phase transitions in two-dimensional systems.J. Phys. C: Solid State Phys.6, 1181, DOI: 10.1088/0022-3719/6/7/010 (1973)

    doi:10.1088/0022-3719/6/7/010 1973

  4. [8]

    Nelson, D. R. & Halperin, B. I. Dislocation-mediated melting in two dimensions.Phys. Rev. B19, 2457–2484, DOI: 10.1103/PhysRevB.19.2457 (1979)

    doi:10.1103/physrevb.19.2457 1979

  5. [9]

    Young, A. P. Melting and the vector coulomb gas in two dimensions.Phys. Rev. B19, 1855–1866, DOI: 10.1103/PhysRevB. 19.1855 (1979). 10.Nelson, D. R.Defects and geometry in condensed matter physics(Cambridge University Press, 2002). 11.Wilczek, F. Quantum time crystals.Phys. Rev. Lett.109, 160401, DOI: 10.1103/PhysRevLett.109.160401 (2012)

    doi:10.1103/physrevb 1979

  6. [12]

    Time crystals: Can diamagnetic currents drive a charge density wave into rotation?Europhys

    Nozières, P. Time crystals: Can diamagnetic currents drive a charge density wave into rotation?Europhys. Lett.103, 57008 (2013)

    2013

  7. [13]

    Impossibility of spontaneously rotating time crystals: A no-go theorem.Phys

    Bruno, P. Impossibility of spontaneously rotating time crystals: A no-go theorem.Phys. Rev. Lett.111, 070402, DOI: 10.1103/PhysRevLett.111.070402 (2013)

    doi:10.1103/physrevlett.111.070402 2013

  8. [14]

    & Oshikawa, M

    Watanabe, H. & Oshikawa, M. Absence of quantum time crystals.Phys. Rev. Lett.114, 251603, DOI: 10.1103/PhysRevLett. 114.251603 (2015)

    doi:10.1103/physrevlett 2015

  9. [15]

    P.et al.Colloquium: Quantum and classical discrete time crystals.Rev

    Zaletel, M. P.et al.Colloquium: Quantum and classical discrete time crystals.Rev. Mod. Phys.95, 031001, DOI: 10.1103/RevModPhys.95.031001 (2023)

    doi:10.1103/revmodphys.95.031001 2023

  10. [16]

    & Zakrzewski, J

    Sacha, K. & Zakrzewski, J. Time crystals: a review.Reports on Prog. Phys.81, 016401, DOI: 10.1088/1361-6633/aa8b38 (2017). 17.Khemani, V ., Moessner, R. & Sondhi, S. A brief history of time crystals.arXiv preprint arXiv:1910.10745(2019). 18.Zhang, J.et al.Observation of a discrete time crystal.Nature543, 217–220 (2017). 19.Kongkhambut, P.et al.Observation...

    Pith/arXiv arXiv 2017

  11. [20]

    Choi, S.et al.Observation of discrete time-crystalline order in a disordered dipolar many-body system.Nature543, 221–225 (2017)

    2017

  12. [21]

    Rovny, J., Blum, R. L. & Barrett, S. E. Observation of discrete-time-crystal signatures in an ordered dipolar many-body system.Phys. Rev. Lett.120, 180603, DOI: 10.1103/PhysRevLett.120.180603 (2018)

    work page doi:10.1103/physrevlett.120.180603 2018

  13. [22]

    Smits, J., Liao, L., Stoof, H. T. C. & van der Straten, P. Observation of a space-time crystal in a superfluid quantum gas. Phys. Rev. Lett.121, 185301, DOI: 10.1103/PhysRevLett.121.185301 (2018). 12/37

    doi:10.1103/physrevlett.121.185301 2018

  14. [23]

    Autti, S., Eltsov, V . B. & V olovik, G. E. Observation of a time quasicrystal and its transition to a superfluid time crystal. Phys. Rev. Lett.120, 215301, DOI: 10.1103/PhysRevLett.120.215301 (2018)

    work page doi:10.1103/physrevlett.120.215301 2018

  15. [24]

    Phys.22, 085001, DOI: 10.1088/1367-2630/ab9fbe (2020)

    O’Sullivan, J.et al.Signatures of discrete time crystalline order in dissipative spin ensembles.New J. Phys.22, 085001, DOI: 10.1088/1367-2630/ab9fbe (2020). 25.Kyprianidis, A.et al.Observation of a prethermal discrete time crystal.Science372, 1192–1196 (2021)

    work page doi:10.1088/1367-2630/ab9fbe 2020

  16. [26]

    https://www.science.org/doi/pdf/10.1126/science.abk0603

    Randall, J.et al.Many-body–localized discrete time crystal with a programmable spin-based quantum simulator.Science 374, 1474–1478, DOI: 10.1126/science.abk0603 (2021). https://www.science.org/doi/pdf/10.1126/science.abk0603. 27.Mi, X.et al.Time-crystalline eigenstate order on a quantum processor.Nature601, 531–536 (2022)

    doi:10.1126/science.abk0603 2021

  17. [28]

    & Zhang, X

    Chen, Y .-H. & Zhang, X. Realization of an inherent time crystal in a dissipative many-body system.Nat. Commun.14, 6161 (2023). 29.Wu, X.et al.Dissipative time crystal in a strongly interacting rydberg gas.Nat. Phys.20, 1389–1394 (2024). 30.Greilich, A.et al.Robust continuous time crystal in an electron–nuclear spin system.Nat. Phys.20, 631–636 (2024)

    2023

  18. [31]

    Science384, 995–1000, DOI: 10.1126/science.adn7087 (2024)

    Carraro-Haddad, I.et al.Solid-state continuous time crystal in a polariton condensate with a built-in mechanical clock. Science384, 995–1000, DOI: 10.1126/science.adn7087 (2024). https://www.science.org/doi/pdf/10.1126/science.adn7087

    doi:10.1126/science.adn7087 2024

  19. [32]

    & Smalyukh, I

    Zhao, H. & Smalyukh, I. I. Space-time crystals from particle-like topological solitons.Nat. Mater.24, 1802–1811 (2025)

    2025

  20. [33]

    & Smalyukh, I

    Zhao, H., Zhang, R. & Smalyukh, I. I. Emergent discrete space-time crystal of majorana-like quasiparticles in chiral liquid crystals.arXiv preprint arXiv:2507.16977(2025)

    Pith/arXiv arXiv 2025

  21. [34]

    C., Elliott, L

    Morrell, M. C., Elliott, L. & Grier, D. G. Nonreciprocal wave-mediated interactions power a classical time crystal.Phys. Rev. Lett.136, 057201, DOI: 10.1103/zjzk-t81n (2026). 35.Shapere, A. & Wilczek, F. Classical time crystals.Phys. Rev. Lett.109, 160402, DOI: 10.1103/PhysRevLett.109.160402 (2012)

    work page doi:10.1103/zjzk-t81n 2026

  22. [36]

    Dai, J., Niemi, A. J. & Peng, X. Classical hamiltonian time crystals–general theory and simple examples.New J. Phys.22, 085006, DOI: 10.1088/1367-2630/aba8d3 (2020)

    work page doi:10.1088/1367-2630/aba8d3 2020

  23. [37]

    L., Oscity, M., Eichler, A., Zilberberg, O

    Heugel, T. L., Oscity, M., Eichler, A., Zilberberg, O. & Chitra, R. Classical many-body time crystals.Phys. Rev. Lett.123, 124301, DOI: 10.1103/PhysRevLett.123.124301 (2019)

    work page doi:10.1103/physrevlett.123.124301 2019

  24. [38]

    Y ., Nayak, C., Balents, L

    Yao, N. Y ., Nayak, C., Balents, L. & Zaletel, M. P. Classical discrete time crystals.Nat. Phys.16, 438–447, DOI: 10.1038/s41567-019-0782-3 (2020)

    work page doi:10.1038/s41567-019-0782-3 2020

  25. [39]

    Liu, T., Ou, J.-Y ., MacDonald, K. F. & Zheludev, N. I. Photonic metamaterial analogue of a continuous time crystal.Nat. Phys.19, 986–991 (2023)

    2023

  26. [40]

    Li, T.et al.Space-time crystals of trapped ions.Phys. Rev. Lett.109, 163001, DOI: 10.1103/PhysRevLett.109.163001 (2012)

    work page doi:10.1103/physrevlett.109.163001 2012

  27. [41]

    Träger, N.et al.Real-space observation of magnon interaction with driven space-time crystals.Phys. Rev. Lett.126, 057201, DOI: 10.1103/PhysRevLett.126.057201 (2021)

    work page doi:10.1103/physrevlett.126.057201 2021

  28. [42]

    Xu, S. & Wu, C. Space-time crystal and space-time group.Phys. Rev. Lett.120, 096401, DOI: 10.1103/PhysRevLett.120. 096401 (2018)

    doi:10.1103/physrevlett.120 2018

  29. [43]

    & Cai, Z

    Yue, M., Yang, X. & Cai, Z. Thermal melting of discrete time crystals: A dynamical phase transition induced by thermal fluctuations.Phys. Rev. B105, L100303, DOI: 10.1103/PhysRevB.105.L100303 (2022)

    work page doi:10.1103/physrevb.105.l100303 2022

  30. [44]

    & Zhang, H

    Yang, P., Baggioli, M., Cai, Z., Tian, Y . & Zhang, H. Holographic dissipative spacetime supersolids.Phys. Rev. Lett.131, 221601, DOI: 10.1103/PhysRevLett.131.221601 (2023). 45.Tan, T. H.et al.Odd dynamics of living chiral crystals.Nature607, 287–293 (2022). 46.Bililign, E. S.et al.Motile dislocations knead odd crystals into whorls.Nat. Phys.18, 212–218 (...

    work page doi:10.1103/physrevlett.131.221601 2023

  31. [48]

    & Vitelli, V

    Fruchart, M., Scheibner, C. & Vitelli, V . Odd viscosity and odd elasticity.Annu. Rev. Condens. Matter Phys.14, 471–510 (2023)

    2023

  32. [49]

    t., Wittkowski, R

    Huang, Z.-F., Vrugt, M. t., Wittkowski, R. & Löwen, H. Anomalous grain dynamics and grain locomotion of odd crystals. Proc. Natl. Acad. Sci.122, e2511350122 (2025). 50.Chen, Y . & Zhang, J. Anomalous flocking in nonpolar granular brownian vibrators.Nat. Commun.15, 6032 (2024). 13/37

    2025

  33. [51]

    & Zhang, J

    Zheng, Z., Jiang, C., Chen, Y ., Baggioli, M. & Zhang, J. Topological signatures of collective dynamics and turbulent-like energy cascades in apolar active granular matter.Proc. Natl. Acad. Sci.123, e2510873123, DOI: 10.1073/pnas.2510873123 (2026). https://www.pnas.org/doi/pdf/10.1073/pnas.2510873123

    doi:10.1073/pnas.2510873123 2026

  34. [52]

    & Ramaswamy, S

    Toner, J., Tu, Y . & Ramaswamy, S. Hydrodynamics and phases of flocks.Annals Phys.318, 170–244, DOI: https: //doi.org/10.1016/j.aop.2005.04.011 (2005). Special Issue

    doi:10.1016/j.aop.2005.04.011 2005

  35. [53]

    Falk, M. L. & Langer, J. S. Dynamics of viscoplastic deformation in amorphous solids.Phys. Rev. E57, 7192–7205, DOI: 10.1103/PhysRevE.57.7192 (1998)

    doi:10.1103/physreve.57.7192 1998

  36. [54]

    Engel, M.et al.Hard-disk equation of state: First-order liquid-hexatic transition in two dimensions with three simulation methods.Phys. Rev. E87, 042134, DOI: 10.1103/PhysRevE.87.042134 (2013)

    doi:10.1103/physreve.87.042134 2013

  37. [55]

    Kapfer, S. C. & Krauth, W. Two-dimensional melting: From liquid-hexatic coexistence to continuous transitions.Phys. Rev. Lett.114, 035702, DOI: 10.1103/PhysRevLett.114.035702 (2015)

    doi:10.1103/physrevlett.114.035702 2015

  38. [56]

    L., Abbott, J

    Thorneywork, A. L., Abbott, J. L., Aarts, D. G. A. L. & Dullens, R. P. A. Two-dimensional melting of colloidal hard spheres.Phys. Rev. Lett.118, 158001, DOI: 10.1103/PhysRevLett.118.158001 (2017)

    work page doi:10.1103/physrevlett.118.158001 2017

  39. [57]

    Bernard, E. P. & Krauth, W. Two-step melting in two dimensions: First-order liquid-hexatic transition.Phys. Rev. Lett. 107, 155704, DOI: 10.1103/PhysRevLett.107.155704 (2011)

    doi:10.1103/physrevlett.107.155704 2011

  40. [58]

    Digregorio, P.et al.Full phase diagram of active brownian disks: From melting to motility-induced phase separation.Phys. Rev. Lett.121, 098003, DOI: 10.1103/PhysRevLett.121.098003 (2018)

    doi:10.1103/physrevlett.121.098003 2018

  41. [59]

    F., Gonnella, G

    Digregorio, P., Levis, D., Cugliandolo, L. F., Gonnella, G. & Pagonabarraga, I. Unified analysis of topological defects in 2d systems of active and passive disks.Soft Matter18, 566–591 (2022)

    2022

  42. [60]

    & Zhang, J

    Jiang, C., Zheng, Z., Chen, Y ., Baggioli, M. & Zhang, J. Experimental observation of gapped shear waves and liquid-like to gas-like dynamical crossover in active granular matter.Commun. Phys.8, 82 (2025)

    2025

  43. [61]

    liquid-gas

    Brazhkin, V . V .et al.“liquid-gas” transition in the supercritical region: Fundamental changes in the particle dynamics. Phys. Rev. Lett.111, 145901, DOI: 10.1103/PhysRevLett.111.145901 (2013). 14/37 Figure 1.Macroscopic continuous spacetime crystal.A.Experimental setup: Particles are placed on an aluminum-alloy plate and confined within a circular bound...

    doi:10.1103/physrevlett.111.145901 2013