Skip to content

Collective Modes

A collective mode is a coherent pattern of motion of many microscopic degrees of freedom that behaves as one dynamical coordinate or response eigenchannel.

The defining feature is organization, not merely particle number. Instead of tracking every microscopic fluctuation δxj(t)\delta x_j(t) separately, one identifies a pattern

Qλq(t)=∑jeλj(q)×δxj(t),\begin{aligned} Q_{\lambda\mathbf q}(t) ={}& \sum_j e_{\lambda j}(\mathbf q) \\ &\times \delta x_j(t), \end{aligned}

whose amplitude evolves approximately independently over a stated range of wavelength, frequency, temperature, and perturbation strength. The label λ\lambda identifies a branch or polarization, and q\mathbf q identifies momentum when translation symmetry is available.

Examples include:

  • compression waves in a fluid;
  • lattice vibrations in a solid;
  • coherent spin precession in an ordered magnet;
  • plasma oscillations of charge density;
  • phase and amplitude oscillations of an order parameter;
  • in-phase and out-of-phase motion in multicomponent systems;
  • diffusive relaxation of a conserved density.

Some collective modes are sharp propagating excitations whose quanta are quasiparticles. Others are overdamped or purely relaxational. A collective mode therefore need not be gapless, ballistic, long lived, or even oscillatory.

Three complementary statements organize the subject:

one coherentcollective coordinate,one homogeneous solutionof linearized dynamics,one pole or eigenchannelof a response matrix.\begin{gathered} \text{one coherent} \\ \text{collective coordinate}, \\ \text{one homogeneous solution} \\ \text{of linearized dynamics}, \\ \text{one pole or eigenchannel} \\ \text{of a response matrix}. \end{gathered}

They become equivalent only after the variables, approximation, boundary conditions, and response convention have been specified.

This page owns the general concept and diagnostic framework of collective modes:

  • collective coordinates and polarization vectors;
  • conservative, damped, and relaxational linearized dynamics;
  • generalized normal-mode and response-matrix eigenproblems;
  • propagating, diffusive, gapped, acoustic, phase, amplitude, and relative modes;
  • mode visibility, oscillator strength, participation, and source dependence;
  • hybridization, avoided crossings, continua, and damping;
  • finite-size and thermodynamic-limit identification;
  • the distinction among collective mode, quasiparticle, normal mode, coherent state, and exact eigenstate;
  • cross-system examples connecting sound, spin waves, plasmons, and order-parameter oscillations.

Neighboring pages retain narrower ownership:

Polarons Preview specializes the framework to a mobile particle dressed by environmental modes and correlations. Lifetime and Spectral Weight develops mode-width, residue, propagation, and validity tests. Luttinger Liquid Preview develops the one-dimensional case in which collective fields replace a finite-residue particle pole. This page supplies their shared mode language without duplicating those derivations.

A microscopic description may involve:

{xj,pj,nj,Sj,ψj,…},\left\{ x_j, p_j, n_j, \mathbf S_j, \psi_j, \ldots \right\},

for a large number of sites, particles, orbitals, or field components. Exact dynamics couples these quantities through the many-body Hamiltonian.

Even when the underlying equations are local, an arbitrary microscopic disturbance generally excites many frequencies and dephases into a complicated pattern.

A collective description selects combinations

Xa(q,t)=∑jfaj(q)Oj(t)X_a(\mathbf q,t) = \sum_j f_{aj}(\mathbf q) O_j(t)

that capture the slow, coherent, or strongly responding sector. Typical XaX_a include:

  • number density;
  • momentum density;
  • energy density;
  • magnetization;
  • lattice displacement;
  • condensate phase;
  • order-parameter amplitude;
  • relative density or phase between components.

The index aa labels coupled variables before their dynamics is diagonalized. A mode is an eigencombination

Qλ=∑aeλaXa.Q_\lambda = \sum_a e_{\lambda a}X_a.

The components eλae_{\lambda a} form the mode polarization or composition.

Collectivity is scale dependent. A sound wave is a coherent density and velocity field at wavelengths long compared with microscopic collision or lattice scales. At atomic resolution, the same disturbance is a complicated superposition of particle motion.

A trustworthy claim states:

collective variablesand background statewavelength and frequency windowdynamical approximationand observable channel.\begin{gathered} \text{collective variables} \\ \text{and background state} \\ \text{wavelength and frequency window} \\ \text{dynamical approximation} \\ \text{and observable channel}. \end{gathered}

There is no basis-independent label attached to an excitation at every resolution.

A normal mode is an eigenpattern of a linearized dynamical problem. A collective mode emphasizes that the pattern organizes many microscopic degrees of freedom or macroscopic fields.

Every harmonic phonon is a normal mode and a collective mode. A two-level atom has a transition frequency, but calling that isolated internal transition a collective mode would usually add no useful information.

The terms overlap strongly, but they answer different questions:

TermMain question
normal modedoes the linearized problem diagonalize into this pattern?
collective modewhich many-body coordinate moves coherently?
quasiparticlecan the excitation be treated as a particle-like object?
resonancedoes the response have a finite-lifetime pole or enhancement?
exact eigenstateis the state stationary under the exact Hamiltonian?

A phonon, magnon, or underdamped plasmon can be both:

collective origin+particle-like quantum.\begin{gathered} \text{collective origin} \\ + \\ \text{particle-like quantum}. \end{gathered}

A diffusive density mode is collective but generally not a ballistic quasiparticle. A Landau fermionic quasiparticle is particle-like but is not primarily described as a coherent oscillation of a macroscopic coordinate.

The categories are overlapping, not mutually exclusive.

A mode is a dynamical degree of freedom or branch. A coherent state is one possible quantum state of an approximately harmonic bosonic mode.

One quantum of a phonon mode is already collective because its displacement pattern extends over many sites. A classical elastic wave corresponds instead to a state with a large coherent occupation of that mode. Large occupation is not required for collectivity.

In a finite closed system, exact eigenstates have real discrete energies. A linearized collective mode may be:

  • an exact eigenexcitation of a quadratic Hamiltonian;
  • an approximate eigenstate concentrated in a narrow spectral window;
  • a resonance that decays into other exact states;
  • a hydrodynamic pole defined only after a thermodynamic and long-wavelength limit.

Calling a mode “exact” requires a stronger statement than identifying a branch in an approximation.

Collective mode versus many-body entanglement

Section titled “Collective mode versus many-body entanglement”

Collective motion does not by itself specify an entanglement pattern. A classical normal mode of coupled oscillators is collective. A quantum one-mode excitation can be entangled in a site partition, but the amount depends on the state and partition.

Conversely, a highly entangled state need not possess a sharp collective mode. Collectivity is a dynamical organization, not an entanglement measure.

Choose a stationary or slowly varying background Xa(0)X_a^{(0)} and write

Xa=Xa(0)+δXa.X_a = X_a^{(0)} + \delta X_a.

The background may be:

  • an equilibrium fluid;
  • a crystal;
  • an ordered magnet;
  • a condensate;
  • a mean-field saddle;
  • a thermal state characterized by conserved densities.

Linearization keeps terms first order in δXa\delta X_a. This step is controlled only when the perturbation remains small in the variables used.

For real coordinates with inertia, a common equation is

∑b[Mab∂t2+Kab(q)]δXb=0.\sum_b \left[ M_{ab}\partial_t^2 + K_{ab}(\mathbf q) \right] \delta X_b = 0.

With

δXb(t)=ebe−iωt,\delta X_b(t) = e_b e^{-i\omega t},

one obtains the generalized eigenproblem

∑bKab(q)eb=ω2∑bMabeb.\sum_b K_{ab}(\mathbf q)e_b = \omega^2 \sum_b M_{ab}e_b.

If MM is positive and KK is positive semidefinite, the stable frequencies are real and nonnegative.

The mode condition is

det⁡[K(q)−ω2M]=0.\det \left[ K(\mathbf q) -\omega^2M \right] = 0.

Not every mode has Newton-like second-order dynamics. Canonically conjugate fields, spin precession, kinetic equations, and quantum amplitudes often obey

i∂tδX=L(q)δX.i\partial_t \delta\mathbf X = \mathcal L(\mathbf q) \delta\mathbf X.

Then

L(q)eλ=ωλ(q)eλ.\mathcal L(\mathbf q)\mathbf e_\lambda = \omega_\lambda(\mathbf q) \mathbf e_\lambda.

The operator L\mathcal L need not be Hermitian in the ordinary Euclidean metric. Bosonic Bogoliubov problems use an indefinite particle–hole metric, while dissipative problems are genuinely non-Hermitian.

A useful schematic interpolation is

∑b[−ω2Mab−iωΓab+Kab]eb=0.\sum_b \left[ -\omega^2M_{ab} -i\omega\Gamma_{ab} +K_{ab} \right] e_b = 0.

Here Γ\Gamma represents damping only within the stated effective description. It may arise from collisions, decay into omitted modes, disorder, viscosity, or a bath.

With the Fourier convention

ei(q⋅r−ωt),e^{i(\mathbf q\cdot\mathbf r-\omega t)},

stable damped poles satisfy

Im⁡ω<0.\operatorname{Im}\omega < 0.

A linearized eigenfrequency can reveal that the chosen background is unstable:

  • ω2<0\omega^2<0 in a conservative second-order problem;
  • Im⁡ω>0\operatorname{Im}\omega>0 in the retarded convention;
  • a negative stiffness or compressibility;
  • a zero mode associated with a flat direction.

An exponentially growing solution is a collective instability, not a stable excitation branch.

Two equal masses already show the essential change of variables. Consider

L=m2(x˙12+x˙22)−k02(x12+x22)−kc2(x1−x2)2.\begin{aligned} L ={}& \frac{m}{2} \left( \dot x_1^2+\dot x_2^2 \right) \\ &- \frac{k_0}{2} \left( x_1^2+x_2^2 \right) \\ &- \frac{k_c}{2} \left( x_1-x_2 \right)^2. \end{aligned}

Define

Q+=x1+x22,Q−=x1−x22.Q_+ = \frac{x_1+x_2}{\sqrt2}, \qquad Q_- = \frac{x_1-x_2}{\sqrt2}.

The Lagrangian becomes

L=m2(Q˙+2+Q˙−2)−k02Q+2−k0+2kc2Q−2.\begin{aligned} L ={}& \frac{m}{2} \left( \dot Q_+^2+\dot Q_-^2 \right) \\ &- \frac{k_0}{2}Q_+^2 \\ &- \frac{k_0+2k_c}{2}Q_-^2. \end{aligned}

Therefore

ω+=k0m,ω−=k0+2kcm.\omega_+ = \sqrt{\frac{k_0}{m}}, \qquad \omega_- = \sqrt{ \frac{k_0+2k_c}{m} }.

The ++ mode moves in phase and does not stretch the coupling spring. The −- mode moves out of phase and pays the additional stiffness 2kc2k_c.

This example teaches four durable lessons:

  1. a mode is a pattern, not one coordinate;
  2. eigenvectors and frequencies are separate pieces of information;
  3. the source can couple differently to the two patterns;
  4. in-phase and relative modes recur in multiband, multisublattice, and multicomponent quantum systems.

After mass normalization, a stable mode may have

Hλ=12(Pλ2+ωλ2Qλ2).H_\lambda = \frac12 \left( P_\lambda^2 + \omega_\lambda^2Q_\lambda^2 \right).

Introduce

Qλ=ℏ2ωλ(bλ+bλ†),Pλ=−iℏωλ2(bλ−bλ†).\begin{aligned} Q_\lambda ={}& \sqrt{ \frac{\hbar}{2\omega_\lambda} } \left( b_\lambda+b_\lambda^\dagger \right), \\ P_\lambda ={}& -i \sqrt{ \frac{\hbar\omega_\lambda}{2} } \left( b_\lambda-b_\lambda^\dagger \right). \end{aligned}

Then

Hλ=ℏωλ(bλ†bλ+12).H_\lambda = \hbar\omega_\lambda \left( b_\lambda^\dagger b_\lambda + \frac12 \right).

The operator bλ†b_\lambda^\dagger creates one quantum of the entire pattern.

The harmonic algebra is exact only for an exactly quadratic canonical mode. In an interacting many-body system:

  • mode operators can contain nonlinear corrections;
  • different modes interact;
  • mode number may not be conserved;
  • high occupation can invalidate the linearization;
  • a branch can decay into other branches or continua.

Bosonic quasiparticle operators are therefore often leading effective variables rather than exact microscopic creation operators.

For a normalized extended mode, a single quantum often gives each microscopic coordinate an amplitude of order N−1/2N^{-1/2}. The total mode energy remains of order one because NN small contributions add coherently.

This is the useful heuristic:

In a collective quantum, many degrees of freedom move a little in a correlated pattern.

It is a heuristic, not a universal definition. Localized collective modes, boundary modes, and strongly inhomogeneous systems require more careful language.

Let sources fbf_b couple to operators XbX_b:

δH(t)=−∑bfb(t)Xb.\delta H(t) = -\sum_b f_b(t)X_b.

Linear response gives

δ⟨Xa⟩=∑bχabRfb.\delta\langle X_a\rangle = \sum_b \chi_{ab}^{\mathrm R}f_b.

For translation-invariant systems, χR\chi^{\mathrm R} depends on (q,ω)(\mathbf q,\omega).

The source vector chooses how the system is driven. The detector vector chooses which response is read out. A mode can exist while remaining dark in a particular source–detector channel.

For Gaussian fluctuations or a linearized effective theory, write

∑bKabRδXb=fa.\sum_b \mathcal K_{ab}^{\mathrm R} \delta X_b = f_a.

When no contact or vertex convention changes the relation,

χR=(KR)−1.\chi^{\mathrm R} = \left( \mathcal K^{\mathrm R} \right)^{-1}.

Collective modes solve the homogeneous equation

KReλ=0,\mathcal K^{\mathrm R} \mathbf e_\lambda = 0,

so their complex frequencies satisfy

det⁡KR(q,ω)=0.\det \mathcal K^{\mathrm R}(\mathbf q,\omega) = 0.

This determinant condition is necessary but not sufficient for visibility. The numerator, source projection, detector projection, and possible cancellations must also be checked.

Near a simple positive-frequency pole,

χabR(q,ω)≃Rab(λ)(q)ω−Ωλ(q)+iκλ(q)+negative-frequency pole+χab,regR.\begin{aligned} \chi_{ab}^{\mathrm R} (\mathbf q,\omega) \simeq{}& \frac{ R_{ab}^{(\lambda)}(\mathbf q) }{ \omega-\Omega_\lambda(\mathbf q) +i\kappa_\lambda(\mathbf q) } \\ &+ \text{negative-frequency pole} \\ &+ \chi_{ab,\mathrm{reg}}^{\mathrm R}. \end{aligned}

For a well-isolated mode, the residue matrix often factorizes:

Rab(λ)∝uλavλb∗,R_{ab}^{(\lambda)} \propto u_{\lambda a} v_{\lambda b}^*,

where right and left mode vectors can differ in a damped non-Hermitian problem.

The residue matrix encodes polarization and visibility. The pole position encodes frequency and damping. Neither alone is a complete mode identification.

For a Hermitian coordinate in a conservative system, positive and negative frequencies occur together. A schematic oscillator response is

χR(ω)=1Ω02−(ω+i0+)2.\chi^{\mathrm R}(\omega) = \frac{1}{ \Omega_0^2 - (\omega+i0^+)^2 }.

Its two poles reflect one real oscillatory coordinate, not two independent positive-energy species.

Every exact transition in a finite isolated system contributes a real-axis delta line. That fact does not make every transition a collective mode.

A finite-size collective identification should show some combination of:

  • a coherent branch across momentum or another control parameter;
  • concentration of weight in a physically motivated collective operator;
  • stable polarization or symmetry quantum numbers;
  • scaling toward a thermodynamic pole or continuum feature;
  • agreement with a linearized effective description.

If several variables mix, the eigenvector

eλ=(eλ1,eλ2,…)T\mathbf e_\lambda = \left( e_{\lambda1}, e_{\lambda2}, \ldots \right)^{\mathsf T}

states which variables oscillate and with what relative amplitude and phase.

Examples include:

  • longitudinal versus transverse displacement;
  • density versus spin density;
  • in-phase versus out-of-phase component motion;
  • order-parameter amplitude versus phase;
  • charge-like versus neutral combinations.

The polarization can change continuously along a branch and exchange character near an avoided crossing.

For a detector

Od=∑adaXa,O_{\mathbf d} = \sum_a d_aX_a,

the observed pole weight depends on

d†R(λ)f,\mathbf d^\dagger R^{(\lambda)} \mathbf f,

where f\mathbf f is the source vector. A mode with nonzero residue matrix can be invisible when symmetry or destructive interference makes this projection vanish.

“No peak in this channel” is not equivalent to “no mode exists.”

For a normalized discrete pattern

∑j∣eλj∣2=1,\sum_j \lvert e_{\lambda j}\rvert^2 = 1,

one possible participation measure is

Pλ=1∑j∣eλj∣4.P_\lambda = \frac{1}{ \sum_j \lvert e_{\lambda j}\rvert^4 }.

It gives Pλ=1P_\lambda=1 for weight on one coordinate and Pλ=NP_\lambda=N for equal magnitude on NN coordinates.

This number is useful but not invariant:

  • it depends on the basis and normalization;
  • internal components and nonorthogonal coordinates require a metric;
  • a long-wavelength mode can be spatially extended but dark;
  • an edge or defect mode can be collective within a localized region.

Participation is evidence about spatial organization, not a universal definition of collectivity.

For

Oq=∑j=1Ne−iq⋅rjOj,O_{\mathbf q} = \sum_{j=1}^{N} e^{-i\mathbf q\cdot\mathbf r_j} O_j,

spectral weight can scale with NN. If instead one uses Oq/NO_{\mathbf q}/\sqrt N, the same coherent branch can carry order-one weight.

Claims about “macroscopic oscillator strength” must therefore state the operator normalization.

Interactions or mechanical couplings make motion of one degree of freedom alter the forces on others. The resulting feedback can synchronize a pattern and supply a restoring force.

Examples:

  • elastic bonds restore lattice displacement;
  • compressibility restores density compression;
  • exchange restores spin twists;
  • Coulomb fields restore charge separation;
  • condensation energy restores order-parameter amplitude.

Conserved densities cannot relax locally without transport. If

∂tn+∇⋅j=0,\partial_t n + \nabla\cdot\mathbf j = 0,

then long-wavelength density fluctuations are slow because gradients are small. Conservation therefore produces hydrodynamic modes even when no broken symmetry is present.

When an exact continuous global symmetry is spontaneously broken, slow variations along the degenerate order-parameter manifold can produce Goldstone modes.

The full assumptions and counting rule belong to Goldstone Modes in Many-Body Systems. Not every collective mode is a Goldstone mode:

  • plasmons can be gapped;
  • amplitude modes are generally gapped;
  • diffusion follows conservation rather than broken symmetry;
  • optical phonons need not soften at zero momentum.

A fluctuation can generate a field that acts back on the complete system. Long-range Coulomb interaction is the canonical example:

δn⟶δϕ⟶collective restoring force.\delta n \longrightarrow \delta\phi \longrightarrow \text{collective restoring force}.

The feedback can move a neutral acoustic branch to a nonzero plasma frequency.

Boundaries, traps, finite samples, and interfaces create global shape, sloshing, breathing, surface, and edge modes. Their frequencies depend on geometry and boundary conditions as well as bulk constitutive data.

A center-of-mass oscillation can be highly collective while remaining protected from internal interactions by a separation theorem. “Collective” does not mean “strongly renormalized.”

Collective variableTypical modeLong-wavelength character
mass or number densitysoundacoustic and propagating
conserved density without inertiadiffusionrelaxational
charge densityplasmonoften gapped in three dimensions
lattice displacementphononacoustic or optical
spin orientationmagnon or spin wavelinear or quadratic, depending on dynamics
condensate or order-parameter phasephase modeoften Goldstone-like when neutral
order-parameter magnitudeamplitude modegenerally gapped and decay sensitive
relative phase of componentsLeggett-type modeusually gapped
relative density of componentsspin-density or counterflow modeacoustic, gapped, or diffusive
shape of a confined cloudbreathing or quadrupole modegeometry dependent

This table classifies by the primary collective coordinate. Real modes can mix several entries.

Combine the continuity equation

∂tn+∇⋅j=0\partial_t n + \nabla\cdot\mathbf j = 0

with Fick’s constitutive law

j=−D∇n.\mathbf j = -D\nabla n.

Linearizing around a uniform state gives

∂tδn=D∇2δn.\partial_t\delta n = D\nabla^2\delta n.

For

δn∝ei(q⋅r−ωt),\delta n \propto e^{i(\mathbf q\cdot\mathbf r-\omega t)},

the mode is

ω(q)=−iDq2.\omega(\mathbf q) = -iDq^2.

It is collective because it governs coherent long-wavelength relaxation of a conserved density. It is not an oscillatory wave and generally does not define a ballistic quasiparticle.

As q→0q\to0,

∣ω∣∼Dq2→0.\lvert\omega\rvert \sim Dq^2 \to 0.

The slow rate is enforced by conservation: a nearly uniform excess can relax only by transporting density over a distance of order q−1q^{-1}.

The diffusion constant is not fixed by conservation alone. It depends on transport coefficients and thermodynamic susceptibilities.

Charge, energy, momentum, and component numbers can mix. Linearized hydrodynamics then produces a matrix

−iω δn=−q2D δn.-i\omega\, \delta\mathbf n = -q^2 \mathsf D\, \delta\mathbf n.

Its eigenvalues determine diffusive combinations. Off-diagonal thermoelectric or spin–charge couplings rotate the hydrodynamic eigenvectors away from the original microscopic labels.

For a neutral fluid of particles with mass mm, linearized number conservation gives

∂tδn+n0∇⋅v=0.\partial_t\delta n + n_0\nabla\cdot\mathbf v = 0.

Neglecting viscosity, the Euler equation is

mn0∂tv=−∇δP.mn_0\partial_t\mathbf v = -\nabla\delta P.

For an isentropic disturbance,

δP=(∂P∂n) ⁣sδn.\delta P = \left( \frac{\partial P}{\partial n} \right)_{\!s} \delta n.

Eliminating v\mathbf v gives

∂t2δn=cs2∇2δn,\partial_t^2\delta n = c_s^2\nabla^2\delta n,

with

cs2=1m(∂P∂n) ⁣s.c_s^2 = \frac{1}{m} \left( \frac{\partial P}{\partial n} \right)_{\!s}.

The two propagating roots are

ω±(q)=±csq.\omega_\pm(\mathbf q) = \pm c_sq.

The density and longitudinal velocity are not separate modes here. They are components of one sound eigenvector.

Viscosity and thermal transport broaden the poles:

ω±(q)=±csq−iΓsq2+⋯ .\omega_\pm(\mathbf q) = \pm c_sq -i\Gamma_sq^2 + \cdots.

At sufficiently small qq,

∣Im⁡ω∣∣Re⁡ω∣∼Γscsq→0,\frac{ \lvert\operatorname{Im}\omega\rvert }{ \lvert\operatorname{Re}\omega\rvert } \sim \frac{\Gamma_s}{c_s}q \to 0,

so hydrodynamic sound can become asymptotically sharp even in an interacting finite-temperature fluid.

The same density channel can support different collective regimes.

For collision time τcoll\tau_{\mathrm{coll}}:

regimeωτcollhydrodynamic≪1collisionless≫1\begin{array}{c|c} \text{regime} & \omega\tau_{\mathrm{coll}} \\ \hline \text{hydrodynamic} & \ll1 \\ \text{collisionless} & \gg1 \end{array}

First sound is controlled by local thermodynamics and transport. Zero sound is a collisionless collective oscillation sustained by self-consistent quasiparticle interactions.

Fermi Liquid Theory Preview owns the Landau kinetic equation and zero-sound condition. The important lesson here is that a mode name is incomplete without its dynamical regime.

Multicomponent and superfluid systems can support:

  • in-phase density sound;
  • counterflow or spin-density sound;
  • first and second sound;
  • entrainment-mixed branches;
  • additional diffusive thermal or concentration modes.

Counting conserved quantities alone is not enough. One must include constitutive relations, inertia, broken symmetries, and couplings among the variables.

Superfluidity in Condensed Matter owns the helium and neutral-superfluid evidence audit for first sound, second sound, and roton branches; this page retains the generic eigenmode taxonomy.

For magnetization M\mathbf M, a schematic nondissipative equation is

∂tM=−γM×Heff.\partial_t\mathbf M = -\gamma \mathbf M \times \mathbf H_{\mathrm{eff}}.

Linearization around an ordered state couples transverse components. The resulting spin wave describes coherent precession of many local moments or spin densities.

Ferromagnetic and antiferromagnetic order can break similar spin symmetries while producing different low-energy dynamics:

ωFM∝q2,ωAFM∝q.\omega_{\mathrm{FM}} \propto q^2, \qquad \omega_{\mathrm{AFM}} \propto q.

The difference reflects canonical structure and background charge densities, not merely the number of broken generators.

Goldstone Modes in Many-Body Systems owns the counting and type-A/type-B explanation. Heisenberg Model owns model-specific one-magnon formulas. Magnons owns the Holstein–Primakoff quantization, ferromagnet–antiferromagnet comparison, and quasiparticle limits.

In an itinerant system, a spin collective mode can coexist with a particle–hole continuum. Its fate depends on:

  • spin-rotation symmetry;
  • exchange enhancement;
  • momentum;
  • available decay phase space;
  • disorder and spin–orbit coupling.

A sharp transverse branch outside a continuum can become Landau damped after entering it.

For a complex order parameter,

Ψ=(Ψ0+h)eiθ,\Psi = \left( \Psi_0+h \right) e^{i\theta},

the fluctuation hh changes the magnitude and θ\theta changes the phase.

Geometrically:

  • hh is radial motion;
  • θ\theta is tangential motion along the order-parameter manifold.

A schematic quadratic Lagrangian is

L=χθ2(∂tθ)2−ρs2(∇θ)2+χh2(∂th)2−ρh2(∇h)2−Δh22h2+λh∂tθ.\begin{aligned} \mathcal L ={}& \frac{\chi_\theta}{2} \left( \partial_t\theta \right)^2 - \frac{\rho_s}{2} \left( \nabla\theta \right)^2 \\ &+ \frac{\chi_h}{2} \left( \partial_th \right)^2 - \frac{\rho_h}{2} \left( \nabla h \right)^2 \\ &- \frac{\Delta_h^2}{2}h^2 + \lambda h\partial_t\theta. \end{aligned}

When λ=0\lambda=0, phase and amplitude separate at quadratic order. The neutral phase branch has

ωθ2=ρsχθq2,\omega_\theta^2 = \frac{\rho_s}{\chi_\theta} q^2,

while the amplitude branch has

ωh2=Δh2+ρhq2χh.\omega_h^2 = \frac{ \Delta_h^2+\rho_hq^2 }{ \chi_h }.

Particle–hole asymmetry, multiple order parameters, explicit symmetry breaking, and nonequilibrium driving can mix amplitude and phase.

Then “phase mode” and “amplitude mode” describe dominant polarization, not exact pure coordinates. The correct object is the eigenvector of the complete fluctuation kernel.

A radial order-parameter fluctuation is often called a Higgs mode by analogy. The analogy is useful only when:

  • an amplitude-like eigenchannel is well defined;
  • the mode is sufficiently isolated from decay continua;
  • its symmetry and selection rules are stated;
  • the measured feature cannot be explained by a threshold or pair-breaking edge alone.

Not every broad bump near twice a gap is a sharp amplitude collective mode.

In a neutral superfluid, phase stiffness can produce an acoustic Goldstone mode. In a charged superconductor, long-range electromagnetic coupling reorganizes the longitudinal phase response and moves spectral weight toward a plasma scale.

This is not ordinary explicit breaking of the global symmetry by a small source. It is dynamical coupling to gauge and charge-density fields. BCS Mean-Field Theory and Goldstone Modes in Many-Body Systems retain the detailed ownership.

Neutral compression versus charged oscillation

Section titled “Neutral compression versus charged oscillation”

For a short-range neutral fluid, long-wavelength density motion is usually acoustic:

ω∼cq.\omega \sim cq.

For a three-dimensional electron fluid with unscreened Coulomb interaction, a uniform displacement of electrons relative to the positive background creates a macroscopic electric restoring field.

The long-wavelength scale is

ωp2=n0e2mϵ0\omega_p^2 = \frac{ n_0e^2 }{ m\epsilon_0 }

in the simplest continuum SI convention.

Thus

lim⁡q→0ωpl(q)=ωp≠0.\lim_{q\to0} \omega_{\mathrm{pl}}(q) = \omega_p \ne 0.

The precise dispersion changes with dimension, screening, lattice band mass, dielectric environment, and geometry. Plasmons Preview owns those charge-mode cases.

The Bohm–Pines construction made the distinction vivid: long-range Coulomb motion can be reorganized into collective plasma coordinates, while residual short-range screened interactions act among particle-like degrees of freedom.

Modern response language expresses the same idea through zeros of a dielectric function or poles of a screened density response.

A charge mode can decay into particle–hole excitations when its frequency and momentum overlap their continuum and the matrix element is nonzero.

The mode then acquires a complex pole:

ωpl⟶Ωpl−iκpl.\omega_{\mathrm{pl}} \longrightarrow \Omega_{\mathrm{pl}} -i\kappa_{\mathrm{pl}}.

Entering a continuum does not always erase the mode immediately, but it removes any automatic guarantee of a narrow Lorentzian.

Suppose two bare coordinates have inverse response

K(ω)=(ω−ω1−g−gω−ω2).\mathcal K(\omega) = \begin{pmatrix} \omega-\omega_1 & -g \\ -g & \omega-\omega_2 \end{pmatrix}.

The collective frequencies solve

det⁡K=0,\det\mathcal K = 0,

giving

ω±=ω1+ω22±(ω1−ω22)2+g2.\begin{aligned} \omega_\pm ={}& \frac{\omega_1+\omega_2}{2} \\ &\pm \sqrt{ \left( \frac{\omega_1-\omega_2}{2} \right)^2 +g^2 }. \end{aligned}

At the uncoupled crossing ω1=ω2\omega_1=\omega_2,

ω+−ω−=2∣g∣.\omega_+-\omega_- = 2\lvert g\rvert.

This is an avoided crossing.

Far to one side of the crossing, the upper branch may resemble coordinate 1 and the lower branch coordinate 2. Far to the other side, their characters exchange.

A reliable mode label should therefore track:

  • frequency;
  • eigenvector overlap;
  • symmetry;
  • source and detector weight;
  • continuity across parameters.

Sorting only by energy can mislabel branches near hybridization.

For two identical components, a uniform source may couple only to

Q+=X1+X22,Q_+ = \frac{X_1+X_2}{\sqrt2},

while

Q−=X1−X22Q_- = \frac{X_1-X_2}{\sqrt2}

is dark in that channel.

The relative mode still exists. A component-selective or antisymmetric source is required to reveal it.

If one coordinate couples to many excitations, its inverse response acquires a self-energy:

KR=K0R−ΣR(ω).\mathcal K^{\mathrm R} = \mathcal K_0^{\mathrm R} - \Sigma^{\mathrm R}(\omega).

Rapid frequency dependence or an imaginary part of ΣR\Sigma^{\mathrm R} can:

  • shift the mode;
  • broaden it;
  • distort the line shape;
  • transfer weight to a continuum;
  • split a feature into several resonances;
  • remove an isolated pole.

An avoided crossing describes two identifiable modes. A threshold cusp or branch-cut feature is a different analytic structure.

An acoustic branch satisfies

ω(q)→0\omega(\mathbf q) \to 0

as q\mathbf q approaches a specified soft momentum.

A gapped branch satisfies

ω(q)→Δ/ℏ>0.\omega(\mathbf q) \to \Delta/\hbar > 0.

“Optical” often means that sublattices or components move against one another and that the zone-center frequency is nonzero. It does not mean the mode is necessarily excited by visible light.

Common long-wavelength forms include

ω=cq+⋯\omega = cq +\cdots

and

ω=Dsq2+⋯ .\omega = D_sq^2 +\cdots.

The exponent reflects symmetry, conservation, canonical structure, dimensionality, and interactions. It should be derived rather than guessed from the word “collective.”

A mode softens when

Ωλ(Q)→0\Omega_\lambda(\mathbf Q) \to 0

as a control parameter approaches an instability or phase transition.

Softening can indicate:

  • loss of a restoring stiffness;
  • onset of spatial order at Q\mathbf Q;
  • restoration or breaking of a symmetry;
  • structural instability;
  • critical slowing down.

The mode may also become strongly damped. At a critical point, a relaxational pole or broad scaling continuum can replace a sharp oscillator.

A branch can remain gapped but develop a local minimum at finite momentum. Such a minimum signals favorable correlations or an incipient instability, but its interpretation is model dependent.

One should not call every finite-qq minimum a roton without stating the system and the evidence.

Consider

Q¨+2κQ˙+Ω02Q=F(t).\ddot Q + 2\kappa\dot Q + \Omega_0^2Q = F(t).

The retarded response is proportional to

χR(ω)=1Ω02−ω2−2iκω.\chi^{\mathrm R}(\omega) = \frac{1}{ \Omega_0^2 - \omega^2 - 2i\kappa\omega }.

The poles are

ω±=−iκ±Ω02−κ2.\omega_\pm = -i\kappa \pm \sqrt{ \Omega_0^2-\kappa^2 }.

For

κ<Ω0,\kappa < \Omega_0,

the response oscillates while decaying. A useful mode frequency is

Ωd=Ω02−κ2.\Omega_d = \sqrt{ \Omega_0^2-\kappa^2 }.

The quality factor in the weak-damping convention is

Qqual≃Ω02κ.Q_{\mathrm{qual}} \simeq \frac{\Omega_0}{2\kappa}.

For

κ>Ω0,\kappa > \Omega_0,

both poles lie on the negative imaginary axis. The response relaxes without completing an oscillation.

The collective coordinate remains meaningful, but particle-like language becomes poor.

At

κ=Ω0,\kappa = \Omega_0,

the two poles coincide in the elementary model. Real many-body systems can have additional memory kernels and frequency-dependent damping, so the simple threshold need not remain exact.

Collective modes can decay or broaden through:

  • anharmonic mode–mode scattering;
  • particle–hole creation;
  • quasiparticle collisions;
  • disorder and boundary scattering;
  • viscosity and thermal conduction;
  • radiation into electromagnetic or substrate modes;
  • explicit symmetry breaking;
  • coupling to an external environment;
  • inhomogeneous broadening.

Different mechanisms have different momentum, temperature, and frequency dependence.

Landau damping is collisionless transfer of collective motion into resonant microscopic excitations. Collisional damping comes from scattering processes that redistribute occupations.

Both broaden response, but they are physically and parametrically distinct.

If a positive-frequency pole is written

ω⋆=Ω−iκ,\omega_\star = \Omega-i\kappa,

then the amplitude decays as e−κte^{-\kappa t}. A Lorentzian in angular frequency has HWHM κ\kappa and FWHM 2κ2\kappa under the simplest isolated-pole convention.

Spectral Functions owns the full conversion among angular-frequency width, energy width, amplitude decay, and population lifetime.

Many microscopic displacements forming one collective coordinate, a sharp branch entering a continuum, and propagating versus relaxational time dependence.

Three views of a collective mode. A mode coordinate combines many microscopic fluctuations with polarization eλje_{\lambda j}. In response space, an underdamped branch Ωλ(q)\Omega_\lambda(q) can remain sharp outside a continuum and broaden after entering it. In time, collective dynamics may be oscillatory, with poles near ±Ω−iκ\pm\Omega-i\kappa, or relaxational, with a pole such as −iDq2-iDq^2.

For a collective operator OqO_{\mathbf q} at zero temperature,

SO(q,ω)∝∑n∣⟨n∣Oq∣0⟩∣2×δ(ω−ωn0).\begin{aligned} S_O(\mathbf q,\omega) \propto{}& \sum_n \left| \langle n| O_{\mathbf q} |0\rangle \right|^2 \\ &\times \delta \left( \omega-\omega_{n0} \right). \end{aligned}

A collective mode appears when substantial operator weight organizes into a coherent branch or controlled pole.

The exact eigenstates remain the foundation. “Collective mode” is an interpretation of their energies, matrix elements, scaling, and dynamical organization.

If one state exhausts most of the relevant spectral weight,

SO(q,ω)≃ZO(q)δ(ω−Ω(q))S_O(\mathbf q,\omega) \simeq Z_O(\mathbf q) \delta \left( \omega-\Omega(\mathbf q) \right)

in an ideal stable limit.

The single-mode approximation can be powerful, but its quality must be checked against:

  • omitted weight;
  • multiparticle continua;
  • sum rules;
  • finite-size level spacing;
  • operator dependence.

A density collective trial state has the schematic form

∣q⟩∝ρq∣0⟩.|\mathbf q\rangle \propto \rho_{\mathbf q}|0\rangle.

Its variational energy connects static structure to an excitation estimate. This does not prove exact single-mode saturation.

Sum Rules owns the exact moment derivation and Feynman bound.

A collective operator can couple mainly to a continuum:

SO=ScontS_O = S_{\mathrm{cont}}

with no isolated pole. The operator is collective, but the spectrum need not contain a sharp collective mode.

This distinction matters near quantum critical points, in one-dimensional fractionalized systems, and in strongly damped metals.

Collective Modes and Particle–Hole Continua

Section titled “Collective Modes and Particle–Hole Continua”

In a fermionic reference state, a density-like collective excitation can have the leading form

∣λ,q⟩∼∑kϕλq(k)×ck+q†ck∣Ψ0⟩+⋯ .\begin{aligned} |\lambda,\mathbf q\rangle \sim{}& \sum_{\mathbf k} \phi_{\lambda\mathbf q}(\mathbf k) \\ &\times c_{\mathbf k+\mathbf q}^\dagger c_{\mathbf k} |\Psi_0\rangle + \cdots. \end{aligned}

The coherent superposition across many k\mathbf k values distinguishes a collective eigenvector from one selected particle–hole configuration.

The set of individual particle–hole energies

ω=ϵk+q−ϵkℏ\omega = \frac{ \epsilon_{\mathbf k+\mathbf q} - \epsilon_{\mathbf k} }{ \hbar }

fills a continuum as k\mathbf k varies.

Interactions can pull or organize a pole outside that continuum. Inside it, resonant decay often broadens the mode.

There is no universal threshold such as “more than half the particles participate.” A mode can be collective because:

  • its eigenvector is a coherent superposition;
  • its frequency is set by a self-consistent field;
  • its residue scales coherently;
  • it is governed by a hydrodynamic variable;
  • it is protected by symmetry or conservation.

The relevant criterion depends on the problem.

In principle, collective modes are read from exact retarded correlators or structure factors. In practice, exact spectra are available only in special models or finite numerical systems.

Linearizing a time-dependent mean-field equation around a stationary saddle produces small-oscillation modes. Examples include:

  • time-dependent Hartree or Hartree–Fock;
  • Gross–Pitaevskii and Bogoliubov–de Gennes equations;
  • time-dependent density-functional approximations;
  • linearized Landau kinetic theory.

The mode solution is only as controlled as the background and closure.

RPA solves a self-consistent response problem. In a scalar density channel,

χRPA=χ01−Vχ0.\chi_{\mathrm{RPA}} = \frac{ \chi_0 }{ 1-V\chi_0 }.

A collective pole can satisfy

1−V(q)χ0(q,ω)=0.1 - V(\mathbf q) \chi_0(\mathbf q,\omega) = 0.

The numerator and imaginary part of χ0\chi_0 still determine visibility and damping. Random Phase Approximation owns signs, screening, Lindhard functions, sum rules, and validity.

For fields Φa\Phi_a, expand an effective action about a saddle:

Seff=S0+12∑abδΦaKabδΦb+⋯ .S_{\mathrm{eff}} = S_0 + \frac12 \sum_{ab} \delta\Phi_a \mathcal K_{ab} \delta\Phi_b + \cdots.

After the correct real-time continuation, zeros of det⁡KR\det\mathcal K^{\mathrm R} give collective poles.

An imaginary-time Hessian alone does not determine retarded damping without analytic continuation.

At long wavelength and low frequency, conservation laws and constitutive relations can determine pole structure without solving microscopic dynamics.

Hydrodynamics predicts the functional form of slow modes. Microscopic theory supplies thermodynamic derivatives and transport coefficients.

Commutators generate exact spectral moments:

⟨[O†,[H,O]]⟩.\langle [O^\dagger,[H,O]] \rangle.

Moment constraints can bound or estimate collective energies, but a finite number of moments cannot by itself distinguish a narrow pole from every possible continuum.

Exact diagonalization, tensor networks, quantum Monte Carlo with analytic continuation, and real-time simulation can reveal candidate modes.

Necessary controls include:

  • system-size scaling;
  • frequency broadening or time-window scaling;
  • momentum resolution;
  • operator normalization;
  • boundary conditions;
  • analytic-continuation uncertainty;
  • sum-rule closure;
  • comparison across several channels.

Neutron, x-ray, electron-energy-loss, and Bragg scattering access momentum- and frequency-resolved correlations.

The measured intensity is schematically

I(q,ω)∼f†[−Im⁡χR(q,ω)]f∗R(q,ω)+background,\begin{aligned} I(\mathbf q,\omega) \sim{}& \mathbf f^\dagger \left[ -\operatorname{Im} \chi^{\mathrm R} (\mathbf q,\omega) \right] \mathbf f \\ &\ast \mathcal R(\mathbf q,\omega) + \text{background}, \end{aligned}

with probe vertices f\mathbf f and resolution kernel R\mathcal R.

Photons carry little crystal momentum in ordinary optical geometries, so these probes emphasize near-zone-center modes and symmetry-selected channels.

A mode can be Raman active, infrared active, both, or neither depending on symmetry and effective charge.

A short perturbation can launch coherent oscillations:

δO(t)∼Ae−κtcos⁡(Ωt+φ).\delta O(t) \sim A e^{-\kappa t} \cos \left( \Omega t+\varphi \right).

The observed frequency and damping can identify a mode, but nonlinear driving, heating, inhomogeneity, and multiple beating frequencies must be controlled.

Collective modes affect:

  • thermal conductivity;
  • sound attenuation;
  • optical conductivity;
  • heat capacity;
  • drag and viscosity;
  • critical relaxation.

These integrated observables constrain a mode description but usually do not determine its complete dispersion and polarization.

A strong experimental or numerical identification combines:

  1. a reproducible branch or pole;
  2. symmetry-consistent polarization;
  3. source and detector selection rules;
  4. controlled intrinsic width;
  5. spectral-weight and sum-rule accounting;
  6. finite-size or resolution scaling;
  7. agreement with an independently measured thermodynamic or transport parameter;
  8. a method whose validity covers the claimed regime.

One peak at one momentum is not yet a collective-mode theory.

Specify the equilibrium, ordered, metastable, or driven state around which fluctuations are defined.

State what oscillates:

δn,δM,u,θ,h,…\delta n, \quad \delta\mathbf M, \quad \mathbf u, \quad \theta, \quad h, \quad \ldots

Declare:

  • hydrodynamic or collisionless;
  • conservative or dissipative;
  • isolated or open;
  • linear or nonlinear;
  • finite or thermodynamic.

Write the linearized equations or inverse response matrix. Do not infer a mode solely from terminology.

Report both

ωλ(q)\omega_\lambda(\mathbf q)

and

eλ(q).\mathbf e_\lambda(\mathbf q).

Check whether the chosen source and detector couple to the eigenvector.

Identify decay thresholds and overlapping branches.

Compare

∣Im⁡ωλ∣\lvert\operatorname{Im}\omega_\lambda\rvert

with the relevant real frequency, branch separation, and experimental resolution.

Vary:

  • size;
  • momentum;
  • broadening;
  • boundary conditions;
  • temperature;
  • control parameters.

Give the wavelength, frequency, amplitude, and occupation range over which the mode description is intended.

Calling every many-body excitation collective

Section titled “Calling every many-body excitation collective”

All exact excitations belong to a many-body Hilbert space. Collectivity requires a coherent coordinate, eigenvector, self-consistent field, hydrodynamic variable, or comparable organizing structure.

A single quantum of an extended normal mode is collective. Macroscopic occupation creates a classical wave but is not the definition.

Plasmons, optical phonons, amplitude modes, and relative-phase modes can be gapped.

Diffusion and overdamped order-parameter relaxation are collective. Sharpness determines quasiparticle usefulness, not collectivity alone.

A finite-system single-particle transition is a pole. Collective identification also requires operator structure, polarization, scaling, or a coherent branch.

A zero of a response denominator can be dark in the chosen channel or canceled by a vanishing numerator. Compute the residue.

Using a scalar response for coupled variables

Section titled “Using a scalar response for coupled variables”

Amplitude, phase, density, current, sublattice, and component channels can mix. Diagonalizing one diagonal entry can miss the true eigenmodes.

Near an avoided crossing, mode character transfers between energy-ordered branches. Track eigenvectors and spectral weight.

A threshold can produce a strong peak or cusp without an isolated pole. Test line shape, analytic structure, and resolution dependence.

Calling every amplitude feature a Higgs mode

Section titled “Calling every amplitude feature a Higgs mode”

Establish an amplitude-like eigenvector, symmetry channel, pole or controlled resonance, and separation from pair-breaking or multiparticle continua.

A neutral phase mode and a charged longitudinal mode can have qualitatively different long-wavelength spectra.

The distinction depends on ωτcoll\omega\tau_{\mathrm{coll}}. A formula from one regime should not be applied in the other without a crossover treatment.

RPA can reveal real collective physics while remaining an approximation with missing exchange, vertex, local-field, or strong-correlation effects.

Participation depends on basis, metric, and normalization. Use it with physical operator and scaling information.

Exercise 1: In-phase and out-of-phase oscillations

Section titled “Exercise 1: In-phase and out-of-phase oscillations”

For the coupled-oscillator Lagrangian

L=m2(x˙12+x˙22)−k02(x12+x22)−kc2(x1−x2)2,\begin{aligned} L ={}& \frac{m}{2} \left( \dot x_1^2+\dot x_2^2 \right) \\ &- \frac{k_0}{2} \left( x_1^2+x_2^2 \right) \\ &- \frac{k_c}{2} \left( x_1-x_2 \right)^2, \end{aligned}

derive the normal coordinates and frequencies. Which mode couples to a uniform force F(x1+x2)F(x_1+x_2)?

Solution

Use

Q±=x1±x22.Q_\pm = \frac{x_1\pm x_2}{\sqrt2}.

The inverse transformation is

x1=Q++Q−2,x2=Q+−Q−2.x_1 = \frac{Q_++Q_-}{\sqrt2}, \qquad x_2 = \frac{Q_+-Q_-}{\sqrt2}.

Then

x12+x22=Q+2+Q−2x_1^2+x_2^2 = Q_+^2+Q_-^2

and

(x1−x2)2=2Q−2.\left( x_1-x_2 \right)^2 = 2Q_-^2.

Therefore

L=m2(Q˙+2+Q˙−2)−k02Q+2−k0+2kc2Q−2.\begin{aligned} L ={}& \frac{m}{2} \left( \dot Q_+^2+\dot Q_-^2 \right) \\ &- \frac{k_0}{2}Q_+^2 \\ &- \frac{k_0+2k_c}{2}Q_-^2. \end{aligned}

The frequencies are

ω+=k0m,ω−=k0+2kcm.\omega_+ = \sqrt{\frac{k_0}{m}}, \qquad \omega_- = \sqrt{ \frac{k_0+2k_c}{m} }.

The uniform force couples through

F(x1+x2)=2FQ+.F(x_1+x_2) = \sqrt2FQ_+.

It excites only the in-phase mode. The out-of-phase mode exists but is dark to a perfectly uniform symmetric source.

Exercise 2: A denominator zero can be dark

Section titled “Exercise 2: A denominator zero can be dark”

Suppose the response in the Q±Q_\pm basis is

χ±R(ω)=1Ω±2−(ω+i0+)2.\chi_\pm^{\mathrm R}(\omega) = \frac{1}{ \Omega_\pm^2 - (\omega+i0^+)^2 }.

A source and detector both couple only to O=X1+X2O=X_1+X_2. Show which pole appears in χOOR\chi_{OO}^{\mathrm R}.

Solution

Since

O=X1+X2=2Q+,O = X_1+X_2 = \sqrt2Q_+,

the source and detector have no projection onto Q−Q_-. Therefore

χOOR=2χ+R.\chi_{OO}^{\mathrm R} = 2\chi_+^{\mathrm R}.

The poles at ω=±Ω+\omega=\pm\Omega_+ appear, while the poles at ±Ω−\pm\Omega_- have zero residue in this channel.

The antisymmetric mode is not absent from the system. It requires an antisymmetric source or detector such as X1−X2X_1-X_2.

For

Q¨+2κQ˙+Ω02Q=0,\ddot Q + 2\kappa\dot Q + \Omega_0^2Q = 0,

classify the motion for:

  1. κ=Ω0/10\kappa=\Omega_0/10;
  2. κ=Ω0\kappa=\Omega_0;
  3. κ=2Ω0\kappa=2\Omega_0.

Find the poles in each case.

Solution

The poles are

ω±=−iκ±Ω02−κ2.\omega_\pm = -i\kappa \pm \sqrt{ \Omega_0^2-\kappa^2 }.

For κ=Ω0/10\kappa=\Omega_0/10,

ω±=−iΩ010±9910Ω0.\omega_\pm = -\frac{i\Omega_0}{10} \pm \frac{\sqrt{99}}{10} \Omega_0.

The poles have nonzero real parts and the mode is underdamped.

For κ=Ω0\kappa=\Omega_0,

ω+=ω−=−iΩ0.\omega_+ = \omega_- = -i\Omega_0.

This is critical damping in the elementary model.

For κ=2Ω0\kappa=2\Omega_0,

ω±=−2iΩ0±i3 Ω0.\omega_\pm = -2i\Omega_0 \pm i\sqrt3\,\Omega_0.

Both poles lie on the negative imaginary axis:

ω±=−i(2∓3)Ω0.\omega_\pm = -i \left( 2\mp\sqrt3 \right) \Omega_0.

The response is overdamped and relaxes as a sum of two exponentials.

Starting from

∂tn+∇⋅j=0\partial_t n + \nabla\cdot\mathbf j = 0

and

j=−D∇n,\mathbf j = -D\nabla n,

derive the dispersion and explain why the mode is collective but not a ballistic quasiparticle.

Solution

Substitution gives

∂tn=D∇2n.\partial_t n = D\nabla^2n.

For

δn=δn0ei(q⋅r−ωt),\delta n = \delta n_0 e^{i(\mathbf q\cdot\mathbf r-\omega t)},

one finds

−iω=−Dq2,-i\omega = -Dq^2,

so

ω=−iDq2.\omega = -iDq^2.

The mode is collective because it describes coherent long-wavelength relaxation of the conserved density field. Its pole has no real oscillation frequency, so there is no ballistic wave packet with a group velocity. Particle-like quasiparticle language is therefore generally inappropriate.

Use

∂tδn+n0∇⋅v=0\partial_t\delta n + n_0\nabla\cdot\mathbf v = 0

and

mn0∂tv=−∇[(∂P∂n) ⁣sδn]mn_0\partial_t\mathbf v = -\nabla \left[ \left( \frac{\partial P}{\partial n} \right)_{\!s} \delta n \right]

to derive the sound speed.

Solution

Differentiate continuity with respect to time:

∂t2δn+n0∇⋅∂tv=0.\partial_t^2\delta n + n_0\nabla\cdot \partial_t\mathbf v = 0.

From the Euler equation,

∂tv=−1mn0(∂P∂n) ⁣s∇δn.\partial_t\mathbf v = - \frac{1}{mn_0} \left( \frac{\partial P}{\partial n} \right)_{\!s} \nabla\delta n.

Substitution gives

∂t2δn−1m(∂P∂n) ⁣s∇2δn=0.\partial_t^2\delta n - \frac{1}{m} \left( \frac{\partial P}{\partial n} \right)_{\!s} \nabla^2\delta n = 0.

Thus

∂t2δn=cs2∇2δn\partial_t^2\delta n = c_s^2\nabla^2\delta n

with

cs2=1m(∂P∂n) ⁣s.c_s^2 = \frac{1}{m} \left( \frac{\partial P}{\partial n} \right)_{\!s}.

Plane waves have

ω±=±csq.\omega_\pm = \pm c_sq.

Exercise 6: Participation is basis dependent

Section titled “Exercise 6: Participation is basis dependent”

In a site basis, let

ej=1Ne_j = \frac{1}{\sqrt N}

for every jj. Compute the participation ratio. Then transform to a basis whose first vector is exactly e\mathbf e. What participation ratio does the same mode have there?

Solution

In the site basis,

∑j∣ej∣4=N(1N2)=1N.\sum_j \lvert e_j\rvert^4 = N \left( \frac{1}{N^2} \right) = \frac1N.

Therefore

P=11/N=N.P = \frac{1}{1/N} = N.

Now choose an orthonormal basis with first basis vector u1=e\mathbf u_1=\mathbf e. The same mode has components

(1,0,…,0).\left( 1,0,\ldots,0 \right).

In that basis,

P=1.P = 1.

The physics did not change. The example proves that participation ratio is basis dependent and must be interpreted in a physically motivated local or component basis.

Exercise 7: Avoided crossing and character exchange

Section titled “Exercise 7: Avoided crossing and character exchange”

For

K=(ω−ω1−g−gω−ω2),\mathcal K = \begin{pmatrix} \omega-\omega_1 & -g \\ -g & \omega-\omega_2 \end{pmatrix},

derive the two frequencies. At ω1=ω2=ω0\omega_1=\omega_2=\omega_0, find normalized eigenvectors.

Solution

The determinant condition is

(ω−ω1)(ω−ω2)−g2=0.\left( \omega-\omega_1 \right) \left( \omega-\omega_2 \right) - g^2 = 0.

Solving the quadratic gives

ω±=ω1+ω22±(ω1−ω22)2+g2.\begin{aligned} \omega_\pm ={}& \frac{\omega_1+\omega_2}{2} \\ &\pm \sqrt{ \left( \frac{\omega_1-\omega_2}{2} \right)^2 +g^2 }. \end{aligned}

At degeneracy,

ω±=ω0±∣g∣.\omega_\pm = \omega_0 \pm \lvert g\rvert.

For g>0g>0, the normalized eigenvectors can be chosen as

e+=12(1−1),e−=12(11).\mathbf e_+ = \frac{1}{\sqrt2} \begin{pmatrix} 1\\-1 \end{pmatrix}, \qquad \mathbf e_- = \frac{1}{\sqrt2} \begin{pmatrix} 1\\1 \end{pmatrix}.

The splitting is 2g2g. Away from the crossing, the eigenvectors rotate continuously between the original coordinate characters.

Exercise 8: Audit a claimed collective mode

Section titled “Exercise 8: Audit a claimed collective mode”

A numerical calculation on one finite lattice shows a broad maximum in a density structure factor at one momentum. The plot uses artificial Lorentzian broadening η=0.15J\eta=0.15J. The fitted FWHM is 0.18J0.18J, and the feature is labeled a long-lived collective mode.

List the missing checks.

Solution

The fitted width is comparable to the artificial broadening, so it does not establish an intrinsic lifetime. A broad maximum at one momentum can be:

  • one broadened discrete line;
  • several unresolved levels;
  • a continuum threshold;
  • a density-of-states enhancement;
  • the finite-size precursor of a resonance;
  • a genuine collective branch.

A responsible analysis should:

  1. repeat the calculation for several η\eta values;
  2. inspect raw finite-size energies and matrix elements;
  3. scale system size and boundary conditions;
  4. trace the feature across momentum;
  5. compare with the particle–hole or multiparticle continuum;
  6. identify its polarization or dominant operator eigenchannel;
  7. check spectral moments and total weight;
  8. compare density response with another appropriate channel;
  9. test a linearized effective or hydrodynamic prediction;
  10. show that an intrinsic width remains after resolution and finite-size effects are controlled.

Until then, “candidate density-response feature” is more accurate than “long-lived collective mode.”

  • A collective mode is a coherent dynamical pattern or response eigenchannel built from many microscopic degrees of freedom.
  • Collectivity is defined relative to variables, background, scale, approximation, and probe.
  • A mode can be identified as a collective coordinate, a homogeneous solution of linearized dynamics, or a pole of a response matrix.
  • Frequency, polarization, residue, damping, and regime are independent pieces of mode data.
  • Normal mode, collective mode, quasiparticle, resonance, coherent state, and exact eigenstate are related but distinct concepts.
  • One quantum can be collective; macroscopic occupation is not required.
  • Conservation laws produce hydrodynamic modes such as diffusion and sound.
  • Broken symmetries can enforce Goldstone modes, but many collective modes are gapped or unrelated to symmetry breaking.
  • Long-range Coulomb feedback can convert a neutral acoustic phase or density mode into a gapped plasma oscillation.
  • Coupled channels hybridize, avoid crossings, exchange character, and can be bright or dark depending on the probe.
  • A continuum can damp or dissolve a mode; a strong threshold feature is not automatically an isolated pole.
  • Participation ratio is useful but basis dependent.
  • A convincing identification combines dispersion, polarization, residue, intrinsic width, sum rules, scaling, and a controlled dynamical model.