---
title: 'Anisotropic Dicke Model: Key Aspects'
url: https://www.emergentmind.com/topics/anisotropic-dicke-model
type: topic
---

# Anisotropic Dicke Model: Key Aspects

Searching arXiv for recent and foundational papers on the anisotropic Dicke model to ground the article in the provided literature.
The anisotropic Dicke model is a collective light–matter Hamiltonian in which a single bosonic mode is coupled to an ensemble of \(N=2j\) two-level systems with unequal rotating-wave and counter-rotating-wave amplitudes. In its most common form, the anisotropy is the independent weighting of the processes \(a^\dagger J_-+aJ_+\) and \(a^\dagger J_++aJ_-\), so that the model interpolates between excitation-conserving \(U(1)\) limits and the parity-symmetric \(\mathbb Z_2\) Dicke regime. This asymmetry underlies the model’s ground-state normal-to-superradiant transition, its excited-state and ergodic structures, its driven and open-system generalizations, and its use as a platform for critical metrology, collective-mode spectroscopy, and quantum transport [1610.06751], [2304.07857], [2406.06301].

## 1. Hamiltonian, notation, and relation to the standard Dicke model

A widely used definition of the anisotropic Dicke model is
$$
H=\omega a^\dagger a+\omega_0 J_z+\frac{g_1}{\sqrt{2j}}\left(a^\dagger J_-+aJ_+\right)+\frac{g_2}{\sqrt{2j}}\left(a^\dagger J_++aJ_-\right),
$$
where \(a,a^\dagger\) are bosonic operators, \(J_z,J_\pm\) are collective spin operators, \(\omega\) is the boson frequency, \(\omega_0\) is the atomic splitting, and \(g_1,g_2\) are the two independent couplings [1610.06751], [2008.10253], [2304.07857], [2505.21169]. The model becomes the standard Dicke model when \(g_1=g_2\), and it reaches integrable limits when either coupling vanishes [1610.06751], [2008.10253].

Different subliteratures use different anisotropy conventions. In the Dicke–Stark formulation, the rotating-wave coupling is \(g\) and the counter-rotating-wave coupling is \(g\tau\), so the isotropic limit is \(\tau=1\) [2405.19776]. In the driven Rydberg-array realization, the normalized anisotropy parameter is
$$
\alpha=\frac{g_{\rm CRW}}{g_{\rm RW}+g_{\rm CRW}},
$$
with \(g_{\rm RW}=(1-\alpha)g\) and \(g_{\rm CRW}=\alpha g\), giving \(\alpha=0\) for pure rotating-wave coupling, \(\alpha=1/2\) for the isotropic Dicke model, and \(\alpha=1\) for pure counter-rotating-wave coupling [2511.22230]. In the matter-interacting model, \(\xi=0\) is the Tavis–Cummings limit and \(\xi=1\) is the Dicke limit [2406.09446]. The symbol \(\gamma\) is not universal: it denotes the ratio \(\lambda_1/\lambda_2\) in the quantum-geometric analysis, but the anisotropy factor multiplying counter-rotating terms in the transport study [2406.06301], [2603.29180]. This suggests that the model is most robustly identified by the structural inequality of rotating-wave and counter-rotating-wave couplings, rather than by any single parameter symbol.

| Convention | Definition | Special values |
|---|---|---|
| \(g_1,g_2\) | independent RW and CRW couplings | \(g_1=g_2\): Dicke; \(g_1=0\) or \(g_2=0\): integrable |
| \(\tau\) | RW \(=g\), CRW \(=g\tau\) | \(\tau=1\): isotropic |
| \(\alpha\) | \(g_{\rm RW}=(1-\alpha)g\), \(g_{\rm CRW}=\alpha g\) | \(\alpha=0\): RW only; \(\alpha=1/2\): isotropic; \(\alpha=1\): CRW only |
| \(\xi\) | coefficient of the CRW channel in the matter-interacting model | \(\xi=0\): TC; \(\xi=1\): Dicke |

A recurrent benchmark for anisotropic studies is the standard isotropic Dicke model. The finite-\(N\) catastrophe-theory, symmetry-adapted-state, and fidelity analyses developed for that case are methodologically important, but that benchmark uses a single coupling multiplying co- and counter-rotating terms equally, and no anisotropy parameter appears [1211.6692].

## 2. Symmetry, integrable limits, and the basic phase structure

The generic anisotropic Dicke model has a parity symmetry
$$
\Pi=\exp\!\big(i\pi[a^\dagger a+J_z+j]\big),\qquad [H,\Pi]=0,
$$
and many finite-size studies restrict the analysis to the positive-parity sector, which contains the ground state in the parameter ranges considered [1610.06751], [2008.10253], [2304.07857]. In the generalized Dicke–Ising realization, the same discrete symmetry is written as
$$
P=(-1)^{a^\dagger a}\prod_{i=1}^{N}\sigma_i^z,
$$
while at \(\alpha=0\) the model acquires an enhanced \(U(1)\) symmetry with conserved total excitation number [2511.22230]. In the generalized Dicke model with couplings \(g\) and \(J\), the same distinction is phrased as a crossover between a generic \(\mathbb Z_2\) theory and the \(U(1)\) Tavis–Cummings and anti-Tavis–Cummings limits [1909.10252].

In the thermodynamic limit, the normal-to-superradiant quantum phase transition occurs at
$$
g_1+g_2=\sqrt{\omega\omega_0},
$$
with the normal phase for \(g_1+g_2<\sqrt{\omega\omega_0}\) and the superradiant phase for \(g_1+g_2>\sqrt{\omega\omega_0}\) [1610.06751], [2008.10253], [2304.07857]. For the common resonant choice \(\omega=\omega_0=1\), the critical line becomes \(g_1+g_2=1\), and along the Dicke line \(g_1=g_2=g\) the critical point is \(g_c=0.5\) [2008.10253]. The order parameter is usually the boson occupation: in the normal phase \(\langle a^\dagger a\rangle/j=0\), whereas in the superradiant phase \(\langle a^\dagger a\rangle/j=\mathcal O(1)\) [1610.06751], [2008.10253].

The integrable limits are structurally important. When \(g_2=0\), the model reduces to the rotating-wave approximation, conserving the total excitation number \(n+m+j\); when \(g_1=0\), it maps to the \(g_2=0\) case by a spin rotation and \(\omega_0\to-\omega_0\) [1610.06751]. These limits organize both the model’s symmetry classification and its later ergodic and transport behavior.

## 3. Critical phenomena, finite temperature, and finite-size structure

The ground-state quantum phase transition retains the Dicke-type normal-to-superradiant structure, but the anisotropic model supports additional critical phenomena. In the finite-size spectral analysis, the ground state is localized in the normal phase with inverse participation ratio near unity, while in the superradiant phase it becomes multifractal with \(P_{\rm gs}\sim \sqrt{N_D}\), so the transition is described as a localized-to-multifractal transition of the ground state [2304.07857]. The same work identifies an excited-state quantum phase transition through two characteristic energies, \(E_{\text{lower}}\) and \(E_{\text{upper}}\), extracted from jumps in von Neumann entanglement entropy; the states between these energies form a central band distinct from lower and upper spectral bands [2304.07857].

At finite temperature and finite \(N\), the sharp superradiant transition broadens into a fluctuational critical region of width
$$
\Delta=\sqrt{\frac{\omega T}{N}},
$$
and the internal structure of that region depends on whether the model is effectively in a \(U(1)\) or \(\mathbb Z_2\) regime [1909.10252]. In the Tavis–Cummings-like sector, the universal relative fluctuation is
$$
r_{\rm TC}=\frac{\pi}{2}-1\approx 0.57080,
$$
whereas in the generalized Dicke regime it is
$$
r_{\rm GD}=4\frac{\Gamma^2(5/4)}{\Gamma^2(3/4)}-1\approx 1.18844.
$$
The same study finds distinct temperature scalings of momentum squeezing,
$$
\delta p_{\rm TC}\propto T^{1/4},\qquad \delta p_{\rm GD}\propto T^{1/2},
$$
and identifies a minimal temperature scale \(T^*\) below which the zero-mode finite-temperature description breaks down and the zero-temperature quantum-critical regime takes over [1909.10252]. A complementary finite-temperature treatment gives the thermal critical line
$$
T_c=\left(\frac{\omega_0}{2\omega}\right)\frac{1}{\tanh^{-1}\!\left(\frac{\omega\omega_0}{(g_1+g_2)^2}\right)},
$$
valid in the superradiant regime, with mutual information serving as a numerical indicator of the thermal transition [2304.07857].

A major generalization is the anisotropic Dicke–Stark model with a nonlinear Stark term and an explicit \(A^2\)-type term,
$$
H = H_{AD}+\frac{U}{N}a^\dagger a J_z + D(a^\dagger+a)^2,
$$
where \(H_{AD}\) contains unequal rotating-wave and counter-rotating-wave couplings [2405.19776]. In that setting, the standard isotropic Dicke model remains subject to the no-go theorem induced by the \(A^2\) term, but sufficiently strong anisotropy and/or Stark coupling can restore equilibrium superradiant phase transitions at both zero and finite temperature. The zero-temperature critical coupling is
$$
g_{c0}^{A^2}=\sqrt{\frac{\Delta\omega-\frac{\Delta U}{2}}{(\tau+1)^2-4\kappa}},
$$
which requires \((\tau+1)^2-4\kappa>0\) and \(U/\omega<2\); in particular, a finite critical coupling exists for \(\tau>2\sqrt{\kappa}-1>1\) under the TRK constraint [2405.19776]. The model is reported to remain in the same universality class as the conventional Dicke model, with
$$
\epsilon \propto |g-g_{c0}^{A^2}|^{1/2},\qquad
N_{\rm ph}\propto |g-g_{c0}^{A^2}|^{-1},\qquad
\Delta x\propto |g-g_{c0}^{A^2}|^{-1/4},
$$
and finite-size scaling controlled by \(\nu=2/3\) [2405.19776].

## 4. Spectral statistics, ergodicity, and dynamical diagnostics

A central question in the anisotropic Dicke model is whether non-ergodic behavior is tied to the normal-to-superradiant transition. Level statistics, the average ratio of consecutive level spacings, and OTOC-like probes show that the ergodic-to-non-ergodic transition is instead controlled by proximity to the integrable axes \(g_1=0\) and \(g_2=0\) [1610.06751]. The relevant spacing benchmarks are Poisson,
$$
P(s)=e^{-s},\qquad \langle r\rangle=2\ln2-1\approx 0.386,
$$
for non-ergodic behavior, and Wigner–Dyson/GOE,
$$
P(s)=\frac{\pi}{2}s\,e^{-\pi s^2/4},\qquad \langle r\rangle\approx 0.5307,
$$
for ergodic behavior [1610.06751]. That work explicitly argues that there is no intrinsic relation between the ergodic–non-ergodic transition and the precursors of the normal–superradiant quantum phase transition [1610.06751].

A later spectral study refines this picture by showing that the excited-state quantum phase transition and the ergodic-to-non-ergodic transition are closely related. In the central spectral band between \(E_{\text{lower}}\) and \(E_{\text{upper}}\), the level-spacing ratio and the participation ratio reveal a crossover from non-ergodic extended or multifractal behavior to ergodic delocalized behavior as the couplings increase [2304.07857]. Read together, these results suggest a distinction between two statements: the ergodic transition is not intrinsically linked to the ground-state normal-to-superradiant boundary, but it is closely organized by excited-state structure in the superradiant regime.

The out-of-time-ordered correlator provides a dynamical diagnostic of both chaos and ground-state criticality. For the anisotropic Dicke model at zero temperature,
$$
F(t)=\frac{1}{2}\left(\langle V^\dagger(0)W^\dagger(t)V(0)W(t)\rangle_\beta+\text{h.c.}\right),
$$
and the residue OTOC
$$
R(t)=1-\frac{F(t)}{F(0)}
$$
distinguishes the normal and superradiant phases: its long-time value is smaller in the normal phase and larger in the superradiant phase, and scanning it over the \((g_1,g_2)\) plane reconstructs the phase boundary \(g_1+g_2=1\) for \(\omega=\omega_0=1\) [2008.10253]. Finite-size analysis shows that the OTOC signature sharpens with increasing \(N\), whereas increasing temperature progressively blurs the boundary [2008.10253].

The model also supports genuinely dynamical reclassifications of the superradiant regime. One thermodynamic-limit effective description argues that the conventional superradiant phase is split by hidden exceptional points into three hierarchic subphases, \(\mathrm{SP}_1\), \(\mathrm{SP}_2\), and \(\mathrm{SP}_3\), characterized by harmonic-oscillator, anti-harmonic-oscillator, and inverted-harmonic-oscillator sectors in the effective Hamiltonian [2505.21169]. In that picture, the Loschmidt echo of the trivial initial state distinguishes fully oscillatory, mixed oscillatory-decaying, and fully decaying quench dynamics, and finite-\(N\) exact diagonalization for \(N=100\) is reported to confirm the predicted structure [2505.21169].

## 5. Collective modes, quantum geometry, and critical metrology

The anisotropic Dicke model has become a natural setting for quantum-geometric analyses because the normal-to-superradiant transition closes an excitation gap and makes the ground state highly parameter-sensitive. The quantum geometric tensor is
$$
\mathcal Q_{\mu\nu}=\langle \partial_\mu \psi|\partial_\nu \psi\rangle-\langle \partial_\mu \psi|\psi\rangle\langle \psi|\partial_\nu \psi\rangle,
$$
with quantum metric \(\mathcal G_{\mu\nu}=\mathrm{Re}\,\mathcal Q_{\mu\nu}\) and Berry curvature \(\mathcal F_{\mu\nu}=2\,\mathrm{Im}\,\mathcal Q_{\mu\nu}\) [2406.06301]. For a nondegenerate eigenstate,
$$
\mathcal Q_{\mu\nu}^{(n)}=\sum_{m\neq n}\frac{\langle n|\partial_\mu H|m\rangle\langle m|\partial_\nu H|n\rangle}{(E_m-E_n)^2},
$$
so singular response near criticality follows directly from gap closing [2406.06301]. The same work defines
$$
g=\frac{\lambda_1+\lambda_2}{\sqrt{\omega\Omega}},\qquad \gamma=\frac{\lambda_1}{\lambda_2},
$$
with \(g_c=1\), and distinguishes two classical limits. In the classical spin limit \(j\to\infty\), the anisotropy persists strongly and the rotating-wave coupling is more favorable for driving the transition; in the classical oscillator limit \(\eta=\Omega/\omega\to\infty\), the contributions of the rotating and counter-rotating terms become symmetric in the projected low-energy theory [2406.06301]. The finite-scale behavior resembles the classical-spin case, and the interplay among anisotropy ratio, spin length, and frequency ratio is reported to enhance the critical response without trade-off [2406.06301].

The low-lying collective spectrum becomes richer when collective matter interactions are included. In the anisotropic Dicke model with
$$
\frac{1}{N}\Big(\eta_xJ_x^2+\eta_yJ_y^2+\eta_zJ_z^2\Big),
$$
the coherent-state energy surface and the Holstein–Primakoff expansion reveal phase and amplitude polariton modes, denoted \(\epsilon_-\) and \(\epsilon_+\), together with shifted critical couplings \(\gamma_{\xi x}^c\) and \(\gamma_{\xi y}^c\) [2406.09446]. Without matter interactions, the Tavis–Cummings limit \(\xi=0\) has a Goldstone phase mode in the superradiant phase, whereas the Dicke limit \(\xi=1\) has a gapped, roton-like phase mode [2406.09446]. With matter interactions, the normal phase can be deformed, the critical couplings shift differently in the \(x\)- and \(y\)-interaction sectors, and the phase mode can be suppressed in a finite coupling window if the two critical lines are sufficiently separated [2406.09446]. The same study derives geometric phases for photon-number and spin contours,
$$
\Gamma_n=2\pi\langle\Psi_0|\hat n|\Psi_0\rangle,\qquad
\Gamma_m=2\pi\langle\Psi_0|\hat J_z-j\mathbb I|\Psi_0\rangle,
$$
which are singular at the superradiant transition and can also detect a first-order transition induced when \(\Delta\eta_x=\omega_0\) [2406.09446].

## 6. Implementations, driven extensions, and nonequilibrium transport

Rydberg-cavity platforms have become a prominent route to realizing anisotropic Dicke physics. In a generalized Dicke–Ising model for driven Rydberg arrays coupled to an optical cavity, the normalized anisotropy parameter \(\alpha\) tunes the balance between rotating-wave and counter-rotating-wave couplings, while nearest-neighbor Rydberg blockade introduces an Ising interaction [2511.22230]. The phase diagram contains a normal phase, a superradiant phase, a checkerboard Solid-\(1/2\) phase, and a superradiant solid in which superradiance coexists with translational symmetry breaking [2511.22230]. In that system, the normal-to-superradiant transition and the Solid-\(1/2\)-to-SRS transition are second-order, whereas the Solid-\(1/2\)-to-SR and SRS-to-SR transitions are first-order for any normalized anisotropy [2511.22230]. Mean-field phase boundaries are independent of \(\alpha\), but finite-size parity behavior, photon-number scaling, and the slight favoring of the SRS phase over the Solid-\(1/2\) state by stronger counter-rotating-wave fluctuations are explicitly anisotropy-dependent [2511.22230].

A complementary proposal uses periodic modulation of a pumping laser in a cavity-coupled Rydberg array to Floquet-engineer an anisotropic Dicke model with dipole-dipole interactions,
$$
H_{\mathrm{ADM}}=\widetilde\omega_C a^\dagger a+\widetilde\omega_A\sum_j n_j+V\sum_{\langle i,j\rangle}n_i n_j+\frac{\Omega}{\sqrt N}\sum_j\Big[(1-\alpha)(\sigma_j^+ a+a^\dagger\sigma_j^-)+\alpha(\sigma_j^+ a^\dagger+a\sigma_j^-)\Big],
$$
where the effective couplings are renormalized by Bessel functions [2503.13949]. The key result is tunability of the counter-rotating to rotating ratio from zero to infinity, enabling access to the pure rotating-wave limit, intermediate anisotropic regimes, and the pure counter-rotating-wave limit [2503.13949]. In the adiabatic state-preparation example, increasing \(\alpha\) converts parity-protected same-sector crossings into avoided crossings, opening a finite gap and improving the preparation of superradiant and superradiant-solid states [2503.13949].

The anisotropic Dicke model also shows sharply different responses to periodic and quasiperiodic driving. Under periodic square-wave modulation of the couplings, the dynamics reaches a plateau not followed by heating on the studied timescales, whereas under Thue–Morse and Fibonacci quasiperiodic drives a prethermal plateau is followed by heating to the infinite-temperature state [2306.16809]. The heating time scales as a stretched exponential in the Thue–Morse case,
$$
\log\tau^{*}=1.55\sqrt{\omega_d}-0.695,
$$
and exponentially in the Fibonacci case,
$$
\log\tau^{*}=0.125\omega_d-0.39
$$
[2306.16809]. The periodic drive also shifts the quantum critical boundary and extends the normal phase to larger couplings [2306.16809].

Nonequilibrium heat transport adds another layer of anisotropic physics. In a reduced anisotropic Dicke system coupled to separate bosonic thermal reservoirs for the qubit and photon sectors, a dressed master equation is used to treat strong qubit–photon coupling in the eigenbasis of the interacting Hamiltonian [2603.29180]. The central transport result is nonmonotonic in anisotropy: at weak coupling the current depends only weakly on the anisotropy factor, at moderate coupling increasing anisotropy enhances the heat current, and at strong coupling large anisotropy strongly suppresses it through multiphoton dressing [2603.29180]. Increasing the number of qubits amplifies both the peaks and the valleys of the heat-flow characteristics, and the limiting cases \(\gamma=0\) and \(\gamma=1\) admit analytical thermodynamic-limit heat currents with a cotunneling transport interpretation that serve as upper boundaries for finite-size currents [2603.29180]. The same study reports that increasing temperature bias and anisotropy enhances thermal rectification, with favorable regimes yielding rectification around \(0.45\)–\(0.5\) [2603.29180].

A distinct but related branch of the subject is the dissipative anisotropic two-photon Dicke model, in which anisotropy and cavity loss generate pole-flip transitions, localized fixed points associated with spectral collapse, Hopf bifurcations, period-doubling cascades, chaotic dynamics, and strong phase coexistence [2207.06019]. Although this is not the one-photon anisotropic Dicke model, it shows how anisotropic weighting of rotating and counter-rotating channels systematically enlarges the dynamical phase portrait of Dicke-type light–matter systems [2207.06019].

Source: https://www.emergentmind.com/topics/anisotropic-dicke-model