Arxiv:2007.08348V3

Dynamical transitions and critical behavior between discrete time crystal phases 1 1,2, Xiaoqin Yang and Zi Cai ∗ 1Wilczek Quantum Center and Key Laboratory of Artificial Structures and Quantum Control, School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 200240, China 2Shanghai Research Center for Quantum Sciences, Shanghai 201315, China In equilibrium physics, spontaneous symmetry breaking and elementary excitation are two con- cepts closely related with each other: the symmetry and its spontaneous breaking not only control the dynamics and spectrum of elementary excitations, but also determine their underlying structures. In this paper, based on an exactly solvable model, we propose a phase ramping protocol to study an excitation-like behavior of a non-equilibrium quantum matter: a discrete time crystal phase with spontaneous temporal translational symmetry breaking. It is shown that slow ramping could induce a dynamical transition between two Z2 symmetry breaking time crystal phases in time domain, which can be considered as a temporal analogue of the soliton excitation spacially sandwiched by two degenerate charge density wave states in polyacetylene. By tuning the ramping rate, we observe a critical value at which point the transition duration diverges, resembling the critical slowing down phenomenon in nonequilibrium statistic physics. We also discuss the effect of stochastic sequences of such phase ramping processes and its implication to the stability of the discrete time crystal phase against noisy perturbations. Usually, the ground state energy of an equilibrium sys- question is how to generalize such an idea into the non- tem does not have much to do with its observable behav- equilibirum quantum matters with intriguing SSB absent ior. What’s physically important are the properties of in equilibrium physics (e.g. DTTS breaking). low-lying excited states, which are likely to be excited ow- In this paper, we address this issue by studying a ing to weak external fields or relatively low temperatures. driven 1D interacting bosonic model that can manifest For example, the thermal and elastic properties of a solid a sub-harmonic response for physical quantities, a signa- are determined by a few number of lattice wave excita- ture of a discrete time crystal (DTC) phase. We impose tions known as phonons[1]. Non-equilibrium quantum a time-dependent perturbation on top of the periodical matter fundamentally differs from its equilibrium coun- driving, which transiently breaks the original time trans- terparts, and has received considerable interest in vari- lational symmetry, and we then monitor the response of ous fields ranging from ultracold atoms[2] to solid state the physical observable. It is shown that a slow perturba- physics[3] over the past decade. However, compared to tion may induce a dynamical transition between two Z2 equilibrium cases, much less is known about the “excita- symmetry breaking DTC phases, as shown in Fig.1. By tion” of non-equilibrium quantum matter, even its defini- tuning the ramping velocity, one can observe a critical tion may be questionable, let alone its relationship with point, at which the transition duration diverges, a remi- fundamental properties of non-equilibrium quantum mat- niscence of the critical slowing down phenomenon[28, 29]. ter, e.g. the symmetries and their spontaneous breaking. Finally, we discuss the effect of multiple random phase ramping processes and its implication to the stability of As a prototypical example of nonequilibrium quantum the discrete time crystal phase against noisy perturba- matter, the time crystal (TC) phase has allowed new pos- tions. sibilities for the spontaneous symmetry breaking (SSB) Model and method: We consider a 1D hard-core paradigm[4], and have attracted considerable interest in bosonic model with an infinite long-range interaction. its different forms[5–22]. Such a intriguing state, despite The Hamiltonian reads as follows: being proven to be forbidden in equilibrium [23, 24], has been experimentally realized in non-equilibrium settings V (t) i j H = J (ˆb†ˆbi + h.c) ( 1) − nînˆj (1) − i +1 − L − with periodic driving [25, 26]. Its physical observables Xi Xij arXiv:2007.08348v3 [cond-mat.quant-gas] 28 Dec 2020 develop persistent oscillations whose periods are an in- teger multiple of the Hamiltonian period, thus sponta- where J is the single-particle nearest-neighbor(NN) hop- ˆ ˆ neously break the discrete temporal translational symme- ping amplitude,n î = bi†bi is the particle number oper- try(DTTS). In equilibrium systems, the SSB is closely re- ator at site i. V (t) is the strength of the all-to-all in- lated to the elementary excitation: it does not only affect teraction, which is time-dependent but does not decay the dynamics and spectrum of an elementary excitation, with distance. L is the system size of the 1D lattice, 1 but also determines its structure. For example, for a one- and the prefactor L in front of the interaction terms of dimensional (1D) system, a spontaneous discrete spatial Eq.(1) guarantees that the total interacting energy lin- translational (e.g. Z2) symmetry breaking allows certain early scales with the system size. In a bipartite lattice soliton-like excitation: a topological defect sandwiched (e.g, a 1D lattice as in our case), Eq.(1) indicates that by two Z2 symmetry breaking phases[27]. A profound the interaction between a pair of bosons are attractive 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î (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î Ψ(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.

Load more