Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2605.01969

Ergodic and Discrete Time Crystal Phases in Periodically Kicked Many-Body Quantum Systems: An Analytical Study

Vijay Kumar, Dibyendu Roy

cond-mat.stat-mechnlin.CDquant-phSubmitted May 3, 2026• Announced May 5, 2026• 0 citations
arXivHTMLPDF

Abstract

We analytically study the time evolution of the expectation values of observables in periodically kicked many-body quantum systems. Starting from an initial state, we compute both the transient and the long-time properties of the observables. Our derivation explains the criteria and the mechanism that lead to the infinite-temperature statistical average of observables at long times, irrespective of the initial state. When the criteria are violated, the observables oscillate with time. These oscillations are subharmonic and robust to small perturbations, suggesting the emergence of a discrete time crystal phase. We demonstrate these features explicitly in periodically kicked nonintegrable spin chains. For a spin chain with two kicks per cycle, we show that the kicked chain can exhibit an ergodic or a discrete-time crystal phase for the same kicking strengths, depending on the initial state preparation. We complement our time-evolution study of observables with the spectral form factor of these kicked models.

Related papers

ERGODIC PROPERTIES OF QUANTUM SPIN CHAINS: KICKED TRANSVERSE ISING MODEL

Tomai Prosen

The problem of calculation of time-correlation functions and relaxation phenomena in interacting quantum many-body systems or quantized fields is put into the framework of ergodic theory and dynamical systems. In this spirit we discuss spectral analysis of the adjoint propagator in a suitable Hilbert space of quantum observables in Heisenberg picture as an alternative approach to characterize many-body dynamics. We illustrate our approach by working out the time-autocorrelation functions at infinite temperature of a large class of observabIes for a spin s = 1/2 Ising chain in a periodically kicking transverse magnetic field, and show that they decay to their asymptotic values as ex: t3/ 2 General understanding of time-correlation functions and relaxation phenomena of interacting quantum many-body systems is a challenging problem with only very partial existing solutions. Situation is much better if one considers the motion of a one-body or few-body systems where the classical ergodic theory (see e.g. [1]) describes a variety of different types of motion, ranging from quasi-periodic integrable motion characterized by a discrete spectrum of the (Liouville) evolution operator, through KAM quasi-integrability, to ergodicity, mixing and chaos characterized by a continuous evolution spectrum. Mixing, which is equivalent of saying that correlation functions of an arbitrary pair of observables (A,B) decay in time, limt-HXJ ((A(t)B) (A)(B)) = 0, is the necessary dynamical property needed to justify the relaxation to equilibrium (micro/grand/canonical) state and the laws of statistical mechanics, such as the fluctuation-dissipation theorem and transport laws. Mixing in quantum many-body system has been investigated by numerical experiment only quite recently [2]. In another numerical experiment [3], a suggestive possibility of dynamical phase transition from non-ergodic/non-mixing to ergodic/mixing dynamics in a non-integrable periodically kicked chain of interacting spinless fermions (kicked t-V model) has been proposed. General absence of mixing in completely integrable quantum many-body systems (such as those solvable by Bethe ansatz), being the consequence of an infinite sequence of conservation 317 J. Bonca et al. (eds.), Open Problems in Strongly Correlated Electron Systems. 317-322. © 2001 Kluwer Academic Publishers. 318 OPEN PROBLEMS IN STRONGLY CORRELATED ELECTRON SYSTEMS laws, is directly responsible for finite temperature anomalous transport properties of such systems [4, 5]. Quantum mixing property is defined (here) relative with respect to a certain (Hilbert) (sub)space 6 of observables and is most conveniently discussed in Heisenberg picture where the equations of motion for quantum observable A(t) are addressed. Our approach may be interpreted also as a kind ofgeneralization ofLee's recurrence-relations methodof constructing continued-fraction representation of time-correlation functions [6]. After a short outline of the main ideas we will illustrate the method by describing a recent calculation [8] of the time-correlation functions for Ising s = 1/2 spin chain periodically kicked with transverse magnetic field. Quantum dynamics in Heisenberg picture is governed by the Heisenberg equation (d/dt)A(t) = i( adH(t))A(t), where (adH)A:= [H, A1 =HA AH. Our aim is to include description of dynamics of periodically time-dependent (e.g. kicked) quantum systems so we will adopt Floquet formalism and integrate equations of motion over one period of time T (for autonomous systems, T may be arbitrary, e.g. infinitesimally small), defining the so called quantum Floquet map or one-step propagatorU := T exp(-i J; dtH(t) )which induces the adjoint Floquet map U ad acting on the (sub)space 6 of observables A(t +T) = U ad A(t), U A:= utAU. We should require the subspace 6 to be invariat under dynamics, namely A E 6 :::} U ad A E 6. Linear space 6 is turned into the Hilbert space by introducing an innner product of two extensive observables in thermodynamic limit (size L -+ 00) as (AlB) = (AtB), (A):= lim .!.. tr A £--+00 L tr 1 (1) (2) which has a meaning of a canonical average at infinite temperature. 1 We note that the adjoint Floquet map U ad is a unitary operator in 6 w.r.t. inner product (1), or equivalently, that ad H is selfadjoint. Thus a pair (U ad , 6) defines a unitary evolution to which we apply standard terminology and results of spectral theory of dynamical systems (see e.g. part III of ref. [1]). To any observable A E 6 we assign spectral measure dI-'A(fJ), where the spectral parameter lies on a unit circle () E 8 1 == [0, 21f), which is directly related to the time-autocorrelation function of A (t) (AtA(mT)) = (AI[U 1 A) = / dI-'A(fJ)e • So, ergodic properties of the unitary evolution are dual to its spectral properties: IThis is the simplest choice which makes sense also for time-dependent systems. For autonomous systems one may work with more general finite-temperature canonical inner product (see e.g. [6]). Ergodic Properties ofQuantum Spin Chains: Kicked Transverse Ising Model 319 • Quantum dynamics is ergodic in 6 iff 1 is (at most) simple (nondegenerate) eigenvalue of Uad (with eigenvector 1 if 1 E 6). Then limT-too ~ JOT dt {(A(t)IB) (All)(IIB)} = O. • Quantum dynamics is mixing in 6 iflimt-too(A(t)IB) = (All)(IIB) for arbitrary A, B E 6. Then U ad has an absolutely continuous spectrum, except possibly for a simple eigenvalue 1. • If 1 is a multiply degenerate eigenvalue of U ad then the corresponding eigenvectors Qn are the conservation laws of 6, U ad Qn = Qn Let Qn be a complete set of orthonormalized (QnIQm) = dnm eigenvectors with eigenvalue 1. The averaged time-autocorrelation function is the spectral weight at point 1 1 jT DA:= lim -T dt(A(t)IA) = L. I(QnIA )1 • T-too 2 -T n (3) For example, if A = J is the particle/spin current D J is proportional to the charge/spin stiffness at infinite temperature and D J > 0 indicates ballistic transport. Even when the set of known conservation laws is incomplete, the above relation gives rigorous lower bound on DA [7, 5]. Next we propose a simple but nontrivial application of the above picture. We consider aHilbert space 6 ofobservables over infinite spin 1/2 chains described by Pauli operators a}, s E {x, y, z}, j E Z, which is spanned by two infinite sequences of selfadjoint observables Un and Vn,nEZ, namely n ~ 1, n =0, n <-1 ,

•Jan 1, 2011•0 citations

Discrete time crystals in many-body quantum chaos

Pekik Nurwantoro, Raditya Weda Bomantara, Jiangbin Gong

Discrete time crystals (DTCs) are phases of matter characterized by the presence of an observable evolving with $nT$ periodicity under a $T$-periodic Hamiltonian, where $n>1$ is an integer insensitive to small parameter variations. In particular, DTCs with $n=2$ have been extensively studied in periodically quenched and kicked spin systems in recent years. In this paper, we study the emergence of DTCs in a many-body system whose semiclassical mean-field dynamics is nonintegrable, using a rather simple model depicting a harmonically driven spin chain. We advocate to first employ a semiclassical approximation to arrive at a mean-field Hamiltonian and then identify the parameter regime at which DTCs exist, with standard tools borrowed from studies of classical chaos. Specifically, we seek symmetric-breaking solutions by examining the stable islands on the Poincar\'e surface of section of the mean-field Hamiltonian. We then turn to the actual many-body quantum system, evaluate the stroboscopic dynamics of the total magnetization in the full quantum limit, and verify the existence of DTCs. Our effective and straightforward approach indicates that in general DTCs are one natural aspect of many-body quantum chaos with mixed classical phase space structure. Our approach can also be applied to general time-periodic systems, which is thus promising for finding DTCs with $n>2$ and opening possibilities for exploring DTCs properties beyond their time-translational breaking features.

arXiv:1909.04337•Physical review. B./Physical review. B•Dec 23, 2019•24 citations

Ergodic properties of a generic nonintegrable quantum many-body system in the thermodynamic limit

Tomaž Prosen

We study a generic but simple nonintegrable quantum many-body system of locally interacting particles, namely, a kicked-parameter (t,V) model of spinless fermions on a one-dimensional lattice (equivalent to a kicked Heisenberg XX-Z chain of 1/2 spins). The statistical properties of the dynamics (quantum ergodicity and quantum mixing) and the nature of quantum transport in the thermodynamic limit are considered as the kick parameters (which control the degree of nonintegrability) are varied. We find and demonstrate ballistic transport and nonergodic, nonmixing dynamics (implying infinite conductivity at all temperatures) in the integrable regime of zero or very small kick parameters, and more generally and importantly, also in the nonintegrable regime of intermediate values of kicked parameters, whereas only for sufficiently large kick parameters do we recover quantum ergodicity and mixing implying normal (diffusive) transport. We propose an order parameter (charge stiffness D) which controls the phase transition from nonmixing and nonergodic dynamics (ordered phase, D>0) to mixing and ergodic dynamics (disordered phase, D=0) in the thermodynamic limit. Furthermore, we find exponential decay of time correlation functions in the regime of mixing dynamics. The results are obtained consistently within three different numerical and analytical approaches: (i) time evolution of a finite system and direct computation of time correlation functions, (ii) full diagonalization of finite systems and statistical analysis of stationary data, and (iii) algebraic construction of quantum invariants of motion of an infinite system, in particular the time-averaged observables.

arXiv:cond-mat/9808150•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Oct 1, 1999•117 citations

Prolonging a discrete time crystal by quantum-classical feedback

Gonzalo Camacho, Benedikt Fauseweh

Nonequilibrium phases of quantum matter featuring time crystalline eigenstate order have been realized recently on noisy intermediate-scale quantum (NISQ) devices. While ideal quantum time crystals exhibit collective subharmonic oscillations and spatiotemporal long-range order persisting for infinite times, the decoherence time of current NISQ devices sets a natural limit to the survival of these phases, restricting their observation to a shallow quantum circuit. Here we propose a time-periodic scheme that leverages quantum-classical feedback protocols in subregions of the system to enhance a time crystal signal significantly exceeding the decoherence time of the device. As a case of study, we demonstrate the survival of the many-body localized discrete time crystal phase in the one-dimensional periodically kicked Ising model, accounting for decoherence of the system with an environment. Based on classical simulation of quantum circuit realizations we find that this approach is suitable for implementation on existing quantum hardware and presents a prospective path to simulate complex quantum many-body dynamics that transcend the low depth limit of current digital quantum computers. Published by the American Physical Society 2024

arXiv:2309.02151•Physical Review Research•Jul 23, 2024•5 citations

Realization of a discrete time crystal on 57 qubits of a quantum computer

Philipp Frey, Stephan Rachel

Unconventional dynamical phases that violate ergodicity have been a subject of extensive research in recent years. A periodically driven system is naively expected to lose all memory of its initial state due to thermalization, yet this can be avoided in the presence of many-body localization. A discrete time crystal represents a driven system whose local observables spontaneously break time translation symmetry and retain memory of the initial state indefinitely. Here, we report the observation of a discrete time crystal on a chain consisting of 57 superconducting qubits on a state-of-the-art quantum computer. We probe random initial states and compare the cases of vanishing and finite disorder to distinguish many-body localization from prethermal dynamics. We further report results on the dynamical phase transition between the discrete time crystal and a thermal regime, which is observed via critical fluctuations in the system's subharmonic frequency response and a substantial speedup of spin depolarization.

Science Advances•Mar 2, 2022•150 citations

Self-ordered Time Crystals: Periodic Temporal Order Under Quasiperiodic Driving

Sayan Choudhury, W. Vincent Liu

A discrete time is a remarkable non-equilibrium phase of matter characterized by persistent sub-harmonic response to a periodic drive. Motivated by the question of whether such time-crystalline order can persist when the drive becomes aperiodic, we investigate the dynamics of a Lipkin-Meshkov-Glick model under quasiperiodic kicking. Intriguingly, this infinite-range-interacting spin chain can exhibit long-lived periodic oscillations when the kicking amplitudes are drawn from the Thue-Morse sequence (TMS). We dub this phase a ``self-ordered time crystal (SOTC), and demonstrate that our model hosts at least two qualitatively distinct prethermal SOTC phases. These SOTCs are robust to various perturbations, and they originate from the interplay of long-range interactions and the recursive structure of the TMS. Our results suggest that quasiperiodic driving protocols can provide a promising route for realizing novel non-equilibrium phases of matter in long-range interacting systems.

•Sep 11, 2021•0 citations
We analytically study the time evolution of the expectation values of observables in periodically kicked many-body quantum systems. Starting from an initial state, we compute both the transient and the long-time properties of the observables. Our derivation explains the criteria and the mechanism that lead to the infinite-temperature statistical average of observables at long times, irrespective of the initial state. When the criteria are violated, the observables oscillate with time. These oscillations are subharmonic and robust to small perturbations, suggesting the emergence of a discrete time crystal phase. We demonstrate these features explicitly in periodically kicked nonintegrable spin chains. For a spin chain with two kicks per cycle, we show that the kicked chain can exhibit an ergodic or a discrete-time crystal phase for the same kicking strengths, depending on the initial state preparation. We complement our time-evolution study of observables with the spectral form factor of these kicked models.
The problem of calculation of time-correlation functions and relaxation phenomena in interacting quantum many-body systems or quantized fields is put into the framework of ergodic theory and dynamical systems. In this spirit we discuss spectral analysis of the adjoint propagator in a suitable Hilbert space of quantum observables in Heisenberg picture as an alternative approach to characterize many-body dynamics. We illustrate our approach by working out the time-autocorrelation functions at infinite temperature of a large class of observabIes for a spin s = 1/2 Ising chain in a periodically kicking transverse magnetic field, and show that they decay to their asymptotic values as ex: t3/ 2 General understanding of time-correlation functions and relaxation phenomena of interacting quantum many-body systems is a challenging problem with only very partial existing solutions. Situation is much better if one considers the motion of a one-body or few-body systems where the classical ergodic theory (see e.g. [1]) describes a variety of different types of motion, ranging from quasi-periodic integrable motion characterized by a discrete spectrum of the (Liouville) evolution operator, through KAM quasi-integrability, to ergodicity, mixing and chaos characterized by a continuous evolution spectrum. Mixing, which is equivalent of saying that correlation functions of an arbitrary pair of observables (A,B) decay in time, limt-HXJ ((A(t)B) (A)(B)) = 0, is the necessary dynamical property needed to justify the relaxation to equilibrium (micro/grand/canonical) state and the laws of statistical mechanics, such as the fluctuation-dissipation theorem and transport laws. Mixing in quantum many-body system has been investigated by numerical experiment only quite recently [2]. In another numerical experiment [3], a suggestive possibility of dynamical phase transition from non-ergodic/non-mixing to ergodic/mixing dynamics in a non-integrable periodically kicked chain of interacting spinless fermions (kicked t-V model) has been proposed. General absence of mixing in completely integrable quantum many-body systems (such as those solvable by Bethe ansatz), being the consequence of an infinite sequence of conservation 317 J. Bonca et al. (eds.), Open Problems in Strongly Correlated Electron Systems. 317-322. © 2001 Kluwer Academic Publishers. 318 OPEN PROBLEMS IN STRONGLY CORRELATED ELECTRON SYSTEMS laws, is directly responsible for finite temperature anomalous transport properties of such systems [4, 5]. Quantum mixing property is defined (here) relative with respect to a certain (Hilbert) (sub)space 6 of observables and is most conveniently discussed in Heisenberg picture where the equations of motion for quantum observable A(t) are addressed. Our approach may be interpreted also as a kind ofgeneralization ofLee's recurrence-relations methodof constructing continued-fraction representation of time-correlation functions [6]. After a short outline of the main ideas we will illustrate the method by describing a recent calculation [8] of the time-correlation functions for Ising s = 1/2 spin chain periodically kicked with transverse magnetic field. Quantum dynamics in Heisenberg picture is governed by the Heisenberg equation (d/dt)A(t) = i( adH(t))A(t), where (adH)A:= [H, A1 =HA AH. Our aim is to include description of dynamics of periodically time-dependent (e.g. kicked) quantum systems so we will adopt Floquet formalism and integrate equations of motion over one period of time T (for autonomous systems, T may be arbitrary, e.g. infinitesimally small), defining the so called quantum Floquet map or one-step propagatorU := T exp(-i J; dtH(t) )which induces the adjoint Floquet map U ad acting on the (sub)space 6 of observables A(t +T) = U ad A(t), U A:= utAU. We should require the subspace 6 to be invariat under dynamics, namely A E 6 :::} U ad A E 6. Linear space 6 is turned into the Hilbert space by introducing an innner product of two extensive observables in thermodynamic limit (size L -+ 00) as (AlB) = (AtB), (A):= lim .!.. tr A £--+00 L tr 1 (1) (2) which has a meaning of a canonical average at infinite temperature. 1 We note that the adjoint Floquet map U ad is a unitary operator in 6 w.r.t. inner product (1), or equivalently, that ad H is selfadjoint. Thus a pair (U ad , 6) defines a unitary evolution to which we apply standard terminology and results of spectral theory of dynamical systems (see e.g. part III of ref. [1]). To any observable A E 6 we assign spectral measure dI-'A(fJ), where the spectral parameter lies on a unit circle () E 8 1 == [0, 21f), which is directly related to the time-autocorrelation function of A (t) (AtA(mT)) = (AI[U 1 A) = / dI-'A(fJ)e • So, ergodic properties of the unitary evolution are dual to its spectral properties: IThis is the simplest choice which makes sense also for time-dependent systems. For autonomous systems one may work with more general finite-temperature canonical inner product (see e.g. [6]). Ergodic Properties ofQuantum Spin Chains: Kicked Transverse Ising Model 319 • Quantum dynamics is ergodic in 6 iff 1 is (at most) simple (nondegenerate) eigenvalue of Uad (with eigenvector 1 if 1 E 6). Then limT-too ~ JOT dt {(A(t)IB) (All)(IIB)} = O. • Quantum dynamics is mixing in 6 iflimt-too(A(t)IB) = (All)(IIB) for arbitrary A, B E 6. Then U ad has an absolutely continuous spectrum, except possibly for a simple eigenvalue 1. • If 1 is a multiply degenerate eigenvalue of U ad then the corresponding eigenvectors Qn are the conservation laws of 6, U ad Qn = Qn Let Qn be a complete set of orthonormalized (QnIQm) = dnm eigenvectors with eigenvalue 1. The averaged time-autocorrelation function is the spectral weight at point 1 1 jT DA:= lim -T dt(A(t)IA) = L. I(QnIA )1 • T-too 2 -T n (3) For example, if A = J is the particle/spin current D J is proportional to the charge/spin stiffness at infinite temperature and D J > 0 indicates ballistic transport. Even when the set of known conservation laws is incomplete, the above relation gives rigorous lower bound on DA [7, 5]. Next we propose a simple but nontrivial application of the above picture. We consider aHilbert space 6 ofobservables over infinite spin 1/2 chains described by Pauli operators a}, s E {x, y, z}, j E Z, which is spanned by two infinite sequences of selfadjoint observables Un and Vn,nEZ, namely n ~ 1, n =0, n <-1 ,
Discrete time crystals (DTCs) are phases of matter characterized by the presence of an observable evolving with nTnTnT periodicity under a TTT-periodic Hamiltonian, where n>1n>1n>1 is an integer insensitive to small parameter variations. In particular, DTCs with n=2n=2n=2 have been extensively studied in periodically quenched and kicked spin systems in recent years. In this paper, we study the emergence of DTCs in a many-body system whose semiclassical mean-field dynamics is nonintegrable, using a rather simple model depicting a harmonically driven spin chain. We advocate to first employ a semiclassical approximation to arrive at a mean-field Hamiltonian and then identify the parameter regime at which DTCs exist, with standard tools borrowed from studies of classical chaos. Specifically, we seek symmetric-breaking solutions by examining the stable islands on the Poincar\'e surface of section of the mean-field Hamiltonian. We then turn to the actual many-body quantum system, evaluate the stroboscopic dynamics of the total magnetization in the full quantum limit, and verify the existence of DTCs. Our effective and straightforward approach indicates that in general DTCs are one natural aspect of many-body quantum chaos with mixed classical phase space structure. Our approach can also be applied to general time-periodic systems, which is thus promising for finding DTCs with n>2n>2n>2 and opening possibilities for exploring DTCs properties beyond their time-translational breaking features.
We study a generic but simple nonintegrable quantum many-body system of locally interacting particles, namely, a kicked-parameter (t,V) model of spinless fermions on a one-dimensional lattice (equivalent to a kicked Heisenberg XX-Z chain of 1/2 spins). The statistical properties of the dynamics (quantum ergodicity and quantum mixing) and the nature of quantum transport in the thermodynamic limit are considered as the kick parameters (which control the degree of nonintegrability) are varied. We find and demonstrate ballistic transport and nonergodic, nonmixing dynamics (implying infinite conductivity at all temperatures) in the integrable regime of zero or very small kick parameters, and more generally and importantly, also in the nonintegrable regime of intermediate values of kicked parameters, whereas only for sufficiently large kick parameters do we recover quantum ergodicity and mixing implying normal (diffusive) transport. We propose an order parameter (charge stiffness D) which controls the phase transition from nonmixing and nonergodic dynamics (ordered phase, D>0) to mixing and ergodic dynamics (disordered phase, D=0) in the thermodynamic limit. Furthermore, we find exponential decay of time correlation functions in the regime of mixing dynamics. The results are obtained consistently within three different numerical and analytical approaches: (i) time evolution of a finite system and direct computation of time correlation functions, (ii) full diagonalization of finite systems and statistical analysis of stationary data, and (iii) algebraic construction of quantum invariants of motion of an infinite system, in particular the time-averaged observables.
Nonequilibrium phases of quantum matter featuring time crystalline eigenstate order have been realized recently on noisy intermediate-scale quantum (NISQ) devices. While ideal quantum time crystals exhibit collective subharmonic oscillations and spatiotemporal long-range order persisting for infinite times, the decoherence time of current NISQ devices sets a natural limit to the survival of these phases, restricting their observation to a shallow quantum circuit. Here we propose a time-periodic scheme that leverages quantum-classical feedback protocols in subregions of the system to enhance a time crystal signal significantly exceeding the decoherence time of the device. As a case of study, we demonstrate the survival of the many-body localized discrete time crystal phase in the one-dimensional periodically kicked Ising model, accounting for decoherence of the system with an environment. Based on classical simulation of quantum circuit realizations we find that this approach is suitable for implementation on existing quantum hardware and presents a prospective path to simulate complex quantum many-body dynamics that transcend the low depth limit of current digital quantum computers. Published by the American Physical Society 2024
Unconventional dynamical phases that violate ergodicity have been a subject of extensive research in recent years. A periodically driven system is naively expected to lose all memory of its initial state due to thermalization, yet this can be avoided in the presence of many-body localization. A discrete time crystal represents a driven system whose local observables spontaneously break time translation symmetry and retain memory of the initial state indefinitely. Here, we report the observation of a discrete time crystal on a chain consisting of 57 superconducting qubits on a state-of-the-art quantum computer. We probe random initial states and compare the cases of vanishing and finite disorder to distinguish many-body localization from prethermal dynamics. We further report results on the dynamical phase transition between the discrete time crystal and a thermal regime, which is observed via critical fluctuations in the system's subharmonic frequency response and a substantial speedup of spin depolarization.
A discrete time is a remarkable non-equilibrium phase of matter characterized by persistent sub-harmonic response to a periodic drive. Motivated by the question of whether such time-crystalline order can persist when the drive becomes aperiodic, we investigate the dynamics of a Lipkin-Meshkov-Glick model under quasiperiodic kicking. Intriguingly, this infinite-range-interacting spin chain can exhibit long-lived periodic oscillations when the kicking amplitudes are drawn from the Thue-Morse sequence (TMS). We dub this phase a ``self-ordered time crystal (SOTC), and demonstrate that our model hosts at least two qualitatively distinct prethermal SOTC phases. These SOTCs are robust to various perturbations, and they originate from the interplay of long-range interactions and the recursive structure of the TMS. Our results suggest that quasiperiodic driving protocols can provide a promising route for realizing novel non-equilibrium phases of matter in long-range interacting systems.