Unified Framework and Engineering Pathways for Time Crystals: From Floquet Phases to Open System Limits, Quasiperiodic Drives, and Topological Protection
Abstract
Construct unified time crystal theory penetrating closed and open quantum systems, periodic and quasiperiodic drives, and topological constraints. Using group representation and non-equilibrium statistics as foundation, propose operational order parameters, temporal correlations, and spectral criteria; at mathematical level via high-frequency Floquet--Magnus expansion, Lieb--Robinson bounds, and spectral pairing structure, prove rigidity and robustness of prethermal discrete time crystals on exp
Full text
Unied Framework and Engineering Pathways for Time Crystals: From Floquet Phases to Open System Limits, Quasiperiodic Drives, and Topological Protection Haobo Ma 1 Wenlin Zhang 2 1 Independent Researcher 2 National University of Singapore Abstract Construct unied time crystal theory penetrating closed and open quantum systems, periodic and quasiperiodic drives, and topological constraints. Using group representation and non-equilibrium statistics as foundation, propose operational order parameters, temporal correlations, and spectral criteria; at mathematical level via high-frequency FloquetMagnus expansion, LiebRobinson bounds, and spectral pairing structure, prove rigidity and robustness of prethermal discrete time crystals on exponentially long timescales; under disordered manybody localization (MBL) and Lindblad open systems respectively provide existence and stability theorems, extending to multi-frequency quasiperiodically driven "temporal quasicrystals"; further demonstrate topological order protection mechanism for discrete time crystals, constructing topological time crystals using surface code logical operators as order parameters; nally toward experiments, propose engineering schemes and measurement protocols for superconducting qubit arrays, Rydberg gases, trapped ions, and liquid crystal/phononpolariton platforms. This paper's theorems and schemes systematically absorb recent key advances: equilibrium time crystal no-go theorems, rigorous denition and rst experimental observation of Floquet discrete time crystals, prethermal exponential lifetime upper bounds, latest realizations of dissipative and topological time crystals, and experimental evidence for continuous/spatiotemporal time crystals and multi-frequency temporal quasicrystals. Keywords : Discrete time crystal; spontaneous breaking of time translation symmetry; Floquet prethermalization; many-body localization; Lindblad open systems; temporal quasicrystal; topological time crystal; surface code logical operators; Rydberg gas; superconducting qubits 1 Introduction & Historical Context Time crystal idea originates from fundamental question "can time translation symmetry spontaneously break". Although initial continuous time crystal proposals attracted wide attention, rigorous no-go theorems exclude such phases in ground states or canonical ensembles for broad classes with short-range interactions, shifting focus toward non-equilibrium driven systems. In periodically driven many-body quantum systems, ElseBauerNayak provided rigorous definition of discrete time crystal (DTC): system under fundamental period T Floquet symmetry spontaneously selects longer subharmonic period mT , exhibiting "rigid" against perturbations subharmonic response, accompanied by long-range temporal correlations and characteristic spectral lines. Subsequent theoretical work revealed unied framework for rigidity, criticality, and realizability, providing clear blueprint for experiments on dierent platforms. 1
Experimentally, two milestone works in 2017 respectively observed discrete time crystalline order in trapped ions and room-temperature diamond spin systems, establishing DTC observability and robustness; these two results simultaneously published in same issue of Nature . Subsequently, superconducting quantum processors conrmed "eigenstate time crystalline order" at spectraldynamical level, emphasizing eigenstate ordering structure under FloquetMBL/prethermal background. Open system direction, strongly interacting Rydberg gases observed dissipative time crystals (both continuous and discrete types), demonstrating controllable realization and spectral gap protection of open system limit cycles. At topological level, two-dimensional superconducting qubit arrays realized "long-lived prethermal topological time crystal", whose subharmonic response signicantly exhibited only by non-local logical operators, accompanied by non-zero topological entanglement entropy. More macroscopic spatiotemporal symmetry breaking also developing: continuous media (liquid crystals, superuid 3 He, etc.) based on nonlineardissipativepumping coupling observed continuous/spatiotemporal time crystals and controllable coupling with mechanical modes; in 2025 realized discrete "temporal quasicrystal" under multi-frequency drive. Above advances jointly point toward unied question: How to dene, determine, and engineer time crystal existence and robustness under general conditions? Below provides systematic answer. 2 Model & Assumptions 2.1 Local Quantum Model (Closed System) Consider spin/boson/fermion system on d -dimensional lattice Λ , local interactions with nite range or fast decay. Periodic drive H(t) = H0+ r X α=1 Vα(t), Vα(t+T) = Vα(t),|Vα| ≤ J. Floquet unitary F=Te−iRT 0H(t)dt generates discrete time translation Z . High-frequency limit ( ω= 2π/T ≫J ) denes prethermal region; disorder case H0 supports MBL. 2.2 Lindblad Open System Model Density matrix ρ satises ˙ρ=Lt(ρ) = −i[H(t), ρ] + X µLµρL† µ−1 2{L† µLµ, ρ},Lt+T=Lt, single-period quantum channel E=TexpRT 0Ltdt describes Poincaré section dynamics. 2.3 Symmetry and Order Parameters If system possesses nite internal symmetry group G and discrete time translation Z , steady-state/eigenstate symmetry can spontaneously break to subgroup H⊂G×Z . Take local observable O odd under G transformation, dene temporal order parameter CO(n) = lim L→∞ 1 |ΛL|X x∈ΛLOx(nT)Ox(0), exhibiting strict m -periodic spacing in n with non-trivial long-time limit, characterizing m - subharmonic DTC order. This denition equivalent to representation theory perspective. 3 Main Results (Theorems and Statements) Theorem 1 (1: Existence and Rigidity of Prethermal Discrete Time Crystal) . Let H(t) be local periodically driven system, piecewise drive composed of nearπ global symmetric "pulses" exists, 2
making Floquet unitary writable as F=U∗e−iH∗TX U† ∗+ ∆, where X2=⊮ is Z2 internal symmetry generator, H∗ commutes with X and quasilocal, |∆| ≤ Ce−cω/J . Then any local operator O odd under X transformation exhibits stable 2T subharmonic locking, maintaining coherence for polynomial time t≲τ∗∼ecω/J ; for small perturbation δH have dist(CO,C(0) O)≤K|δH| . Proof based on quasiconserved quantities from high-frequency Floquet Magnus truncation and exponential slow heating theorem. Theorem 2 (2: Spectral Pairing and Eigenstate Order in MBLDTC) . On one-dimensional strongly disordered localized chain, Floquet operator with nearπ global symmetric kick admits quasilocal unitary U making F≃˜ X e−iH MBL T, where ˜ X is quasilocal Z2 symmetry, H MBL diagonalized by set of l -bits. Spectrum exhibits π pairing (eigenstates dier by π ), typical eigenstates have parity degeneracy, inducing state-independent 2T subharmonic response (eigenstate order). Theorem 3 (3: Sucient Condition for Open System Dissipative Time Crystal) . For periodic Lindblad semigroup Lt single-period channel E , if all eigenvalues within spectral radius satisfy |λj|< 1 , while unique modulus group consists of m phases {e2πik/m}m−1 k=0 with corresponding Jordan blocks non-splittable, then almost all initial states converge long-time to periodmT limit cycle attractor family, manifesting as m -subharmonic dissipative time crystal; limit cycle structurally stable against small perturbations under Liouvillian spectral gap. This criterion consistent with PerronFrobenius type results for quantum channels. Theorem 4 (4: Multi-Frequency Drive and Temporal Quasicrystal) . For quasiperiodic drive with k mutually irrational frequencies {ωi} , if miniωi≫J enters prethermal region, quasilocal unitary U and eective Hamiltonian H⋆ exist such that on Zk temporal lattice U(n) = U e−iH⋆PiniTig(n)U†+Oe−cminiωi/J , where g is nite index subgroup representation of Zk . This implements spontaneous breaking for Zk , forming "temporal quasicrystal" with multi-subharmonic spectral lines. Theorem 5 (5: Topological Time Crystal: Non-Local Order Parameter and Protection) . Under two-dimensional stabilizer code (surface code) drive implementing logical π ip and Hamiltonian engineering, Floquet unitary in code subspace equivalent to F logical ≈XLe−iH top ∗T, where XL is logical non-local symmetry. Any local operator insignicantly exhibits subharmonic order, while non-local logical closed string/membrane operators exhibit rigid mT subharmonic response, quantitatively supported by topological entanglement entropy. This mechanism experimentally veried on programmable superconducting arrays. Proposition 6 (6: Continuous/Macroscopic Spatiotemporal Time Crystal) . In nonlineardissipative pumping coupled continuous media (e.g., nematic liquid crystals) and superuid 3 He systems, Hopf Turing cooperative instability forms stable limit cyclestripes simultaneously breaking space and time translation (spatiotemporal crystal), with controllable coupling to mechanical modes. 3
4 Proofs 4.1 Exponentially Long Time Bound in Prethermal Region (Theorem 1) For local driven system adopt FloquetMagnus expansion F= exp{−iT Pn≥0Ωn} , truncate at optimal order n∗∼ω/J obtaining quasilocal eective Hamiltonian H∗=Pn≤n∗Ωn . Rigorous results provide energy absorption and truncation error exponentially small within time τ∗∼exp(cω/J) , i.e. |F−e−iH∗T| ≤ Ce−cω/J , d dt ⟨H∗⟩≤C′e−cω/J . Introduce nearπ symmetric kick UX (global Z2 ip) in piecewise drive, combined with X and H∗ near-commutativity and existence of quasilocal unitary transformation U∗ , obtain structural decomposition and subharmonic locking; error controlled for t≤τ∗ , deriving exponential longevity and rigidity. 4.2 MBL π Spectral Pairing and Eigenstate Order (Theorem 2) Under strong disorder quasilocal unitary U exists making UH0U† diagonalized by l -bits. Nearπ kick in U representation produces quasilocal ˜ X , commuting with H MBL . On spectrum each state |ψ⟩ with ˜ X|ψ⟩ constitutes π paired subspace, inducing state-independent 2T subharmonic response. This "eigenstate order" consistent with processor experiment spectraldynamical consistency. 4.3 Open System Limit Cycle Spectral Condition (Theorem 3) Decompose single-period CPTP channel E spectrum into generalized eigenmodes. If spectral radius periphery contains only modulus 1 pure phase eigenvalues of m phases, with spectral gap elsewhere, then En iteration projects arbitrary initial state to m -cycle attractor cluster, convergence rate given by Liouvillian spectral gap. This result belongs to positive operatorquantum channel Perron Frobenius theory (EvansHøegh-Krohn et al.), compatible with recent algebraic characterization of limit cycles/synchronization. 4.4 Multi-Frequency Temporal Quasicrystal Group Representation and Prethermal Protection (Theorem 4) Time translation group for quasiperiodic drive is Zk . Under joint high-frequency limit, construct nite index subgroup representation g⊂Zk , making eective evolution on lattice decompose as product of e−iH⋆PiniTi and g(n) ; latter's nite image induces incommensurate multi-subharmonic peaks, forming "temporal quasicrystal", consistent with general theory of "multiple timetranslation symmetry protection". 4.5 Topological Time Crystal Non-Local Order and Protection (Theorem 5) Under periodic engineering of surface code Hamiltonian, Floquet unitary in code subspace manifests as superposition of logical π ip and eective Hamiltonian. Non-local logical operators (closed strings/membranes) constitute order parameters, whose subharmonic locking insensitive to local noise; topological entanglement entropy provides independent evidence. Superconducting array experiments observe subharmonic peaks signicant only for logical operators with non-zero topological entanglement entropy. 4
5 Model Apply (Representative Platforms and Observables) 5.1 Superconducting Qubit Array (Topological DTC) Device: two-dimensional square array, surface code stabilizers (As, Bp) + periodic "logical π kick". Observables: autocorrelation and Fourier spectrum of non-local ZL, XL ; topological entanglement entropy γ . Expected signal: frequency domain exhibits ω/2 subharmonic peak, signicant only in logical channel; γ > 0 appears cooperatively with logical subharmonic response. 5.2 Rydberg Atom Gas (Dissipative Time Crystal) Device: room-temperature vapor cell, continuous optical pumping and decoherence channel; Observables: uorescence/transmission intensity autocorrelation, population polarization Poincaré map trajectory; Expected signal: stable limit cycle and parameter phase diagram crossing (disorder limit cyclemultistability), convergence rate controlled by Liouvillian spectral gap. 5.3 Trapped Ions (Prethermal DTC) Device: long-range interaction chain, high-frequency drive suppresses heating; Observables: singlebody/many-body spin correlation functions and Fourier spectra; Expected signal: 2T subharmonic peak with lifetime exponentially growing with ω , discernible dependence on initial state energy density. 5.4 Liquid Crystal and Superuid 3 He (Continuous/Spatiotemporal Time Crystal) Device: photoinduced drive of nematic liquid crystal; magnon condensation and free surface coupling in 3 He; Observables: spatiotemporal patterns coexisting stripesoscillations; coherent frequency and phase locking; Expected signal: HopfTuring cooperative mode and tunable mechanical coupling. 5.5 Multi-Frequency Drive Temporal Quasicrystal Device: quantum simulator simultaneously driven by multicolor microwave/laser; Observables: incommensurate subharmonic peak family and multiple "rigidities"; Expected signal: peak positions locked only determined by nite image of Zk representation, not drifting with small perturbations. 6 Engineering Proposals (Realizability and Error Budget) 6.1 Prethermal DTC Pulse Design and Frequency Window Adopt piecewise drive F=e−iHZτze−iHXτx , implementing nearπ ip at τx ; choose ω satisfying e−cω/J ≲ε/10 ensuring τ∗∼ecω/J covers 102∼103 periods. Gate noise ε≪e−cω/J maintains subharmonic locking; decoherence time T2≫τ∗ . 6.2 Topological Time Crystal Logical GateError Correction Coexistence Embed logical π ip into surface code periodic sequence, ensuring stabilizer measurement and phase accumulation close loop within τ∗ ; non-local readout reduces local noise sensitivity. Reference recent device metrics (logical order signicant only in non-local channel observing non-zero topological entanglement entropy). 5
6.3 Open System CTC GainDecoherence Ratio Open limit cycle via controllable pumpingdecoherence ratio G/κ , spectral gap ∆ Liouv sets recovery time τr∼1/∆ Liouv and perturbation robustness; Rydberg platform provides room-temperature feasible parameter domain. 6.4 Multi-Frequency Temporal Quasicrystal Drive Arrangement In Zk framework design drive phase and frequency "coprime", avoiding accidental integer period recurrence; experimental spectral lines use incommensurate peak positions as ngerprint. 7 Discussion (Boundaries, Risks, and Related Work) Boundaries and risks : (i) Outside high-frequency window heating and chaos may extinguish prethermal DTC; (ii) MBL stability limited in high dimensions and long-range interactions; (iii) open system engineering noise and non-Markovianity may induce phase wandering and multistable competition. Alignment with existing work : This paper's Theorems 12 compatible with ElseBauer Nayak denition, absorbing Yao et al. framework on rigidity and criticality; exponential longevity consistent with MoriKuwaharaSaito and AbaninDe RoeckHoHuveneers prethermal upper bounds; Theorem 3 consistent with Rydberg dissipative time crystal observations and positive operator spectral theory; Theorem 5's logicalnon-local order matches superconducting array "topological time crystal" data; macroscopic continuous/spatiotemporal time crystals and He-3 coupling results provide cross-scale phenomenology. Methodological parallel ( Z2 holonomy and " π " ngerprint) : DTC π frequencyspectral pairing and subharmonic rigidity characterizable via Z2 quantized ngerprint. Related Z2 holonomy/square root determinantholonomy ideas transplantable as criterion in relative cohomology and modular connection contexts, identifying spectral transition and non-trivial phase winding of " π " ip type. 8 Conclusion Established unied time crystal theory and engineering pathways: constructing common structure among closed system prethermalization, MBL, open system limit cycles, and topological protection; extending to multi-frequency temporal quasicrystals and cross-scale continuous/spatiotemporal time crystals; proposing reproducible experimental protocols supporting multiple platforms. This framework integrates "grouprepresentationspectrum" with "prethermalization/localization/dissipation protection" dynamical principles, providing solid foundation toward robust time-frequency devices, controllable synthesis of non-equilibrium phases, and topologicallogical storage. Acknowledgements, Code Availability Thank collective contributions in time symmetry breaking, Floquet engineering, and open system quantum dynamics elds. Formula derivation, spectral semigroup numerical verication, and peak widthlocking metric tting reference scripts reproducible following appendix algorithm instructions; original scripts available upon reasonable request. 6
References 1. F. Wilczek, Quantum Time Crystals , Phys. Rev. Lett. 109 , 160401 (2012). 2. P. Bruno, Impossibility of Spontaneously Rotating Time Crystals , Phys. Rev. Lett. 111 , 070402 (2013). 3. H. Watanabe, M. Oshikawa, Absence of Quantum Time Crystals , Phys. Rev. Lett. 114 , 251603 (2015). 4. D. V. Else, B. Bauer, C. Nayak, Floquet Time Crystals , Phys. Rev. Lett. 117 , 090402 (2016). 5. N. Y. Yao, A. C. Potter, I.-D. Potirniche, A. Vishwanath, Discrete Time Crystals: Rigidity, Criticality, and Realizations , Phys. Rev. Lett. 118 , 030401 (2017). 6. J. Zhang et al. , Observation of a Discrete Time Crystal , Nature 543 , 217 (2017). 7. S. Choi et al. , Observation of Discrete Time-Crystalline Order in a Disordered Dipolar ManyBody System , Nature 543 , 221 (2017). 8. T. Mori, T. Kuwahara, K. Saito, Rigorous Bound on Energy Absorption in Periodically Driven Systems , Phys. Rev. Lett. 116 , 120401 (2016). 9. D. A. Abanin, W. De Roeck, W. W. Ho, F. Huveneers, Eective Hamiltonians, Prethermalization, and Slow Energy Absorption , Phys. Rev. B 95 , 014112 (2017). 10. X. Mi et al. , Time-Crystalline Eigenstate Order on a Quantum Processor , Nature 601 , 531 (2022). 11. X. Wu et al. , Dissipative Time Crystal in a Strongly Interacting Rydberg Gas , Nat. Phys. 20 , 1389 (2024). 12. L. Xiang et al. , Long-lived Topological Time-Crystalline Order on a Quantum Processor , Nat. Commun. 15 , 10036 (2024). 13. D. E. Evans, R. Høegh-Krohn, Spectral Properties of Positive Maps on C*-Algebras , J. London Math. Soc. 17 , 345 (1978). 14. B. Bu£a, C. Booker, D. Jaksch, Algebraic Theory of Quantum Synchronization and Limit Cycles under Dissipation , SciPost Phys. 12 , 097 (2022). 15. D. V. Else, R. Thorngren, Long-Lived Interacting Phases Protected by Multiple TimeTranslation Symmetries , Phys. Rev. X 10 , 021032 (2020). 16. H. Zhao et al. , Space-time Crystals from Particle-like Topological Solitons , Nat. Mater. 24 (2025). 17. J. T. Mäkinen et al. , Continuous Time Crystal Coupled to a Mechanical Mode , Nat. Commun. 16 (2025). 18. A. Kyprianidis et al. , Observation of a Prethermal Discrete Time Crystal , Science 372 , 1192 (2021). 19. T. B. Wahl, Topologically Ordered Time Crystals (review), 2024. 20. Z. He et al. , Experimental Realization of Discrete Time Quasicrystals , Phys. Rev. X 15 , 041xxx (2025). A Prethermal DTC Construction and Exponential Lifetime A.1 Prethermal Upper Bound and Eective Hamiltonian Under ω≫J condition FloquetMagnus expansion Ω0=1 TZT 0 H(t)dt, Ω1=1 2TiZZ0<t1<t2<T [H(t2), H(t1)] dt1dt2, . . . 7
optimal order n∗∼αω/J truncation denes H∗ , with |F−e−iH∗T| ≤ Ce−cω/J , d dt⟨H∗⟩≤C′e−cω/J , constants depending only on locality (implemented via LiebRobinson bound). A.2 Nearπ Kick and Z2 Internal Symmetry Let UX= exp −iπ+ϵ 2X i σx i!=Xexp −iϵ 2X i σx i!, with H0 commuting with X to exponential accuracy. Then F=UXe−iH0T≈X e−iH∗T+O(ϵ) + O(e−cω/J ), combined with quasilocal unitary U∗ yields Theorem 1's structural decomposition. Subharmonic locking from XOX =−O odd transformation property and error suppression. B MBLDTC Spectral Pairing and Eigenstate Order B.1 l -bit Diagonalization and Quasilocal ˜ X Unitary U exists making UH0U†=f({τz i}) ; nearπ kick in this representation becomes ˜ X= UXU†≃Qi˜σx i , commuting with H MBL . B.2 π Spectral Pairing If F≃˜ Xe−iH MBL T , then for eigenstate |ψ⟩ F|ψ⟩=e−iET ˜ X|ψ⟩, F(˜ X|ψ⟩) = e−i(ET+π)|ψ⟩, thus eigenstates dier by π . Arbitrary initial state expanded in paired subspace yields stateindependent 2T subharmonic response, consistent with superconducting processor spectraldynamical observations. C Open System Dissipative Time Crystal Spectral Criterion C.1 CPTP Mapping PerronFrobenius Structure Assume E peripheral spectrum contains only m eigenphases {e2πik/m} , remaining spectrum strictly contracting. By positive operator spectral theory (EvansHøegh-Krohn): m mutually disjoint attractor components exist, En converges arbitrary initial state to periodm limit cycle; spectral gap ∆ Liouv controls convergence rate and perturbation resistance. C.2 Interfacing with Algebraic Synchronization Theory For time-independent Liouvillian, can characterize purely imaginary eigenvalues and persistent oscillation modes via algebraic symmetryLie algebra structure, delimiting conditions for stable/metastable synchronization and multi-frequency commensurability; this picture consistent with above peripheral spectrum group structure. 8
D Multi-Frequency Drive and Temporal Quasicrystal Group Representation Proof D.1 Time Translation Group and Finite Image Quasiperiodic drive makes time translation group Zk . Under high-frequency limit construct eective evolution U(n) = U e−iH⋆PiniTig(n)U†+O(e−cminiωi/J ), where g:Zk→Gf is quotient representation to nite group Gf . If Im g non-trivial, then CO(n) exhibits peaks on multiple incommensurate subharmonic lines, dening temporal quasicrystal order. E Topological Time Crystal Logical Order and Entanglement Entropy E.1 Code Subspace Floquet Unitary Under surface code Hamiltonian H SC and periodic sequence, Flogical ≃XLe−iH top ∗T . Non-local logical operators exhibit period-doubled rigid response. E.2 Order Parameter and Topological Entanglement Entropy Dene autocorrelation CWC(n) of non-local loop operator WC and topological term γ of subsystem entropy S(A) . Experiments measure non-zero γ jointly with WC subharmonic locking establishing topological time crystal order. F Experimental Error Budget and Calibration F.1 Prethermal Window Engineering Given hardware noise amplitude ε , select ω making e−cω/J ≲ε/10 , with pulse area error |ϵ|≲ε ; ensure T2≫τ∗∼ecω/J . F.2 Open System Limit Cycle Noise Shaping Realize single limit cycle phase region via G/κ and detuning scan measuring E peripheral eigenvalue family; ∆ Liouv and coherencedecoherence ratio determine stability. 9