Hawkes Process

The defining mechanism of a Hawkes process is its self-exciting property: the arrival of an event instantaneously increases the conditional intensity, raising the probability of future events in the near term.

• Conditional Intensity Function: The rate of event arrival depends explicitly on the history of past events.
• Branching Structure: Each event can be thought of as an 'immigrant' (spontaneous) or an 'offspring' triggered by a previous event.
• Clustering: This feedback loop naturally generates the volatility clustering observed in real markets without requiring an external regime switch.

The influence of a past event on the current intensity is governed by a kernel function, which dictates how the excitation decays over time.

• Exponential Kernel: The most common choice, offering Markovian properties and computational tractability. Decay follows α * exp(-β * t).
• Power-Law Kernel: Captures long-memory effects where excitation decays more slowly, aligning with empirical evidence of long-range dependence in financial time series.
• Sum of Exponentials: A flexible approximation that allows for multi-timescale excitation, capturing both immediate reflexive reactions and slower diffusion of information.

The branching ratio (n) is the expected number of offspring events generated by a single parent event. It is the key parameter determining system stability.

• Sub-critical (n < 1): The process is stationary. Event clusters eventually die out. This is the standard regime for modeling financial order flow.
• Critical (n = 1): The process is at the edge of stability, with clusters persisting indefinitely on average.
• Super-critical (n > 1): The intensity explodes exponentially. This is non-stationary and typically avoided unless modeling a flash crash or panic event.

The Mutually Exciting Hawkes Process extends the framework to multiple event types, where events of one type can excite or inhibit events of another.

• Cross-Excitation: A buy market order increasing the intensity of subsequent buy orders, capturing herding behavior.
• Bid-Ask Symmetry: Modeling the excitation between orders on opposite sides of the Limit Order Book to capture spread dynamics.
• Excitation Matrix: A matrix of parameters where entry α_ij quantifies the excitation from event type j to event type i, allowing for rich causal modeling of market microstructure.

Parameters are typically estimated using Maximum Likelihood Estimation (MLE). The log-likelihood for a Hawkes process has a closed-form expression that separates the endogenous and exogenous components.

• Compensator: The integral of the intensity function, representing the expected number of events.
• Computational Complexity: Naive MLE is O(N²) due to the double summation over all event pairs, but recursive formulations for exponential kernels reduce this to O(N).
• Regularization: L1 or Lasso penalties are often applied to the excitation matrix in multivariate models to enforce sparsity and prevent overfitting.

Model validation relies on the random time change theorem. If the model is correctly specified, transforming the event times by the integrated conditional intensity yields a unit-rate Poisson process.

• Residual Process: The transformed inter-arrival times should be i.i.d. standard exponential.
• QQ-Plots: Quantile-quantile plots of residuals against the exponential distribution visually diagnose misfit.
• Compensator Plot: The cumulative number of events plotted against the compensator should lie along a 45-degree line if the model is well-calibrated.

The core mechanism where each event arrival increases the conditional intensity of future events. In market microstructure, a large trade triggers a cascade of subsequent activity. The intensity function λ(t) = μ + Σ g(t - t_i) captures this, where μ is the baseline rate and g is the excitation kernel—typically an exponential decay. This creates realistic volatility clustering without requiring external regime-switching mechanisms.

A critical stability parameter representing the expected number of child events generated by a single parent event. Mathematically, it is the integral of the excitation kernel. A ratio < 1 ensures the process is stationary and non-explosive; a ratio approaching 1 creates long memory and critical behavior. In adversarial market simulation, tuning this parameter controls the persistence of order flow cascades and prevents unrealistic infinite feedback loops.

Extends the univariate model to capture cross-excitation between dimensions. The intensity of process i becomes λ_i(t) = μ_i + Σ_j Σ g_{ij}(t - t_k^j). This models how buy orders excite sell orders and vice versa, creating realistic bid-ask bounce. The excitation matrix G encodes the mutually exciting network structure, essential for simulating correlated order book dynamics across multiple price levels or assets.

The most common parametric form for the excitation function: g(t) = α · e^{-βt}. The parameter α controls jump size (immediate intensity increase) while β controls decay rate (how quickly excitation dissipates). This kernel yields Markovian properties, enabling efficient simulation via Ogata's thinning algorithm. Alternatives include power-law kernels for long-memory processes observed in inter-trade durations.

The standard simulation method for Hawkes processes. It uses a rejection sampling approach: propose candidate event times from a homogeneous Poisson process with intensity exceeding the true conditional intensity, then accept or reject based on the ratio. Key steps:

• Maintain an upper bound M(t) ≥ λ(t)
• Generate candidate from Poisson(M(t))
• Accept with probability λ(t)/M(t) This ensures exact simulation without discretization error.

The primary method for calibrating Hawkes parameters to historical order flow data. The log-likelihood for observed event times {t_i} on [0, T] is: L = Σ log λ(t_i) - ∫₀ᵀ λ(s) ds. Optimization is typically performed via gradient descent or EM algorithms. For adversarial simulation, MLE-fitted parameters ensure synthetic markets replicate the statistical signatures of real exchange microstructure.