---
title: Dynamical Structure Factors in Quantum Systems
url: https://www.emergentmind.com/topics/dynamical-structure-factors-dsfs
type: topic
---

# Dynamical Structure Factors in Quantum Systems

The dynamical structure factor (DSF) is a fundamental two-point correlation function that encodes the quantum and thermal dynamics of density (or other local operators) in many-body systems. It is central in diverse contexts ranging from inelastic neutron and X-ray scattering to cold-atom Bragg spectroscopy and quantum simulation. The DSF accesses the spectral content and momentum structure of operator fluctuations and collective excitations, providing direct links between microscopic theory and experimental observables.

## 1. Formal Definition and Physical Significance

In a translationally invariant quantum lattice system, the DSF for the density operator is defined as
$$
S(q, \omega) = \int_{-\infty}^{\infty} dt \, e^{i\omega t} \langle \rho_q(t) \, \rho_{-q}(0) \rangle,
$$
where $\rho_q(t) = \sum_j e^{-iqj} n_j(t)$ is the spatial Fourier transform of the number operator, and the average is taken at thermal equilibrium (or on the ground state at $T=0$). This form captures the propagation of collective excitations at wavevector $q$ and frequency $\omega$ [2409.07030].

The DSF more generally extends to any local operator $O_j$, with
$$
S_O(q, \omega) = \int_{-\infty}^{\infty} dt \, e^{i\omega t} \langle O_q^\dagger(t) \, O_q(0) \rangle
$$
and admits the Lehmann spectral representation
$$
S_O(q, \omega) = \sum_f |\langle\psi_f | O_q | \psi_0\rangle|^2 \delta(\omega - (E_f - E_0)),
$$
where $|\psi_0\rangle$ is the ground state and $|\psi_f\rangle$ are excited states [1912.06076]. DSFs provide direct access to the excitation spectra, collective mode structure, and many-body response of the system.

In experimental terms, $S(q, \omega)$ is the central quantity accessed by inelastic scattering, e.g., the differential cross section of neutron or X-ray scattering is proportional to the DSF. In cold-atom systems, DSFs quantify the excitability of the system by weak density probes, encoding signals such as phonons, gaps, and spinon continua.

## 2. Measurement Theory: From Projective to Weak Measurements

Conventional measurement of unequal-time, two-point correlators underlying the DSF is fundamentally limited by projective (strong) measurement: a strong measurement of $n_j$ collapses the wavefunction and destroys the quantum coherence required to access $⟨n_j(t_2) n_k(t_1)⟩$ for $t_2 \neq t_1$. Therefore, standard protocols relying on projective detection are insufficient for reconstructing the DSF [2409.07030].

To circumvent this, weak measurement protocols have been developed, notably in cold-atom experiments. The approach is structured as follows:

- **Weak measurement protocol:** Two brief, spatially resolved weak (homodyne) measurements are performed, separated in time by $\Delta t$. Each measurement yields a noisy outcome,
  $$
  n_{j, t} = \langle n_j(t) \rangle + m_{j, t}/(2\Gamma^{1/2}),
  $$
  with $\Gamma \ll 1$ the measurement strength and $m_{j, t}$ Gaussian white noise.

- **Ensemble cross-correlation:** Repeatedly sampling the noisy outputs, the noise-averaged cross-correlation
  $$
  C_w(j, j'; \Delta t) = \langle n_{j,0} n_{j',\Delta t} \rangle_\mathrm{noise}
  $$
  reconstructs $\mathrm{Re} \langle n_j(0) n_{j'}(\Delta t)\rangle$ to leading order in $\Gamma$.

- **Spatial and temporal Fourier analysis:** By Fourier transforming both in space and time,
  $$
  S(q, \omega) = \sum_{\Delta j} \int_0^\infty d\Delta t\, e^{i[\omega\Delta t - q\Delta j]} G_{\Delta j}(\Delta t),
  $$
  with $G_{\Delta j}(\Delta t)$ the spatially averaged cross-correlation, one recovers the DSF [2409.07030].

This protocol enables direct, minimally invasive measurement of $S(q, \omega)$ in quantum gases without the need for explicit external driving at $(q,\omega)$.

## 3. Computational and Experimental Protocols

Several operational routes to DSF estimation in quantum simulators and experiments have emerged:

- **Weak measurement in cold atoms (homodyne/PCI):** Time-separated weak phase-contrast imaging of site populations, as detailed above, allows direct DSF estimation in systems modeled by, e.g., the Bose–Hubbard Hamiltonian [2409.07030].

- **Quantum simulation protocols:** For generic lattice models, a Ramsey-style sequence can map unequal-time correlators to single-time observables. For example, in spin systems, applying a local $\pi/4$ pulse and evolving under the Hamiltonian, measurement of single-spin observables yields retarded correlators whose spatial and temporal Fourier transform reconstructs the DSF [1912.06076].

- **Simulation and scaling considerations:**
  - Statistical errors in weak measurement scale as $\Gamma^{-1}$; systematic back-action errors scale $\propto\Gamma^{1/2}$. An experimentally practical compromise has been demonstrated at $\Gamma\sim0.1$ for sample sizes on the order of tens.
  - The frequency and momentum resolution depend on the sampling interval ($\Delta t$) and total observation time ($T$), as well as numerical aperture (NA) and depth-of-field for imaging systems.
  - In practical cold-atom experiments, phase-contrast imaging with sufficient NA captures long-wavelength modes with percent-level accuracy, with artifacts confined to high-$q$ components [2409.07030].

- **Matrix product state (MPS) and classical simulations:** Numerical verification employs MPS-based sampling of measurement protocols, where ground-state preparation and time evolution are performed on, e.g., the 1D Bose–Hubbard model. The reconstruction matches exact DSF results both in the gapless superfluid and gapped Mott-insulating regimes [2409.07030].

## 4. Theoretical Structure and Exact Results

The formal theoretical structure of DSFs reveals rich connections to many-body physics:

- **Sum rules:** The DSF satisfies the $f$-sum rule,
  $$
  \int_{-\infty}^{\infty} d\omega\, \omega\, S(q, \omega) = N \frac{q^2}{2m},
  $$
  ensuring consistency with total spectral weight from fluctuation–dissipation relations [1912.06076].

- **Spectral properties:** The DSF encodes the spectrum of elementary excitations, including phonon branches, Mott gaps, multi-magnon continua, and signatures of fractionalization or topology, depending on the model. In the 1D Bose–Hubbard model, the DSF directly displays the transition from a gapless linear phonon branch in the superfluid to a finite Mott gap in the insulating regime [2409.07030].

- **Computational complexity:** For general local Hamiltonians, estimation of the DSF is BQP-hard—classically intractable—since accurate evaluation of unequal-time correlators would allow simulation of arbitrary quantum circuits. This computational hardness is inherited by DSFs under any experimental or numerically practical approximation within polynomial precision [1912.06076].

- **Practical quantum advantage:** Quantum simulators operating with moderate system sizes ($N\sim50$–70), commonly accessed in noisy intermediate-scale quantum (NISQ) hardware, already explore DSFs beyond the reach of classical diagonalization or Krylov subspace approaches [1912.06076].

## 5. Specific Applications, Error Analysis, and Experimental Implementation

The DSF has been computed and measured in a variety of contexts:

- **Cold-atom systems:** The DSF measured from PCI/weak-measurement protocols exhibits quantitative agreement with exact MPS correlators for the 1D Bose–Hubbard model, including the identification of phononic and Mott-insulating signatures [2409.07030].

- **Error sources and trade-offs:**
  - Small $\Gamma$ improves linearity but requires more trajectories for statistical convergence; large $\Gamma$ increases backaction errors. An empirical optimum is near $\Gamma \approx 0.1$ with $\sim50$ trajectories.
  - Limited imaging NA primarily affects high-$q$ response; phonon features and low-$q$ response remain robust, even with finite depth-of-field and aberrations [2409.07030].

- **Temporal and spatial resolution:** Adequate frequency resolution ($\Delta\omega$) requires long observation windows, while momentum resolution is determined by the Fourier-space coverage of the imaging system.

- **Broader platforms:** Analogous protocols apply to trapped ions (e.g., Raman addressing and fluorescence imaging), Rydberg atom arrays (site-addressed pulses and van-der-Waals evolution), and superconducting qubits (microwave pulse sequences and dispersive readout), all yielding DSFs via similar Fourier analysis of measured correlators [1912.06076].

## 6. Physical Interpretation and Impact

Measurement of $S(q, \omega)$ provides direct access to the fundamental dynamical properties:

- **Dynamical response and excitation spectrum:** Peaks and gaps in $S(q, \omega)$ map onto elementary collective modes (e.g., phonons, rotons, magnons) and many-body spectral features (e.g., Mott gaps, spinon continua).
- **Probing of phase transitions:** The DSF provides clear spectroscopic fingerprints for phase transitions, such as closing of the phonon gap at the superfluid–Mott-insulator transition or the appearance of critical signatures and mode softening.
- **Benchmarking and quantum simulation:** The DSF is a stringent benchmark for both theoretical models and experimental platforms, enabling verification of quantum advantage and the study of real-time many-body dynamics in regimes inaccessible to classical computation.

The development of weak measurement protocols, along with scalable quantum-simulation-based measurement, now offers a minimally invasive, all–measurement-based route to the full dynamical response function, positioning the DSF as a central observable in quantum many-body physics [2409.07030, 1912.06076].

---

**References:**
- "Dynamical Structure Factor from Weak Measurements" [2409.07030]
- "Dynamical structure factors of dynamical quantum simulators" [1912.06076]

Source: https://www.emergentmind.com/topics/dynamical-structure-factors-dsfs