arXiv:2007.08348v3 [cond-mat.quant-gas] 28 Dec 2020 inmyb usinbe e ln t eainhpwith relationship ter, mat- its quantum non-equilibrium alone of let properties fundamental questionable, be defini- may its even “excita- tion matter, the quantum about to non-equilibrium known of compared tion” is However, less state much decade. cases, solid equilibrium past to the atoms[2] vari- over ultracold in physics[3] from interest ranging considerable fields received ous coun- has equilibrium and quantum its terparts, from Non-equilibrium excita- differs wave fundamentally lattice phonons[1]. matter of as number known few solid tions a a by of determined properties elastic are and thermal temperatures. low the relatively example, or For of fields ow- external excited properties weak be to to the ing likely are are which states, important behav- excited observable physically low-lying its What’s with do to ior. much have not does tem iesoa 1)sse,asotnosdsrt spatial discrete ( spontaneous one- translational a a for system, example, (1D) For dimensional excitation, structure. its elementary determines an also affect of but only spectrum not and does dynamics it the excitation: re- elementary closely the is to SSB lated the systems, sponta- equilibrium thus symme- In translational try(DTTS). period, temporal discrete in- Hamiltonian the break an the neously observables are of physical periods multiple Its whose teger oscillations 26]. persistent [25, develop driving settings has periodic non-equilibrium 24], with [23, in equilibrium realized in experimentally forbidden been despite be state, to intriguing proven a being Such in forms[5–22]. interest different considerable (SSB) its attracted breaking have symmetry and spontaneous paradigm[4], pos- the new allowed for has phase sibilities (TC) crystal time the matter, oio-ieectto:atplgcldfc sandwiched defect two topological by a excitation: soliton-like sal,tegon tt nryo neulbimsys- equilibrium an of energy state ground the Usually, sapooyia xml fnnqiiru quantum nonequilibrium of example prototypical a As yaia rniin n rtclbhvo ewe discr between behavior critical and transitions Dynamical e.g. Z h ymtisadtersotnosbreaking. spontaneous their and symmetries the 2 ymtybekn hss2] profound A phases[27]. breaking symmetry n ol nueadnmcltasto ewe two between p transition br dynamical we a symmetry model, q induce translational solvable could non-equilibrium temporal ing exactly a spontaneous an of with behavior on phase excitation-like al based an but paper, study excitations, s this to elementary its In and of symmetry spectrum tures. the and other: dynamics each the with related closely cepts iceetm rsa hs gis os perturbations. and noisy processes against ramping phase transiti phase crystal the such time statis point discrete of nonequilibrium which sequences in stochastic at phenomenon of value o down critical p analogue slowing a in critical temporal states observe wave a we density as charge rate, considered degenerate two be by can sandwiched which domain, time e.g. 1 neulbimpyis pnaeu ymtybekn and breaking symmetry spontaneous physics, equilibrium In colo hsc n srnm,Saga ioTn Univers Tong Jiao Shanghai Astronomy, and Physics of School ice unu etradKyLbrtr fAtfiilStru Artificial of Laboratory Key and Center Quantum Wilczek Z 2 ymtybekn loscertain allows breaking symmetry ) 2 hnhiRsac etrfrQatmSine,Saga 20 Shanghai Sciences, Quantum for Center Research Shanghai ioi Yang Xiaoqin 1 n iCai Zi and eato,wihi iedpnetbtde o decay prefactor the not and does but distance. time-dependent with is which teraction, tra site at ator al clswt h ytmsz.I iatt lattice bipartite a In size. lin- system energy the ( interacting with total scales the early that guarantees Eq.(1) h neato ewe aro oosaeattractive are bosons of pair a between interaction the igapiue ˆ amplitude, ping where usini o ognrlz uha dait h non- ( the physics into absent equilibrium SSB idea in intriguing an with such matters quantum generalize equilibirum to how is question h iceetm rsa hs gis os perturba- noisy against phase of crystal stability tions. the time phase to discrete random implication the multiple its of and effect processes ramping the discuss 29]. we remi- phenomenon[28, down Finally, a slowing critical diverges, critical a the duration of transition observe niscence the can which By one at velocity, Fig.1. point, in ramping shown the as tuning phases, DTC breaking symmetry ooi oe iha nnt ogrneinteraction. long-range follows: as infinite reads an Hamiltonian The with model bosonic inmyidc yaia rniinbtentwo between transition dynamical perturba- a slow induce a that may of shown tion response is It the observable. monitor physical then the we trans- and time periodical symmetry, original the lational the of breaks impose transiently We top which on driving, phase. (DTC) perturbation crystal time-dependent time a discrete signa- a a manifest quantities, of can physical ture for that response model sub-harmonic bosonic a interacting 1D driven e.g, Z H oe n method: and Model nti ae,w drs hsiseb tdiga studying by issue this address we paper, this In 2 = ymtybekn iecytlpae in phases crystal time breaking symmetry Dltiea norcs) q()idctsthat indicates Eq.(1) case), our in as lattice 1D a J atmmte:adsrt iecrystal time discrete a matter: uantum − i hsc.W lodsusteeffect the discuss also We physics. tic stesnl-atcenaetniho(N hop- nearest-neighbor(NN) single-particle the is ndrto iegs eebigthe resembling diverges, duration on J otnosbekn o nycontrol only not breaking pontaneous aig ti hw htso ramp- slow that shown is It eaking. t mlcto otesaiiyo the of stability the to implication its ,2, 1, laeyee ytnn h ramping the tuning By olyacetylene. odtrieterudryn struc- underlying their determine so X lmnayectto r w con- two are excitation elementary i i ooeapaermigprotocol ramping phase a ropose ∗ . h oio xiainspacially excitation soliton the f ( t,Saga 020 China 200240, Shanghai ity, trsadQatmControl, Quantum and ctures ˆ b V i † L ˆ b n ( i i t +1 L 1 stesse ieo h Dlattice, 1D the of size system the is stesrnt ftealt-l in- all-to-all the of strength the is ) = nfoto h neato em of terms interaction the of front in + t iecytlphases crystal time ete 35 China 1315, ˆ b e.g. h.c i † ˆ b ecnie Dhard-core 1D a consider We i ) stepril ubroper- number particle the is TSbreaking). DTTS − V L ( t ) X ij ( − 1) i − j n ˆ i n ˆ j (1) Z 2 2
(repulsive) if they are located in the same (different) sub- lattice. Such an infinitely long-range interaction has been realized in recent high-finesse cavity experiments[30, 31] by coupling bosons to the cavity vacuum mode whose period is twice of that of the optical lattice. The total particle number N is conserved in our system, and we fo- cus on the case of half-filling (N = L/2) throughout this paper. In the equilibrium case (V (t) = V0), the ground state of the 1D Ham.(1) is always a Mott-insulator with a charge-density-wave (CDW) order for arbitrary positive V0. Owing to the all-to-all coupling of the interaction in Ham.(1), in the thermodynamical limit the mean- field method is not an approximation but an exact FIG. 1: (Color online). Schematic diagram of (a) a period method. In particular, it provides an exact descrip- doubling dynamics in the presence of periodical driving and tion of both the equilibrium properties and real-time two degenerate TC phases and (b)the phase ramping protocol dynamics of the system[32, 33] (see also Supplementary in our model and the dynamical transition induced by it. Material(SM)[34]). The mean-field method allows us to decouple the all-to-all interaction by introducing an aux- iliary staggered field which is self-consistently determined (see SM[34]). The dynamical behavior is even richer and during the evolution. The Ham.(1) can be expressed as: more interesting when we introduce periodical driving into Ham.(1), e.g. V (t > 0) = Vf + δ cos2πt with δ i H¯ (t)= J [b†bi + h.c]+ m(t)V (t) ( 1) nˆi (2) the driving amplitude. We observe that[34] depending − i +1 − Xi Xi on the different choices of the initial states and driving amplitude δ, the long-time dynamics could exhibit a pe- 1 i where m(t) = Ψ(t) ( 1) nˆi Ψ(t) and Ψ(t) is the h | L − | i | i riodic oscillation with a frequency that is either identi- wavefunction of the system at time t. The time evolu- cal to or independent of the driving frequency; the for- tion under the Ham.(2) can be solved exactly by perform- mer can be considered as a synchronization phenomenon ing the Jordan-Wigner transformation to transform 1D which has been observed in periodically-driven integrable hard-core bosons into spinless fermions, and the Ham.(2) systems[36]. In addition, it can also exhibit quasi- transform into a non-interacting fermionic model. periodic oscillations with a multi-period structure[34]. We choose the periodic boundary condition (PBC), Most interesting dynamics can be observed in the in- which allows us to perform the Fourier transformation, termediate driving regime (see Fig.2 b), where m(t) ex- after which the fermionic Hamiltonian turns to hibits a persistent oscillation whose period is twice that
[c† c† ] εk m(t) ck of the external driving period: a signature of DTC[37] H¯ (t)= k k+π (3) m(t) εk+π ck+π that the DTTS in Ham.(1) H(t) = H(t + T ) has been Xk spontaneously broken (m(t) = m(t +2T ) = m(t + T ) 6 where the summation is over the momentum in the first with T = 1 is the period of driving). This phenomena π π is rooted in the non-linearity of the self-consistent mean- Brillouin zone of Ham.(2) (k [ 2 , 2 ]), and ck (ck† ) denotes the annihilation(creation)∈ − operator of the spin- field equation of motion. However, unlike the period dou- less fermion. εk = 2J cos k, thus εk = εk+π. Eq.(3) bling phenomena in the non-linear classical systems (e.g.a indicates that the dynamics− of the system− can be consid- driven-dissipative pendulum[38]) or the “dissipative TC” ered as a collective behavior of different k modes, each of in open quantum systems[13], our system is dissipation- which is a two-level quantum system subjected to time- less without entropy generation. Ham.(1) is free from dependent field that was self-consistently determined as disorder, hence it is the integrability[39] rather than the 1 m(t)= Ψk(t) c† ck π Ψk(t) . many-body localization[9] that prevents our model from L kh | k + | i DiscreteP time crystal. Despite the triviality of the being heated to an infinite temperature state. ground state phase diagram, the system can exhibit rich Dynamical transition between Z2 symmetry breaking dynamical behavior in the presence of a time-dependent DTC phases: Despite the richness of the dynamics behav- V (t). For instance, in a quantum quench protocol, one ior, here we will study neither the global non-equilibrium can start from a ground state of Ham.(1) with V (t =0)= phase diagram of our model, nor the mathematical origin Vi, suddenly change it to a different value V (t> 0) = Vf behind this DTC phase. Instead, we will use this exactly and let the system evolve under this new Hamiltonian. solvable model as a starting point to study a non-trivial It has been shown that the long-time dynamics of this dynamical behavior of the DTC phase, which can be con- model can exhibit either persistent oscillations or ther- sidered as a temporal analogy of the soliton excitation malization depending on the choice of initial states[35] in Su-Schrieffer-Heeger (SSH) model[27]. Owing to the 3
4 (a)
1 (a) 1.0
3
10
3 (t) V(t) n S
D T =0.5
0
0.5
0.5 (b) T =8.0
0
2
2
10
c
T -T 0.2 0.4
0 0 0.0 m(t)
0.0
(b)
transition duration T =0.5 0.5 (c)
0 1.0
T =8.0
0 (t) n
0.0 S m(t)
T =2.4
0
0.5
T =2.45
0 T =2.40
0.4 (d) 0
T =2.45
0
0.0 m(t)
0.0
-1
[J ] 200 400 t 600 800
-0.4
-1
[J ]
1500 1505 1510 t 1515 FIG. 3: (Color online). Dynamics of the relative displacement of the peak positions Sn(t) with (a) fast and slow and (b) in- FIG. 2: (Color online). (a)Periodical driving V (t > 0) = termediate ramping rates close to the dynamical transition Vf +δ cos 2πt with a period T0 = 1. (b) Time crystal dynamics c point T0 = 2.41. The inset of Fig.2 (a) shows the transition with a period 2 in the absence of ramping. (c) Long-time dy- duration close to the dynamical critical point. Other param- namics of m(t) with slow and fast ramping, which correspond eters are chosen to be the same as those in Fig.2 to two “degenerate” TC phases. (d)Long-time dynamics of m(t) with an intermediate ramping rate close to the dynami- cal transition point. The parameters are chosen as Vi = 10J, V = 3J, δ = 0.8J and L = 5000. f ters (Vf = 3J, δ = 0.8J) except the duration of ramp- ing T0. In the limit of T0 = 0 the periodic driving is abruptly changed by a phase of 2π at t = ti, thus the spontaneous breaking of the DTTS, the DTC phase are Hamiltonian is identical to that without ramping, so is supposed to be two-fold “degenerate”, each of which has the long-time dynamics. For a rapid ramping, the system a period 2 and can be connected to the other one by shift- exhibits a similar long-time dynamics, but with a weaker ing a half-period (1) along the temporal direction. Mo- amplitude, as shown in Fig.2 (c). In the opposite limit tivated by the soliton excitation 1D CDW system which of slow ramping, the peak positions of m(t) are pinned separates these two Z2 symmetry breaking “degenerate” to those of V (t). Thus, after the ramping, the phase of states in space[27], to realize such an object, one needs to V (t) is pushed forward by 2π, so is m(t). However, due impose a temporal perturbation on the system that tran- to the period doubling feature of the DTC, m(t) is only siently breaks the original time translational symmetry. shifted by a half-period, and thus falls into the other de- Here, we introduce an additional linear ramping of the generate state that differs from the one before the ramp- phase on top of the periodic driving: ing, although V (t) remains unchanged. Thus, by slowly ramping the driving phase, one can induce a dynamical V (t)= V + δ cos2π[t + φ(t)], (4) f transition between two “degenerate” DTC phases. On t ti the contrary, for a fast ramping, the system cannot follow where φ(t) = −T0 for t [ti,ti + T0], and φ(t) = 0 oth- ∈ V (t) “adiabatically”, thus finally relaxes to a DTC phase erwise. ti is the initial time of the ramping, and T0 is its duration, after which the external driving accumulates an similar to the original one. These two distinct dynamical additional 2π phase thus the Hamiltonian is identical the behaviors of slow and fast ramping indicate a transition one before the ramping. We assume prior to the ramping between them. From Fig.2 (d), we can find that this transition occurs suddenly at T 2.41(1) and there is (t
(a) tial decay of the DTC order parameter in time, which 0.1 indicates the instability of the DTC phase against such
0.0 noisy perturbations. On the contrary, the high-frequency
-0.1 fluctuation will not induce the dynamical transition, thus 1.0 preserve the long-range DTC order (Fig.8 (b)). It agrees
0.8
0.6 with our result obtained from a single ramping process: (t)
n
0.4
S there is a frequency threshold for a phase ramping pro- 0.2 tocol to induce a transition between the DTC phases.
0.0 Notice that this kind of particular telegraph frequency 4
(b)
2 fluctuation is different from those encountered in the real-
0
istic environment, thus an interesting question is whether
-2 our conclusion can apply to DTC phases with more re- -4 alistic noises, which usually perturb the systems contin- 0.8 uously thus can not be considered as rare events. An-
0.6 (t)
n
0.4 other important question is the effect of the thermal fluc- S
0.2 tuations induced by coupling the system to a thermal
0.0
-1
[J ] bath. Different from the noisy systems considered above, 0 200 400 t 600 800 a thermal bath does not only introduce noisy fluctua-
FIG. 4: (Color online). Dynamics of Sn(t) in the presence tions, but also dissipations. Since our original Ham.(1) of two typical stochastic sequences of phase ramping with is a quantum many-body model, its dynamics in the pres- α = 0.01, where the driving frequency randomly fluctuates ence of a thermal bath is beyond the scope of this work within the regime (a) |δω| ∈ [0.05, 0.2] and (b) |δω| ∈ [1, 5]. and will be left for the future. In general, dissipation tends to balance the noise effect, thus makes the DTC Stochastic sequences of phase ramping and its implica- phases more stable compared to the cases with only noisy tion to noisy systems. Up to now, we focus on the effect of driving. a single phase ramping protocol on the DTC phases, and Experimental realization and detection: The proposed find a dynamical transition between two DTC phases can model can be realized in an experimental setup simi- only be induced by a sufficiently slow ramping. What’s lar with to the one in Ref.[30]. In terms of detections, its implication to the stability of a DTC phase (e.g. its the CDW order parameter can be directly measured us- robustness against thermal fluctuation or external noise)? ing the superlattice band-mapping technique[40], or in- To address this issue, we perform a stochastic sequence directly through the heterodyne detection[30, 31]. For of phase ramping processes, whose initial time and du- an optical lattice with J 400Hz, one can estimate that ration are chosen randomly. Since the proposed phase the typical oscillation period≈ is in the order of a mag- ramping protocol in Eq.(4) equals to a shift of the driv- nitude of 10ms, which is a measurable time scale in the ing frequency from ω = 1(2π) to ω′ = [1 + δω](2π) with current experimental setup. The major obstacle for the 1 δω = T0 , a stochastic sequence of phase ramping pro- experimental observation is the dephasing mechanisms cesses with random durations (T0) can be considered as caused by various experimental imperfections including a telegraph-like fluctuation of the driving frequency (see the dissipation (particle loss) or decoherence effect, ther- Fig.8 for an instance), which can be considered as a con- mal fluctuation and the extra heating effect generated sequence of the perturbation from external noises. by our driving protocols, all of which will lead to the 5 damping of oscillations. Thus, the major experimental Supplemental material challenge is to control the dephasing rate and make it considerably smaller than the typical energy scale of the Hamiltonian, which enables us to observe and analyze VALIDITY OF THE TIME-DEPENDENT the dynamics of persistent oscillations before they damp. SELF-CONSISTENT MEAN-FIELD METHOD Conclusion and outlook: In conclusion, we propose a general protocol to realize a transition between Z2 sym- It is well known that for a model with infinitely long- metry breaking DTC phases, a temporal analogue of the range interaction, the mean-field method provides an ex- soliton excitation in 1D CDW system. A dynamical crit- act treatment for the equilibrium properties in the ther- ical point has been observed, at which the transition du- modynamic limit. In the following, we will prove that ration diverges. Future developments will include the this is also the case for the dynamical problems. analysis of the effect of the time-dependent perturba- We denote F is the Fock basis of hard-core bosons (e.g. tion on other time crystal phases with different trans- 0110 11 ), and divide the evolution period [0,t] into M | ··· i lational symmetry breaking. For instance, besides the slices with ∆t = t/M. In the limit of M , the path → ∞ DTC phase, we also observe a time crystal phase whose integral expression of the propagators can be written as: period is incommensurate with that of the driving[34], 0 i∆tH(t1) 1 1 i∆tH(tM ) M = F e− F F e− F dubbed “time quasicrystal”[41]. By making an anal- K h | | ih | · · · | i F1, XFM−1 ogy to the incommensurate CDW states, we expect that { ··· } the dynamical response of this time quasicrystal may be (5) j significantly different from that of the DTC. Further- where F is a Fock basis at the j-th time slice (j = | i j j more, in the absence of periodic driving, our model could 0, M), with Fj F F a unit matrix. We fur- ··· { } | ih | manifest a persistent oscillation in its quantum quench ther divide ourP Hamiltonian as H = Tˆ + Vˆ where dynamics[35], which represents a TC phase with broken Tˆ = J (ˆb†ˆbi + h.c), and the interacting energy can − i i +1 symmetry with respect to continuous temporal transla- be rewrittenP as: tions. It is interesting to study the dynamical response of 1 i 2 such continuous TC phase to time-dependent perturba- Vˆ = LV (t)[ ( 1) nˆi] . (6) L − tion, which should differ from excitation in the dissipative Xi time crystal[42] since our system is disspationless. Acknowledgments.—We acknowledge support by the By performing Suzuki-Trotter decomposition, one can National Key Research and Development Program of obtain that China (Grant No. 2016YFA0302001), NSFC of China i∆tH i∆tTˆ i∆tVˆ 2 e− = e− e− + (∆t ) (7) (Grant No. 11674221, No.11574200), Shanghai Mu- O nicipal Science and Technology Major Project (Grant Therefore, in the limit of ∆t 0, the propagator in No.2019SHZDZX01) and the Program Professor of Spe- Eq.(5) can be expressed as: → cial Appointment (Eastern Scholar) at Shanghai Institu- tions of Higher Learning. M 1 − i∆tV j (tj ) = Tj,j e− {F } (8) K +1 F1, XFM−1 jY=1 { ··· }
j i∆tTˆ j+1 where Tj,j+1 = F e− F , and V Fj (tj ) = h | | j i2 { } j i∆tVˆ j i∆tLV (ti) F j F e− F = e− S{ } , where F = h1 | i| i S{ } j L i( 1) ni ni F . − |{ }∈{ } PThe quadratic part in Eq.(8) can now be decoupled using the Hubbard-Stratonovic transformation by intro- ducing auxillary fields mj with j =1,M 1 as: { } − M 1 j − j 2 j j dm i∆tLV (tj )[(m ) 2m F ] = Tj,j+1e− − S{ } K Z jY=1 XFj { } N (9) where is a normalization factor to keep the identity of the GaussianN integral. In the thermodynamic limit L , Eq.(9) indicates that the saddle point approximation→ j ∞δ [m ] j j 2 j j Lδmj = 0 with [m ] = V (tj )[(m ) 2m F ], be- comes exact, whichL means that the auxillary− S{ fields} mj are given by the CDW order parameter of the state{ as} 6
mj = Fj . As a consequence, For the time evolution ) f( 0.5 (a) of a stateS{ ψ} in each infinitesimal time step (∆t 0), 0.02 the evolution| i operator eiH(t)δt are equal to the one→ af- 0.01 m(t) ¯ iH(t,m(t))∆t 0.00 ter mean-field approximation e with a time- 0 1 2 dependent CDW order parameters m(t), which is deter- 0.4 mined self-consistently during the time evolution by the ) ( saddle point method. f (b)
0.2 m(t)
0.5 0.1
0.0
0.0 0.5 1.0 1.5 2.0
V =0.75J
0
0.0
V =10J
0
0.4 )
(c) f(
0.01 m(t)
0.2
m(t)
0.00
0 1 2
0.2
0.0 ) f(
-1 (d)
0.2
[J ] m(t) 0 20 40 60 80 t
0.5
0.0
0.0 0.5 1.0 1.5 2.0 FIG. 5: (Color online). Quench dynamics of M(t) calculated by the self-consistent mean-field method with different Vi and 0.0 −5 the parameters Vf = 3J, L = 5000, ∆t = 10 , δ = 0;
-1
[J ]
200 220 240 t
FIG. 6: (Color online). Dependence of periodically driven QUANTUM QUENCH DYNAMICS WITHOUT dynamics of M(t) on driving amplitudes [(a) δ = 0.3J; (b) DRIVING δ = 0.8J; (c) δ = 8J;] and initial states [(d) Vi = 0.75J, δ = 0.8J]. Other parameters are chosen as Vi = 10J for − − (a)-(c), L = 5000 and ∆t = 10 5J 1 for (a)-(d). In this section, we study quantum quench dynamics of the model (δ = 0) by choosing an initial state as the ground state of H with V (t) = V , and suddenly i In Fig.(7), we plot the long-time dynamics of m(t) with change the interaction strength to a different value V f different driving amplitude δ and initial states V . For and let the system evolve under this new Hamiltonian. i Fig.7 (a)-(c), we fixed the initial state as the ground state The long-time dynamics of this model has been analyzed of the Hamiltonian with V = 10J, while for Fig.7 (d), in Ref.[1], here we will outline the mean-field results of i we change it to V =0.75J. As shown in Fig.7 (a), for a Ref.[1]. We plot the time evolution of the CDW order i weak periodical driving, the long-time dynamics becomes parameter m(t) = 1 ( 1)i n starting from two dif- L i i a persistent quasi-periodic oscillation with a multi-period ferent initial states in Fig.5,− fromh i which we can find two P structure, which can be reflected in its Fourier spectrum distinct dynamical behaviors even though the evolution f(ω) = 1 dteiωtm(t). As shown in the inset of Fig.7 Hamiltonian are the same (V = 3J). For the quench T f (a), the Fourier spectrum f(ω) exhibits sharp peaks at dynamics from the initial state with V = 10J, m(t) will R i different characteristic frequency, while the dominant one converge to a constant after sufficiently long time, while locates at the place exactly the same with the external for the other case with V = 0.75J, it persistently oscil- i driving frequency (ω = 1), which indicates m(t) is syn- lates with a period that spontaneously emerges during chronous with the external periodic driving. For an in- the quench dynamics. termediate driving amplitude, for instance δ = 0.8J as shown in Fig.7 (b), the peak position of f(ω) has been shifted from 1 to 0.5, which indicates a period-doubling PERIODICALLY DRIVEN DYNAMICS time crystal phase as we studied in the main text. When we further increase δ to a strongly driven regime (e.g. The periodically driven dynamical behavior in our δ = 8J as shown in Fig.7 c), the frequency spectrum model is very rich, while in the main text we only focus on exhibits a broad distribution, a signature of a chaotic the time crystal phase. Here we list several typical long- dynamics without a dominant time scale. If we start time behaviors, even though it is difficult to enumerate from an different initial state, for instance, the ground all the possibilities. state of Vi = 0.75J as shown in Fig.7 (d), we can find 7
)
t T =4J
0 (
0.5
L=500 T =J (a)
0
L=1000 5 T =0.2J
0
L=2000 m(t)
0.4
0
-1
[J ] 0 500 1000 1500 t
-1
[J ] 0 10 t 20
-3 0.5 (b) t=10 ˜ -4 FIG. 8: (Color online). Dynamics of θ(t) predicted by the
t=10
-5 Ham.(11) in the presence of fast, intermediate and slow ramp- t=10 ing with parameters α = γ = J, and θ˜ ∈ [− π , 3π ] with PBC.
m(t) 2 2
0.0 DYNAMICAL TRANSITIONS IN OTHER DTC MODELS
-1
[J ] 950 960 970 t 980 990 The proposed driving protocol is general in the sense that it can induce similar dynamical transitions between FIG. 7: (Color online). (a) Quench dynamics of M(t) with the DTC phases in other models. As an example, we −5 different system size L and the parameters δ =0, ∆t = 10 ; consider a classical nonlinear oscillator with the Hamil- (b) Periodically driven dynamics with different ∆t and δ = tonian: 0.8J, L = 5000; Vi = 10J and Vf = 3J for (a) and (b) H(t)= αI γI2 + V˜ (t)I cos2 θ (10) − m(t) exhibits an persistent oscillation with a frequency where we assume p = √I sin θ, q = √I cos θ with p/q ω = 0.46, which has nothing to do with the external the momentum/coordinate of the oscillator. The non- 2 ˜ driving frequency. linearity of the oscillator is encoded in I . V (t) = cos[2π(t + φ(t))] with φ(t) obeys the same phase ramp- ing protocol proposed above. In the absence of ramping (φ(t) = 0), it is known that such a driven oscillator could FINITE SIZE EFFECT AND CONVERGENCE OF support period doubling solution q(t) cos πt in cer- THE NUMERICAL RESULTS tain parameter regime[22]. By imposing∼ phase ramping on top of the period doubling dynamics, we assume that q(t) cos[πt+θ˜(t)], where θ˜(t) represents the phase devi- In our numerical simulations, we choose the lattice size ∼ as L = 5000. Here we will show that within our simu- ation from the period doubling solution, which resembles 1 lation period (0