- 1State Key Laboratory of Low Dimensional Quantum Physics, Department of Physics, Tsinghua University, Beijing, China
- 2Frontier Science Center for Quantum Information, Beijing, China
Understanding phases of matter is of both fundamental and practical importance. Prior to the widespread appreciation and acceptance of topological order, the paradigm of spontaneous symmetry breaking, formulated along the Landau–Ginzburg–Wilson (LGW) dogma, is central to understanding phases associated with order parameters of distinct symmetries and transitions between phases. This work proposes to identify ground-state phases of the quantum many-body system in terms of time order, which is operationally defined by the appearance of the non-trivial temporal structure in the two-time auto-correlation function of a symmetry operator (order parameter) while the system approaches thermodynamic limit. As a special case, the (symmetry protected) time crystalline order phase detects continuous time crystal (CTC). We originally discover the physical meaning of CTC’s characteristic period and amplitude. Time order phase diagrams for spin-1 atomic Bose–Einstein condensate (BEC) and quantum Rabi model are fully worked out. In addition to time-crystalline order, the intriguing phase of time-functional order is discussed in two non-Hermitian interacting spin models.
1 Introduction
A consistent theme for studying the many-body system, particularly in condensed matter physics, concerns the classification of phases and their associated phase transitions [20, 52, 68]. In the celebrated Landau–Ginzburg–Wilson (LGW) paradigm [35, 70], spontaneous symmetry breaking plays a central role with order parameters characterizing different phases of matter possessing respective broken symmetries. Other schemes for classifying phases as well as their associated transitions are, however, beyond the Landau–Ginzburg–Wilson paradigm, which are by now well accepted since first established decades ago [53, 63, 64]. For example, topological order, which classifies the gapped quantum many-body system, constitutes a topical research direction [63, 64, 66, 67]. Our current understanding categorizes gapped systems into gapped liquid phases [74] and gapped non-liquid phases, with the former broadly including phases of topological order [63, 64], symmetry-enriched topological order [9, 12, 25, 65], and symmetry-protected trivial order [10, 11, 23], while the recently discussed fracton phases [55, 56, 60] belong to the latter of gapped non-liquid phases.
Temporal properties of phases are also worthy of investigations as exemplified by many recent studies [41, 50, 69]. For instance, time crystal (TC) or perpetual temporal dependence in a many-body ground state that breaks spontaneously time translation symmetry (TTS) constitutes an exciting new phenomenon. First proposed by Wilczek [69] for quantum systems and followed by Shapere and Wilczek [54] for classical systems in 2012, TC in their original sense is unfortunately ruled out by Bruno’s no-go theorem the following year [3, 42]. Watanabe and Oshikawa (WO) reformulate the idea of quantum TC, and present a refined no-go theorem for many-body systems without too long-range interactions [62]. Continued efforts are directed at searching for continuous time crystal in open systems [4, 5, 30] and classical driven-diffusive systems [28]. Most recent efforts on this topic are directed toward non-equilibrium discrete/Floquet TC breaking discrete TTS [16, 17, 31, 51, 61, 72], particularly in systems with disorder that facilitate many-body localizations [61, 72], in addition to clean systems [19, 27, 39, 49]. Ongoing studies are further extended to open systems with Floquet driving in the presence of dissipation [14, 21, 22, 37, 46], with experimental investigations reported for a variety of systems [2, 13, 43, 47, 48, 57, 75]. A recent study addresses TC and its associated physics along the imaginary time axis [6].
We introduce time order in this work as the essential element for a new perspective to identify and categorize quantum many-body phases, based on different ground-state temporal patterns. Each quantum many-body Hamiltonian
We will adopt the WO definition of CTC based on two-time auto-correlation function of an operator. It was first outlined in the now famous no-go theorem work [62], and it establishes a general and rigorous subtype of CTC. Recently, Kozin and Kyriienko claim to have realized such a genuine ground-state CTC in a multi-spin model with long-range interaction [33], buttressing much confidence to the search for exotic CTCs. The operational definition that we introduced for time order encompasses WO CTC as one type of time order phases. We will also explore and elaborate a variety of possible exotic phases.
2 Results
2.1 Time Order
We argue that ground-state temporal properties of a quantum many-body system can be used to characterize or classify its phases. Hence, the concept of time order can be introduced analogous to an order parameter by bestowing it in the non-trivial temporal dependence. To exemplify the essence of the associated physics, we shall present an operational definition for time order and accordingly work out the exhaustive list of all allowed phases. According to the WO proposal [62], a witness to CTC is the following two-time (or unequal time) auto-correlation function (with respect to the ground state):
for operator
with
If f(t) is time periodic in the thermodynamic limit, the system is in a state of CTC. This can be reformulated into an explicit operational protocol by introducing a twisted vector. For a quantum many-body system with energy eigen-state
The orthonormal set of eigen-wavefunctions
where ηj ≡|aj|2 denotes weights of the ground-state twisted vector, η0 the corresponding ground-state weight, and ηj (with j > 0) the excited-state weight.
When the coarse-grained order parameter
with
Given an order parameter
If f(t) is a constant, the time dependence will be trivial. However, a subtlety appears when f(t) is vanishingly small with respect to the system size. Since what we are after is the system’s explicit temporal behavior or time dependence, it is easily washed out to f(t) = 0 by a vanishing norm of the twisted vector. Such a difficulty can be mitigated by multiplying system volume V, that is, using the twisted vector |v⟩ → V|v⟩, to check if the correlation for the bulk order parameter F(t) ≡ V2f(t) exhibits temporal dependence, or vanishes as follows:
When f(t) = 0 but F(t) remains a periodic function, the system can still be considered a CTC. Such a remedy surprisingly captures the essence of generalized CTC of Ref. [40].
The analysis presented above can be directly extended to excited states [59]. It is also straightforwardly applicable to non-Hermitian systems, as long as a plausible “ground state” can be identified, for example, by requiring its eigen-energy to possess the largest imaginary part or the smallest norm. Denoting the imaginary part of energy eigen-value Ei as Im(Ei), a prefactor
Therefore, quantum many-body phases can be classified according to time order. The two-time auto-correlation function-based complete operational procedure for classifying time order thus extends the definition of WO CTC as provided in Ref. [62]. Our central results can be simply stated as follows: if f(t) exhibits non-trivial time dependence, then time order exists. If f(t) = 0 but F(t) displays non-trivial time dependence instead, then generalized time order exists.
More specifically, if f(t) = const. is non-zero, the system exhibits time trivial order. The same applies when f(t) = 0 and F(t) = const. For all other situations, non-trivial time order prevails. A complete classification for all time order ground-state phases is shown in Table 1, according to the temporal behaviors of their auto-correlation functions f(t) or F(t). As discussed in Section 4, the above discussion and classification on time order can be extended to finite temperature systems as well.
The operational procedure outlined previously presents a straightforward approach for detecting time order, albeit with reference to an order parameter operator. Hence, more appropriately, this approach should be called order parameter assisted time order or symmetry-based (or -protected) time order to emphasize its reference to symmetry order parameter of a quantum many-body system. The twisted vector facilitates easy calculations to distinguish between different time order phases from time trivial ones, as we illustrate in the following text in terms of a few concrete examples. It is reasonable to expect that transitions between different time order phases can occur, reminiscent of phase transitions in the LGW spontaneous symmetry breaking paradigm.
2.2 Time Order Phase in a Spin-1 Atomic Condensate
A spin-1 atomic Bose–Einstein condensate (BEC) under single spatial mode approximation (SMA) [36, 44, 73] is described by the following Hamiltonian:
where
The validity of this model is well established based on extensive theoretical [8, 24, 71, 76] and experimental [1, 7, 38, 45] studies of spinor BEC over the years. The fractional population in spin states |1, 1⟩ and |1, − 1⟩,
We will concentrate on the zero magnetization Fz = 0 subspace and employ exact diagonalization (ED) to calculate eigen-states. p = 0 is assumed since Fz is conserved. Figure 1 illustrates the system’s complete time order phase diagram. For ferromagnetic interaction c2 < 0 as with 87Rb atoms, the critical quadratic Zeeman shift q/|c2| = 2 splits the whole region into the time trivial order (TT) phase for smaller q that observes TTS and the generalized time crystalline (gTC) order phase for q/|c2| > 2, where TTS is spontaneously broken. The latter (gTC phase) is found to coincide with the polar phase [71]. Limited by available computation resources, the system sizes we explored with ED remain moderate which prevent us from mapping out the finer details in the immediate neighborhood of q = 2|c2|. Further elaboration of time order properties in this region is therefore needed. On the other hand, for antiferromagnetic interaction c2 > 0 with 23Na atoms, we find q = 0 separates TT phase from gTC order. We note here that q = 2|c2| is the second-order quantum phase transition (QPT) critical point between the polar phase and the broken-axisymmetry phase of the ferromagnetic spin-1 BEC, while q = 0 corresponds to the first-order QPT critical point for antiferromagnetic interaction.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/847409/fphy-10-847409-HTML/image_m/fphy-10-847409-g001.jpg)
FIGURE 1. Time order phase diagram for spin-1 atomic BEC, where TT and gTC, respectively, denote time trivial and generalized time crystalline order. The region of (hashed) line segments surrounding c2 = 0 for the non-interacting system is to be excluded.
More detailed discussions including the dependence of time order phases on system size, possible approaches to detect them, and extension to thermal state phases can be found in Section 4.
2.3 Time Order Phase Diagram for Quantum Rabi Model
As a second example, we consider time order phases of the quantum Rabi model described by the Hamiltonian as follows:
where
It is known that the aforementioned model exhibits a QPT to a superradiant state, despite its simplicity [29]. The transition occurs at the critical point gc ≡ 1, with the dimensionless parameter
whose low-energy eigen-states are
whose eigen-states are given by
For this model, the scaled average cavity photon number
For g < 1, we find
respectively, where η0 = sinh4(rnp) and η2 = cosh2(rnp) sinh2(rnp). For g > 1, we obtain the following equation:
The time order phase diagram is shown in Figure 2. When g < 1, the system ground state corresponds to a generalized time crystalline order phase, while the system exhibits time trivial order when g > 1. Despites such a simple model composed of a two-level system and a bosonic field mode, the ground state of the quantum Rabi model displays an intriguing temporal phase structure accompanied by a finite-component quantum phase transition.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/847409/fphy-10-847409-HTML/image_m/fphy-10-847409-g002.jpg)
FIGURE 2. Time order phase diagram for the quantum Rabi model, where TT and gTC, respectively, denote time trivial and generalized time crystalline order.
2.4 Non-Hermitian Many-Body Interaction Model
Finally, we consider two effective models with many-body spin–spin interaction and non-Hermitian effects. The first is described by the Hamiltonian:
with two σy operators at sites i and j in a string of otherwise σx N-body spin interaction. 1/[N(N − 1)] is the normalization factor, λ is the spin interaction strength, and γ represents an effective dissipation rate. λ > 0 and γ ≥ 0 are both real numbers.
We observe that the Greenberger–Horne–Zeilinger (GHZ) states
correspond to two non-degenerate system eigen-states with eigen-energies ± (λ + iγ)/2. The spectra of this model system are bounded inside the circle of radius
An appropriate order parameter operator in this case becomes the average magnetization
A second non-Hermitian model Hamiltonian is given by
where [⋅] denotes the integer part, σN+1 ≡σ1 corresponds to the periodic boundary condition, and λ and γ are spin-string interaction strength and dissipation strength, respectively, as in the previous model, and are real numbers. This Hamiltonian contains [(N + 1)/2]-body interaction terms and supports GHZ state |G+⟩ as a non-degenerate excited state [18] with eigen-energy ϵ+ = − N. The other two eigen-states of concern are
The eigen-energies for |Ψ(±)⟩ are given by
At γ = 0, the above non-Hermitian Hamiltonian Eq. 17 reduces to a Hermitian one, whose ground state |Ψ0⟩ corresponds to the one with smaller ϵ from |Ψ(−)⟩ and |Ψ(+)⟩, or
Therefore, we directly obtain the following equation:
When λ ≠ 0 and γ ≠ 0, the system exists in time functional order phase and again results from the non-Hermitian Hamiltonian. When λ ≠ 0 but γ = 0, the auto-correlation function reduces to
as for a genuine time crystal of the WO type exhibiting time crystalline order. When λ = 0 and 0 < |γ| ≤ 1, we find
The system ground state again exhibits time-crystalline order. When λ = 0 and |γ| > 1, we obtain
by choosing
the ground state reduces to the time trivial order phase.
The above two non-Hermitian models represent direct generalizations of the Hermitian system as considered in Refs. [18, 33]. While slightly more complicated, they remain sufficiently simple for compact analytical treatment, thus helping to reveal interesting and clear physical meanings of the underlying time order.
2.5 Some Remarks About Continuous Time Crystal
According to the WO no-go theorem [62], f(t) for the ground state or the Gibbs ensemble of a general many-body Hamiltonian whose interactions are not-too-long ranged exhibits no temporal dependence; hence, it belongs to time trivial order according to our classification scheme. At first sight, this seems to sweep many important models of condensed matter physics into the same boring class of the time trivial order phase. However, it remains to explore, for instance, many-body systems with more than two-body (or k-body) interactions, or non-Hermitian systems, which might support the existence of CTC. Inspired by the recent results on CTC [33], we believe more time crystalline phases will be uncovered and further understanding will be gained in the future.
The two-time auto-correlation function f(t) measures the CTC phase as in earlier studies [33, 62], while both auto-correlation functions f(t) and F(t) together define different time-order phases we propose in this work. The absence of the local temporal behavior f(t) = 0 does not imply the absence of any temporal property in the bulk, when F(t) could have various temporal behaviors. Based on this, our operational definitions of time order are developed. We also extend the scope of CTC to include both TC order and gTC order phases. This distinction between f(t) and F(t) gives more insights into quantum many-body phases.
As emphasized earlier, continuous time crystal results from spontaneously breaking continuous time translation symmetry. Due to the genuine time periodicity contained in CTC, it might be possible to explore and design new types of clocks based on macroscopic many-body systems, as the time period is directly related to energy spectra, and whose physical meaning is clearly the same as for atomic clock states. Furthermore, they are not affected by finite size effect in contrast to periodicity in DTC.
3 Discussion and Conclusion
While ground-state phases of a quantum many-body system are mostly classified with their Hamiltonian based on the following two paradigms: LGW symmetry breaking order parameter or topological order, this work proposes to study phases from time dimension using time order or more specifically with the proposed symmetry-based time order. Compared to the recent progress and understanding gained for topological order [66, 67], one could try to develop a framework for entanglement-based time order instead of the symmetry-based time order we employ here in this study. Quantum entanglement in a many-body system is responsible for topological order, whose origin lies at the tensor product structure of the quantum many-body Hilbert space
Through time order, one focuses on temporal structure of the evolution operator
In conclusion, understanding the phases of matter constitutes a cornerstone of contemporary physics. Capitalizing on the concept of CTC for the many-body ground state with perpetual time dependence, this study argues that information from time domain can be employed to classify the quantum phase as well, which provides a new perspective toward the understanding of ground-state time dependence, significantly beyond existing studies on CTC. We introduce time order, provide its operational definition in terms of two-time auto-correlation function of an appropriate symmetry order operator, bestow physical meaning to characteristic frequencies and amplitudes of the correlation function, and present a complete classification of time order phases. Time order phase diagrams for a spin-1 BEC system and the quantum Rabi model are fully worked out. Interesting time order phases in non-Hermitian spin models with multi-body interaction are presented. In addition to the time crystalline order which already attracts broad attention from its studies in terms of CTC, other phases we identify, for example, time quasi-crystalline order and time functional order, represent exciting new possibilities.
4 Methods and Calculation Details
Here, in this section, we provide more supporting material for our main results and related details for the aforementioned presentation. It is organized as follows: in Section 4.1, we extend the discussion of time order to finite temperature; in Section 4.2, we present calculation details related to the spin-1 atomic BEC example considered and give intriguing results for finite temperature scenario in spin-1 BEC; in Section 4.3,as a more straightforward approach to understand numerical results, we present a variational approach for treating the polar ground state of a spin-1 BEC. Finally, we give the details about ground state and eigen-energy calculation in the non-Hermitian quantum many-body model with multi-body interaction in Section 4.4.
4.1 Time Order at Finite Temperature
At finite temperature T, excited states will be populated, which can be taken into account with the Gibbs ensemble
where |vk⟩ is the eigen-state twisted vector for |ψk⟩ and cjk is its associated weight. Analogously, for the non-Hermitian case, we find
where
It is easily noted that f(t) at finite temperature contains contributions from all eigen-states of the quantum many-body system
4.2 Time Order in a Spin-1 Atomic BEC
For typical interaction parameters of a spin-1 BEC (e.g., of ground state 87Rb or 23Na atoms) in a tight trap, spin domain formation is energetically suppressed when the atom number is not too large as spin-dependent interaction strength is much weaker than spin-independent interaction [26, 32, 36, 58]. This facilitates a single-spatial-mode approximation (SMA) by assuming all spin states share the same spatial wave function ϕ(r), which effectively decouples the spatial degrees of freedom from the spin and results in the Hamiltonian [36, 44] in Eq. 6 for the model many-body system, where
As discussed in the main text, a suitable order parameter for this model system is
where bj ≡|aj|2 is the weight of the ground-state twisted vector,
where Ai = N⟨ψi|v⟩, Bj ≡|Aj|2 is the weight of the enlarged ground state twisted vector, and
Our study below is for the zero magnetization Fz = 0 subspace and employs exact diagonalization (ED) to calculate eigen-states as well as eigen-energies. The overall time order phase diagram for spin-1 BEC is shown in Figure 1. For ferromagnetic interaction c2 < 0, the critical quadratic Zeeman energy q/|c2| = 2 splits the whole region into the time trivial order (TT) phase for smaller q that observes TTS, and the generalized time crystalline (gTC) order phase for q/|c2| > 2 where TTS is spontaneously broken. The latter (gTC phase) is found to coincide with the ground-state polar phase. The available computation resource limits the calculation to a finite system size, which prevents us from mapping out the exact details in the immediate neighborhood of q = 2, where further elaboration is needed for its time order properties. On the other hand, for antiferromagnetic interactions, we find q = 0 separates TT phase from gTC order.
In Figure 3, the weights for the ground state as well as for the low-lying excited states are shown as functions of q for a typical system size of N = 10 000. Only the ground-state weight b0 is non-vanishing in the q < 2 (q < 0) region for ferromagnetic (antiferromagnetic) interactions, but total weight b is zero in the q > 2 (q > 0) region for ferromagnetic (antiferromagnetic) interaction, which prompts us to examine further the enlarged weights Bi corresponding to the bulk order parameter. For ground and the first excited states, the volume enlarged weights B0,1 are found to be non-vanishing, although both decrease as q increases and grow with N as q approaches q = 2 (q = 0) for ferromagnetic (antiferromagnetic) interaction. However, as mentioned above, limited to a system size of N = 10 000 by computation resource in the ED calculation, we cannot exactly map out the behavior near q = 2 (q = 0) for ferromagnetic (antiferromagnetic) interaction. This consequently leaves empty for q in region [2.0, 2.02] ([0, 0.01]) for ferromagnetic (antiferromagnetic) interaction.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/847409/fphy-10-847409-HTML/image_m/fphy-10-847409-g003.jpg)
FIGURE 3. Weights of ground-state twisted vector in the ground and low-lying excited states as functions of q at system size N = 10 000. The upper panel is for ferromagnetic interaction, where weights bi for q < 2 are shown in (A), while weights Bi for q > 2 are shown in (B). The lower panel is for antiferromagnetic interaction, where weights bi for q < 0 are shown in (C), while weights Bi for q > 0 are shown in (D).
The dependence on system size N is clearly revealed by Figure 4, with the enlarged weights in the gTC regime attaining fixed values as the system approaches thermodynamic limit (N → ∞). In regions away from q = 2 (q = 0) for ferromagnetic (antiferromagnetic) interaction, ED numerics can always approach thermodynamic limit, except for the immediate neighborhood near q = 2 (q = 0), where we infer with confidence the tendencies to divergence of the weights B0,1 as q approaches q = 2 (q = 0).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/847409/fphy-10-847409-HTML/image_m/fphy-10-847409-g004.jpg)
FIGURE 4. Weights of ground-state twisted vector in the ground and low-lying excited states as functions of system size N at q = 2.1 for ferromagnetic interaction (A) and at q = 0.2 for antiferromagnetic interaction (B).
The time evolution of two-time auto-correlation function F(t) is plotted in Figures 5A,C for ferromagnetic and for antiferromagnetic interactions, while Figures 5B,D display energy gaps between ground and the first excited states as a function of q for ferromagnetic and antiferromagnetic interactions, respectively, at a system size of N = 5000. The behavior of F(t) is quantitatively consistent with that of the weights Bi(q) (i = 0, 1) shown in Figure 3 and the energy gap ϵ1 − ϵ0 shown in Figures 5B,D.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/847409/fphy-10-847409-HTML/image_m/fphy-10-847409-g005.jpg)
FIGURE 5. F(t) for different q as a function of time t. The solid and dotted lines correspond to Re(F) and Im(F), respectively. The red, green, and blue lines correspond to q = 2.5, q = 3, and q = 5, respectively, for ferromagnetic interaction (A). The red, green, and blue lines correspond to q = 0.7, q = 1, and q = 3, respectively for antiferromagnetic interaction (C). The energy gap between ground and the first excited state ϵ1 − ϵ0 as a function of q for ferromagnetic (B) and antiferromagnetic interactions (D), at system size N = 5000.
At finite temperature, excited states come into play by also contributing to the correlation function. We find the gTC order hosted in the polar phase persists for both ferromagnetic and antiferromagnetic interactions. The corresponding time evolution and Fourier transform of F(t) are shown in Figure 6, calculated for N = 500 at a temperature of β ≡ 1/T = 1. The Fourier transform is performed for Re(F) over t = [0, 1000] with the zero frequency (DC) component subtracted or for Im(F). The upper (lower) panel corresponds to ferromagnetic (antiferromagnetic) interaction at q = 3 (q = 2). For ferromagnetic interaction, two distinct frequency components are clearly identified for q = 3, associated with the two different energy level gaps. The beautiful beat pattern for F(t) would appear, while we only show the short time behavior in Figure 6 (a). Thus, the gTC phase remains at a finite temperature. Moreover, we also find a generalized time quasi-crystalline order phase assuming the two frequencies are incommensurate, by fine-tuning their corresponding energy gaps such that the relation Δ1/Δ2 = m1/m2 with m1 and m2 being co-primes is not satisfied. The gTC phase at finite temperature here is robust which is in contrast to the melting behavior of continuous time crystal (CTC) shown in Ref. [33].
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/847409/fphy-10-847409-HTML/image_m/fphy-10-847409-g006.jpg)
FIGURE 6. F(t) as a function of time t at q = 3 for ferromagnetic interaction (A) and q = 2 for antiferromagnetic interaction (C). The red and blue solid lines, respectively, correspond to Re(F) and Im(F). The Fourier transform spectrum
Finally, we hope to address the critical question about how could this time order, sort of a perpetual time dependence, can be observed. We note the bulk two-time auto-correlation function introduced
In Figure 7A, we show the behavior of oscillation amplitude for F(t)/Nsum(0). The time dependence of F(t)/Nsum(0) at q = 2.5 for ferromagnetic interaction is shown in Figure 7B.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/847409/fphy-10-847409-HTML/image_m/fphy-10-847409-g007.jpg)
FIGURE 7. (A) Bi/Nsum(0) as a function of q for ferromagnetic interaction. (B) F(t)/Nsum(0) as a function of time t at q = 2.5 for ferromagnetic interaction. The red and blue solid lines correspond to the real and imaginary part of F(t)/Nsum(0), respectively.
4.3 A Variational Polar State for Ferromagnetic Spin-1 BEC
One might naively expect that nothing particularly interesting could happen in the polar phase of a ferromagnetic spin-1 BEC, where essentially all atoms reside in the single particle state |1, 0⟩. Nevertheless, due to the competition between spin exchange interaction c2 and quadratic Zeeman shift q, the ground state of our system differs from |N1 = 0, N0 = N, N−1 = 0⟩, which can be affirmed based on a simple variational analytical calculation given in this section.
We use the number-state basis |N1, N0, N−1⟩ ≡|[N], M, k⟩, where
We see the extreme value (the minimum) of E is reached when cos(ϕ) = ±1, that is, for a real variational parameter a, which will be assumed from now on. This gives the following equation:
with
which determines the locations for the extreme values.
and the corresponding extreme values are as follows:
In the thermodynamic limit N → ∞, they reduce, respectively, to
For ferromagnetic interaction (c2 < 0), E− assumes the minimum, which corresponds to the ground state
Despite the vanishing order parameter nsum in the polar phase (here the gTC order phase from the time order perspective), the enlarged quantity Nsum retains a finite value. Hence, the physics we present here clearly belongs to the realm of quantum effects, beyond the reach of mean-field theory.
4.4 The Non-Hermitian Spin Model With Multi-Body Interaction
The non-Hermitian quantum many-body model Hamiltonian is
with
where [⋅] denotes the integral part, σN+1 ≡σ1, λ, and γ are spin-string interaction strength and dissipation strength, respectively. λ and γ are both real numbers. i is the imaginary unit. σx,y,z are Pauli operators. N is the qubit number of the system. The Hamiltonian has the [(N + 1)/2]-body interaction term and supports the GHZ state |G+⟩ as a non-degenerate excited state.
First, the Greenberger–Horne–Zeilinger (GHZ) states are denoted as follows:
and
where
We immediately know that |G±⟩ is the degenerate ground state of the ferromagnetic Ising Hamiltonian
The action of
Then we know the two eigen-states of
where α1,2 are the undetermined coefficients. By substituting into the Schrödinger equation
We obtain the eigen-energy
If γ = 0, the Hamiltonian is Hermitian, and we have the ground state |Ψ0⟩ ≡|Ψ(−)⟩ with energy
For the GHZ state |G+⟩, we can know it is a non-degenerate excited state with energy ϵ+ = − N, for
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author Contributions
T-CG proposed and conducted the research, LY supervised the research, and T-CG, LY discussed the results and wrote the manuscript.
Funding
This work is supported by the National Key R&D Program of China (grant no. 2018YFA0306504), the National Natural Science Foundation of China (NSFC) (grants nos. 11 654001 and U1930201), and the Key-Area Research and Development Program of GuangDong Province (grant no. 2019B030330001).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
T-CG thanks Qi Liu and Ming Xue for their helpful discussion about spinor Bose–Einstein condensate.
References
1. Anquez M, Robbins BA, Bharath HM, Boguslawski M, Hoang TM, Chapman MS Quantum Kibble-Zurek Mechanism in a Spin-1 Bose-Einstein Condensate. Phys Rev Lett (2016) 116:155301. doi:10.1103/physrevlett.116.155301
2. Autti S, Eltsov VB, Volovik GE Observation of a Time Quasicrystal and its Transition to a Superfluid Time crystal. Phys Rev Lett (2018) 120:215301. doi:10.1103/PhysRevLett.120.215301
3. Bruno P Impossibility of Spontaneously Rotating Time Crystals: A No-Go Theorem. Phys Rev Lett (2013) 111:070402. doi:10.1103/PhysRevLett.111.070402
4. Buča B, Tindall J, Jaksch D Non-stationary Coherent Quantum many-body Dynamics through Dissipation. Nat Commun (2019) 10:1730. doi:10.1038/s41467-019-09757-y
5. Buča B, Jaksch D Dissipation Induced Nonstationarity in a Quantum Gas. Phys Rev Lett (2019) 123:260401. doi:10.1103/PhysRevLett.123.260401
6. Cai Z, Huang Y, Vincent Liu W Imaginary Time crystal of thermal Quantum Matter. Chin Phys. Lett. (2020) 37:050503. doi:10.1088/0256-307x/37/5/050503
7. Chang M-S, Hamley CD, Barrett MD, Sauer JA, Fortier KM, Zhang W, et al. Observation of Spinor Dynamics in Optically TrappedRb87Bose-Einstein Condensates. Phys Rev Lett (2004) 92:140403. doi:10.1103/physrevlett.92.140403
8. Chang M-S, Qin Q, Zhang W, You L, Chapman MS Coherent Spinor Dynamics in a Spin-1 Bose Condensate. Nat Phys (2005) 1:111–6. doi:10.1038/nphys153
9. Chen X, Burnell FJ, Vishwanath A, Fidkowski L Anomalous Symmetry Fractionalization and Surface Topological Order. Phys Rev X (2015) 5:041013. doi:10.1103/PhysRevX.5.041013
10. Chen X, Gu Z-C, Liu Z-X, Wen X-G Symmetry Protected Topological Orders and the Group Cohomology of Their Symmetry Group. Phys Rev B (2013) 87:155114. doi:10.1103/PhysRevB.87.155114
11. Chen X, Liu Z-X, Wen X-G Two-dimensional Symmetry-Protected Topological Orders and Their Protected Gapless Edge Excitations. Phys Rev B (2011) 84:235141. doi:10.1103/PhysRevB.84.235141
12. Cheng M, Gu Z-C, Jiang S, Qi Y Exactly Solvable Models for Symmetry-Enriched Topological Phases. Phys Rev B (2017) 96:115107. doi:10.1103/PhysRevB.96.115107
13. Choi S, Choi J, Landig R, Kucsko G, Zhou H, Isoya J, et al. Observation of Discrete Time-Crystalline Order in a Disordered Dipolar many-body System. Nature (2017) 543:221–5. doi:10.1038/nature21426
14. Cosme JG, Skulte J, Mathey L Time Crystals in a Shaken Atom-Cavity System. Phys Rev A (2019) 100:053615. doi:10.1103/PhysRevA.100.053615
15. Damski B, Zurek WH Dynamics of a Quantum Phase Transition in a Ferromagnetic Bose-Einstein Condensate. Phys Rev Lett (2007) 99:130402. doi:10.1103/PhysRevLett.99.130402
16. Else DV, Bauer B, Nayak C. Floquet Time Crystals. Phys Rev Lett (2016) 117:090402. doi:10.1103/PhysRevLett.117.090402
17. Else DV, Monroe C, Nayak C, Yao NY Discrete Time Crystals. Annu Rev Condens Matter Phys (2020) 11:467–99. doi:10.1146/annurev-conmatphys-031119-050658
18. Facchi P, Florio G, Pascazio S, Pepe FV Greenberger-horne-zeilinger States and Few-Body Hamiltonians. Phys Rev Lett (2011) 107:260502. doi:10.1103/PhysRevLett.107.260502
19. Fan C-h., Rossini D, Zhang H-X, Wu J-H, Artoni M, La Rocca GC Discrete Time crystal in a Finite Chain of Rydberg Atoms without Disorder. Phys Rev A (2020) 101:013417. doi:10.1103/PhysRevA.101.013417
20. Fradkin E Field Theories of Condensed Matter Physics. 2nd ed. Oxford, UK: Oxford University Press (2004).
21. Gambetta FM, Carollo F, Marcuzzi M, Garrahan JP, Lesanovsky I Discrete Time Crystals in the Absence of Manifest Symmetries or Disorder in Open Quantum Systems. Phys Rev Lett (2019) 122:015701. doi:10.1103/PhysRevLett.122.015701
22. Gong Z, Hamazaki R, Ueda M Discrete Time-Crystalline Order in Cavity and Circuit Qed Systems. Phys Rev Lett (2018) 120:040404. doi:10.1103/PhysRevLett.120.040404
23. Gu Z-C, Wen X-G Tensor-entanglement-filtering Renormalization Approach and Symmetry-Protected Topological Order. Phys Rev B (2009) 80:155131. doi:10.1103/PhysRevB.80.155131
24. Guzman J, Jo G-B, Wenz AN, Murch KW, Thomas CK, Stamper-Kurn DM Long-time-scale Dynamics of Spin Textures in a Degenerate F = 187rb Spinor Bose Gas. Phys Rev A (2011) 84:063625. doi:10.1103/physreva.84.063625
25. Heinrich C, Burnell F, Fidkowski L, Levin M Symmetry-enriched String Nets: Exactly Solvable Models for Set Phases. Phys Rev B (2016) 94:235136. doi:10.1103/PhysRevB.94.235136
26. Ho T-L Spinor Bose Condensates in Optical Traps. Phys Rev Lett (1998) 81:742–5. doi:10.1103/PhysRevLett.81.742
27. Huang B, Wu Y-H, Liu WV Clean Floquet Time Crystals: Models and Realizations in Cold Atoms. Phys Rev Lett (2018) 120:110603. doi:10.1103/PhysRevLett.120.110603
28. Hurtado-Gutiérrez R, Carollo F, Pérez-Espigares C, Hurtado PI Building Continuous Time Crystals from Rare Events. Phys Rev Lett (2020) 125:160601. doi:10.1103/PhysRevLett.125.160601
29. Hwang M-J, Puebla R, Plenio MB Quantum Phase Transition and Universal Dynamics in the Rabi Model. Phys Rev Lett (2015) 115:180404. doi:10.1103/PhysRevLett.115.180404
30. Iemini F, Russomanno A, Keeling J, Schirò M, Dalmonte M, Fazio R Boundary Time Crystals. Phys Rev Lett (2018) 121:035301. doi:10.1103/PhysRevLett.121.035301
31. Khemani V, Lazarides A, Moessner R, Sondhi SL Phase Structure of Driven Quantum Systems. Phys Rev Lett (2016) 116:250401. doi:10.1103/PhysRevLett.116.250401
32. Koashi M, Ueda M Exact Eigenstates and Magnetic Response of Spin-1 and Spin-2 Bose-Einstein Condensates. Phys Rev Lett (2000) 84:1066–9. doi:10.1103/PhysRevLett.84.1066
33. Kozin VK, Kyriienko O Quantum Time Crystals from Hamiltonians with Long-Range Interactions. Phys Rev Lett (2019) 123:210602. doi:10.1103/PhysRevLett.123.210602
34. Lamacraft A Quantum Quenches in a Spinor Condensate. Phys Rev Lett (2007) 98:160404. doi:10.1103/physrevlett.98.160404
36. Law CK, Pu H, Bigelow NP Quantum Spins Mixing in Spinor Bose-Einstein Condensates. Phys Rev Lett (1998) 81:5257–61. doi:10.1103/PhysRevLett.81.5257
37. Lazarides A, Roy S, Piazza F, Moessner R Time Crystallinity in Dissipative Floquet Systems. Phys Rev Res (2020) 2:022002. doi:10.1103/PhysRevResearch.2.022002
38. Luo XY, Zou YQ, Wu LN, Liu Q, Han MF, Tey MK, et al. Deterministic Entanglement Generation from Driving through Quantum Phase Transitions. Science (2017) 355:620–3. doi:10.1126/science.aag1106
39. Machado F, Else DV, Kahanamoku-Meyer GD, Nayak C, Yao NY Long-range Prethermal Phases of Nonequilibrium Matter. Phys Rev X (2020) 10:011043. doi:10.1103/PhysRevX.10.011043
40. Medenjak M, Buča B, Jaksch D Isolated Heisenberg Magnet as a Quantum Time crystal. Phys Rev B (2020) 102:041117. doi:10.1103/PhysRevB.102.041117
41. Mierzejewski M, Giergiel K, Sacha K Many-body Localization Caused by Temporal Disorder. Phys Rev B (2017) 96:140201. doi:10.1103/PhysRevB.96.140201
42. Nozières P Time Crystals: Can Diamagnetic Currents Drive a Charge Density Wave into Rotation? Europhysics Lett (2013) 103:57008. doi:10.1209/0295-5075/103/57008
43. Pal S, Nishad N, Mahesh TS, Sreejith GJ Temporal Order in Periodically Driven Spins in star-shaped Clusters. Phys Rev Lett (2018) 120:180602. doi:10.1103/PhysRevLett.120.180602
44. Pu H, Law CK, Raghavan S, Eberly JH, Bigelow NP Spin-mixing Dynamics of a Spinor Bose-Einstein Condensate. Phys Rev A (1999) 60:1463–70. doi:10.1103/PhysRevA.60.1463
45. Qiu LY, Liang HY, Yang YB, Yang HX, Tian T, Xu Y, et al. Observation of Generalized Kibble-Zurek Mechanism across a First-Order Quantum Phase Transition in a Spinor Condensate. Sci Adv (2020) 6:eaba7292. doi:10.1126/sciadv.aba7292
46. Riera-Campeny A, Moreno-Cardoner M, Sanpera A Time Crystallinity in Open Quantum Systems. Quantum (2020) 4:270. doi:10.22331/q-2020-05-25-270
47. Rovny J, Blum RL, Barrett SE P31 NMR Study of Discrete Time-Crystalline Signatures in an Ordered crystal of Ammonium Dihydrogen Phosphate. Phys Rev B (2018) 97:184301. doi:10.1103/PhysRevB.97.184301
48. Rovny J, Blum RL, Barrett SE Observation of Discrete-Time-crystal Signatures in an Ordered Dipolar many-body System. Phys Rev Lett (2018) 120:180603. doi:10.1103/PhysRevLett.120.180603
49. Russomanno A, Iemini F, Dalmonte M, Fazio R Floquet Time crystal in the Lipkin-Meshkov-Glick Model. Phys Rev B (2017) 95:214307. doi:10.1103/PhysRevB.95.214307
50. Sacha K Anderson Localization and mott Insulator Phase in the Time Domain. Sci Rep (2015) 5:10787. doi:10.1038/srep10787
51. Sacha K Modeling Spontaneous Breaking of Time-Translation Symmetry. Phys Rev A (2015) 91:033617. doi:10.1103/physreva.91.033617
53. Senthil T, Vishwanath A, Balents L, Sachdev S, Fisher MPA Deconfined Quantum Critical Points. Science (2004) 303:1490–4. doi:10.1126/science.1091806
54. Shapere A, Wilczek F Classical Time Crystals. Phys Rev Lett (2012) 109:160402. doi:10.1103/physrevlett.109.160402
55. Shirley W, Slagle K, Chen X Universal Entanglement Signatures of Foliated Fracton Phases. Scipost Phys (2019) 6:15. doi:10.21468/SciPostPhys.6.1.015
56. Shirley W, Slagle K, Wang Z, Chen X Fracton Models on General Three-Dimensional Manifolds. Phys Rev X (2018) 8:031051. doi:10.1103/PhysRevX.8.031051
57. Smits J, Liao L, Stoof HTC, van der Straten P Observation of a Space-Time crystal in a Superfluid Quantum Gas. Phys Rev Lett (2018) 121:185301. doi:10.1103/PhysRevLett.121.185301
58. Stamper-Kurn DM, Ueda M Spinor Bose Gases: Symmetries, Magnetism, and Quantum Dynamics. Rev Mod Phys (2013) 85:1191–244. doi:10.1103/RevModPhys.85.1191
59. Syrwid A, Zakrzewski J, Sacha K Time crystal Behavior of Excited Eigenstates. Phys Rev Lett (2017) 119:250602. doi:10.1103/PhysRevLett.119.250602
60. Vijay S, Haah J, Fu L Fracton Topological Order, Generalized Lattice Gauge Theory, and Duality. Phys Rev B (2016) 94:235157. doi:10.1103/PhysRevB.94.235157
61. von Keyserlingk CW, Khemani V, Sondhi SL Absolute Stability and Spatiotemporal Long-Range Order in Floquet Systems. Phys Rev B (2016) 94:085112. doi:10.1103/PhysRevB.94.085112
62. Watanabe H, Oshikawa M Absence of Quantum Time Crystals. Phys Rev Lett (2015) 114:251603. doi:10.1103/physrevlett.114.251603
63. Wen XG Vacuum Degeneracy of Chiral Spin States in Compactified Space. Phys Rev B (1989) 40:7387–90. doi:10.1103/PhysRevB.40.7387
64. Wen XG Topological Orders in Rigid States. Int J Mod Phys B (1990) 04:239–71. doi:10.1142/S0217979290000139
65. Wen XG Quantum Orders and Symmetric Spin Liquids. Phys Rev B (2002) 65:165113. doi:10.1103/physrevb.65.165113
66. Wen XG Colloquium : Zoo of Quantum-Topological Phases of Matter. Rev Mod Phys (2017) 89:041004. doi:10.1103/RevModPhys.89.041004
67. Wen XG. Choreographed Entanglement Dances: Topological States of Quantum Matter. Science (2019) 363:eaal3099. doi:10.1126/science.aal3099
69. Wilczek F Quantum Time Crystals. Phys Rev Lett (2012) 109:160401. doi:10.1103/PhysRevLett.109.160401
70. Wilson K, Kogut J The Renormalization Group and the ε Expansion. Phys Rep (1974) 12:75–199. doi:10.1016/0370-1573(74)90023-4
71. Xue M, Yin S, You L Universal Driven Critical Dynamics across a Quantum Phase Transition in Ferromagnetic Spinor Atomic Bose-Einstein Condensates. Phys Rev A (2018) 98:013619. doi:10.1103/PhysRevA.98.013619
72. Yao NY, Potter AC, Potirniche ID, Vishwanath A Discrete Time Crystals: Rigidity, Criticality, and Realizations. Phys Rev Lett (2017) 118:030401. doi:10.1103/PhysRevLett.118.030401
73. Yi S, Müstecaplıoğlu ÖE, Sun CP, You L Single-mode Approximation in a Spinor-1 Atomic Condensate. Phys Rev A (2002) 66:011601. doi:10.1103/PhysRevA.66.011601
74. Zeng B, Wen XG Gapped Quantum Liquids and Topological Order, Stochastic Local Transformations and Emergence of Unitarity. Phys Rev B (2015) 91:125121. doi:10.1103/PhysRevB.91.125121
75. Zhang J, Hess PW, Kyprianidis A, Becker P, Lee A, Smith J, et al. Observation of a Discrete Time crystal. Nature (2017) 543:217–20. doi:10.1038/nature21413
Keywords: time order, time crystal, quantum phase, Bose–Einstein condensate, non-Hermitian many-body physics, fully connected model, exotic phase
Citation: Guo T- and You L (2022) Quantum Phases of Time Order in Many-Body Ground States. Front. Phys. 10:847409. doi: 10.3389/fphy.2022.847409
Received: 02 January 2022; Accepted: 24 January 2022;
Published: 22 March 2022.
Edited by:
Weibin Li, University of Nottingham, United KingdomReviewed by:
Yongqiang Li, National University of Defense Technology, ChinaYunbo Zhang, Zhejiang Sci-Tech University, China
Copyright © 2022 Guo and You. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Tie-Cheng Guo, gtc16@tsinghua.org.cn; Li You, lyou@mail.tsinghua.edu.cn