Time crystals in the driven transverse field Ising model under quasiperiodic modulation

We investigate the transverse field Ising model subject to a two-step periodic driving protocol and quasiperiodic modulation of the Ising couplings. Analytical results on the phase boundaries associated with Majorana edge modes and numerical results on the localization of single-particle excitations are presented. The implication of a region with fully localized domain-wall-like excitations in the parameter space is eigenstate order and exact spectral pairing of Floquet eigenstates, based on which we conclude the existence of time crystals. We also examine various correlation functions of the time crystal phase numerically, in support of its existence.