---
title: Nonlinear Self-Excited Hawkes Processes
url: https://www.emergentmind.com/topics/nonlinear-self-excited-hawkes-processes
type: topic
---

# Nonlinear Self-Excited Hawkes Processes

A nonlinear self-excited Hawkes process is a class of point processes in which the conditional intensity at time \( t \) is a nonlinear, history-dependent functional of previous events. These processes generalize classical linear Hawkes models by allowing nonlinear or state-dependent "link" functions, richer memory kernels, excitation and inhibition, and complex interactions in multivariate and marked extensions. Nonlinear self-excited Hawkes processes capture a broad array of phenomena in fields such as neuroscience, finance, seismology, and network science, where event clustering, feedback, and power-law statistics are prevalent.

## 1. Formal Definition and Existence Theory

Let \( \{N_t\}_{t\ge0} \) denote a simple point process adapted to its natural filtration. The canonical univariate nonlinear Hawkes process is specified by a link (rate) function \( \phi:\mathbb{R}_+\to\mathbb{R}_+ \) and an excitation kernel \( h:[0,\infty)\to[0,\infty) \), yielding the conditional intensity

\[
\lambda_t = \phi\left(\int_0^{t-} h(t-s)\,dN_s\right),
\]
with the process satisfying the compensator condition \( \mathbb{E}[N(a,b]|\mathcal{F}_a] = \mathbb{E}[\int_a^b \lambda_s\,ds|\mathcal{F}_a] \) for all \( a < b \).

For existence and uniqueness, standard assumptions are:
- \( h \) is nonnegative, càdlàg, decreasing, and integrable (\( \|h\|_1 < \infty \));
- \( \phi \) is nonnegative, increasing, and Lipschitz with constant \( a \) such that \( a \|h\|_1 < 1 \).

Under these, the process is stationary, ergodic, and nonexplosive [1204.1067][1304.7531].

Extensions include:
- Multivariate versions, with intensity for each \( k \) given by
  \[
  \lambda^k_t = \phi_k\left(\nu_k + \sum_{l=1}^K \int_{-\infty}^{t^-} h_{lk}(t-s) dN^l_s\right),
  \]
  supporting both excitation (\( h_{lk} > 0 \)) and inhibition (\( h_{lk} < 0 \)) [2103.17164][1707.04928].
- Marked and path-dependent processes, where the conditional intensity is determined by the history of event times and marks, potentially with highly nonlinear structure [2505.22659].

Recent results relax the global Lipschitz requirement on \( \phi \), establishing existence under continuity and subcritical linear growth: \( \Phi(x) \leq C + Lx \) with \( L\|h\|_1 < 1 \) [2212.11660].

Processes with superlinear activations (e.g., polynomial \( \phi(x) = a x^k,\,k>1 \)) are explosive: there is no stationary solution, and events may accumulate in finite time [2212.11660].

## 2. Regimes, Dynamical Properties, and Phase Transitions

Nonlinear self-excited Hawkes processes display rich dynamical regimes determined by growth rates of \( \phi \) and the integrated kernel \( m_1 = \|h\|_1 \):

- **Sublinear:** \( \lim_{z\to\infty} \phi(z)/z = 0 \). The process is stationary and exhibits standard LLN, CLT, and LDP.
- **Subcritical linear:** \( \phi(z) \sim z \) for large \( z \), \( m_1 < 1 \). The process is stationary. Mean event rate is \( \mu = \phi(0)/(1-m_1) \).
- **Critical:** \( \phi(z) \sim z \), \( m_1 = 1 \). Stationarity fails, with quantities such as \( N_t \) growing polynomially.
- **Supercritical:** \( \phi(z) \sim z, m_1 > 1 \) or \( \phi(z)/z \to \infty \). Exponential or explosive growth, no nontrivial stationary law.
- **Explosive:** Sufficiently fast-growing \( \phi \) (e.g., superlinear) cause clustering to coalesce events within finite time [1304.7531][2212.11660].

In multivariate and high-dimensional cases, excitation and inhibition create further complexity but analogous spectral-radius conditions ensure stationarity [1707.04928][2103.17164].

## 3. Limit Theorems: Gaussian Fluctuations and Large Deviations

Under regularity and stability, the stationary nonlinear Hawkes process exhibits principled asymptotic behavior:

- **Functional Central Limit Theorem (FCLT):** For large observation windows \( T \),
  \[
  X^T_t = \frac{N_{Tt}-\mu T t}{\sqrt{T}},\qquad 0\le t\le1
  \]
  converges weakly in Skorokhod topology to a Gaussian process with variance
  \[
  \sigma^2 = \operatorname{Var}(N[0,1]) + 2\sum_{j=1}^\infty \operatorname{Cov}(N[0,1], N[j,j+1]).
  \]
  In the linear case, this reduces to \( \sigma^2 = \nu/(1-\|h\|_1)^3 \) [1204.1067].

- **Law of the Iterated Logarithm (Strassen's Principle):** Sample paths exhibit oscillations bounded by
  \[
  \limsup_{T\to\infty} \frac{N_{Tt}-\mu\,Tt}{\sqrt{2T\log\log T}} = \sigma \sqrt{t}\quad \text{in } C([0,1])
  \]
  [1204.1067].

- **Large Deviation Principles (LDP):** Under broad conditions (sublinear or subcritical regimes), the empirical event rate \( N_t / t \) satisfies an LDP at speed \( t \) with rate function obtained via a variational/Legendre transform of the log moment generating function [1108.2432][1304.7531][1702.05852]. In Markovian cases (kernel as sum of exponentials), explicit variational formulas are available.

- **Moderate Deviations:** Scaling between CLT and LDP, moderate deviations are governed by explicit quadratic rate functions [1702.05852].

- **Non-Markovian Field and Master Equation Approaches:** Infinite-dimensional Markov embeddings yield closed master equations and functional Hamilton–Jacobi PDEs, allowing analysis of multifractality and tail behaviors [2001.01197][2102.00242].

## 4. Power-Law Distributions and Nonlinear Mechanisms

A salient property of broad classes of nonlinear Hawkes processes is the emergence of power-law (heavy-tailed) intensity and event count distributions in stationarity.

If the intensity map \( g(T) \) is fast-accelerating (e.g., \( g(T) = \lambda_0 \exp(\beta T) \) or \( T^n \) with \( n > 2 \)), and marks are two-sided with nonpositive mean and sufficiently fast decreasing tails, the stationary intensity distribution satisfies:
\[
P_{\text{st}}(\lambda) \sim \lambda^{-2} \quad \text{(Zipf's law, mean-zero marks)}.
\]
More generally, \( P_{\text{st}}(\lambda) \sim \lambda^{-1-\mu} \), with the exponent depending on distributional parameters [2102.00242].

This mechanism fundamentally differs from branching-process-driven power laws, arising instead from the nonlinear feedback plus heavy-tailed, sign-ambiguous marks. The same scaling appears in bursty phenomena in seismology, finance, and critical networks.

## 5. Inference: Bayesian, Frequentist, and Nonparametric Approaches

Inference in nonlinear Hawkes models is challenged by the loss of linear structure and the lack of branching representations. Recent methods include:

- **Bayesian Nonparametric Estimation:** Hierarchical priors on baseline rates and kernels (e.g., splines, Gaussian Processes), with MCMC or variational inference. Posterior contraction rates are established under mild entropy and prior-mass conditions. Granger-causality (network) estimation is consistent [2103.17164][2105.09618].
- **Graph Recovery:** Edges in the interaction graph are estimated via spike-and-slab or variable selection priors, with guaranteed consistency and control of false discovery rate [2103.17164].
- **Neural and GP Surrogates:** Feed-forward neural networks (NNNH) or nonparametric GP surrogates model kernels and intensity maps directly from data, enabling flexible capture of both excitation and inhibition [2303.03073][2105.09618].
- **Likelihood Optimization:** For processes with Markovian structure (e.g., exponential kernels), the log-likelihood and gradients have efficient recursions, facilitating scalable MLE or variational methods [2507.22867][2106.04844].
- **State-Augmentation and Thinning Algorithms:** For simulation and inference, thinning (Ogata) schemes are adapted to nonlinear cases, sometimes with path-dependent or latent state augmentations [2505.22659][2106.04844].

## 6. Applications, Empirical Features, and Extensions

Nonlinear self-excited Hawkes processes have been successfully deployed in diverse scientific domains:

- **Neuroscience:** Modeling spike trains with excitation and inhibition, revealing mechanisms of synchrony and variable-length memory corresponding to neuronal refractoriness [2507.22867][1707.04928].
- **Finance:** Interpreting clustering of trades, volatility bursts, and power-law returns via nonlinear mechanism and fat-tailed intensity statistics [2102.00242][2001.01197].
- **Seismology:** ETAS-type models use broad nonlinear Hawkes frameworks to model aftershock distributions and multifractality [2001.01197][2102.00242].
- **Network Science:** Path-dependent nonlinear marked Hawkes models describe time-varying, feedback-driven network growth, including social contact networks with higher-order influence mechanisms [2505.22659].
- **High-Dimensional Data:** Scalable inference (e.g., neural/deep models, sparsity-aware Bayesian schemes) supports applications to tens or hundreds of interacting streams [2303.03073][2103.17164][1707.04928].

Extensive simulation and empirical studies support the empirical validity and flexibility of nonlinear Hawkes processes in these domains, with parameter recovery and causal graph estimation evaluated via held-out likelihood, cross-validation, and hypothesis testing [2507.22867][2303.03073][2103.17164].

## 7. Open Problems and Extensions

Key open directions and recent advances include:
- Relaxation of stability and regularity conditions, notably via the Palm space and Markov-chain perspectives, admitting more general nonlinearities [2212.11660].
- Generalization to variable memory, network-dependent memories, and marked processes [2507.22867][2505.22659].
- Quantitative understanding of explosive regimes and scaling laws near criticality, including transcritical bifurcations in the critical manifold [2212.11660][2001.01197].
- Full characterization of weak dependence, concentration inequalities, and non-asymptotic error in high-dimensional settings [1707.04928].
- Unified field-theoretic treatments allowing extensions to spatially distributed, multifractal, or path-dependent nonlinear settings [2001.01197].

These advances continue to expand the theoretical reach and practical impact of nonlinear self-excited Hawkes processes in modern applied probability and stochastic modeling.

Source: https://www.emergentmind.com/topics/nonlinear-self-excited-hawkes-processes