Probing Quantum Thermalization of a Disordered Dipolar Spin Ensemble with Discrete Time-Crystalline Order

Joonhee Choi, Hengyun Zhou, Soonwon Choi, Renate Landig, Wen Wei Ho, Junichi Isoya, Fedor Jelezko, Shinobu Onoda, Hitoshi Sumiya, Dmitry A. Abanin, and Mikhail D. Lukin

Phys. Rev. Lett. 122, 043603 – Published 1 February 2019

Abstract

We investigate thermalization dynamics of a driven dipolar many-body quantum system through the stability of discrete time crystalline order. Using periodic driving of electronic spin impurities in diamond, we realize different types of interactions between spins and demonstrate experimentally that the interplay of disorder, driving, and interactions leads to several qualitatively distinct regimes of thermalization. For short driving periods, the observed dynamics are well described by an effective Hamiltonian which sensitively depends on interaction details. For long driving periods, the system becomes susceptible to energy exchange with the driving field and eventually enters a universal thermalizing regime, where the dynamics can be described by interaction-induced dephasing of individual spins. Our analysis reveals important differences between thermalization of long-range Ising and other dipolar spin models.

Received 26 June 2018
Revised 19 November 2018

DOI:https://doi.org/10.1103/PhysRevLett.122.043603

Physics Subject Headings (PhySH)

Long-range interactions Nonequilibrium statistical mechanics Quantum phase transitions Spin dynamics Spin relaxation Spontaneous symmetry breaking

Floquet systems Nitrogen vacancy centers in diamond

Atomic, Molecular & OpticalCondensed Matter, Materials & Applied PhysicsStatistical Physics & ThermodynamicsQuantum Information, Science & Technology

Authors & Affiliations

Joonhee Choi^1,2,*, Hengyun Zhou^1,*, Soonwon Choi¹, Renate Landig¹, Wen Wei Ho¹, Junichi Isoya³, Fedor Jelezko⁴, Shinobu Onoda⁵, Hitoshi Sumiya⁶, Dmitry A. Abanin⁷, and Mikhail D. Lukin^1,†

¹Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA
²School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, USA
³Faculty of Pure and Applied Sciences, University of Tsukuba, Tsukuba, Ibaraki 305-8573 Japan
⁴Institut für Quantenoptik, Universität Ulm, 89081 Ulm, Germany
⁵Takasaki Advanced Radiation Research Institute, National Institutes for Quantum and Radiological Science and Technology, 1233 Watanuki, Takasaki, Gunma 370-1292, Japan
⁶Sumitomo Electric Industries Ltd., Itami, Hyougo, 664-0016, Japan
⁷Department of Theoretical Physics, University of Geneva, 1211 Geneva, Switzerland

^*These authors contributed equally to this work
^†lukin@physics.harvard.edu

Article Text (Subscription Required)

Click to Expand

Supplemental Material (Subscription Required)

Click to Expand

References (Subscription Required)

Click to Expand

Issue

Vol. 122, Iss. 4 — 1 February 2019

Reuse & Permissions

Access Options

Article Available via CHORUS

Download Accepted Manuscript

Author publication services for translation and copyediting assistance advertisement

Images

Figure 1
Experimental system and observation of DTC order. (a) Periodically driven, interacting NV centers. During each Floquet period $T$ (dotted box), NV centers interact for time $τ_{1}$ , then experience pulsed microwave rotations of duration $τ_{2}$ , at resonant frequencies $ω_{1}$ or $ω_{2}$ . After $n$ Floquet cycles, the population difference between $| 0 ⟩$ and $| - 1 ⟩$ is measured. (b) Distinct Floquet time evolutions realized: $Z_{2}$ Ising, with only Ising interactions, and $Z_{2}$ and $Z_{3}$ , with both Ising and spin-exchange interactions. (c)–(f) Representative time traces of the normalized spin polarization $P (n T)$ and Fourier spectra $| S (ν) |^{2}$ of the $Z_{3}$ DTC order at a perturbation $ε / π = 0.06$ with $T = 70 ns$ (c),(d) and 130 ns (e),(f). In (c),(e), blue, gray, and green points correspond to $P (t)$ at $t \equiv T, 2 T, 3 T$ (mod $3 T$ ), respectively.
Reuse & Permissions
Figure 2
Phase diagram of DTC order for short drive periods. (a)–(c) Phase diagrams in semi-log scale. The phase boundary (markers) is identified as a crystalline fraction of 10%. The dotted line indicates the linear phase boundary predicted by a self-consistent mean-field analysis [52]. Shaded areas denote a universal dephasing regime corresponding to Markovian thermalization where $T > T^{*}$ (see Fig. 4). In (a), the dashed line represents the theoretical prediction from Ref. [25]. Errorbars denote 95% confidence intervals of the phase boundary [52]. (d) Short- $T$ phase diagram in linear scale [markers as in (a)–(c)]. (e) Bloch sphere illustrating the screening effect of spin-exchange interactions. $h_{z}$ and $h_{y}$ are mean fields arising from Ising and spin-exchange interactions, respectively, and $ε / T$ is the perturbing field. The black arrow corresponds to the mean field solution $| ψ_{MF} ⟩$ .
Reuse & Permissions
Figure 3
Long interaction time behavior of DTC order. (a),(b) Temporal decay of $Z_{3}$ DTC order $| S (ν = 1 / 3) |^{2}$ , for different $ε$ , as a function of sweeping window position $n_{sweep}$ . Dashed lines denote the noise floor. (c) Late-time decay rate $Γ$ as a function of $ε$ , with phenomenological quadratic fit. Each data point results from an average over simple exponential fits of $| S (ν = 1 / 3) |^{2}$ starting from $n_{sweep} = 15 - 20$ to exclude initial transients resulting from a short-time, non-universal dephasing. Error bars denote the statistical error of the fit results. The arrow indicates the mean-field phase boundary.
Reuse & Permissions
Figure 4
Universal thermalization dynamics for different Floquet Hamiltonians. (a) Exponents of the stretched exponential fits versus $T$ . Data points denote the average over $β$ values extracted at different $ε$ , error bars are the standard deviation of the mean. Lines denote fits to extract the saturation timescale $T^{*}$ (arrows), identified where $\bar{β} = 0.9$ . (b) Late-time decay rate $Γ$ as a function of $ε$ (after a global offset $Γ_{0}$ has been subtracted, see [52]) at $T = 3.5$ , 3.5, 2.3, and $0.3 μ s$ for the $Z_{2}$ -Ising, $Z_{2}$ , $Z_{3}$ , and two-group $Z_{3}$ cases, respectively [markers as in (a)]. Solid line indicates a dephasing model fit, predicting $Γ = ε^{2} / 2 + Γ_{0}$ . Error bars as in Fig. 3. Cross markers show numerical results of an infinite-range, interacting spin system with both Ising and spin-exchange interactions [52]. Error bars in the numerics represent the standard deviation of decay rate distributions for different realizations.
Reuse & Permissions

Physical Review Letters