Phase and Amplitude Modes in the Anisotropic Dicke Model with Matter Interactions

Phase and amplitude modes, also called polariton modes, are emergent phenomena that manifest across diverse physical systems, from condensed matter and particle physics to quantum optics. We study their behavior in an anisotropic Dicke model that includes collective matter interactions. We study the low-lying spectrum in the thermodynamic limit via the Holstein–Primakoff transformation and contrast the results with the semi-classical energy surface obtained via coherent states. We also explore the geometric phase for both boson and spin contours in the parameter space as a function of the phases in the system. We unveil novel phenomena due to the unique critical features provided by the interplay between the anisotropy and matter interactions. We expect our results to serve the observation of phase and amplitude modes in current quantum information platforms.


I. INTRODUCTION
Amplitude and phase modes are collective excitations of an order parameter usually studied in connection to a continuous symmetry breaking.For the U (1) symmetry, they are known as the massive Anderson-Higgs (AH) [1,2] and the massless Nambu-Goldstone (NG) modes [3,4], respectively, and are present in several systems, from condensed matter [5][6][7][8][9], quantum antiferromagnets [10][11][12] to particle physics [13][14][15].Phase and amplitude modes help to describe critical phenomena such as quantum phase transitions (QPT) [3], i.e., a parametric sudden change of the ground-state's properties of a quantum system [16,17].The paradigmatic case is the quartic potential V (α) = rα 2 + gα 4 expressed in terms of the scalar complex parameter α, that possesses the U (1) continuous symmetry [16].The parameter space is separated into two phases: the normal or non-ordered phase (r, g > 0), where the potential has the form of a single well with a minimum at α = 0, and the ordered phase (r < 0), where the potential has a Mexican hat form.The ground state is a degeneration of states with |α| > 0. In the ordered phase, fluctuations of the parameter quantify the two possible excitations: 1) the phase mode (NG), which occurs at constant amplitude and has vanishing excitation energy, thus massless; 2) the amplitude mode (AH), which has constant phase and possess finite excitation energy, thus being massive.In general, both modes are hard to observe experimentally.The NG is particularly sensitive to symmetry deviations, as it gains finite excitation energy [18,19].Likewise, the AH mode is challenging to observe because of its massive character.Ultracold atoms setups are particularly suitable for their observation and control due to their high tunability [19][20][21][22][23][24].
Under the rotating-wave approximation, the Dicke model becomes the Tavis-Cummings (TC) model [77], which is integrable in a classical sense, conserving the total number of excitations.This is reflected in the semiclassical limit, where a Mexican hat potential emerges in the superradiant phase [78].In this symmetry-broken phase, many states cluster together to form a quasicontinuous band of energies, where the phase can take any value in [0, 2π).Hence, there is no cost of energy to displace the state or change the order parameter along the potential's minima.Because of this, the TC model exhibits both a gapless phase mode and a gaped amplitude mode.Instead, the Dicke model has only a Z 2 symmetry.As a result, it has a doubly degenerate spectrum in the superradiant phase corresponding to a double-well semi-classical potential [78].The pair of arXiv:2406.09446v1[quant-ph] 12 Jun 2024 (symmetry-breaking) ground-states correspond to those at the bottom of a potential, reflecting two possible values for the phase (0 and π).Hence, the phase mode becomes gaped (Higgs-like) and corresponds to what has been called a roton-type mode [79][80][81] in analogy to the roton excitations in superfluid helium.The phase and amplitude modes in the Dicke model are called polariton modes [82], as well.A polariton is a hybrid quantum state emerging from the strong light-matter interactions [83].Because the phase and amplitude modes in the Dicke model represent the collective light-matter quantum state at low energy, they correspond to effective upper-and lower-polaritons [84,85].The differences between the Dicke and TC models can be identified by studying the gap between the ground and first excited states.Its has been studied for finite-size [19,86], and in the thermodynamical limit [35,87].Also, the phase and amplitude modes have been observed for large lightmatter coupling [81,88].A way to do it is by exploiting their singular behavior at the superradiant QPT [18].The symmetry-breaking effect in the Dicke model has been measured experimentally [27,70], and employed for quantum sensing [89].
The so-called anisotropic, generalized or unbalanced Dicke model [18,38,78,[90][91][92] is a useful tool to explore the parametric change of the gap and understanding the polariton modes as one passes from a continuous to a discrete symmetry, where the phase model can gain mass.There, one assumes a relative strength between the rotating (excitation conserving) and counter-rotating terms, allowing to tune the Hamiltonian between the TC and Dicke limits.Both the equilibrium [93][94][95][96] and non-equilibrium properties [72,[97][98][99][100] of the anisotropic Dicke model has been studied, exhibiting a rich phase space and the presence of novel critical phenomena.Phase modes have been found in other extensions of the Dicke model, such as generalizations of the spinboson coupling [101], Jahn-Teller-Dicke models [102], two-mode Dicke models [88,[103][104][105], Jaynes-Cummings lattices [106], the three-level Dicke model [107,108], and also including matter-matter interactions [19].
In light of experimental progress in the combination of strong matter-matter interactions and ultra-strong light-matter systems, such as the recent realization of the Dicke model in a solid-state system with tunable spin-magnon interactions [109][110][111], in this work, we investigate the effects of collective matter interactions over the phase and amplitude or polariton modes in the anisotropic Dicke model.Collective matter interactions mediated by light introduce unique behavior in the Dicke model such as first-order QPT [112][113][114][115][116][117][118], parametric shifts over the superradiant QPT [54,119,120], new quantum phases [117,[121][122][123], and novel perspectives on quantum chaos [117,124], such as the amplification of regularity-to-chaos transition [125].
To achieve our end, we employ the standard approach of applying the Holstein-Primakoff (HP) transformation to obtain a low-energy approximation of the excitation modes around the minima [126].Given the limitations of the HP transformation [127], we interpret the results under the light of the semiclassical corresponding energy surfaces that can be obtained via coherent states [122].Additionally, we calculate the geometric or Berry phase around the critical points and relate its behavior to the phase and amplitude modes [128,129].The geometric phase possess a non-analytic behavior along a QPT, so it has been employed as a signature of criticality [130][131][132].It has been observed experimentally in Heisenberg chains [133] and qubit systems [134,135].It has been studied for the standard Dicke model [136][137][138] and extensions such as including two-photon procesess [139], impurity coupling [140], and collective matter interactions in the z direction [141].
The article is organized as follows.In Sec.II, we present the anisotropic Dicke Hamiltonian with matter interactions and details over its classical limit.Next, in Sec.III, we obtain exact expressions for the low-lying polariton branches using the Holstein-Primakoff approximation.In Sec.IV, we discuss the gap behavior as a function of the different phases in the system.We calculate the geometrical phase for arbitrary circulations generated by the photon number or the relative population operator in Sec.V. Finally, in Sec.VI, we present our conclusions and perspectives.

II. ANISOTROPIC DICKE HAMILTONIAN WITH MATTER INTERACTIONS
The anisotropic Dicke Hamiltonian including collective matter interactions is given by The first two terms denote the non-interacting parts of the Hamiltonian, where ω is the boson frequency, ω 0 is the qubit energy splitting, â † (â) is the creation (annihilation) boson operator, Ĵz is the relative population of excited qubits, given that Ĵz,x,y are pseudo-spin operators representing the collective degrees of freedom of the set of N qubits, which follow the rules of the su(2)-algebra.The third term is the spin-boson interaction, with γ the coupling strength and ξ the anisotropy parameter, that allows tuning the relative strength of the counter-rotating term, so ξ = 0 and ξ = 1 are the TC and Dicke limits, respectively.The last term accounts for the matter interactions, where η i with i = x, y, z are the couplings in each direction, whose meaning depends on the particular setup realizing the Dicke model, including Josephson dynamics [115,142,143], atomic dipolar couplings [119,144], optomechanical setups [145,146], or interactions between superconducting qubits [54,56,57,147,148].A standard approach is regarding η z as the strength of collective onsite interaction and η x (η y ) as a strength of the collective inter-qubit interactions.It is worth emphasizing that the relevant interacting parameters are ∆η zµ = η z − η µ , with µ = x, y, given that the pseudospin length is conserved, and one can express one direction in terms of the others Ĵ2 z = j(j + 1) Î − Ĵ2 x − Ĵ2 y [122].The Hamiltonian in Eq. 1 commutes with the squared pseudospin length operator, [ Ĥξ , Ĵ2 ], so the collective ground-state lies in the totally symmetric subspace that corresponds to j = N/2, and the Hilbert space of the system is effectively reduced to N + 1 states.As a result, a classical corresponding Hamiltonian can be obtained by mapping the boson and collective pseudospin degrees of freedom to the classical limit via coherent states.To do so, one takes the value of the quantum Hamiltonian over the tensor product of Glauber |z⟩ and Bloch |w⟩ coherent states as trial states [38,40,78,149,150], where |0⟩ and |j, −j⟩ are the boson and pseudo-spin vacuum states, respectively [151].In canonical classical variables for the boson (q, p), z = j/2 (q + ip) and the pseudospin (j z , ϕ), w = (1 + j z )/(1 − j z )e −iϕ the classical corresponding Hamiltonian reads [122] A standard method to determine quantum phases in the system is to analyze the extreme points of the classical energy surface in Eq. 2, i.e., finding fixed points of the corresponding Hamilton equations because they can be related to critical changes in the systems.Hence, absolute minima signal ground-state QPT and other extreme points to excited-state QPT [78,93,117].In the following, we will focus only on the absolute minima of the energy surfaces, as we are interested in polariton modes.We employ the classical energy surfaces to interpret the behavior of these modes across the parameter space.We picture the surfaces in a set of variables (u, v) = arccos(−j z )(cos ϕ, sin ϕ) that, by using Hamilton equations to eliminate the bosonic quadratures q and p, allows one to visualize the fixed points in the collective pseudospin space alone [78,122].In this picture, the energy surface reads as where f ξ± = γ 2 /γ c ξ± , and γ c ξ± = √ ωω 0 /(1 ± ξ), are the critical couplings for the superradiant-(±) phases as detailed below.
The anisotropy introduces new quantum phases to the standard Dicke model that matter-matter interactions then modify.For zero matter interactions (η I = 0), the Dicke (ξ = 1) and TC (ξ = 0) models exhibit only two phases, the normal (non-ordered) and the superradiant (ordered), separated by the critical coupling γ c ξ+ .The anisotropy modifies the energy landscape such that two different superimposed phases appear, the superradiant-(+) and superradiant-(−) phases, the latter with a critical coupling γ c ξ− = [44,93,94].These phases are also called the electric and magnetic superradiant, with a different coupling combination for the counter and counterrotating terms [18,92].Each phase has its extreme points, but only one predominates as the ground-state, depending on which of the two terms, the rotating or the counter-rotating, is stronger.For our case of study ξ ∈ [0, 1], the superradiant-x phase defines the groundstate.

III. PHASE AND AMPLITUDE MODES IN THE ANISOTROPIC DICKE MODEL
The Holstein-Primakoff transformation (HPT) [164] provides a quadratic approximation to the lower energy modes of Eq. 1 in the thermodynamic limit (N → ∞).After diagonalizing the quadratic, approximated Hamiltonian, the resulting branches correspond to the polariton modes or mean-field low-lying excitations of the photon-matter quantum superposition [18,80].The HPT reads where b ( b † ) is an annhilation (creation) boson operator, such that [ b, b † ] = I.First, by substituting Eqs. 4 one obtains Next, to study the excitation modes we employ a mean field approximation by displacing both â and b bosonic operators where α, β ∈ C are scaled displacements.The square roots in the HPT become Next, we make use of the thermodynamic limit N → ∞.
Here, the customary assumption is that ⟨ b †b ⟩ ≪ 2j, so By keeping only terms up to quadratic order, the Hamiltonian becomes A. Deformed Normal Phases The deformed normal phases are characterized by a zero expectation value of photon and matter excitations.Thus, by setting α, β → 0, we obtain Next, one writes the Hamiltonian in terms of the boson quadratures [126] By substituting Eqs. 10, the Hamiltonian is now represented as follows: where We considering the following dimensionless frequencies containing the matter interaction strengths ωzx = 1 − ∆η zx ω 0 , and ωzy = 1 − ∆η zy ω 0 .
Next, two rotations are performed to eliminate the cross terms xŷ and px py .We will apply a rotation with angle θ 1 for the variables x and ŷ, and another with angle θ 2 for px and py.
x ŷ cos θ px py The following conditions arise: and tan 2θ 2 = 2f where f ξx = (γ/γ ξx ) 2 , and the critical coupling indicating the onset of superradiance is Substituting these conditions yields a two-mode decoupled Hamiltonian, determining the energies ϵ N 1± that represent the low-lying excitation modes.
with energies are where g ξy = (γ/γ ξy ) 2 , with the critical coupling of the superradiant-y behavior.It is convenient to rewrite Eq. 17 in a second quantized form in order to obtain a decoupled Hamiltonian with new bosonic excitations âi ( where , and ω N+ = ϵ N 1+ /ϵ N 2+ .Thus, we obtain the quantum oscillation modes and the low-energy spectrum.It is essential to highlight that these expressions depend on the ∆η xz and ∆η yz parameters, allowing us to modulate the behavior of the modes both in the normal phase and in the superradiant phase as we will see later on.Finally, the Hamiltonian reads with The phase ϵ N − and amplitude ϵ N + modes in the deformed normal phases become B. Superradiant phases Now, we consider the case where α, β ̸ = 0. To obtain a quadratic Hamiltonian, we need to eliminate the linear terms in ĉ, ĉ † , and d, d † from Eq. III.Hence, we obtain the following conditions determining the values of α and β, ωα − γ Solving the system of equations leads to where The critical coupling γ c ξx in Eq. 16 for the existence of superradiant behavior is determined by µ x = 1, because for µ x > 1 α and β in Eq. 28 become complex.Writing the resulting Hamiltonian in terms of Eq. 28, we have: where we have defined a set of reduced variables Just as we did in the deformed normal phases, we will proceed to diagonalize the Hamiltonian in Eq.III B, by transforming it into quadratures and then applying rotations that decoupled them.This transformation reads so the Hamiltonian is expressed as where Next, we consider the following variable changes Again, we apply two rotations: ϕ 1 for x and y, and ϕ 2 for p x and p y .
x ŷ cos ϕ px py We obtain the following conditions for the angles ϕ 1 and ϕ 2 from eliminating the crossed terms q1q2 and p1p2, .
The Hamiltonian reduces to: We decouple the Hamiltonian using quadratures.
where, as in the deformed normal case, Finally, the Hamiltonian reads in second quantized form as with with the phase ϵ S − and amplitude ϵ S + modes being C. Deformed phase Here, we comment on the deformed phase in the Dicke limit (ξ = 1) that suppresses the superradiant-y phase [117,122].It exists independently of the spinboson coupling γ, while ∆η zy ≤ ω 0 , and it is characterized by zero expectation value of the photon number but is not a normal phase.This can be recognized in the energy surface by a π rotation of the extreme points.If we take ξ = 1 in Eqs.18 and 40, g ξy → 0, and we get Because ϵ S,N ± = ϵ S,N 2± ϵ S,N 1± , the phase mode becomes vanishing, effectively suppressed, and the amplitude mode is imaginary for ∆η zx ≥ 1, so it is not defined in the deformed phase.

IV. ROLE OF MATTER-MATTER INTERACTIONS
Next, we explore the behavior of the phase and amplitude modes as a function of the matter-matter interaction parameters for the Dicke (ξ = 1), TC (ξ = 0) and an anisotropic case, particularly ξ = 0.5 as an example.

A. Absence of matter-matter interactions
As shown in Fig. 1, first, we revisit the polariton modes without matter interactions.In Fig. 1 (a1) we show the polariton modes in the Tavis-Cummings limit.There, the lower curve corresponds to the phase mode ϵ N − , while the upper one to the amplitude mode ϵ N + .In the normal phase, the parametric evolution of both modes increases proportional to the light-matter coupling.Their energy gap is just the Rabi splitting 2γ.This agrees with the overall behavior of the energy surface, which adopts a spherical well shape, a reflection of the conserved U (1) symmetry [see Fig. 1 (a2)].The dashed, vertical line marks the critical coupling γ 0+ , separating the normal from the superradiant phases.There, the energy surface experiences a sudden change, taking the shape of a Mexican hat (see [Fig. 1 (a3-a4)].Because of the onset of degeneration, we have spontaneous symmetry breaking.Once we enter the superradiant phase, the phase mode ϵ S − becomes a Goldstone mode with zero energy, depicted on the energy surface as the minimum energy ring of the Mexican hat [22,23,88,106,165,166].Instead, the amplitude mode ϵ S + depends on the brim of the Mexican hat.As it grows, it requires increasingly higher energy costs than the phase modes, which remain constant.Note that, in this case γ 0x = γ 0y = γ 0− , so there is no presence of the superradiant-(−) effects.In polariton terms, the two modes represent a maximal light-matter hybridization in resonance (ω = ω 0 ).The amplitude mode, being massive, can be identified as an upper-polariton branch that contains the most matter content, while the phase mode, being massless, becomes the lower-polariton branch with the most photonic content.Here, entering the superradiant phase implies a vanishment of the lower-polariton energy.
The energy spectra on Fig. 1 (b1) and (c1) correspond to the standard Dicke model (ξ = 1) and the anisotropic case when ξ = 0.5, respectively.Both cases exhibit quite similar behavior.In the normal phase, the behavior resembles the TC limit because the system is in the stronginteracting regime.At the critical coupling γ c ξx = γ c ξ+ , however, the phase mode tends to zero, while the amplitude mode to a fixed value given by (1 , which in the absence of interactions becomes and for the standard Dicke model is [80].When one passes the QPT, the phase mode gains a finite energy, becoming the rotonic mode [23,88,165], which converges to the finite value for larger couplings.This value comes from the energy cost to pass between the two degenerate minima (whose phases are 0 and π) observed in Figs. 1 (b3) and (b4).The amplitude mode behaves quite similarly to the TC since the center of the energy surface presents an unstable point and tends quadratically as a function of γ, ϵ γ≫1 + = (1 + ξ) 2 γ 2 /ω in the absence of matter interactions, which has the value 4γ 2 /ω in the standard Dicke model [80].It has been shown that increasing the anisotropy decreases the amplitude mode gap till one obtains the Goldstone mode.This can be seen from Eq. 48, where the limit ξ → 0 leads to ω returning to the non-interacting value of photon field for the isotropic model [80].This critical feature has been proposed as an experimental signature for detecting the Goldstone mode [18,22].We observe this feature in Fig. 1 (c1), and can be explained in terms of the change of the energy surface that leads to a new set of extremal points related to the superradiant-(−) phase that reduces the energy of the phase mode [122].In this case, the presence of the counter-rotating terms that don't conserve the excitation number is reflected in the finite energy of the phase mode as a lower-polariton mode where light gains "mass" close to the QPT and then becomes a massless but highly correlated state.In polariton terms, the lower-polariton energy becomes gaped.We notice, however, that the criticality of these cases is not affected by the superradiant-(−) extreme points because, in the absence of interactions, they cannot be shifted, as it happens below.For the Dicke model γ c 1− → ∞, and for the anisotropic case with ξ = 0.5, the critical coupling γ c 0.5− = 3γ c 0.5+ is deep in the superradiant phase and not shown in Fig. 1 (c1).

B. Presence of matter-matter interactions
Next, we explore the effect of matter interactions over the polariton modes.In general, the z-interactions alone have the role of shifting both the critical coupling and the energy spectrum.To simplify the analysis we will focus on the role of the x and y interactions via the relevant parameters ∆η zx and ∆η zy , respectively.
In Fig. 2 we show a representative value of the interactions given by η y = 0.9ω 0 (∆η zy = −0.9ω0 ), keeping η x = 0.Because of the interactions, in the normal phase, the energy surface widens along the u direction for all ξ, so we call it the deformed normal phase [122] [See Figs. 2 (a2), (b2), and (c2)].In this case, the deformation does not substantially affect the phase and amplitude modes in the normal phase, except for an energy gap of the am-plitude mode for zero spin-boson coupling given by ϵ N + = ω 0 ωzx ωzy .(49) and representing the cost to climb the well in presence of interactions.Next, we cross the superradiant QPT.
For the three cases under study (ξ = 0, 0.5, 1), the amplitude mode is not significantly changed by the interactions compared to the case without.However, this is not the case for the phase mode in the TC limit (ξ = 0).
As the system enters the superradiant phase, the excitation number is not conserved anymore because of the deformation, so we do not have the Mexican hat potential anymore.As shown in Fig. 2 (a1), the phase mode gains energy, becoming rotonic.As we increase coupling, the onset of new extreme points at higher energies contributes to somehow restoring the original symmetry of the Mexican hat potential in the absence of interactions and reducing the mode's energy.Thus, as the spinboson coupling increases, the phase mode slowly tends to a Goldstone-like phase mode.As a result, interactions produce a "variable" mass roton mode.In both the anisotropic and Dicke limit cases, the results are similar to the case without interactions.After transitioning to the superradiant phase, two stable degenerate minimum points form, so there is no mechanism capable of decreasing the energy of the phase mode and eventually converging to a general energy value given by Eq. 48.
The situation is different for x-interactions, as shown in Fig. 2 for a representative value η x = 0.9 (∆η zx = −0.9).We observe similar dynamics for the upper and lower-polariton modes as in the previous case.However, new phenomena emerge from the interplay between the anisotropy and matter interactions, i.e., the ability of the latter to shift the position of the superradiant-(−) phenomena.Again, the normal phase becomes deformed for every ξ, but now with a different orientation, alongside the v direction corresponding to the change from y-to xinteractions [122] [See Figs. 3 (a2), (b2), and (c2)].While the overall behavior of the amplitude mode remains the same, in the TC limit [Fig.3 (a1)], we observe that the phase mode vanishes in the normal phase before the onset of the critical coupling.This is not a Goldstone mode but a suppression due to matter-matter interactions.As it can be seen from Eqs. 18, two conditions make the phase mode vanish, either f ξx = 1 or g ξy = 1.As it can be shown in Fig. 3 (a1), the relative position of the critical coupling γ c 0x with respect to γ c 0y is shifted thanks to η x ̸ = 0, creating a coupling regime where the phase mode cannot exist, although the amplitude mode does.This can be seen from the relation between the two critical couplings and it occurs when ∆γ < 0. This unique effect results from the influence of the matter interactions over the competence between the superradiant-+ and superradiant-y phases resulting from the anisotropy.This behavior does not occur in the Dicke and anisotropic cases [See Figs. 3 (a1) and (c1)], provided the onset of the superradiant-y behavior occurs deep in the superradiantx one.However, for the anisotropic case, the interplay between the anisotropy and the matter interactions might shift the critical coupling γ c ξy = to suppress the mode.Beyond that, the Dicke and anisotropic cases, the behavior of the polariton modes remain qualitatively similar to previous cases, except for the π rotation of the energy surfaces as one crosses the superradiant QPT, as it can be seen in Figs. 3 (b2) and (b3), as well as in Figs. 3 (c2) and (c3).

V. GEOMETRIC PHASE IN PRESENCE OF MATTER INTERACTIONS
Next, we investigate the influence of matter interactions over the Berry or geometric phase induced by the bosonic field and the collective pseudospin over both polariton modes.The geometric phase is a topological quantum phase that the eigenstates of a Hamiltonian acquire when one goes adiabatically over a close path in the system's parameter space, accounting for the Hilbert space geometry [128,129].The geometric phase has already been studied for the Dicke model [136,137], extensions [139,140], and under the presence of matter interactions in the z-direction [141].Mainly because it becomes singular at criticality, it has been shown to serve as a signature for the superradiant QPT.Also, matter interactions could modify its scaling behavior [141].
First, for comparison, we obtain the approximated ground-state energy without excitations.According to Eq. 9, the scaled energy reads It agrees to that calculated using coherent states as trial states [122].Then, by substituting Eq. 28 into Eq.51, we obtain the ground-state energy for both the normal and superradiant phases where |Ψ 0 ⟩ is the ground-state of the system, and µ x (∆η zx , ξ) = (∆η zx /ω 0 + f ξ+ ) −1 , with f ξ+ = γ 2 /γ 2 ξ+ .The ground-state energy signals the superradiant QPT modified by the anisotropy and matter interactions at γ c ξx (µ x = 1).
To calculate the geometric phase, we follow the standard procedure [137,138,[167][168][169][170], where one considers a unitary transformation producing a cyclic trajectory over the parameter space.It may be an arbitrary adiabatic circulation generated by the photon number [139] or over the Bloch sphere [170].First, let us consider an adiabatic circulation generated by the photon number It provides an additional phase to the bosonic annihilation (creation) operator â → âe iϕn (â † → â † e −iϕn ), effectively displacing them.Next, we introduce a timedependent unitary transformation, where the parameter ϕ n (t) = ω n t varies adiabatically in the interval [0, 2π) with angular frequency ω n .The geometric phase is calculated in terms of the circuit integral of the Berry connection where ⟨Ψ 0 |n|Ψ 0 ⟩ = |α| 2 is just the expectation value of the photon number within the HPA.Hence, the geometric phase reads For ξ = 1, we recover the results of the standard [137] and modified Dicke models [141].Similarly, if now we We notice that, for these contours, both the Γ n and Γ m geometric phases depend only in the ratio f ξx , as from Eq. , µ −1 x = ωzx /ω 0 (f ξx − 1) + 1.As a result, the behavior is qualitatively independent from ξ and η y .That it does not depend on η y is because the extreme points in the energy surface belonging to the superradiant-y occur higher in energy, far away from the ground-state energy.Instead, the independence from ξ is a feature that captures the universal behavior of the superradiant QPT.In Fig. 4, we show the results for both geometric phases Γ n and Γ m and compare them to the scaled ground-state energy.We observe that the geometric phase identifies the position of the superradiant QPT, as expected.The geometric phases are sensitive to the sign of η x increasing with η x .Furthermore, we observe that the geometric phase also helps to identify the first-order phase transition occurring at ∆η x = ω 0 , where the geometric phase becomes zero independently of γ in both cases and then turns negative.

VI. DISCUSSION AND CONCLUSIONS
We have derived analytic expressions for low-energy excitations of the anisotropic Dicke model in the presence of matter interactions by employing the Holstein-Primakoff transformation and neglecting higher-order terms in the thermodynamic limit.We explored the modes in the normal (ordered) phase, deformed by matter interactions, and the superradiant (ordered) phase from a unified point of view, whose interpretation is supported by the analysis of energy surfaces obtained using the corresponding classical Dicke model via coherent states [122].
n N e f a M 7 J / h y S P 3 I + v g E y L I 5 G < / l a t e x i t > (b) < l a t e x i t s h a 1 _ b a s e 6 4 = " J h / 0 z I G o P 8 9 2 l A T v H L s g q Beyond the shifting of critical couplings and groundstate energy associated to matter-matter interactons, we find unique effects from combining them with the anisotropy.In most cases, amplitude, massive, or upper-polariton mode behave qualitatively the same, just changing energetically to reflect the deformations of the energy surface, so interactions increase the gap to the phase mode, creating an energy cost given by ω 0 ωzx ωzy .Novel behavior appears mostly over the ground-state mode, i.e., the phase or lower-polariton mode.For the Dicke model, it becomes a roton mode that, for larger couplings, tends to the non-interacting photon energy ω and decreases in energy with the anisotropy.Without interactions it eventually transforms into the Goldstone mode.Instead, the TC limit exhibits a Goldstone mode in the absence of matter-matter interactions that becomes rotonic as the latter breaks the U (2) symmetry, leaving conserved only the Z 2 one.
A major result is that tuning the combination of matter interactions in the z and y directions can reduce the rotonic mode energy in the TC Hamiltonian, producing a "variable" mass roton mode, or a lower-polariton that starts losing matter content with interactions.This is understood in terms of the energy surface, as the interactions influence on higher energies becomes less relevant to the ground-state, thus restoring to some extent the Goldstone-like behavior.Second, we discover that the Goldstone mode can be suppressed by tuning the interactions in the z and x directions by shifting the relative position between the critical couplings of the superradiantx and -y phases.This suppression does not affect the amplitude mode.
We anticipate these features could be detected and employed as an experimental signature to observe the Goldstone mode in novel setups, such as spin-magnon systems, where the interactions are tunable [109][110][111], as well as in others within the broad range of setups where the Dicke model can be realized, from quantum optics and atomic physics to condensed matter.
We also calculated the geometric phase for both boson and pseudospin contours, which equals the expectation value of the photon number and excited pseudospin in the phase mode within this approximation.We find that the role of the z and x-interactions is to shift the phase.Still, it is independent of the y interactions and the anisotropy parameter ξ, as it is a low-energy approach capable only of capturing the universal behavior of the superradiant-x transition.Furthermore, the geometric phase becomes singular at the superradiant QPT, serving as a signature of criticality, where we recover previous results from Refs.[119,141].It also detects the first-order QPT that results from changing matter-matter interactions by changing sign.
We conjecture that the predominance of the superradiant-x phase will be inverted to the superradiant-y phase in the case of the Anti-Tavis-Cummings model where now the rotating terms become smaller to the counter-rotating ones [18,92,171].Our model can be extended to explain this case but goes beyond the scope of this work.Finally, exploring the effect of the anisotropy and matter interactions over the polariton modes beyond mean-field is left for future work [82,172].

( 56 )
Like in the boson case, this shifts the bosonic annihilation (creation) operator b → be iϕm ( b † → b † e −iϕm ), effectively displacing them.Changing adiabatically in time allows us to calculate the geometric phase as
The anisotropy in the Dicke model that, in this work, tunes the Hamiltonian between the Tavis-Cummings to the standard Dicke models produces a new set of extreme points and quantum phases, the superradiant-(+) and (−) phases, where superradiance is determined by the (+) phase given that it lies lower in energy.Matter interactions serve two purposes: to shift the values of critical couplings and ground-state energy and to produce new extreme points in the energy surface, reflecting the onset of new quantum phases.Particularly, interactions in the x and y directions directly modify the criticality of the (+) and (−) phases, respectively, and shift their relative position in terms of the spin-boson coupling.
FIG. 4. (a) Scaled ground-state energy ϵ0 = E0/2j as a function of the light-matter coupling.(b) Scaled geometric phase Γn considering an adiabatic circulation generated by the photon number.(c) Scaled geometric phase Γm considering an adiabatic circulation generated by the collective pseudoespin z-projection.As indicated in the figures, we consider different values of the matter interactions ∆ηzx with ηz = 0.The case with ∆ηzx = 0 is indicated as a dashed line.All cases are calculated in resonance ω = ω0 = 1.