Springer for Research & Development
Proof of the Absence of Long-Range Temporal Orders in Gibbs States
Haruki Watanabe ORCID: orcid.org/0000-0002-8112-021X1,
Masaki Oshikawa2 &
Tohru Koma3
Journal of Statistical Physics (2020)Cite this article
We address the question whether time translation symmetry can be spontaneously broken in a quantum many-body system. One way of detecting such a symmetry breaking is to examine the time-dependence of a correlation function. If the large-distance behavior of the correlation function exhibits a nontrivial time-dependence in the thermodynamic limit, the system would develop a temporal long-range order, realizing a time crystal. In an earlier publication, we sketched a proof for the absence of such time dependence in the thermal equilibrium described by the Gibbs state (Watanabe and Oshikawa in Phys Rev Lett 114:251603, 2015). Here we present a complete proof and extend the argument to a more general class of stationary states than the Gibbs states.
Time crystals are a newly proposed state of matter that spontaneously breaks the time translation symmetry. The idea of time crystals in the case of the continuous time translation symmetry was first proposed by Wilczek [1], although the validity of the concrete model in this original proposal was soon questioned in Ref. [2]. Then a no-go theorem for a wider but still restricted class of models was presented in Ref. [3]. In a more general setting, the absence of time crystalline orders in the ground state or in the Gibbs state was proven in Ref. [4] without specifying the Hamiltonian but assuming only its locality. These developments triggered further investigation of so-called Floquet time crystals or discrete time crystals in nonequilibrium setting [5,6,7,8,9,10] that break a discrete time translation symmetry into its subgroup. See Refs. [11,12,13] for recent reviews on this topic.
The argument for the no-go theorem at finite temperatures in Ref. [4] was based on the Lieb-Robinson bound [14, 15], which was used to constrain finite time behavior of the correlation function. However, in Ref. [4], Fourier transformation of the correlation function was performed with respect to an infinitely long time, out of the validity of the constraint. This issue was recently pointed out by Ref. [13]. In this work, we present a complete version of the proof without such an issue. Furthermore, we examine the conditions on the density operator to which our argument can be straightforwardly extended. Clarifying these subtleties and settling down the limitations on what the Gibbs state and similar type of stationary states can do should in turn accelerate our exploration of new states that exhibit nontrivial temporal orders.
The Theorem and Its Proof
Setup and Statement
Let us consider a static Hamiltonian \(\hat{H}\) defined on a d-dimensional lattice \(\Lambda \) that is a finite subset of \(\mathbb {Z}^d\). We assume that the Hamiltonian \(\hat{H}\) is written as a sum of local bounded Hamiltonians:
$$\begin{aligned} \hat{H}=\sum _{\mathbf {x}\in \Lambda }\hat{h}_{\mathbf {x}}. \end{aligned}$$
More precisely, we assume that the support of the local Hamiltonian \(\hat{h}_{\mathbf {x}}\) is limited to a finite range \(R_h\) from \(\mathbf {x}\in \Lambda \) and that the operator normFootnote 1 of \(\hat{h}_{\mathbf {x}}\) is bounded by a constant \(N_h\). Both \(R_h\) and \(N_h\) are independent of the position \(\mathbf {x}\in \Lambda \) or the system size \(|\Lambda |\). This setting includes a wide variety of quantum spin systems, fermion systems, and "hard-core" boson systems Footnote 2.
Similarly, we consider observables (not necessarily Hermitian) \(\hat{A}\) and \(\hat{B}\) written as a sum of local observables:
$$\begin{aligned} \hat{A}:=\frac{1}{|\Lambda |}\sum _{\mathbf {x}\in \Lambda } \hat{a}_{\mathbf {x}},\quad \hat{B}:=\frac{1}{|\Lambda |}\sum _{\mathbf {x}\in \Lambda } \hat{b}_{\mathbf {x}}. \end{aligned}$$
The support of \(\hat{a}_{\mathbf {x}}\) and \(\hat{b}_{\mathbf {x}}\) (\(\mathbf {x}\in \Lambda \)) are within a finite range \(R_a\), \(R_b\) from \(\mathbf {x}\) and their operator norm is bounded by constants \(N_a\), \(N_b\), respectively. All of these constants are independent of \(\mathbf {x}\) or \(|\Lambda |\).
We introduce the time evolution of operators for \(t\in \mathbb {R}\) by
$$\begin{aligned} \hat{A}(t):=e^{i\hat{H} t}\hat{A}e^{-i\hat{H} t}. \end{aligned}$$
Our interest is in the time-dependence of the correlation function
$$\begin{aligned} \langle \hat{A}(t)\hat{B}\rangle :=\mathrm{Tr}\left( \hat{A}(t)\hat{B}\hat{\rho }\right) . \end{aligned}$$
Here \(\hat{\rho }\) is the Gibbs state
$$\begin{aligned} \hat{\rho }:=\frac{1}{Z}e^{-\beta \hat{H}} \end{aligned}$$
at the inverse temperature \(\beta \) and \(Z:=\mathrm{Tr}\,e^{-\beta \hat{H}}\) is the partition function. Our claim is that \(\langle \hat{A}(t)\hat{B}\rangle \) is independent of t in the thermodynamic limit \(|\Lambda |\rightarrow +\infty \) [4], i.e.,
$$\begin{aligned} \lim _{|\Lambda |\rightarrow \infty }\left| \langle \hat{A}(t)\hat{B}\rangle -\langle \hat{A}\hat{B}\rangle \right| =0. \end{aligned}$$
Proof for \(\beta >0\)
To prove Eq. (6), it is sufficient to treat the special case \(\hat{B}=\hat{A}^\dagger \):
$$\begin{aligned} \lim _{|\Lambda |\rightarrow \infty }\left| \langle \hat{A}(t)\hat{A}^\dagger \rangle -\langle \hat{A}\hat{A}^\dagger \rangle \right| =0. \end{aligned}$$
This is because \(\langle \hat{A}(t)\hat{B}\rangle \) can be rewritten as
$$\begin{aligned} \langle \hat{A}(t)\hat{B}\rangle&= \frac{1}{2}\Big \langle \big (\hat{A}(t)+\hat{B}^\dagger (t)\big )\big (\hat{A}+\hat{B}^\dagger \big )^\dagger \Big \rangle +\frac{i}{2}\Big \langle \big (\hat{A}(t)+i\hat{B}(t)^\dagger \big )\big (\hat{A}+i\hat{B}^\dagger \big )^\dagger \Big \rangle \nonumber \\&\quad -\frac{1+i}{2}\langle \hat{A}(t)\hat{A}^\dagger \rangle -\frac{1+i}{2}\langle \hat{B}^\dagger (t)\hat{B}\rangle . \end{aligned}$$
Once Eq. (7) is established, it applies to all four correlation functions in the right-hand side and we obtain Eq. (6).
We denote by \(|\Phi _n\rangle \) the eigenstate of the Hamiltonian \(\hat{H}\) with the eigenvalue \(E_n\) (\(n\in \mathbb {N}\)). Using the complete system and writing
$$\begin{aligned} \rho (E_n):=\frac{1}{Z}e^{-\beta E_n}, \end{aligned}$$
we get
$$\begin{aligned} \langle \hat{A}(t)\hat{A}^\dagger \rangle =\sum _{m,n}|\langle \Phi _m|\hat{A}|\Phi _n\rangle |^2\rho (E_m)\,e^{i(E_m-E_n)t}. \end{aligned}$$
We split the summation over m and n into four intervals of \(E_n-E_m\):
$$\begin{aligned} {\left\{ \begin{array}{ll} \text {(i):}&{}2\varepsilon \le E_n-E_m\le K,\\ \text {(ii):}&{}-K\le E_n-E_m\le -2\varepsilon ,\\ \text {(iii):}&{}K<|E_n-E_m|,\\ \text {(iv):}&{}|E_n-E_m|<2\varepsilon , \end{array}\right. } \end{aligned}$$
where \(\varepsilon \) is a small positive number and K is a large positive number. Then the time-dependence of \(\langle \hat{A}(t)\hat{A}^\dagger \rangle \) can be bounded as
$$\begin{aligned}&\left| \langle \hat{A}(t)\hat{A}^\dagger \rangle -\langle \hat{A}\hat{A}^\dagger \rangle \right| \nonumber \\&\quad \le {} 2\sum _{m,n\,:\,2\varepsilon \le E_n-E_m\le K}|\langle \Phi _m|\hat{A}|\Phi _n\rangle |^2\rho (E_m)\nonumber \\&\qquad +2\sum _{m,n\,:\,-K\le E_n-E_m\le -2\varepsilon }|\langle \Phi _m|\hat{A}|\Phi _n\rangle |^2\rho (E_m)\nonumber \\&\qquad +2\sum _{m,n\,:\,K<|E_n-E_m|}|\langle \Phi _m|\hat{A}|\Phi _n\rangle |^2\rho (E_m)\nonumber \\&\qquad +\sum _{m,n\,:\,|E_n-E_m|<2\varepsilon }|\langle \Phi _m|\hat{A}|\Phi _n\rangle |^2\rho (E_m)\left| e^{i(E_m-E_n)t}-1\right| . \end{aligned}$$
In the following, we derive an upper bound for each term in the right hand side one by one. The first two terms will be bounded using the Lieb-Robinson bound and the monotonically decreasing nature of the Boltzmann factor (9). The third term will be evaluated by making use of the large energy difference. Finally, the last term is trivially small because of the time-dependent factor with small energy difference. Plugging these results [Eqs. (30), (33), (38), and (39) below] into the right-hand side of Eq. (12), we get
$$\begin{aligned} \left| \langle \hat{A}(t)\hat{A}^\dagger \rangle -\langle \hat{A}\hat{A}^\dagger \rangle \right|&\le 2\varepsilon +2\varepsilon +\frac{2C}{K^2}+2N_a^2\varepsilon |t|, \end{aligned}$$
where C is a positive constant independent of the system size. Since we can take \(\varepsilon \) to be small and K to be large by choosing a sufficiently large system size \(|\Lambda |\), we obtain the desired result.
The Range (i): \(2\varepsilon \le E_n-E_m\le K\)
Let us start with the contribution from the range \(2\varepsilon \le E_n-E_m\le K\). To this end, we introduce a cutoff function \(\eta ^{+}\in C_0^\infty (\mathbb {R})\) (i.e., an infinitely differentiable function with a compact support ) that satisfies the following conditions: Footnote 3
$$\begin{aligned} {\left\{ \begin{array}{ll} \eta ^{+}(\omega )=1 &{} (2\varepsilon \le \omega \le K),\\ \eta ^{+}(\omega )=0 &{} (\omega \le \varepsilon \text { or }K+\varepsilon \le \omega ),\\ 0\le \eta ^{+}(\omega )\le 1&{} (\text {otherwise}). \end{array}\right. } \end{aligned}$$
The Fourier transform of \(\eta ^+(\omega )\) is given by
$$\begin{aligned} \tilde{\eta }^+(t):=\frac{1}{2\pi }\int _{-\infty }^{+\infty }d\omega \, e^{i\omega t}\eta ^+(\omega ), \end{aligned}$$
which decays faster than any power of t. This can be shown by performing an integration by parts repeatedly:
$$\begin{aligned} \tilde{\eta }^+(t)=\frac{1}{2\pi }\left( \frac{i}{t}\right) ^\ell \int _{-\infty }^{+\infty }d\omega \, e^{i\omega t}\frac{\partial ^\ell }{\partial \omega ^\ell }\eta ^+(\omega )\quad \text { for }t\ne 0, \end{aligned}$$
which implies, for any integer \(\ell \in \mathbb {N}\), that
$$\begin{aligned}&|\tilde{\eta }^+(t)|\le \mathcal {C}_\ell |t|^{-\ell }, \end{aligned}$$
$$\begin{aligned}&\mathcal {C}_\ell :=\frac{1}{2\pi }\int _{-\infty }^{+\infty }d\omega \, \left| \frac{\partial ^\ell }{\partial \omega ^\ell }\eta ^+(\omega )\right| . \end{aligned}$$
We consider a correlation function
$$\begin{aligned} g(t):=\langle [\hat{A}(t),\hat{A}^\dagger ]\rangle =\sum _{m,n}|\langle \Phi _m|\hat{A}|\Phi _n\rangle |^2\left( \rho (E_m)-\rho (E_n)\right) e^{i(E_m-E_n)t}. \end{aligned}$$
On one hand, we have
$$\begin{aligned} \int _{-\infty }^{+\infty } dt\, g(t)\tilde{\eta }^+(t)&=\sum _{m,n}|\langle \Phi _m|\hat{A}|\Phi _n\rangle |^2 \left( \rho (E_m)-\rho (E_n)\right) \int _{-\infty }^{+\infty } dt\,e^{-i(E_n-E_m)t}\tilde{\eta }^+(t)\nonumber \\&=\sum _{m,n}|\langle \Phi _m|\hat{A}|\Phi _n\rangle |^2\left( \rho (E_m)-\rho (E_n)\right) \eta ^+(E_n-E_m)\nonumber \\&\ge \sum _{m,n\,:\,2\varepsilon \le E_n-E_m\le K}|\langle \Phi _m|\hat{A}|\Phi _n\rangle |^2\left( \rho (E_m)-\rho (E_n)\right) \nonumber \\&\ge \sum _{m,n\,:\,2\varepsilon \le E_n-E_m\le K}|\langle \Phi _m|\hat{A}|\Phi _n\rangle |^2\rho (E_m) \Delta . \end{aligned}$$
In passing to the third line, we used \(\rho (E_n)< \rho (E_m)\) when \(E_n> E_m\) and the conditions (14) of \(\eta ^+(\omega )\). In the last line, we defined
$$\begin{aligned} \Delta :=1-\max _{m,n\,:\,2\varepsilon \le E_n-E_m \le K}\frac{\rho (E_n)}{\rho (E_m)}. \end{aligned}$$
For the Gibbs state (9), we have
$$\begin{aligned}&\Delta \ge h(\varepsilon )>0\quad \text {for }\varepsilon >0, \end{aligned}$$
$$\begin{aligned}&h(\varepsilon ):=1-e^{-2\beta \varepsilon }. \end{aligned}$$
On the other hand, we can decompose the integral into two parts as
$$\begin{aligned} \int _{-\infty }^{+\infty } dt\, g(t)\tilde{\eta }^+(t) =\int _{|t|\ge T}dt\, g(t)\tilde{\eta }^+(t)+\int _{-T}^{+T}dt\, g(t)\tilde{\eta }^+(t) \end{aligned}$$
where T is a large positive number. For the first integral in the right-hand side, we use the property Eq. (17) of the function \(\tilde{\eta }^+(t)\) as well as the trivial bound \(|g(t)|\le 2N_a^2\).Footnote 4 For a given function \(\tilde{\eta }^+(t)\) with the parameters \(\varepsilon \) and K, we can find a large T such that
$$\begin{aligned} \left| \int _{|t|\ge T}dt\, g(t)\tilde{\eta }^+(t)\right| \le \varepsilon h(\varepsilon ). \end{aligned}$$
For the second integral, we can use the Lieb-Robinson bound [14, 15], from which we have [4]
$$\begin{aligned} \Big \Vert \big [\hat{A}(t),\hat{A}^\dagger \big ]\Big \Vert \le \frac{C_1+C_2|t|^d}{|\Lambda |} \end{aligned}$$
for system-size-independent constants \(C_1\) and \(C_2\). Thus
$$\begin{aligned} \left| \int _{-T}^{+T}dt\, g(t)\tilde{\eta }^+(t)\right|&\le \int _{-T}^{+T}dt\, \big |\langle [\hat{A}(t),\hat{A}^\dagger ]\rangle \big |\left| \tilde{\eta }^+(t)\right| \nonumber \\&\le \int _{-T}^{+T}dt\, \frac{C_1+C_2|t|^d}{|\Lambda |}\frac{K}{2\pi }=K\frac{(d+1)C_1T+C_2T^{d+1}}{\pi (d+1)|\Lambda |}, \end{aligned}$$
where in the second inequality we used
$$\begin{aligned} \left| \tilde{\eta }^+(t)\right| \le \frac{1}{2\pi }\int _{-\infty }^{+\infty }d\omega \, \eta ^+(\omega )\le \frac{1}{2\pi }\int _{\varepsilon }^{K+\varepsilon }d\omega \, 1=\frac{K}{2\pi }. \end{aligned}$$
Therefore, for any given large K and T, there exists a large volume \(|\Lambda |\) such that
$$\begin{aligned} \left| \int _{-T}^{+T}dt\, g(t)\tilde{\eta }^+(t)\right| \le \varepsilon h(\varepsilon ). \end{aligned}$$
Combining Eqs. (25) and (29) with the bound (20), we get
$$\begin{aligned} \sum _{m,n\,:\,2\varepsilon \le E_n-E_m\le K} |\langle \Phi _m|\hat{A}|\Phi _n\rangle |^2\rho (E_m)\le \frac{2\varepsilon h(\varepsilon )}{\Delta }\le 2\varepsilon . \end{aligned}$$
The Range (ii): \(-K\le E_n-E_m\le -2\varepsilon \)
Similarly, to estimate the contribution from the range \(-K\le E_n-E_m\le -2\varepsilon \), we introduce a cutoff function \(\eta ^{-}\in C_0^\infty (\mathbb {R})\) that satisfies the following conditions:
$$\begin{aligned} {\left\{ \begin{array}{ll} \eta ^{-}(\omega )=-1 &{} (-K\le \omega \le -2\varepsilon ),\\ \eta ^{-}(\omega )=0 &{} (\omega \le -K-\varepsilon \text { or }-\varepsilon \le \omega ),\\ -1\le \eta ^{-}(\omega )\le 0&{} (\text {otherwise}). \end{array}\right. } \end{aligned}$$
This time we have
$$\begin{aligned} \int _{-\infty }^{+\infty } dt\, g(t)\tilde{\eta }^-(t)&=\sum _{m,n}|\langle \Phi _m|\hat{A}|\Phi _n\rangle |^2 \left( \rho (E_m)-\rho (E_n)\right) \int _{-\infty }^{+\infty } dt\,e^{-i(E_n-E_m)t}\tilde{\eta }^-(t)\nonumber \\&=\sum _{m,n}|\langle \Phi _m|\hat{A}|\Phi _n\rangle |^2\left( \rho (E_n)-\rho (E_m)\right) (-\eta ^-(E_n-E_m))\nonumber \\&\ge \sum _{m,n\,:\,-K\le E_n-E_m\le -2\varepsilon }|\langle \Phi _m|\hat{A}|\Phi _n\rangle |^2\left( \rho (E_n)-\rho (E_m)\right) \nonumber \\&\ge \sum _{m,n\,:\,-K\le E_n-E_m\le -2\varepsilon }|\langle \Phi _m|\hat{A}|\Phi _n\rangle |^2\rho (E_n) \Delta \nonumber \\&\ge \sum _{m,n\,:\,-K\le E_n-E_m\le -2\varepsilon }|\langle \Phi _m|\hat{A}|\Phi _n\rangle |^2\rho (E_m) \Delta . \end{aligned}$$
Here, \(\Delta \) is defined in Eq. (21). In the same way as before, we find
$$\begin{aligned} \sum _{m,n\,:\,-K\le E_n-E_m\le -2\varepsilon } |\langle \Phi _m|\hat{A}|\Phi _n\rangle |^2\rho (E_m)\le 2\varepsilon . \end{aligned}$$
The Range (iii): \(K<|E_n-E_m|\)
The third contribution can be easily bounded by using a trick.
$$\begin{aligned} \sum _{m,n\,:\,K<|E_n-E_m|}|\langle \Phi _m|\hat{A}|\Phi _n\rangle |^2\rho (E_m)&\le \frac{1}{K^2}\sum _{m,n}(E_n-E_m)^2|\langle \Phi _m|\hat{A}|\Phi _n\rangle |^2\rho (E_m)\nonumber \\&= \frac{1}{K^2} \langle [\hat{A},\hat{H}][\hat{H},\hat{A}^\dagger ]\rangle \nonumber \\&\le \frac{1}{K^2} \big \Vert [\hat{A},\hat{H}]\big \Vert ^2. \end{aligned}$$
Thanks to the assumed locality of the Hamiltonian, this operator norm can be bounded as
$$\begin{aligned}&\big \Vert [\hat{A},\hat{H}]\big \Vert =\frac{1}{|\Lambda |}\sum _{\mathbf {x},\mathbf {y}\in \Lambda \,:\,|\mathbf {x}-\mathbf {y}|\le R_h+R_a}\big \Vert [\hat{a}_{\mathbf {x}},\hat{h}_{\mathbf {y}}]\big \Vert \le C, \end{aligned}$$
$$\begin{aligned}&C:= 2N_aN_hv(R_a+R_h), \end{aligned}$$
$$\begin{aligned}&v(R):=\frac{1}{|\Lambda |}\sum _{\mathbf {x},\mathbf {y}\in \Lambda \,:\,|\mathbf {x}-\mathbf {y}|\le R}1. \end{aligned}$$
Note that v(R) does not grow with \(|\Lambda |\). Therefore,
$$\begin{aligned} \sum _{m,n\,:\,K<|E_n-E_m|}|\langle \Phi _m|\hat{A}|\Phi _n\rangle |^2\rho (E_m)\le \frac{C}{K^2}. \end{aligned}$$
The Range (iv): \(|E_n-E_m|<\varepsilon \)
Finally, using the fact that \(|e^{ix}-1|=2|\sin \frac{x}{2}|\le |x|\) for any real number x, we get
$$\begin{aligned}&\sum _{m,n\,:\,|E_n-E_m|<2\varepsilon }|\langle \Phi _m|\hat{A}|\Phi _n\rangle |^2\rho (E_m)\left| e^{i(E_m-E_n)t}-1\right| \nonumber \\&\quad \le 2\varepsilon |t|\sum _{m,n:\,|E_n-E_m|<2\varepsilon }|\langle \Phi _m|\hat{A}|\Phi _n\rangle |^2\rho (E_m)\nonumber \\&\quad \le 2\varepsilon |t|\sum _{m,n}|\langle \Phi _m|\hat{A}|\Phi _n\rangle |^2\rho (E_m)=2\varepsilon |t| \langle \hat{A}\hat{A}^\dagger \rangle \nonumber \\&\quad \le 2\Vert \hat{A}\Vert ^2\varepsilon |t| \le 2N_a^2\varepsilon |t|. \end{aligned}$$
This completes the verification of Eq. (13) and hence the proof of Eq. (7).
Proof for \(\beta =0\)
Interestingly, the proof in the previous section does not apply to the infinite temperature (\(\beta =0\)) where \(\rho (E_n)\) in Eq. (9) becomes constant:
$$\begin{aligned} \rho (E_n)=\frac{1}{\mathcal {D}}. \end{aligned}$$
Here \(\mathcal {D}\) is the dimension of the entire Hilbert space. In this special case, however, we can directly prove Eq. (6) using the clustering property of the infinite-temperature state.Footnote 5 Thus the "absence of the time crystals" also holds at the infinite temperature, consistently with the intuition that the infinite temperature is the most disordered limit.
At \(\beta =0\), the equal-time correlation function trivially exhibits the locality
$$\begin{aligned} \langle \hat{a}_{\mathbf {x}} \hat{b}_{\mathbf {y}} \rangle = \langle \hat{a}_{\mathbf {x}} \rangle \langle \hat{b}_{\mathbf {y}} \rangle \end{aligned}$$
if the support of \(\hat{a}_{\mathbf {x}}\) and \(\hat{b}_{\mathbf {y}}\) do not overlap. This implies the clustering property and the absence of any spatial long-range order. However, quantum dynamics is nontrivial even at the infinite temperature (see, for example, Ref. [16] and references therein) and the question of the time crystal is not totally trivial.
To prove Eq. (7) in this setting, let us define
$$\begin{aligned} \delta \hat{A}:=\hat{A}-\langle \hat{A}\rangle =\frac{1}{|\Lambda |}\sum _{\mathbf {x}\in \Lambda } (\hat{a}_{\mathbf {x}}-\langle \hat{a}_{\mathbf {x}}\rangle ). \end{aligned}$$
It follows that
$$\begin{aligned} \Big |\langle \delta \hat{A}(t)\,\delta \hat{A}^\dagger \rangle \Big |&=\Big |\sum _{n,m}|\langle \Phi _n|\delta \hat{A}|\Phi _m\rangle |^2 \rho (E_n)e^{i(E_n-E_m)t}\Big |\nonumber \\&\le \sum _{n,m}|\langle \Phi _n|\delta \hat{A}|\Phi _m\rangle |^2\rho (E_n)=\langle \delta \hat{A}\,\delta \hat{A}^\dagger \rangle \end{aligned}$$
$$\begin{aligned} \Big |\langle \hat{A}(t)\hat{A}^\dagger \rangle -\langle \hat{A}\,\hat{A}^\dagger \rangle \Big |=\Big |\langle \delta \hat{A}(t)\,\delta \hat{A}^\dagger \rangle -\langle \delta \hat{A}\,\delta \hat{A}^\dagger \rangle \Big |\le 2 \langle \delta \hat{A}\,\delta \hat{A}^\dagger \rangle . \end{aligned}$$
The locality in Eq. (41) implies that \(\langle \delta \hat{A}\,\delta \hat{A}^\dagger \rangle \) is inversely proportional to \(|\Lambda |\). Therefore, we obtain Eq. (7) which gives Eq. (6) as explained in Sec. 2.2.1.
The proof for the Gibbs state at a finite temperature in Sec. 2.2 equally applies to other stationary states (\([\hat{\rho },\hat{H}]=0\)) of a local, static Hamiltonian as long as the density operator \(\hat{\rho }\) satisfies the following conditions.
The weight \(\rho (E_n)\) for each eigenstate \(|\Phi _n\rangle \) should be a strictly decreasing function of \(E_n\), i.e., \(\rho (E_n)<\rho (E_m)\) when \(E_n>E_m\) for any n and m.
There exists a smooth function \(h(\varepsilon )\) of \(\varepsilon \) such that the quantity \(\Delta \) defined in Eq. (21) is bounded below as in Eq. (22). Furthermore, \(h(\varepsilon )\) must be independent of the system size \(|\Lambda |\).
This argument can also be straightforwardly modified when the weight \(\rho (E_n)\) is a strictly increasing function of \(E_n\). Therefore, the same theorem holds for Gibbs state with a "negative temperature" (\(\beta <0\)), which is well-defined for a bounded Hamiltonian we discuss here.
Recently, in the interesting paper [17], Huang showed the inequality
$$\begin{aligned} \left| \langle \hat{A}(t)\hat{B}\rangle -\langle \hat{A}\hat{B}\rangle \right| = O(|\Lambda |^{-1}) \end{aligned}$$
for an arbitrary stationary state \(\hat{\rho }\) assuming a sufficiently fast decay of spatial correlation functions but without assuming the locality of the Hamiltonian. Here we comment that, in order to prove the weaker statement Eq. (6), which already implies the absence of long-range temporal order, it is sufficient to assume the "clustering" property of spatial correlation functions. We say a state exhibits clustering if there exists a system-size independent function f(r) with \(\lim _{r\rightarrow \infty }f(r)=0\) such that correlation functions of any local bounded operator obey
$$\begin{aligned} |\langle \delta \hat{a}_{\mathbf {x}}\,\delta \hat{a}_{\mathbf {y}}^\dagger \rangle |\le f(|\mathbf {x}-\mathbf {y}|)\quad \text {for all}\quad \mathbf {x}, \mathbf {y}\in \Lambda . \end{aligned}$$
As we discuss in Appendix A, Eq. (46) readily gives
$$\begin{aligned} \lim _{|\Lambda |\rightarrow \infty }\langle \delta \hat{A}\,\delta \hat{A}^\dagger \rangle =0. \end{aligned}$$
Moreover, Eq. (43) holds as long as \(\rho (E_n)\ge 0\) for all n even when \(\rho (E_n)\) is not constant [17]. Thus, we get Eq. (7) by combining Eqs. (44) and (47).
In conclusion, in order to realize a temporal long-range order in a stationary state, the system has to fulfill both of the following two conditions:
Either the weight \(\rho (E_n)\) breaks some of the above conditions or the Hamiltonian \(\hat{H}\) violates the locality.
The state \(\hat{\rho }\) does not possess the clustering property.
For example, a time-crystalline behavior may be observed in a single eigenstate (except for the ground state) or a micro-canonical ensemble of a local Hamiltonian (e.g. Ref. [18]), and in the ground state of a non-local Hamiltonian (e.g. Ref. [19]).
After we revised our manuscript and updated the version on arXiv, Huang also posted the second version of his paper (Y. Huang, arXiv:1912.01210v2) in which the author made a revision for the main result in the first version [17], and added a new result which is similar to ours in Sec. 3.
The operator norm of an operator \(\hat{O}\) is defined as \({\Vert }\hat{O}{\Vert }:=\text {sup}_{|\psi \rangle , {\Vert }{|\psi \rangle }|{\Vert }>0}{\Vert }\hat{O}{|\psi \rangle }{\Vert }/{\Vert }{|\psi \rangle }{\Vert }\).
For the assumed boundedness, the maximum number of bosons that can occupy a single site must be a finite number independent of the system size.
An example of \(\eta ^{+}(\omega )\) for the range \(\varepsilon \le \omega \le 2\varepsilon \) and \(K\le \omega \le K+\varepsilon \) can be constructed using \(m(x):=\int _{-1}^{x}dy\,e^{-\frac{1}{1-y^2}}\) (\(-1\le x\le +1\)). For example, one can set \(\eta ^{+}(\omega )=m(\frac{2\omega -3\varepsilon }{\varepsilon })/m(+1)\) for \(\varepsilon \le \omega \le 2\varepsilon \).
Here and hereafter, we use the standard properties of the operator norm, such as \({|}{\langle }{\hat{O}}{\rangle }{|}\le {\Vert }{\hat{O}}{\Vert }\), \({\Vert }{\hat{O}}{\Vert }={\Vert }{\hat{O}}^\dagger {\Vert }\), and \({\Vert }{\hat{O}\hat{O}'}{\Vert }\le {\Vert }{\hat{O}}{\Vert }{\Vert }{\hat{O}'}{\Vert }\) for operators \(\hat{O}\) and \(\hat{O}'\).
In the earlier version of the manuscript [20], our proof for \(\beta =0\) was based on the Lieb-Robinson bound. The present simpler proof was informed by Yichen Huang. See also Ref. [17] and Sect. 3.
Wilczek, F.: Quantum time crystals. Phys. Rev. Lett. 109, 160401 (2012)
Bruno, P.: Comment on "quantum time crystals". Phys. Rev. Lett. 110, 118901 (2013)
Bruno, P.: Impossibility of spontaneously rotating time crystals: a no-go theorem. Phys. Rev. Lett. 111, 070402 (2013)
Watanabe, H., Oshikawa, M.: Absence of quantum time crystals. Phys. Rev. Lett. 114, 251603 (2015)
Sacha, K.: Modeling spontaneous breaking of time-translation symmetry. Phys. Rev. A 91, 033617 (2015)
Khemani, V., Lazarides, A., Moessner, R., Sondhi, S.L.: Phase structure of driven quantum systems. Phys. Rev. Lett. 116, 250401 (2016)
Else, D.V., Bauer, B., Nayak, C.: Floquet time crystals. Phys. Rev. Lett. 117, 090402 (2016)
Yao, N.Y., Potter, A.C., Potirniche, I.-D., Vishwanath, A.: Discrete time crystals: rigidity, criticality, and realizations. Phys. Rev. Lett. 118, 030401 (2017)
Choi, S., Choi, J., Landig, R., Kucsko, G., Zhou, Hengyun, Isoya, Junichi, Jelezko, Fedor, Onoda, Shinobu, Sumiya, Hitoshi, Khemani, Vedika, von Keyserlingk, Curt, Yao, Norman Y., Demler, Eugene, Lukin, Mikhail D.: Observation of discrete time-crystalline order in a disordered dipolar many-body system. Nature 543, 221–225 (2017)
Zhang, J., Hess, P.W., Kyprianidis, A., Becker, P., Lee, A., Smith, J., Pagano, G., Potirniche, I.D., Potter, A.C., Vishwanath, A., Yao, N.Y., Monroe, C.: Observation of a discrete time crystal. Nature 543, 217–220 (2017)
Sacha, K., Zakrzewski, J.: Time crystals: a review. Rep. Prog. Phys. 81, 016401 (2017)
Else, D.V., Monroe, C., Nayak, C., Yao, N.Y.: Discrete time crystals. arXiv:1905.13232
Khemani, V., Moessner, R., Sondhi, S.L.: A brief history of time crystals. arXiv:1910.10745
Lieb, E.H., Robinson, D.W.: The finite group velocity of quantum spin systems. Commun. Math. Phys. 28, 251–257 (1972)
Hastings, M.B.: Locality in quantum systems. arXiv:1008.5137
Žnidarič, M.: Spin transport in a one-dimensional anisotropic Heisenberg model. Phys. Rev. Lett. 106, 220601 (2011)
Huang, Y.: Absence of temporal order in states with fast decay of spatial correlations. arXiv:1912.01210v1
Syrwid, A., Zakrzewski, J., Sacha, K.: Time crystal behavior of excited eigenstates. Phys. Rev. Lett. 119, 250602 (2017)
Kozin, V.K., Kyriienko, O.: Quantum time crystals from Hamiltonians with long-range interactions. Phys. Rev. Lett. 123, 210602 (2019)
Watanabe, H., Oshikawa, M., Koma, T.: Absence of temporal order in states with fast decay of spatial correlations. arXiv:1911.12939v1
We thank Hal Tasaki for fruitful discussions. We also thank Shivaji Sondhi for encouraging us to publish the present result, which is based on an earlier unpublished note, and Yichen Huang for informing us of a simplification of the proof for \(\beta =0\). The work of H. W. was supported by JSPS KAKENHI Grant No. JP17K17678 and by JST PRESTO Grant No. JPMJPR18LA. The work of M. O. was supported in part by JSPS KAKENHI Grant No. JP19H01808 and by US National Science Foundation Grant No. NSF PHY-1748958 through Kavli Institute for Theoretical Physics, UC Santa Barbara.
Department of Applied Physics, University of Tokyo, Tokyo, 113-8656, Japan
Haruki Watanabe
Institute for Solid State Physics, University of Tokyo, Kashiwa, 277-8581, Japan
Masaki Oshikawa
Department of Physics, Gakushuin University, Mejiro, Toshima-ku, Tokyo, 171-8588, Japan
Tohru Koma
Search for Haruki Watanabe in:
Search for Masaki Oshikawa in:
Search for Tohru Koma in:
Correspondence to Haruki Watanabe.
Communicated by Hal Tasaki.
Appendix A: Fluctuation and Clustering
In this appendix, we show that fluctuations of normalized macroscopic observables are negligible in the large volume limit when the state possesses the clustering property. The clustering defined in Eq. (46) means that, for any \(\varepsilon >0\), there exists \(R>0\) (independent of \(|\Lambda |\)) such that
$$\begin{aligned} |\langle \delta \hat{a}_{\mathbf {x}}\,\delta \hat{a}_{\mathbf {y}}^\dagger \rangle |<\frac{1}{2}\varepsilon \quad \text {if}\quad |\mathbf {x}-\mathbf {y}|>R. \end{aligned}$$
(A1)
$$\begin{aligned} \langle \delta \hat{A}\,\delta \hat{A}^\dagger \rangle&\le \frac{1}{|\Lambda |^2}\sum _{\mathbf {x},\mathbf {y}\in \Lambda \,:\,|\mathbf {x}-\mathbf {y}|\le R}| \langle \delta \hat{a}_{\mathbf {x}}\,\delta \hat{a}_{\mathbf {y}}^\dagger \rangle |+\frac{1}{|\Lambda |^2}\sum _{\mathbf {x},\mathbf {y}\in \Lambda \,:\,|\mathbf {x}-\mathbf {y}|> R}| \langle \delta \hat{a}_{\mathbf {x}}\,\delta \hat{a}_{\mathbf {y}}^\dagger \rangle |\nonumber \\&\quad \le \frac{2N_a^2v(R)}{|\Lambda |}+\frac{1}{2}\varepsilon , \end{aligned}$$
where v(R) is defined in Eq. (37). Since \(N_a\) is independent of the system size, we can find \(|\Lambda |\) such that \(\frac{2N_a^2v(R)}{|\Lambda |}<\frac{1}{2}\varepsilon \). Therefore, for a sufficiently large system size, we have
$$\begin{aligned} \langle \delta \hat{A}\,\delta \hat{A}^\dagger \rangle <\varepsilon . \end{aligned}$$
This completes the proof of Eq. (47).
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Watanabe, H., Oshikawa, M. & Koma, T. Proof of the Absence of Long-Range Temporal Orders in Gibbs States. J Stat Phys (2020). https://doi.org/10.1007/s10955-019-02471-5
Accepted: 13 December 2019
DOI: https://doi.org/10.1007/s10955-019-02471-5
Not logged in - 3.214.184.250