Quantum Phase Transition in the Coupled-Top Model: From Z2 to U(1) Symmetry Breaking

TL;DR

The paper studies how symmetry breaking in a quantum model leads to different types of quantum phase transitions.

Contribution

The study reveals a novel critical exponent of 1 for U(1) symmetry breaking, distinct from the usual 1/2 for Z2 symmetry.

Findings

The anisotropic model shows a Z2 symmetry breaking leading to second-order quantum phase transitions.

U(1) symmetry breaking results in a vanishing energy gap and a new critical exponent of 1.

The model's rich symmetry structure makes it a key system for studying quantum phase transitions.

Abstract

We investigate the coupled-top model, which describes two large spins interacting along both x and y directions. By tuning coupling strengths along distinct directions, the system exhibits different symmetries, ranging from a discrete Z2 to a continuous U(1) symmetry. The anisotropic coupled-top model displays a discrete Z2 symmetry, and the symmetry breaking induced by strong coupling drives a quantum phase transition from a disordered paramagnetic phase to an ordered ferromagnetic or antiferromagnetic phase. In particular, the isotropic coupled-top model possesses a continuous U(1) symmetry, whose breaking gives rise to the Goldstone mode. The phase boundary can be well captured by the mean-field approach, characterized by the distinct behaviors of the order parameter. Higher-order quantum effects beyond the mean-field contribution can be achieved by mapping the large spins to bosonic operators via the Holstein–Primakoff transformation. For the anisotropic coupled-top model with Z2 symmetry, the energy gap closes, and both quantum fluctuations and entanglement entropy diverge near the critical point, signaling the onset of second-order quantum phase transitions. Strikingly, when U(1) symmetry is broken, the energy gap vanishes beyond the critical point, yielding a novel critical exponent of 1, rather than 1/2 for Z2 symmetry breaking. The rich symmetry structure of the coupled-top model underpins its role as a paradigmatic model for studying quantum phase transitions and exploring associated physical phenomena.

Figures

Full text

1. Introduction

Quantum phase transitions have emerged as one of the central topics in modern physics [1,2,3], offering profound insights into the fundamental nature of quantum many-body systems. These transitions, which occur at zero temperature and are driven by quantum fluctuations rather than classical thermal fluctuations, often lead to dramatic changes in the ground-state properties of materials, such as the emergence of new ordered phases and the appearance of exotic quantum states [2]. Understanding quantum phase transitions not only enriches our theoretical knowledge of quantum physics but also holds great promise for technological applications, such as the critical quantum metrology [4,5,6,7,8,9].

In spin interaction systems, numerous prototype models, such as the transverse-field Ising model (TFIM), XY model, and their generalizations [10,11,12,13,14,15,16,17], have been extensively studied to explore quantum phase transitions in both the ground states and dynamics. Recently, several approaches involving anisotropic coupling, such as the Monte Carlo method, mean-field theory, and feedback-based quantum algorithms, have attracted attention for investigating phase transition phenomena [18,19,20,21]. Specifically, the anisotropic XY model [18], the anisotropic two-dimensional Ising case with triangular or rectangular lattice [19,20], and the anisotropic next-nearest-neighbor Ising model [21] provide insights into the system’s symmetries and critical properties.

When spin- $[eqn]$ s in the TFIM are replaced with large spins, we obtain the coupled-top model (CTM) [22,23,24,25,26]. As a paradigmatic bipartite system, the CTM is employed to study the interaction between two large spins. Persistent investigations of dynamical and statistical behaviors in the CTM, including ergodic behavior, quantum scars, and level spacing distribution [23,24,25], have been conducted. In particular, various types of quantum phase transitions in the CTM and its generalized cases [25,27,28,29,30], such as the dynamical, ground-state, and excited-state phase transitions, have attracted significant attention. The detection of novel physical phenomena can be effectively carried out using the out-of-time-order correlator [31,32], which is experimentally implementable [33,34]. Recently, the triangular coupled-top model has emerged as a new platform for studying unconventional frustrated magnetic behaviors, exhibiting phenomena that merit further investigation [35].

Unlike the TFIM or XY model, which is typically studied on a one-dimensional spin chain or higher-dimensional lattices, the coupled-top model is a bipartite system. This model is more analogous to the finite-component light-matter interaction systems [36,37,38,39,40], undergoing the second-order quantum phase transitions from the disorder to ordered phase [41,42,43,44]. Paradigmatic models, such as Rabi, Dicke, and spin-boson model, allow control over symmetries either $[eqn]$ or U(1) by tuning anisotropic coupling strengths along two orthogonal directions [45,46,47,48,49]. Notably, in systems with either symmetry preservation or breaking, anisotropic matter–light interactions induce novel phase diagrams and quantum triple points [47,48,49], which have received widespread attention. Significantly, these models reduce to the Jaynes–Cummings model [50], Tavis–Cummings model [51], and spin–boson model [52] with rotating wave approximation, accompanied by the symmetry changing from $[eqn]$ to U(1). The latter exhibits distinct critical behaviors related to the Goldstone mode [53,54]. Similarly, the anisotropy in the coupled-top model can be introduced by adjusting coupling constants along x- and y-directions separately. This scheme provides an important platform for both experimental and theoretical control of symmetries, as well as for investigating associated physical properties.

In this paper, we investigate quantum phase transitions of the coupled-top model resulting from $[eqn]$ and U(1) symmetry breaking and the associated quantum criticality. The paper is organized as follows. In Section 2, we define the coupled-top model and explain the $[eqn]$ and U(1) symmetries. In Section 3, we derive the mean-field solutions and corresponding phase diagram. In Section 4, we extend the analysis beyond the mean-field approach, calculating the excitation energy, quantum fluctuations, entanglement entropy, and critical phenomena. In Section 5, we examine the energy gap related to the Goldstone mode for U(1) symmetry. Finally, Section 6 provides the conclusions of this study.

2. The Model and Symmetries

The coupled-top model can be regarded as a generalization of the XY model in a transverse field, with spin- $[eqn]$ replaced by large spins. Its Hamiltonian can be written as

$$[eqn]$$

where $[eqn]$ denotes the d-component $[eqn]$ collective spin operators of magnitude J for the i-th large spin. $[eqn]$ characterizes the strength of the transverse field. $[eqn]$ and $[eqn]$ represent coupling constants between spins in the x- and y-directions, respectively. For clarity, we express the dimensionless Hamiltonian as

$$[eqn]$$

with coupling constants $[eqn]$ and $[eqn]$ .

Previous studies on the coupled-top model generally consider the interaction in x-direction ( $[eqn]$ ) [22,23,24,25,30], which is highly anisotropic and only has $[eqn]$ symmetry. By introducing both $[eqn]$ and $[eqn]$ , Hamiltonian (2) becomes much more flexible. Significantly, the generalized coupled-top model bridges the gap between discrete $[eqn]$ symmetry and continuous U(1) symmetry. We focus on the spontaneous symmetry breaking and critical behaviors associated with quantum phase transitions in the thermodynamic limit ( $[eqn]$ ).

For $[eqn]$ , the coupled-top model is isotropic in the $[eqn]$ plane. We can easily prove that the Hamiltonian is invariant under an arbitrary rotation along z direction, namely, $[eqn]$ , with

$$[eqn]$$

It indicates that the isotropic coupled-top model possesses the U(1) symmetry.

For $[eqn]$ , the coupled-top model becomes anisotropic, which exists a $[eqn]$ symmetry. It manifests the invariance of Hamiltonian under the spin flipping along x- and y-axis, namely, $[eqn]$ , $[eqn]$ , and $[eqn]$ , which corresponds to $[eqn]$ in Equation (3). Hamiltonian (2) satisfies the commutation relation $[eqn]$ , with the parity operator

$$[eqn]$$

For weak coupling without symmetry breaking, the transverse field plays the dominated role, and the coupled-top model tends to be in the paramagnetic phase (PP), with magnetization aligned along the transverse field (z-direction). The preserved symmetry leads to $[eqn]$ and $[eqn]$ . In contrast, strong coupling can break the $[eqn]$ or U(1) symmetry, which leads to the emergence of ferromagnetic phase (FP) or antiferromagnetic phase (AFP). The symmetry breaking is accompanied by the spontaneous magnetization in x- or y-directions, namely, $[eqn]$ or $[eqn]$ . Therefore, $[eqn]$ and $[eqn]$ serve as the order parameter, which can be employed to distinguish different phases.

3. Phase Diagram

To facilitate the analytical study, we first investigate the phase diagram of the coupled-top model in the thermodynamic limit ( $[eqn]$ ). The mean-field approach, where correlations between large spins are neglected, provides an efficient method to distinguish different phases, especially in the thermodynamic limit ( $[eqn]$ ) [55,56]. Within the mean-field framework [57,58], we can construct a trial ground state by a tensor product of Bloch coherent states [59],

$$[eqn]$$

The expectation values of the large spin operators are represented by the corresponding points on the Bloch sphere as

$$[eqn]$$

with $[eqn]$ , which can be employed to derive the averaged energy expectation value of the system,

$$[eqn]$$

The variational parameters $[eqn]$ and $[eqn]$ can be determined by minimizing the energy expectation value given in Equation (7). Applying the variational principle, we can achieve the following equations:

$$[eqn]$$ $$[eqn]$$ $$[eqn]$$ $$[eqn]$$

Solving Equations (8a)–(8d) self-consistently, we can achieve $[eqn]$ and $[eqn]$ , as well as the ground-state energy $[eqn]$ and order parameters $[eqn]$ . According to distinct behaviors of the order parameters, the phase diagram is shown in Figure 1a, and the corresponding phase transition points are identified as

$$[eqn]$$

The quantum phase transitions along the x-direction occur at $[eqn]$ , corresponding to transitions from PP to FP and AFP, respectively. Similar behaviors are observed along the y-direction with respect to $[eqn]$ . The white dashed lines separate the disordered phase (PP) from ordered phases (FP and AFP), while the red solid lines ( $[eqn]$ ) distinguish between FP and AFP. Notably, each point along the red solid lines breaks the continuous U(1) symmetry and exhibits an infinite degeneracy. Overall, we classify the ground states into five phases. To visualize behaviors of the ground states in different phase, we display the energy surface $[eqn]$ [Equation (7)] as a function of $[eqn]$ and $[eqn]$ in Figure 1b–e, whose minimum corresponds to the ground state. A detailed analysis of their physical properties is presented below:

(i) The symmetry-preserving region, where the disordered paramagnetic phase (PP) dominates, occupies the center of the phase diagram, with $[eqn]$ $[eqn]$ . The PP can be characterized by the analytical solutions

$$[eqn]$$

As evidenced in Figure 1b, the numerical energy minimum resides at $[eqn]$ for $[eqn]$ and $[eqn]$ , confirming the theoretical predictions.

(ii) The $[eqn]$ symmetry-breaking region along the x-direction arises when the coupling strength $[eqn]$ exceeds the critical value while maintaining $[eqn]$ . This region can be described by

$$[eqn]$$

Specifically, for $[eqn]$ , the ordered ferromagnetic phase (x-FP) dominates with

$$[eqn]$$

For $[eqn]$ , the system stabilizes in the ordered antiferromagnetic phase (x-AFP) where

$$[eqn]$$

Notably, the two possible values of $[eqn]$ induce a double degeneracy in both the FP and AFP. This is consistent with numerical results shown in Figure 1c, where the energy displays two minima, one with a positive and the other with a negative value, along the $[eqn]$ axis for $[eqn]$ and $[eqn]$ .

(iii) The $[eqn]$ symmetry-breaking region along the y-direction occurs when $[eqn]$ varies and crosses the critical value under the condition $[eqn]$ . This regime is governed by the following:

$$[eqn]$$

Clearly, for $[eqn]$ , the ordered ferromagnetic phase (y-FP) is dominant with

$$[eqn]$$

For $[eqn]$ , the system enters the ordered antiferromagnetic phase (y-AFP) where

$$[eqn]$$

Analogous to the x-direction case, a two-fold degeneracy with dual minima takes place along the y-axis, which is visualized in Figure 1d for $[eqn]$ and $[eqn]$ .

(iv) The U(1) symmetry-breaking region emerges beyond the critical points (indicated by red solid lines in Figure 1 for $[eqn]$ ). This region is analytically defined by:

$$[eqn]$$

Crucially, the infinite degeneracy arises from the continuous freedom of $[eqn]$ , with the following conditions: $[eqn]$ for $[eqn]$ , $[eqn]$ for $[eqn]$ , $[eqn]$ for $[eqn]$ , and $[eqn]$ for $[eqn]$ , respectively. As shown in Figure 1e for $[eqn]$ , the ground states form a circular manifold of energy minima resembling a Mexican hat profile. All spin polarizations align uniformly in the $[eqn]$ -direction while permitting continuous $[eqn]$ variations, a hallmark of continuous symmetry breaking.

Above all, by adjusting coupling strength $[eqn]$ and $[eqn]$ , we achieve precise control over the conservation and breaking of $[eqn]$ and U(1) symmetries. This capability establishes a promising platform for investigating quantum phase transitions and their associated emergent phenomena.

4. Quantum Criticality

Although the mean-field approach successfully produces the phase diagram, certain critical behaviors in the vicinity of the phase transition points remain hidden. This is because these behaviors are associated with higher-order contributions that lie outside the scope of the mean-field approximation. To systematically address these effects, we first apply a unitary transformation [35,57,60,61,62], $[eqn]$ , with

$$[eqn]$$ $$[eqn]$$

where $[eqn]$ and $[eqn]$ are the mean-field solutions from Section 3. This unitary transformation rotates spins sequentially around the z-axis with angles $[eqn]$ followed by rotations about y-axis by angles $[eqn]$ . The spin operators in the transformed Hamiltonian $[eqn]$ can be expressed as

$$[eqn]$$

Subsequently, we introduce the mapping from the spin operators to the bosonic creation and annihilation operators via Holstein–Primakoff transformation [63]. In the thermodynamic limit $[eqn]$ where $[eqn]$ , the transformation is expressed as

$$[eqn]$$ $$[eqn]$$ $$[eqn]$$

where the bosonic operators satisfy $[eqn]$ .

By combining Equations (20) and (21), we can expand the Hamiltonian $[eqn]$ and neglect the terms of order $[eqn]$ with $[eqn]$ due to the limit $[eqn]$ , yielding:

$$[eqn]$$

The first term $[eqn]$ represents the ground-state energy obtained from the mean-field approach. The second term can be eliminated, namely, $[eqn]$ in Equation (A1), as $[eqn]$ and $[eqn]$ satisfy Equations (8a)–(8d). The third term $[eqn]$ has a quadratic form, as expressed in Equation (A2) in Appendix A. It can be further simplified by substituting $[eqn]$ and $[eqn]$ in different phases given in Section 3. For example, by substituting the values of $[eqn]$ and $[eqn]$ for the x- and y-direction symmetry breaking, we can simplify the general form of $[eqn]$ in Equation (A2) as

$$[eqn]$$

where $[eqn]$ . It indicates that the index $[eqn]$ designates the principal symmetry-breaking axis, while $[eqn]$ marks the orthogonal direction with preserved symmetry. The discrete parameter $[eqn]$ serves as a region indicator, corresponding to the regimes $[eqn]$ and $[eqn]$ , respectively. Furthermore, Equation (23) is also suitable to the symmetry-preserving PP region, which corresponds to $[eqn]$ .

To systematically diagonalize the quadratic Hamiltonian $[eqn]$ , we employ the $[eqn]$ representation by constituting the canonical vector $[eqn]$ with $[eqn]$ and $[eqn]$ . $[eqn]$ satisfies the canonical commutation relations, $[eqn]$ , where $[eqn]$ with $[eqn]$ and $[eqn]$ the $[eqn]$ identity and null matrices, respectively. Then, the quadratic Hamiltonian $[eqn]$ can be rewritten as

$$[eqn]$$

with $[eqn]$ and

$$[eqn]$$

Based on the Williamson theorem [64], for the positively defined real matrix H, we can find a symplectic transformation S ( $[eqn]$ ) such that

$$[eqn]$$

In virtue of the symplectic transformation S, we can construct a new vector of canonical operators $[eqn]$ , with which the quadratic Hamiltonian $[eqn]$ can be decoupled into two independent harmonic oscillators, as follows:

$$[eqn]$$

In the symmetry-breaking regions where $[eqn]$ , above parameters are given by $[eqn]$ , $[eqn]$ and $[eqn]$ . In contrast, in the symmetry-preserving PP region where $[eqn]$ , the parameters simplify to $[eqn]$ , $[eqn]$ and $[eqn]$ .

Moreover, using the symplectic transformation S, we can construct $[eqn]$ covariance matrix $[eqn]$ [35,60,64,65], given by $[eqn]$ . The quantum fluctuations $[eqn]$ in position space and $[eqn]$ in momentum space are defined by the diagonal elements of the covariance matrix, namely $[eqn]$ and $[eqn]$ for $[eqn]$ .

Close to distinct quantum phase transition points $[eqn]$ and $[eqn]$ , we can analytically compute the the excitation energy $[eqn]$ , quantum fluctuations $[eqn]$ and $[eqn]$ . The explicit expressions for these quantities in both the symmetry-breaking and the symmetry-preserving paramagnetic phases are derived in Appendix B; see in particular Equations (B3) and (B4).

The entanglement entropy [35,60,65,66,67], serving as a measure of the entanglement between distinct subsystems, is intimately connected to Heisenberg’s uncertainty relation in the context of the quadratic Hamiltonian for interacting bosonic systems. In terms of $[eqn]$ and $[eqn]$ , the entanglement entropy can be calculated by

$$[eqn]$$

As shown in Figure 2, we plot the panels for the excitation energy $[eqn]$ , the quantum fluctuation $[eqn]$ , and the entanglement entropy $[eqn]$ as functions of the coupling strength $[eqn]$ with $[eqn]$ held constant. In panel (a), the lowest excitation energy $[eqn]$ , which corresponds to the energy gap of the system, vanishes at the critical points $[eqn]$ . This behavior indicates the second-order quantum phase transitions from the FP to the PP, and from the PP to the AFP, along the x-direction. Similarly, the quantum fluctuation $[eqn]$ in panel (d) and the entanglement entropy $[eqn]$ in panel (g) both exhibit divergences close to the critical points, but tend to stabilize at constant values away from these points. Clearly, these features provide strong evidence for the existence of quantum phase transitions.

All the quantum phase transitions shown in Figure 2 are associated with the $[eqn]$ symmetry breaking along the x-direction. To better understand the universal quantum criticality, it is essential to investigate the critical behaviors of these phase transitions. In panels (b,c), the ln-ln plots of the lowest excitation energy consistently follow the exponential law $[eqn]$ near critical points in four regions: $[eqn]$ (x-FP), $[eqn]$ (PP), $[eqn]$ (PP), and $[eqn]$ (x-FP). The divergences of the quantum fluctuation in panels (e,f) exhibit exponential scaling relations in each region, where $[eqn]$ . Subsequently, we obtain that $[eqn]$ . Moreover, symmetry-breaking phenomena along the y-direction exhibit analogous mechanisms to those observed in the x-direction, and are consequently omitted from the figure to avoid redundancy.

Through exact analysis in Equations (B3) and (B4) and numerical fitting in Figure 2, we summarize the complete behaviors of the quantum criticality in Table 1. These results are consistent with the phenomena observed in the coupled-top model, as well as in the light–matter interaction systems such as the Rabi and Dicke models, which all exhibit the $[eqn]$ symmetry [41,47,48].

5. Excitation Spectra for U(1) Symmetry

Special attention should be paid to $[eqn]$ , which corresponds to the isotropic coupled-top model with U(1) symmetry. Strong spin–spin interaction can break the U(1) symmetry, with critical points $[eqn]$ . Furthermore, these critical points represent the intersections of three phases: PP, FP, and AFP, as shown in Figure 1a. Each point in the U(1) symmetry-breaking region exhibits infinite degeneracy, associated with the Goldstone mode.

Using the transformations as in the previous section, we can derive a similar form of the quadratic term $[eqn]$ as in Equations (24) and (27). In the U(1) symmetry-breaking regions, the parameters are $[eqn]$ , $[eqn]$ , and $[eqn]$ , whereas in the symmetry-preserving regions, the parameters become $[eqn]$ , $[eqn]$ , and $[eqn]$ . The discrete parameter $[eqn]$ serves as an indicator, corresponding to $[eqn]$ and $[eqn]$ , respectively.

Consequently, for $[eqn]$ , the excitation spectra for U(1) symmetry are given by

$$[eqn]$$

While in symmetry-breaking regions for $[eqn]$ and $[eqn]$ ,

$$[eqn]$$

As shown in Figure 3, two branches of excitation energies, represented by blue and red dashed lines, intersect at $[eqn]$ . In contrast to the dispersions shown in Figure 2a, the energy gap of the lowest branch vanishes beyond the critical points $[eqn]$ , consistent with the infinite degeneracy of the Goldstone mode [48,50,68]. Moreover, near the critical points where $[eqn]$ and $[eqn]$ , the energy gap follows a significantly different exponential relation,

$$[eqn]$$

where the exponent value is 1. This unusual exponent value has also been found in other U(1)-symmetry-breaking systems, such as the Jaynes–Cummings model [50].

The vanishing energy gap and the corresponding distinct exponent highlight the uniqueness of the continuous symmetry in contrast to the discrete symmetry. The coupled-top model provides an ideal platform for investigations of critical behaviors associated with distinct symmetries.

6. Conclusions

In this paper, we investigate the quantum phase transition in the coupled-top model. It describes the interactions between two spin ensembles along both x- and y-directions, which can be regarded as a generalization of the XY model in a transverse field. Notably, by tuning the coupling constants, the coupled-top model bridges the gap between discrete $[eqn]$ symmetry and continuous U(1) symmetry, which provides an important platform to study distinct critical behaviors associated with different symmetry breaking.

Using the mean-field approach, we derive the novel phase diagram in the thermodynamic limit $[eqn]$ , which consists of five phases: the disordered paramagnetic phase (PP), the ordered ferromagnetic phase with symmetry breaking along x- or y-direction (x-FP or y-FP), the ordered antiferromagnetic phase with symmetry breaking along x- or y-direction (x-AFP or y-AFP). When the coupling strengths exceed the critical points $[eqn]$ for $[eqn]$ , the discrete $[eqn]$ symmetry breaking occurs in d-direction ( $[eqn]$ ), depending on the larger one of $[eqn]$ and $[eqn]$ . It results in a double degeneracy with spontaneous magnetization in x- or y-directions. On the other hand, when $[eqn]$ , the continuous U(1) symmetry is spontaneously broken beyond the critical points, leading to the emergence of infinite degeneracy associated with the Goldstone mode.

To investigate quantum criticality and fluctuations beyond the mean-field ansatz, we retain and solve the quadratic term in the transformed Hamiltonian using a sequence of transformations. Regarding the $[eqn]$ symmetry breaking, we present evidence of the second-order quantum phase transitions, specifically from PP to FP and AFP. At the critical points, the excitation energy vanishes, and both quantum fluctuations and entanglement entropy diverge. Furthermore, the analysis of the critical behaviors near the phase transition points confirms the breaking of $[eqn]$ symmetry. In contrast, when the system undergoes U(1) symmetry breaking, the energy gap vanishes beyond the critical points, yielding a novel exponential relation.

Our approach can be easily extended to systems with multiple large spins, where both anisotropy and geometric frustration are of substantial significance. The competition between them is likely to give rise to intricate quantum phase transitions and critical phenomena, rendering it worthy of further exploration.