Power-Law PN-PSD Models Overview

• Power-law PN-PSD models are stochastic processes characterized by a power spectral density that scales as 1/|f|^β over a defined frequency range.
• They leverage frameworks like stochastic differential equations, nonlinear GARCH, and pulse train models to capture burstiness, heavy tails, and long-range correlations.
• These models are essential in fields such as statistical physics, finance, and astronomy for interpreting and synthesizing complex noise phenomena.

Power-law PN-PSD models constitute a class of stochastic processes and statistical models in which the power spectral density (PSD) exhibits power-law scaling over a range of frequencies. These models are central to the analysis and synthesis of nontrivial temporal or spatial fluctuations characterized by long-range correlations, heavy-tailed distributions, and scale-free behavior, and they find application across statistical physics, finance, climatology, and astronomy. In particular, power-law PN-PSD models underpin much of modern analysis of flicker noise ($1/f$ noise), red noise, and systems displaying self-organized criticality or persistent memory.

1. Mathematical Definitions and Fundamental Properties

Let $X(t)$ denote a (zero-mean) stochastic process, either in time (e.g., sampled returns, phase noise, flux) or over a spatial domain. The central defining property of a power-law PN-PSD process is that its power spectral density, $S(f)$, obeys

$S(f) \propto \frac{1}{|f|^\beta}$

over a nontrivial frequency interval $[f_{\min}, f_{\max}]$, where $\beta > 0$ is the spectral index. For white noise, $\beta=0$, while for Brownian (integrated white) noise, $\beta=2$.

Such processes exhibit nontrivial autocorrelation properties ($C(\tau) \sim \tau^{\beta-1}$ for $1<\beta <3$), and their time series often display heavy tails or bursty fluctuations, violating the assumptions of weak memory implicit in classical models.

A key feature is the link between the exponent $X(t)$0, the process variance and stationarity: for $X(t)$1, the integrated power diverges as $X(t)$2, so practical models often introduce explicit low-frequency cutoffs or allow for bending (smoothly broken) power laws that flatten below a break frequency, thereby ensuring physical stationarity (Chakraborty et al., 2020).

2. Stochastic Differential Equations and Scaling Principles

Stochastic differential equations (SDEs) provide a unifying framework for continuous-time power-law PN-PSD models. A general form that yields stationary distributions with power-law tails and a PSD scaling as $X(t)$3 is

$X(t)$4

where $X(t)$5 is the standard Wiener process, $X(t)$6 sets the noise amplitude, $X(t)$7 controls the nonlinearity, and $X(t)$8 determines the PDF tail exponent. The stationary probability density is $X(t)$9 for $S(f)$0. The corresponding spectral index is (Kaulakys et al., 2010, Ruseckas et al., 2014):

$S(f)$1

This relation allows the construction of models covering the full range $S(f)$2 (and beyond for more elaborate SDE forms). For example, for $S(f)$3 one obtains pure $S(f)$4 noise.

Such SDEs naturally generate bursty, intermittent time series and unify disparate phenomena—avalanche statistics, long-range autocorrelations, structural scaling of increments—under a single scaling principle: time rescaling or stretching maps to amplitude rescaling, enforcing scale invariance across the process (Ruseckas et al., 2014).

3. Discrete and Nonlinear GARCH-Based Power-Law Models

Within discrete-time modeling, the Generalized Autoregressive Conditional Heteroskedasticity (GARCH) family offers a parametric route to power-law PN-PSD behavior, particularly relevant to financial volatility series.

• The linear GARCH(1,1) process for variance,

$S(f)$5

produces a heavy-tailed stationary distribution for $S(f)$6, but its PSD scales as $S(f)$7 (Brownian noise), lacking true $S(f)$8 behavior outside the $S(f)$9 case.

• Nonlinear extensions such as “$S(f) \propto \frac{1}{|f|^\beta}$0-GARCH-I” and “$S(f) \propto \frac{1}{|f|^\beta}$1-GARCH-II,” where the variance update is made nonlinear in the previous variance or past shock by introducing powers $S(f) \propto \frac{1}{|f|^\beta}$2, yield, in the continuous-time limit, diffusion processes whose Itô SDEs fall into the nonlinear class detailed above (Kononovicius et al., 2014). The stationary density then has the tail $S(f) \propto \frac{1}{|f|^\beta}$3 and the PSD exponent is

$S(f) \propto \frac{1}{|f|^\beta}$4

For $S(f) \propto \frac{1}{|f|^\beta}$5, pure $S(f) \propto \frac{1}{|f|^\beta}$6 scaling is recovered; larger $S(f) \propto \frac{1}{|f|^\beta}$7 gives $S(f) \propto \frac{1}{|f|^\beta}$8. This shows ARCH-family models can, via nonlinearity, generate both empirical power-law distributions (volatility clustering) and long-memory spectra observed in many financial and other complex systems (Kononovicius et al., 2014).

4. Simulation Algorithms and Practical Construction

Practical synthesis of power-law noise is typically performed in the discrete Fourier domain, exploiting the spectral-shaping filter approach. Starting from white Gaussian noise, a frequency-dependent shaping is applied:

• For a desired PSD exponent $S(f) \propto \frac{1}{|f|^\beta}$9, set the Fourier coefficients as $[f_{\min}, f_{\max}]$0 with $[f_{\min}, f_{\max}]$1, then inverse FFT to the time domain (Ashby, 2011).
• For finite-length processes, circular convolution with appropriately shaped filters (see Wold representation) yields exactly the correct $[f_{\min}, f_{\max}]$2 structure (Kimberk et al., 2022).
• The iterative generation of multicolored (piecewise power-law) noise, fitted to empirically measured Allan deviation profiles, can be done by mapping local time-domain slopes to frequency exponents and assembling a piecewise continuous PSD (Marchi et al., 2023).

From the process statistics, one can derive analytic formulae for the variance, autocovariances, and even the statistics of zero crossings, all in terms of the underlying power-law exponent (Kimberk et al., 2022).

5. Bending Power-Law PSDs and Stationarity Considerations

For processes with steep spectra ($[f_{\min}, f_{\max}]$3), weak stationarity is not preserved, as total power diverges for $[f_{\min}, f_{\max}]$4 (integrated variance explodes). To resolve this, bending (or smoothly broken) power-law PSD models are employed:

$[f_{\min}, f_{\max}]$5

Here, above the break frequency $[f_{\min}, f_{\max}]$6, the spectrum follows power-law scaling ($[f_{\min}, f_{\max}]$7), while below it saturates ($[f_{\min}, f_{\max}]$8 constant), restoring finite total power and statistical stationarity (Chakraborty et al., 2020). This approach is prevalent in modeling long-term variability in AGN lightcurves and similar stochastic astrophysical contexts (Gonzalez-Martin et al., 2012, Diamantopoulos et al., 12 Mar 2026).

The explicit break introduces a timescale (typically linked to physical mechanisms, e.g., viscous disk timescales in AGN), and avoids the convergence pathologies of pure power-law models, while retaining broad flexibility for joint PSD and PDF modeling.

6. Pulse Train, Intermittency, and Other Physical Constructions

An alternative to SDE-based models is the pulse sequence, where the process is expressed as a sum of random pulses (in arrival time, amplitude, and duration). If the pulse duration distribution follows a power law ($[f_{\min}, f_{\max}]$9), the resulting PSD obeys

$\beta > 0$0

asymptotically for large frequencies (Garcia et al., 2016). The statistical properties of the pulse train—variance, skewness, kurtosis—are determined by pulse amplitude statistics and arrival rates but not by $\beta > 0$1 itself (which only sets the spectrum). This shot-noise formalism provides deep connections between phenomenological and physically grounded models; for instance, the mapping between pulse models and SDEs is explicit via time rescaling arguments (Ruseckas et al., 2014).

7. Empirical Applications and Model Diagnostics

Power-law PN-PSD models underpin quantitative analysis of noise in oscillators, clocks, financial time series, X-ray variability in AGN, and more. Experimental diagnostics commonly include:

• Fitting of periodograms/PSD to single or bending power-law models, extracting $\beta > 0$2 or bend parameters (Gonzalez-Martin et al., 2012, Diamantopoulos et al., 12 Mar 2026).
• Computing Allan variance or modified Allan variance, then inferring PSD via analytic inversion or Monte Carlo sampling (Ashby, 2011, Marchi et al., 2023).
• Assessing stationarity by examining the behavior at low frequencies and via tail integrability of variance or autocorrelations (Chakraborty et al., 2020).
• Evaluating intermittency, kurtosis, skewness, and burst statistics to quantify deviation from classical Gaussian or white noise behavior (Garcia et al., 2016).

Power-law PN-PSD models thus provide both a flexible parametric and a process-based theoretical infrastructure for modeling, simulating, and diagnosing nontrivial scale-invariant fluctuations in a wide variety of complex stochastic systems. They bridge microscopic phenomenological mechanisms (nonlinear response, burstiness, aggregation of impacts) and the macroscopic scaling laws observed in empirical data across disciplines (Kononovicius et al., 2014, Kaulakys et al., 2010, Ruseckas et al., 2014).