**Dispersion and density of states explain bands beyond the zone center.** The set of $\omega_s(\mathbf{q})$ values across the Brillouin zone forms the phonon dispersion. Its slopes near the zone center determine acoustic group velocities, while avoided crossings, soft branches, and flat regions reveal coupling, instability, or high density of states. Raman spectroscopy normally sees only a restricted projection of this band structure, not the complete dispersion. Second-order Raman processes create or annihilate two phonons whose wavevectors sum appropriately. They can sample the Brillouin zone and produce overtone or combination bands shaped by joint phonon density of states and matrix elements. Defects, finite size, interfaces, superlattice periodicity, and disorder relax momentum conservation, activating non-zone-center phonons in nominally first-order spectra. A broad band matching a calculated density-of-states maximum is suggestive, but assignment still requires energy, symmetry, excitation, and defect controls. Inelastic neutron and x-ray scattering measure energy versus momentum more directly, with different cross sections and sample requirements. Electron energy-loss and ultrafast methods can access still other regions or populations. First-principles calculations connect these measurements by predicting eigenvectors, dispersion, Raman activities, Born effective charges, dielectric response, and anharmonic couplings. Agreement at one zone-center frequency is insufficient validation of an entire calculated phonon model. |Analysis layer|Primary observable|What it constrains|Common ambiguity|Strongest cross-check| |---|---|---|---|---| |Mode inventory and symmetry|Peak count, polarization, Raman/IR activity|Phase and point-group consistency|Weak, overlapped, or geometry-forbidden modes|Group theory plus polarized Raman and IR| |Frequency and splitting|Peak centers and degeneracy lifting|Stress, temperature, composition, symmetry breaking|Several variables shift the same mode|Multiple modes with independent coefficients| |Line shape and linewidth|FWHM, asymmetry, continuum coupling|Lifetime, disorder, carriers, confinement|Instrument resolution and unresolved components|Resolution standard, temperature series, alternate model| |Intensity and excitation profile|Area versus polarization or laser energy|Tensor elements and resonance coupling|Optical interference, absorption, focus, texture|Response calibration and layered optical model| |Dispersion or two-phonon bands|Energy versus momentum or broad combination structure|Force constants and lattice dynamics|Matrix-element weighting and defects|Neutron/x-ray data or converged first-principles calculation| **Peak position is a state variable with multiple sensitivities.** A measured mode frequency can be expanded locally around a reference state as $$ \Delta\omega_s=\sum_{i,j}\Pi_{sij}\sigma_{ij}+\chi_{sT}\Delta T+\chi_{sc}\Delta c+\chi_{sn}\Delta n_c+\cdots $$ The terms represent stress through phonon deformation potentials, temperature, composition, carrier density, and other relevant variables. They need not be independent: temperature changes strain through thermal expansion, composition changes lattice constant and electronic resonance, and carriers can modify both phonon self-energy and local heating. One peak shift cannot uniquely solve several unknowns. For hydrostatic volume change, a mode Grüneisen parameter is commonly defined by $$ \gamma_s=-\frac{\partial\ln\omega_s}{\partial\ln V} $$ This scalar is useful for hydrostatic or quasiharmonic reasoning but does not replace the full deformation-potential tensor under arbitrary stress. Uniaxial or shear stress can split degenerate modes and rotate their eigenvectors; polarization-resolved spectra then contain more information than a scalar shift. The elastic constants and mechanical boundary condition are required to convert strain to stress. A universal silicon coefficient such as a single number of inverse centimeters per gigapascal is not valid across all wafers and geometries. The observed component depends on crystal orientation, incident and analyzed polarization, stress tensor, phonon deformation potentials, temperature, and sign convention. Calibration must match the substrate orientation, device geometry, and reference state. Similar caution applies to GaN, SiC, diamond, III–V compounds, oxides, and two-dimensional materials. Temperature shifts arise from explicit anharmonic phonon–phonon interactions and implicit thermal expansion. Linewidth usually grows as additional decay channels become populated, but defects, electron–phonon coupling, phase transitions, or changing resonance can create nonmonotonic behavior. A calibrated stage series at low probe power provides empirical $\omega_s(T)$ and linewidth relations for the actual material; literature coefficients are transferable only when composition, stress, carrier density, and measurement conditions are compatible. Composition can generate one-mode, two-mode, or mixed alloy behavior. In SiGe, for example, Si–Si, Si–Ge, and Ge–Ge-like bands contain composition and local-environment information, but strain, clustering, temperature, and resonance influence their frequencies and intensities. Standards or a joint multi-band model are stronger than inserting one band into a universal composition equation. **Line shape contains lifetime information only after deconvolution and model selection.** A damped harmonic oscillator has a susceptibility form such as $$ \chi_s''(\omega)\propto\frac{\Gamma_s\omega}{(\omega_s^2-\omega^2)^2+(\Gamma_s\omega)^2} $$ In a limited spectral region and weak-damping limit, a Lorentzian may approximate the band. Gaussian broadening can represent static inhomogeneity or instrument response, and a Voigt profile combines both phenomenologically. The measured linewidth is a convolution with spectrometer resolution, slit width, pixel sampling, laser linewidth, and spatial variation within the spot. Subtracting widths in quadrature is valid only for compatible line-shape assumptions. An intrinsic lifetime is related to a properly defined homogeneous linewidth, but conventions differ between angular frequency, ordinary frequency, wavenumber, half width, and full width. State the convention before using a relation like lifetime proportional to inverse linewidth. Inhomogeneous strain, composition variation, unresolved isotope components, grain orientations, or temperature gradients broaden a peak without representing a shorter microscopic phonon lifetime. Fano asymmetry occurs when a discrete phonon interferes with a continuum, as in sufficiently carrier-rich semiconductors or electronically resonant systems. The asymmetry parameter, center, and width covary strongly with the background model. A Fano fit can describe an asymmetric peak without proving the continuum’s identity. Carrier-density extraction needs an appropriate coupled dielectric or microscopic model and independent electrical or optical evidence. Polar longitudinal optical phonons can couple to free-carrier plasmons, producing longitudinal optical phonon–plasmon coupled modes. Their frequencies and damping depend on carrier density, effective mass, mobility, dielectric constants, phonon parameters, and geometry. Assigning the strongest longitudinal feature to the uncoupled LO frequency can produce incorrect stress or composition when carrier coupling is substantial. **Disorder and confinement relax momentum selection rather than merely broadening peaks.** When translational coherence is finite, phonons away from the zone center contribute to first-order Raman scattering. A generic phonon-confinement spectrum can be written $$ I(\omega)\propto\int_{BZ}\frac{|C(\mathbf{q};L)|^2}{[\omega-\omega_s(\mathbf{q})]^2+(\Gamma_s/2)^2}\,d\mathbf{q} $$ The weighting $C(\mathbf{q};L)$ depends on the assumed real-space confinement and characteristic size $L$. The result follows the actual dispersion, so confinement may shift and asymmetrically broaden a band rather than simply add a symmetric width. Extracted size is model-dependent and should be validated by microscopy or diffraction. Nanocrystals add surface modes, interface modes, size distributions, stress, ligand coupling, alloying, and dielectric-environment effects. A fitted confinement size can absorb all of these missing variables. Core–shell structures may show distinct core, shell, surface, and interface response whose resonance conditions differ. Multiwavelength Raman and structural characterization help separate spatial regions and electronic selectivity. Amorphous materials lack long-range translational selection and produce broad bands related to vibrational density of states weighted by coupling coefficients. Labels borrowed from crystalline zone-center modes can be useful descriptors but should not imply the same eigenvectors. Nanocrystalline spectra often mix crystalline peaks, grain-boundary response, amorphous background, and stress distributions. Defects can activate otherwise forbidden or finite-wavevector modes and also change electronic resonance. Intensity ratios used as defect metrics are generally regime- and material-specific. At high defect density, coherent domains shrink, peaks overlap, and the same ratio may reverse trend. A calibration must span the intended microstructure and use an independent defect or domain-size measurement. **First-principles calculations require convergence and physical validation.** Density-functional perturbation theory, finite-displacement force constants, and molecular-dynamics correlation methods can produce vibrational spectra. The result depends on structure, exchange-correlation approximation, pseudopotentials, supercell, displacement size, wavevector sampling, plane-wave cutoff, electronic smearing, non-analytic polar correction, and convergence of forces. Imaginary frequencies are usually plotted as negative values and indicate negative curvature of the modeled energy surface, but they can mean a real structural instability or a numerical problem. Acoustic sum-rule violations, insufficient supercell size, loose relaxation, poor sampling, or omitted long-range electrostatics can create spurious soft modes. Re-relaxation and systematic convergence should precede a phase-instability claim. Polar crystals need non-analytic corrections near the zone center to describe longitudinal-optical/transverse-optical splitting. Born effective charges and the high-frequency dielectric tensor enter this correction, and the limiting frequency depends on the direction of approach to the zone center. Comparing a calculation lacking this term with an experimental LO band is not a meaningful validation. Calculated Raman “activity” is not raw experimental intensity. Excitation energy, resonance, temperature population, scattered-frequency factors, polarization, absorption, optical interference, and instrument response transform it. Frequency scaling is also sometimes applied to compensate systematic computational error, but a fitted scale factor can hide mode-dependent deficiencies. Validate eigenvectors and symmetry assignments, not only scaled frequencies. ```flowchart Establish phase, structure, composition, orientation, and measurement geometry -> Predict zone-center irreducible representations and Raman/IR activity -> Acquire calibrated Raman, polarization, temperature, and reference spectra -> Separate substrate, fluorescence, instrument, and second-phase features -> Fit peaks with resolution convolution and alternative line-shape tests -> Compare frequencies, tensors, eigenvectors, and dispersion calculations -> Test stress, temperature, composition, carriers, disorder, and confinement -> Jointly fit multiple modes with independently constrained variables -> Validate conclusions using diffraction, microscopy, transport, or IR data ``` **A production analysis preserves spectra, coordinates, and competing explanations.** Record laser wavelength and power, spot and focus, polarization, objective, spectral resolution, grating, slit, detector, calibration, temperature, sample orientation, atmosphere, acquisition time, baseline, fit window, and line-shape convention. Store raw counts and residuals alongside peak tables. Automated fitting should reject saturation, cosmic rays, low signal, unresolved overlap, parameter-bound hits, and spatially implausible jumps. Report frequency, linewidth convention, integrated area, symmetry label, assignment confidence, and covariance or uncertainty. For derived stress, composition, carrier density, crystallite size, or temperature, report the calibration coefficients, reference state, optical and mechanical model, and uncertainty propagation. A database match should rank candidate phases but must not replace selection-rule and context checks. The strongest conclusion is the narrowest one supported by independent observables. A new peak plus polarization change and an unstable calculated eigenvector can support a structural transition; a single shift may support only a changed lattice state. Multiple modes and orthogonal measurements keep a plausible vibrational story from becoming a false material constant. The durable way to interpret phonon modes is through an eigenvector-symmetry-dispersion-line-shape-state-variable-confinement-calculation-and-identifiability lens.