---
title: Dynamical Power Density Spectrum
url: https://www.emergentmind.com/topics/dynamical-power-density-spectrum
type: topic
---

# Dynamical Power Density Spectrum

Searching arXiv for recent papers on dynamical power density spectra and closely related PSD formulations across fluids, stochastic processes, and dynamical systems.
A dynamical power density spectrum is a spectral characterization of fluctuations generated by an evolving system, but its precise meaning depends on the underlying dynamical context. In the frequency-domain formulation used for stochastic time series, it is the power spectral density of a time-dependent observable and quantifies the amount of fluctuations at a given frequency [2306.00417]. In the setting of semiperiodic pulse trains, it captures dynamical behavior such as oscillations, chaos, and intermittent bursts through the frequency power spectral density of the signal [2106.15904]. In the fluid-mechanical setting of two-dimensional Poiseuille and Couette flow with the van der Waals effect, the corresponding object is not a temporal spectrum but a spatial Fourier power spectrum in streamwise wavenumber, averaged over time, of small fluctuations about laminar stationary states [2603.19190]. Across these formulations, the unifying theme is that a dynamical power density spectrum resolves how fluctuation power is distributed over spectral variables and links that distribution to the governing dynamics, their instabilities, and their characteristic scales.

## 1. Definitions and mathematical formulations

In the standard frequency-domain formulation for a real scalar signal \( \Phi(t) \), the finite-time Fourier transform is
\[
\mathcal{F}_\Phi(\omega) = \int_{-T/2}^{T/2} \Phi(t)\,e^{-i\omega t}\,dt ,
\]
and the frequency power spectral density is
\[
S_\Phi(\omega) = \lim_{T\to\infty} \frac{1}{T}\,\big\langle \big|\mathcal{F}_\Phi(\omega)\big|^2\big\rangle .
\]
For a stationary process, this is equivalent, via the Wiener–Khinchin theorem, to the Fourier transform of the autocorrelation function \(R_\Phi(\tau)\) [2106.15904]. In the Markov-process formulation, for an observable \(z(\mathbf{x})\) measured along a trajectory, the PSD is defined by
\[
S^z(\omega) = \lim_{\tau \to \infty} \frac{1}{\tau} \left( \big\langle |\hat{z}_\tau(\omega)|^2 \big\rangle - \big|\langle \hat{z}_\tau(\omega)\rangle \big|^2 \right),
\]
with
\[
\hat{z}_\tau(\omega) = \int_0^\tau dt \, e^{i \omega t}\, z(\mathbf{x}(t)),
\]
and again
\[
S^z(\omega) = \int_{-\infty}^{\infty} dt \, e^{i \omega t} \, \mathrm{Cov}\big(z(t),z(0)\big)
\]
by Wiener–Khinchin [2306.00417].

A distinct but related formulation appears for stochastic fields governed by linear autonomous stochastic differential equations. For a field \(\phi\) satisfying
\[
(\mathcal{L}\phi)(\mathbf{x},t) = \xi(\mathbf{x},t),
\]
with Gaussian driving \(\xi\), the field covariance in Fourier space is diagonal under homogeneity, and the spectral density becomes
\[
P_\phi(k) = \frac{P_\xi(k)}{|f(k)|^2},
\]
where \(f(k)\) is the Fourier representation of the linear operator \(\mathcal{L}\) [1708.05250]. In that framework, the spectral density is treated as the central object that encodes the dynamics of linear autonomous SDEs.

A further variation arises in finite-time, single-trajectory analysis. For Brownian motion, the single-trajectory PSD for one component is
\[
S_T^{(j)}(f) =\frac{1}{T}\int_0^T\!dt_1\int_0^T\!dt_2\, \cos\big(f(t_1-t_2)\big)\,X_{t_1}^{(j)}X_{t_2}^{(j)},
\]
which is a random functional of a single path rather than an ensemble average [1801.02986]. For active Ornstein–Uhlenbeck particles, the finite-time PSD is defined by
\[
\mu(f,T)=\frac{1}{T}\Big\langle \Big|\int_{0}^{T} X(t)\,e^{-ift}\,dt\Big|^{2}\Big\rangle,
\]
which remains meaningful for both stationary and non-stationary dynamics [2605.09399].

In the fluid case, the term is used in a spatial rather than temporal sense. For a scalar field \(f(x,y,t)\), the procedure is: restrict to a transverse band, average over \(y\), compute a one-dimensional discrete Fourier transform in \(x\), take the modulus, then perform a time average over a late-time window. The resulting spectra \(\langle |\hat f(k)|^2 \rangle\) are spatial power spectra in streamwise wavenumber, averaged over time, and the reported power laws have the form
\[
P_f(k)\sim k^{-\alpha}
\]
[2603.19190]. This usage makes explicit that “dynamical” need not imply temporal Fourier analysis; it may refer instead to the spectrum of time-averaged fluctuations generated by the dynamics.

## 2. Relation to dynamics, covariance, and governing equations

The defining feature of a dynamical power density spectrum is its direct relation to the system’s equations of motion and correlation structure. For stationary stochastic dynamics, the PSD is the Fourier transform of the covariance function, so peaks, plateaus, and power-law tails correspond to characteristic timescales and correlation structures of the dynamics [2605.09399]. In overdamped diffusions and continuous-time Markov jump processes, the low-frequency limit satisfies
\[
S^z(0) = \lim_{\tau \to \infty} \Big( \tau \, \mathrm{Var}(\bar{z}_\tau) \Big),
\qquad
\bar{z}_\tau = \frac{1}{\tau}\int_0^\tau dt\,z(t),
\]
while the high-frequency limit obeys
\[
\lim_{\omega \to \infty} \big(\omega^2 S^z(\omega)\big)
= 2\langle \nabla z \cdot \mathbf{B} \nabla z \rangle_{\mathrm{st}}
\]
[2306.00417]. The spectrum therefore interpolates between long-time fluctuation statistics and short-time increments.

For linear stochastic differential equations, the same connection is expressed through the transfer structure of the dynamics. In the homogeneous, stationary, Gaussian setting, the observable PSD is \(P_\phi(k)=P_\xi(k)/|f(k)|^2\), so the deterministic operator contributes through \(|f(k)|^{-2}\), while the driving appears through \(P_\xi(k)\) [1708.05250]. For the damped oscillator
\[
(\alpha \partial_t^2 + \beta \partial_t + m^2)\,\phi(t) = \xi(t),
\]
with white noise, the PSD becomes
\[
P_\phi(\omega)=\frac{1}{(\gamma - \alpha\omega^2)^2 + (\beta\omega)^2},
\]
making resonance, damping, and high-frequency decay explicit in spectral form [1708.05250].

For mechanical Langevin oscillators, the power spectral density likewise emerges from the dynamical susceptibility. For the simple harmonic Langevin oscillator,
\[
\ddot{x} + \Gamma \dot{x} + \omega_0^2 x = \frac{F_{\mathrm{th}}(t)}{m},
\]
the PSD is
\[
S_x(\omega) = \frac{2 \Gamma k_B T / \pi m}{(\omega^2 - \omega_0^2)^2 + \Gamma^2 \omega^2},
\]
which is the baseline for the paper’s time-dependent and anharmonic generalizations [2305.19260].

In inertial channel flow with the van der Waals effect, the spectral behavior is tied to a specific instability mechanism. The governing variables are density \(\rho\), divergence \(\chi=\nabla\cdot\mathbf{u}\), and vorticity \(\omega=\nabla^\perp\cdot\mathbf{u}\), and the divergence equation contains the nonlinear term \(-2\det(\nabla\mathbf{u})\) together with the van der Waals coupling
\[
\frac{4p_0}{\rho_{HS}}\nabla\cdot\left(\frac{\nabla\rho}{\rho}\right).
\]
These two together generate an “inertial-range” direct cascade in \(\rho\) and \(\chi\), and the resulting power-law spectra are therefore linked to the dynamics of compressive modes rather than to rotational turbulence in the incompressible sense [2603.19190].

## 3. Canonical spectral structures

Several recurrent spectral structures appear across the literature. A strictly periodic signal with period \(p\) has a Dirac-comb spectrum,
\[
S_\Phi(\omega)=\sum_{k=-\infty}^{\infty} W_k\,\delta(\omega-k\Omega_0),
\qquad
\Omega_0=\frac{2\pi}{p},
\]
with weights determined by the harmonic content of the waveform [2106.15904]. In the semiperiodic pulse framework, the full PSD factorizes into a pulse-spectrum envelope and a forcing spectrum determined by amplitude and arrival statistics:
\[
S_\Phi(\omega)=I_2\,\varrho_\varphi(d\omega)\times S_{\rm forcing}(\omega).
\]
Strict periodicity plus nonzero mean amplitude yields a Dirac comb modulated by the pulse spectrum; jitter or renewal waiting times suppress or broaden the comb, leaving mainly the pulse spectrum [2106.15904].

For Brownian motion, the ensemble PSD is
\[
\mu_S^{(j)}(f,\infty)=\frac{4D}{f^2},
\]
so the canonical diffusive signature is a \(1/f^2\) spectrum [1801.02986]. For active Ornstein–Uhlenbeck particles in free space, activity does not alter the Brownian \(f^{-2}\) exponent, but modifies its amplitude and introduces a crossover at the persistence frequency. The exact infinite-time PSD is
\[
\mu(f,\infty) = \frac{4D + 4 D_A }{f^{2}}  -  \frac{2D_A}{f^2 + \tau_A^{-2}},
\]
so both low and high frequencies remain \(f^{-2}\), with different amplitudes separated by \(f\sim \tau_A^{-1}\) [2605.09399].

Under harmonic confinement, richer structures arise. The stationary AOUP PSD is
\[
\mu(f,\infty) = \frac{2D}{f^{2} + \tau_R^{-2}} + \frac{2 v_p ^2 \tau_A^{-1} }{ \big(f^2 + \tau_A^{-2}\big)\big(f^2 + \tau_R^{-2}\big)},
\]
which produces two characteristic signatures absent in both thermal systems and free AOUPs: a two-plateau structure from a double-trapping mechanism due to two noise sources, and a new \(f^{-4}\) spectral scaling associated with transient ballistic motion [2605.09399]. This suggests that spectral exponents can diagnose not only diffusive versus ballistic motion but also the coexistence of multiple relaxation mechanisms.

In trapped nano-oscillator models, oscillatory frequency modulation leads to sideband spectra. For
\[
\omega(t)=\omega_0+\Delta\omega\cos(\Omega t),
\]
the PSD becomes a sum of Lorentzians centered at \(\omega_n=\omega_0+n\Omega\), with weights \(J_n^2(\Delta\omega/\Omega)\), so the power spectrum develops a Floquet-like sideband structure [2305.19260]. By contrast, slow frequency drift broadens the PSD into an almost flat band, while the quadrature PSD remains approximately Lorentzian [2305.19260].

In laminar channel flow with the van der Waals effect, the spectra are power laws in wavenumber. For Poiseuille flow, the paper reports
\[
P_{\rho^2}(k)\sim k^{-7/3},\qquad
P_{\chi^2}(k)\sim k^{-5/3},\qquad
P_{\omega^2}(k)\sim k^{-8/3},
\]
while for Couette flow
\[
P_{\rho^2}(k)\sim k^{-7/3},\qquad
P_{\chi^2}(k)\sim k^{-4/3},\qquad
P_{\omega^2}(k)\sim k^{-7/3}.
\]
The exponents depend both on the variable and on the background shear profile [2603.19190].

## 4. Measurement, estimation, and inference

The practical construction of a dynamical power density spectrum depends strongly on the data model. In the semiperiodic pulse framework, the PSD is computed from the finite-time Fourier transform and interpreted through a filtered point-process model, in which the pulse shape determines the envelope and the arrival statistics determine line structure, broadening, or its disappearance [2106.15904]. This yields a direct statistical interpretation of departures from a spectral Dirac comb.

For steady-state Markov processes, the PSD can be constrained even without a complete dynamical model. The spectrum is bounded by rational functions involving two constants \(C^0\) and \(C^\infty\):
\[
\frac{1}{C^\infty + \omega^2 C^0}
\le
\frac{S^z(\omega)}{2\mathrm{Var}_{\mathrm{st}}(z)}
\le
\frac{1}{\dfrac{1}{C^0} + \dfrac{\omega^2}{C^\infty}},
\]
and in equilibrium these constants are the low- and high-frequency limits of the normalized PSD [2306.00417]. Out of equilibrium, the same bounds hold but \(C^0\) and \(C^\infty\) no longer coincide with the limiting values, which permits intermediate-frequency peaks associated with oscillatory modes [2306.00417].

For linear stochastic fields, spectral density inference can be posed as an inverse problem. In the Information Field Theory formulation, the field covariance is diagonal in Fourier space and the spectral density is parameterized as
\[
P_\phi(k)=\exp\big[\tau(k)+\tan(\delta(k))\big].
\]
The smooth component \(\tau\) captures the background spectral shape, while \(\tan(\delta)\) captures divergent or peaked features. Noisy or incomplete observations are modeled by
\[
d_u = R_u[\phi] + n_u,
\]
and the spectrum is inferred by minimizing a marginal Hamiltonian after analytically marginalizing over the field [1708.05250]. This provides a non-parametric route to reconstructing dynamical spectra from noisy and masked data.

Mechanical-vibration analysis offers a different estimation paradigm. High Order Dynamic Mode Decomposition represents data as
\[
x_k = \sum_{j=1}^{r} \phi_j \lambda_j^k b_j,
\]
or equivalently
\[
x(t) = \sum_{j=1}^{r} \phi_j e^{\mu_j t} b_j,
\qquad
\lambda_j=e^{\mu_j\Delta t},
\]
with \(\mu_j=\alpha_j+i\omega_j\) encoding damping and oscillation frequency [2306.10864]. A Kernel Density Spectrum is then formed by replacing each modal line with a kernel centered at the HODMD frequency, using Gaussian or Lorentzian kernels and weights such as \(A_k^2\) or \(A_kT_k\). This produces a high-resolution, tunable spectral density from identified damped modes rather than FFT bins [2306.10864].

In observational cosmology, the relevant objects are power spectra of galaxy density and momentum fields in redshift space. Weighted density and momentum fields are defined from galaxy catalogues, and the density–momentum cross power spectrum multipoles are estimated alongside auto-density and auto-momentum multipoles. The combined density monopole, momentum monopole, and cross dipole power spectrum are then fitted to infer the growth rate \(f\sigma_8\) [2411.09571]. This suggests that “dynamical” power spectra can refer to spectra of fields that directly encode dynamical growth, not only to explicit time-series observables.

## 5. Variable dependence, mechanisms, and physical interpretation

A recurrent theme is that the spectral form is variable-dependent. In thermodynamic nonequilibrium settings, the PSD of an observable \(z\) is constrained by both the observable and the properties of the system, so different observables can have different limiting behaviors and peak structures [2306.00417]. In biological stochastic dynamical systems linearized at a fixed point, the matrix-valued PSD
\[
\mathbfcal{S}(\omega) = (\mathrm{i}\omega \mathbf{I} + \mathbf{J})^{-1}\,\mathbf{L}\,\mathbf{D}\,\mathbf{L}^\top\,(-\mathrm{i}\omega \mathbf{I} + \mathbf{J})^{-\top}
\]
shows explicitly how each auto- and cross-spectrum depends on the Jacobian, dispersion, and diffusion matrices [2305.19890]. Each entry is a complex rational function of frequency whose denominator encodes the eigenvalues of \(\mathbf{J}\), hence damping rates and oscillatory modes [2305.19890].

In the van der Waals channel-flow problem, variable dependence is especially pronounced. The direct cascade and the power-law spectra are primarily a phenomenon of \(\rho\)–\(\chi\) dynamics driven by the van der Waals term and the nonlinear strain term in the divergence equation, while vorticity acts mainly as a stationary background shear [2603.19190]. When vorticity is pinned to its background state and only density and divergence are evolved, the spectra of \(\rho^2\), \(\chi^2\), and the potential velocity components retain the same exponents as in the full system [2603.19190]. This strongly indicates that the underlying physics of the power spectra reside primarily in the density and velocity divergence variables, and are not directly related to the vorticity of the flow [2603.19190].

In semiperiodic pulse trains, the relevant variable split is between pulse shape and forcing statistics. The final PSD always factorizes into a pulse-spectrum envelope and a forcing spectrum. The pulse shape controls the high-frequency envelope, while amplitude and timing statistics control whether the spectrum exhibits a Dirac comb, broadened peaks, or a smooth background [2106.15904]. For Lorentzian pulses, the envelope decays exponentially, so a smooth exponential PSD need not imply a special “chaos-only” mechanism; it may arise from pulse structure together with broad waiting-time statistics and nearly zero-mean amplitudes [2106.15904].

A similar decomposition underlies active-particle spectra. In free AOUPs, the persistence time \(\tau_A\) sets a crossover frequency but does not change the \(f^{-2}\) exponent. Under confinement, the interplay between persistence \(\tau_A\), trap relaxation \(\tau_R\), and activity strength creates multi-regime spectra, including the \(f^{-4}\) range associated with transient ballistic motion [2605.09399]. This suggests that spectral exponents can separate mechanisms even when mean-square displacement alone does not.

## 6. Nonequilibrium, finite-time effects, and interpretation limits

Dynamical power density spectra are particularly sensitive to nonequilibrium structure. For continuous-time Markov processes, the spectrum at arbitrary frequency is bounded by short- and long-time behavior in equilibrium, but out of equilibrium the constants entering the bounds can no longer be identified with the limiting behavior of the spectrum, allowing for peaks that correspond to oscillations in the dynamics [2306.00417]. A key finite-frequency entropy production bound is
\[
\sigma_{\mathrm{st}}
\ge
\frac{2}{\Delta z^2}
\left[
\omega^2 S^z(\omega) -
\lim_{\omega\to\infty}\big(\omega^2 S^z(\omega)\big)
\right],
\]
so the height of peaks in \(\omega^2 S^z(\omega)\) above the high-frequency plateau is related to dissipation [2306.00417].

Finite observation time introduces additional structure. For free AOUPs, the finite-time PSD develops a low-frequency plateau whose leading behavior is
\[
\mu(0,T)\simeq \frac{2}{3}D_{\rm eff}T^2
\]
for large \(T\), together with high-frequency oscillations whose amplitude decays as \(1/(f^3T)\) [2605.09399]. Under confinement, the short-window plateau instead scales as
\[
\mu(0,T)\simeq \langle x^2\rangle\,T,
\]
and the oscillatory corrections decay as \(1/(f^2T)\) with exponential dependence on \(T/\tau_R\) and \(T/\tau_A\) [2605.09399]. This difference shows that finite-time effects can distinguish non-stationary free motion from stationary confined motion.

Single-trajectory analysis imposes further limitations. For Brownian motion, the single-trajectory PSD has the same \(1/f^2\) scaling as the ensemble PSD in the regime of large \(fT\), but the amplitude is a fluctuating random factor:
\[
\tilde S_T^{(k)}(f)=A^{(k)}\,\mu_S(f,T).
\]
Thus the scaling exponent can be inferred from a single trajectory, but the numerical amplitude, and hence the diffusion coefficient, cannot be reliably extracted from a single PSD because the amplitude remains broadly distributed even as \(T\to\infty\) at fixed nonzero frequency [1801.02986]. This shows that a dynamical power density spectrum may encode robust dynamical class information while leaving absolute transport coefficients poorly determined in finite or single-realization data.

A related interpretive caution appears in interferometric astronomy. When the power spectrum is estimated from reconstructed images rather than from visibilities, the image-based estimator acquires a scale-dependent bias due to incompleteness in baseline coverage, whereas the visibility-based estimator reproduces the true power spectrum [1805.11464]. This suggests that the reliability of a dynamical spectrum depends not only on the definition of the spectral object but also on the estimator and observation operator.

## 7. Broader significance and cross-domain usage

The phrase “dynamical power density spectrum” spans several scientific idioms, but the underlying role is consistent: it is a frequency- or wavenumber-resolved measure of fluctuations generated by a dynamical process, used to identify timescales, modes, transport regimes, and underlying mechanisms. In stochastic thermodynamics it serves as a frequency-resolved variance constrained by relaxation and dissipation [2306.00417]. In semiperiodic nonlinear dynamics it distinguishes periodic, quasi-periodic, and chaotic regimes through Dirac combs, broadened harmonics, or pulse envelopes [2106.15904]. In active matter it diagnoses persistence, confinement, and transient ballistic motion through plateaus and power-law crossovers [2605.09399]. In fluid dynamics it reveals inertial-like cascades in nominally laminar flows and can isolate the variables carrying the instability [2603.19190].

The term also generalizes beyond time series. In 21-cm studies of interstellar turbulence, the relevant “dynamical power density spectrum” is realized as spatial power spectra of density and velocity fields in wavenumber space, from which the energy associated with structures and the type of forcing of the turbulence are inferred [1805.11464]. In cosmology, joint power spectra of density and momentum fields in redshift space form a dynamical spectrum of large-scale structure growth, constraining \(f\sigma_8\) through the density monopole, momentum monopole, and density–momentum cross dipole [2411.09571]. In biological stochastic dynamical systems, PSD matrices at fixed points provide an exact rational representation of the second-order statistics of fluctuations, with coefficients determined by the Jacobian, dispersion, and diffusion matrices [2305.19890].

Taken together, these usages suggest a broad definition: a dynamical power density spectrum is a spectral representation of fluctuation power that is constructed from, and interpreted through, the governing dynamics. The spectral variable may be temporal frequency, spatial wavenumber, or a mixed representation. The observable may be a scalar time series, a vector field, or a set of coupled stochastic degrees of freedom. What makes the spectrum dynamical is not the notation \(P(\omega)\) or \(P(k)\) alone, but the explicit connection between spectral structure and the mechanisms that generate, transport, damp, or organize fluctuations.

Source: https://www.emergentmind.com/topics/dynamical-power-density-spectrum