Fact-checked by Grok 2 weeks ago

Jump diffusion

Jump diffusion is a class of stochastic processes that combines continuous paths modeled by diffusion processes, such as Brownian motion, with discontinuous jumps, often driven by a Poisson process.^[1] These models, a type of Lévy process, have applications across fields including physics, pattern theory, and notably financial mathematics, where they capture both gradual changes and sudden shifts in asset prices due to events like economic shocks.^[2] The foundational jump diffusion model in finance was introduced by economist Robert C. Merton in his 1976 paper, extending the Black-Scholes-Merton framework to accommodate non-Gaussian return distributions with heavy tails and excess kurtosis observed in empirical financial data.^[3] In Merton's model, the asset price process follows a geometric Brownian motion augmented by random jumps, where the jump sizes are lognormally distributed and occur at a constant intensity rate λ via a Poisson process.^[2] Mathematically, the stock price S_t evolves as dS_t / S_t = (r - λ κ) dt + σ dW_t + (J - 1) dN_t, with r as the risk-free rate, σ as volatility, W_t as Wiener process, κ as the expected jump size, and N_t as the Poisson counter.^[1] Jump diffusion models have become essential in derivative pricing, risk management, and volatility forecasting, particularly for options where the Black-Scholes assumptions of continuous paths fail to explain phenomena like volatility smiles or implied volatility skews.^[1] Merton's approach yields semi-closed-form solutions for European options via convolution of Black-Scholes prices over Poisson-distributed jump counts, enhancing tractability for short-maturity instruments.^[2] Extensions, such as Kou's 2002 double-exponential jump diffusion, further refine the model by using asymmetric jump sizes to better fit empirical leptokurtosis in indices like the S&P 500.^[1] Despite their strengths in modeling rare events, these models can suffer from estimation challenges and may underperform pure Lévy process alternatives in capturing infinite activity jumps.^[1]

Mathematical Foundations

Definition

Jump diffusion refers to a class of stochastic processes that model systems exhibiting both continuous, gradual variations and abrupt, discontinuous shifts. These processes hybridize a diffusion component, which captures smooth evolution akin to random walks, with a jump component that introduces sudden changes at random times. The diffusion part is typically driven by Brownian motion, representing incremental fluctuations, while jumps occur as discrete events, often governed by a Poisson process for arrival times and random sizes for magnitudes. This combination allows jump diffusion to represent real-world dynamics where change is not purely continuous nor exclusively stepwise.^[4] In contrast to pure diffusion processes, such as the Wiener process, which feature continuous sample paths with no discontinuities, jump diffusion incorporates finite-activity jumps that create kinks or breaks in the trajectory. Pure diffusion models assume all variability arises from infinitely many small, continuous increments, leading to paths that are nowhere differentiable but always continuous. Jump diffusion extends this by adding a finite number of larger, irregular displacements, enabling the modeling of rare, high-impact events alongside routine diffusion. Similarly, it differs from pure jump processes, exemplified by the compound Poisson process, which lack any continuous element and evolve only through discrete leaps at Poisson-distributed times, resulting in piecewise constant paths between jumps. Intuitively, jump diffusion can describe scenarios like the position of a particle undergoing Brownian motion but occasionally experiencing impulsive forces that cause instantaneous shifts, blending pervasive small-scale diffusion with sporadic large displacements. This structure provides a more flexible framework for stochastic modeling than either pure diffusion or pure jump alternatives, accommodating phenomena with mixed continuity and discontinuity without requiring infinite jump activity. Prerequisite to understanding jump diffusion are basic concepts of stochastic processes, which are families of random variables indexed by time, evolving probabilistically to depict uncertainty in dynamic systems.^[5]

Stochastic Differential Equation

The jump diffusion process X_t is governed by the stochastic differential equation (SDE)

dX_t = \mu(X_t) \, dt + \sigma(X_t) \, dW_t + \int_{\mathbb{R}} F(X_{t-}, z) \, \tilde{N}(dt, dz),

where W_t is a standard Brownian motion, \tilde{N}(dt, dz) = N(dt, dz) - \nu(dz) dt is a compensated Poisson random measure with intensity measure \nu(dz) dt, and F(x, z) specifies the jump size as a function of the pre-jump state x and mark z.^[6] This equation decomposes the process into three components: the drift term \mu(X_t) \, dt captures the deterministic trend; the diffusion term \sigma(X_t) \, dW_t models continuous random fluctuations driven by the Wiener process W_t; and the jump term \int_{\mathbb{R}} F(X_{t-}, z) \, \tilde{N}(dt, dz) accounts for discontinuous changes, where the jump intensity is determined by the Lebesgue measure on the marks integrated against \nu(dz), often with \nu(dz) = \lambda f(z) dz for intensity \lambda > 0 and jump size density f(z) (e.g., normal or log-normal).^[6]^[7] In the specific case of finite-activity jumps, the jump component simplifies to dJ_t = \sum_{i=1}^{N_t} Y_i, where N_t is a Poisson process with rate \lambda and Y_i are i.i.d. jump sizes drawn from distribution f(y), leading to the SDE dX_t = \mu(X_t) \, dt + \sigma(X_t) \, dW_t + dJ_t - \lambda \mathbb{E}[Y] dt to ensure the compensator is martingale.^[8] The general framework for jump diffusions arises from the Lévy-Itô decomposition of a Lévy process L_t, which expresses L_t = b t + \Sigma W_t + \int_0^t \int_{|z|>1} z \, N(ds, dz) + \int_0^t \int_{|z|\leq 1} z \, \tilde{N}(ds, dz), where b is the drift, \Sigma the volatility, N the Poisson random measure, and the integrals separate large and small jumps under the condition \int_{|z|\leq 1} |z|^2 \nu(dz) + \int_{|z|>1} 1 \wedge z^2 \nu(dz) < \infty; jump diffusions correspond to the case where the small-jump integral is absent or absorbed into the diffusion, yielding finite variation jumps.^[7] An example is Kou's double exponential jump diffusion model, where for an asset price S_t, the SDE is

\frac{dS_t}{S_{t-}} = \mu \, dt + \sigma \, dW_t + d\left( \sum_{i=1}^{N_t} (V_i - 1) \right),

with \mu, \sigma constants, N_t a Poisson process of rate \lambda, and \log V_i = Y_i following a double exponential distribution with density f_Y(y) = p \eta_1 e^{-\eta_1 y} \mathbf{1}_{y \geq 0} + q \eta_2 e^{\eta_2 y} \mathbf{1}_{y < 0}, where p + q = 1, \eta_1 > 0, \eta_2 > 0.^[9] The infinitesimal generator \mathcal{L} of the jump diffusion process, which governs the evolution of expectations \mathbb{E}[f(X_t) | X_0 = x] for suitable test functions f, is derived using Itô's formula for jump processes: for small h > 0,

\mathbb{E}[f(X_h) - f(x)] = \mathbb{E}\left[ \int_0^h \left( \mu(X_s) f'(X_s) + \frac{1}{2} \sigma^2(X_s) f''(X_s) \right) ds + \sum_{0 < s \leq h} \left( f(X_s) - f(X_{s-}) - f'(X_{s-}) \Delta X_s \right) \right],

and dividing by h then taking h \to 0 yields

\mathcal{L} f(x) = \mu(x) f'(x) + \frac{1}{2} \sigma^2(x) f''(x) + \lambda \int_{\mathbb{R}} \left[ f(x + y) - f(x) - f'(x) y \mathbf{1}_{|y| < 1} \right] f_Y(y) \, dy,

where the integral term arises from the expected jump contribution, with the indicator ensuring compensation for small jumps in the general Lévy case.^[10]^[6]

Properties and Characteristics

Jump diffusion processes exhibit distinct moment properties arising from their combined continuous and discontinuous components. The first moment, or expected value, is E[X_t] = X_0 + \mu t, where \mu is the total drift coefficient in the compensated SDE. The second moment, or variance, is \mathrm{Var}(X_t) = \sigma^2 t + \lambda t \mathbb{E}[Y^2], reflecting the diffusive variance \sigma^2 t plus the jump-induced variability, provided \mathbb{E}[Y^2] < \infty. Higher moments exist under similar integrability conditions on the Lévy measure \nu, specifically if \int_{|y| > 1} |y|^n \nu(dy) < \infty for the n-th moment.^[7]^[11] The presence of jumps imparts a non-Gaussian nature to the process distributions, leading to fat-tailed marginals and potential skewness. Unlike pure diffusions, which yield normal increments, the jump component introduces discontinuities that generate heavier tails—often with kurtosis exceeding 3—and asymmetry depending on the jump size distribution. For instance, log-normal jump sizes in financial models produce positively skewed returns, while symmetric jumps maintain zero skewness but still elevate tail risks. This non-normality is evident in the characteristic function \mathbb{E}[e^{i z X_t}] = e^{t \psi(z)}, where \psi(z) includes the jump integral \int (e^{i z y} - 1 - i z y \mathbf{1}_{|y| \leq 1}) \nu(dy), deviating from the Gaussian quadratic form.^[7]^[11] As Lévy processes, jump diffusions feature stationary and independent increments, meaning the distribution of X_{t+\Delta} - X_t depends only on \Delta and is independent of past increments. Full process stationarity requires initialization from an invariant measure, which exists for specific parameterizations with mean-reverting dynamics, such as the jump Ornstein-Uhlenbeck process dX_t = -\kappa X_t dt + \sigma dW_t + dJ_t under \kappa > 0 and confining jumps. Ergodicity, implying convergence to this stationary distribution, holds exponentially in such cases when the drift dominates jumps and the Lévy measure satisfies moment conditions, ensuring long-term averaging properties.^[7]^[12]^[11] Simulation of these processes combines numerical schemes for the continuous part with exact sampling for jumps. The diffusion component is approximated via the Euler-Maruyama method: \Delta X_t^{\text{diff}} \approx \mu \Delta t + \sigma \sqrt{\Delta t} Z, where Z \sim \mathcal{N}(0,1). Jumps are generated by simulating a Poisson process for arrival times and drawing independent sizes from the distribution governed by \nu, with thinning applied for state-dependent intensities to accept/reject proposals. For infinite activity cases, small jumps are often truncated and approximated by additional Gaussian noise to manage computational complexity.^[13]^[11] A key distinction lies between finite-activity and infinite-activity jump diffusions. Finite-activity models employ a compound Poisson process, where \int \nu(dy) = \lambda < \infty, resulting in finitely many jumps over any finite interval and paths of finite variation if \int_{|y| \leq 1} |y| \nu(dy) < \infty. Infinite-activity variants, driven by general Lévy measures with \nu(\mathbb{R}) = \infty, involve infinitely many small jumps, leading to paths of infinite variation and better capturing clustering or high-frequency discontinuities, as in variance gamma or normal inverse Gaussian processes.^[7]^[11]

Historical Development

Origins in Physics

The conceptual origins of jump diffusion models trace back to studies in statistical mechanics during the 1950s and 1960s, where researchers extended classical Brownian motion to account for impulsive forces arising from discrete collisions or reorientations in molecular systems.^[14] Early work modeled rotational Brownian motion in liquids using random jump mechanisms to describe large-amplitude reorientations, contrasting with the continuous diffusion assumed in Debye's original framework.^[15] These models incorporated time-fluctuating jump rates to capture non-Markovian effects in dense fluids, providing a foundation for handling intermittent, discontinuous changes in particle orientation driven by impulsive interactions.^[14] A key microscopic foundation emerged from derivations starting with the N-body Liouville equation, particularly in the context of stellar dynamics, where collective particle motions under gravitational interactions lead to effective single-particle descriptions using diffusion approximations for relaxation processes in star clusters.^[16] Subrahmanyan Chandrasekhar's analyses in this area demonstrated how the Liouville equation for multi-particle systems reduces to Fokker-Planck equations incorporating diffusion terms for velocity changes due to encounters, highlighting the role of molecular chaos assumptions in projecting the full N-body dynamics onto lower-order distributions.^[16] A more rigorous microscopic analysis was provided in 2009 by Reguera, Rubí, and Pérez-Madrid, who derived the jump diffusion model directly from the N-body Liouville equation using projection operator techniques akin to Zwanzig's formalism.^[17] Their work assumes overdamped conditions and molecular chaos to obtain a generalized Smoluchowski equation featuring both diffusive and jump terms, with the latter arising from the influence of surrounding particles on a tagged particle's motion in molecular systems.^[17] This derivation confirms the jump term's origin in large-amplitude perturbations beyond small-angle approximations, linking it explicitly to the empirical jump diffusion observed in physical systems.^[17] In physical interpretations, jumps in these models represent rare, discontinuous events such as particle captures in gravitational fields or sudden quantum transitions, which dominate the stochastic evolution when continuous diffusion alone fails to capture intermittent dynamics.^[17] For instance, in molecular contexts, jumps correspond to bond-breaking events enabling rapid reorientation, while in astrophysical settings, they model velocity changes from encounters.^[16] These interpretations underscore jump diffusion's utility for rare events that punctuate otherwise smooth particle trajectories in statistical mechanics.^[17]

Introduction in Finance

Jump diffusion processes were first adapted to financial modeling by Robert C. Merton in his seminal 1976 paper, where he extended the Black-Scholes framework by incorporating discontinuous jumps into the underlying asset's return process to better account for sudden market movements.^[18] This innovation addressed key empirical shortcomings of the pure diffusion model introduced by Black and Scholes in 1973, which assumed continuous price paths via geometric Brownian motion and failed to capture the observed leptokurtic distributions and fat tails in asset returns.^[19] Merton's jump-augmented model, combining a diffusion component with Poisson-driven jumps of log-normal size, enabled more accurate option pricing by explaining phenomena like the volatility smile, where implied volatilities vary with strike prices.^[1] The evolution from continuous diffusion models to jump diffusion gained momentum as empirical evidence highlighted the limitations of Gaussian assumptions in replicating real-world return distributions, particularly the higher probability of extreme events.^[20] Key advancements include David S. Bates' 1996 model, which integrated jumps with stochastic volatility to enhance realism in exchange rate and equity pricing, allowing for correlated jump risks under risk-neutral measures. Similarly, Steven G. Kou's 2002 double exponential jump diffusion model introduced asymmetric jump sizes via a Laplace distribution, facilitating closed-form solutions for European options and better fitting skewness and kurtosis in return data.^[9] The 1987 stock market crash significantly amplified the adoption of jump diffusion in quantitative finance, as it exposed the inability of pure diffusion models to price the crash risk implicit in option markets and the resulting volatility skews.^[1] Post-crash analyses demonstrated that jump models could retrospectively and prospectively capture the heightened tail risks, influencing derivative pricing, risk management, and regulatory frameworks in modern finance.^[20]

Advancements in Pattern Theory

Ulf Grenander developed pattern theory in the 1980s and 1990s as a mathematical framework for representing and analyzing complex configurations in high-dimensional spaces, such as images, where jump diffusions played a central role in modeling discrete changes and continuous deformations within generative models.^[21] In this approach, jump diffusions enabled the synthesis of posterior measures through sequential solutions to jump-diffusion equations of generalized Langevin form, facilitating the creation and annihilation of structural elements in pattern configurations.^[21] This integration allowed pattern theory to address variability in representations by incorporating nuisance parameters and hierarchical structures, providing a unified perspective on generative modeling for real-world signals.^[22] During the 1990s, advancements extended jump-diffusion Markov processes to posterior sampling within Bayesian frameworks, enhancing the ability to infer patterns from noisy observations. For instance, Zhu and Mumford introduced stochastic jump-diffusion processes for tasks like computing medial axes, where jumps handled discrete boundary adjustments and diffusions smoothed continuous paths, improving efficiency in Bayesian estimation. Similarly, Zhu's jump-diffusion method for range image segmentation combined birth-death jumps with anisotropic diffusions to sample from posteriors in cluttered data, demonstrating robust performance in separating foreground from background under uncertainty.^[23] These developments solidified jump diffusions as a cornerstone for scalable inference in pattern theory, paralleling contemporaneous extensions in financial modeling for handling discontinuous asset dynamics. A key advancement in the 2000s involved extending jump-diffusion processes to orthogonal groups, such as SO(3), for pose estimation in object recognition, as proposed by Srivastava, Grenander, Jensen, and Miller. This framework constructed ergodic Markov processes on SO(3)^k—where k is the unknown number of objects—using jumps for object instantiation or removal and diffusions for rotational adjustments, ensuring convergence to posterior expectations via sample path averages.^[24] Such extensions built on Grenander's deformable template theory, enabling precise handling of rigid-body transformations in high-dimensional parameter spaces.^[24] In pattern theory, jump diffusions facilitated the separation of representation, observation, and inference phases, particularly in cluttered scenes, by defining distinct probabilistic measures for each: representations via generative templates, observations through likelihoods, and inference via sampling from posteriors.^[22] This modular structure allowed robust analysis of complex systems by isolating variability in templates from observational noise, a principle central to Grenander's unifying perspective on pattern synthesis and recognition.^[25]

Applications in Physics

Particle and Kinetic Systems

Jump diffusion models are employed to simulate the movement of colloidal particles and ions in physical systems where discrete events, such as binding and unbinding to substrates or other particles, interrupt continuous Brownian motion. In crowded environments, these models account for the effective diffusion of tracer particles by incorporating stochastic jumps that represent temporary immobilization during binding events followed by release, leading to anomalous diffusion profiles that deviate from pure Gaussian behavior. For instance, in solutions with high concentrations of macromolecules or ions, the binding/unbinding kinetics reduces the long-time diffusion coefficient while enhancing short-time caging effects, as derived from many-body interaction frameworks. In plasma physics, jump diffusion processes describe particle trajectories in collisionless or weakly collisional regimes, where continuous drift is punctuated by rare, large-angle deflections due to Coulomb interactions or stochastic scattering. These models capture the transport of charged particles across magnetic fields or in turbulent plasmas, approximating the cumulative effect of multiple small collisions as discrete jumps to improve computational efficiency over traditional Fokker-Planck descriptions. Such approaches are particularly useful for simulating test particle diffusion in high-energy plasmas, where the jump term models impulsive velocity changes from long-range electrostatic forces. Numerical simulations of these systems often rely on Monte Carlo methods that hybridize Langevin dynamics for the diffusive component with Poisson processes for jump events, enabling efficient sampling of rare transitions in multi-particle ensembles. In this scheme, particle positions evolve via overdamped Langevin equations between jumps, while jump occurrences and magnitudes are drawn from a Poisson distribution with state-dependent rates, preserving the underlying stochastic differential equation structure. This combination facilitates scalable simulations of kinetic systems, such as colloidal suspensions or plasma sheaths, by avoiding the stiffness associated with pure event-driven algorithms.

Boltzmann Equation Approximations

Recent developments from 2021 to 2023 have introduced particle schemes that leverage jump-diffusion processes to approximate the collision operators in the Boltzmann equation, particularly for modeling rarefied gas dynamics. These schemes, notably the Gamma-Boltzmann model proposed by Mies, Sadr, and Torrilhon, reformulate the Boltzmann equation as a stochastic jump-diffusion process that captures both the free streaming of particles and their binary collisions through discrete jumps.^[26] In this model, the diffusion component represents the continuous free streaming phase, while the jump term explicitly models binary elastic collisions, with parameters tuned to match the relaxation rates of moments up to the heat fluxes, achieving a Prandtl number of 2/3 for Maxwellian molecules.^[26] This hybrid continuous-discrete approach offers significant advantages over traditional Direct Simulation Monte Carlo (DSMC) methods, especially in low-density flows where collision rates are sparse. Unlike DSMC, which relies on fully stochastic particle interactions and scales poorly with decreasing density due to the need for many particles to resolve rare events, jump-diffusion schemes maintain efficiency by decoupling streaming (via deterministic or diffusive advection) from collisions (via targeted jumps), requiring fewer particles per cell—typically around 100—for comparable accuracy.^[26] Computational tests in benchmark problems like Couette flow and lid-driven cavity flows demonstrate that the Gamma-Boltzmann model converges to DSMC reference solutions with reduced variance and better scaling in Knudsen number regimes relevant to microflows and aerospace applications.^[27] Validation of these approximations relies on rigorous mathematical analysis, including error bounds derived from Wasserstein distances between the empirical measures of particle systems and the true Boltzmann solution. Convergence proofs establish that the particle scheme weakly converges to the Gamma-Boltzmann process, with rates depending on the time step and number of particles, ensuring asymptotic fidelity to the original collision operator in the kinetic theory framework.^[26] These bounds confirm the model's suitability for numerical simulations of the Boltzmann equation in transitional flow regimes, where continuum approximations fail.^[27]

Applications in Finance and Economics

Merton's Jump Diffusion Model

Merton's jump diffusion model, introduced in 1976, extends the geometric Brownian motion framework by incorporating discontinuous jumps to better capture the empirical behavior of asset returns, particularly their fat-tailed distributions. The model posits that the stock price S_t follows a stochastic differential equation (SDE) under the physical measure:

\frac{dS_t}{S_{t-}} = (\alpha - \lambda \kappa) \, dt + \sigma \, dW_t + (J - 1) \, dN_t,

where \alpha is the expected return, \sigma is the volatility of the continuous diffusion component driven by a standard Wiener process W_t, N_t is a Poisson process with constant intensity \lambda representing the average number of jumps per unit time, J is the random jump amplitude following a log-normal distribution such that \ln J \sim \mathcal{N}(\gamma, \delta^2), and \kappa = \mathbb{E}[J - 1] is the compensator ensuring the process is a martingale after adjustment.^[18] This formulation allows the model to generate both small continuous changes and occasional large discontinuous shifts, reflecting sudden market events like news announcements or economic shocks. Parameter estimation in Merton's model typically relies on historical return data, employing methods such as maximum likelihood estimation (MLE) or moment matching to infer the values of \alpha, \sigma, \lambda, \gamma, and \delta. MLE maximizes the likelihood of observed returns under the jump-diffusion dynamics, accounting for the compound Poisson jumps, while moment matching equates theoretical moments (e.g., mean, variance, skewness, kurtosis) of the log-return distribution to sample moments from data. These approaches have been applied successfully to equity returns, with MLE providing efficient estimates when jump arrivals are infrequent.^[28]^[29] Empirically, the model improves upon the Black-Scholes framework by explaining the excess kurtosis observed in stock return distributions, which exhibit fatter tails than the normal distribution assumed in pure diffusion models. Studies on common stocks demonstrate that incorporating jumps reduces model mispricing and better replicates the leptokurtic nature of daily or monthly returns, with estimated jump intensities typically higher for individual equities than for indices, often in the range of 10-50 per year depending on the stock and period.^[30] Despite its advantages, Merton's model has notable limitations, including the assumption of a constant jump risk premium embedded in the drift adjustment, which does not account for time-varying compensation for jump risk, and the lack of mechanisms for volatility clustering, where jumps occur independently without correlation to past volatility regimes. These features can lead to underestimation of tail risks during turbulent periods.

Derivative Pricing and Risk Neutral Valuation

In the context of derivative pricing, jump diffusion models, such as Merton's framework, are adapted to the risk-neutral measure to ensure no-arbitrage valuation. Under this measure, the asset price process follows a stochastic differential equation (SDE) where the drift term is adjusted to r - \lambda \kappa, with r denoting the risk-free rate, \lambda the jump intensity, and \kappa the expected relative jump size, while the diffusion and jump components remain unchanged from the physical measure.^[18] This adjustment compensates for the jump risk premium, allowing the discounted asset price to act as a martingale, which is essential for pricing derivatives like options. Pricing solutions for European options in this setting often leverage the characteristic function of the log-price process, which can be solved analytically for jump diffusion models. The closed-form expression for a European call option, known as Merton's formula, represents the price as an infinite sum of Black-Scholes prices, each weighted by the Poisson probability of n jumps occurring over the option's life:

C = \sum_{n=0}^{\infty} \frac{e^{-\lambda' T} (\lambda' T)^n}{n!} C_{BS}(S_0 e^{n \gamma}, K, r_n, \sigma_n^2 T, T),

where \lambda' = \lambda (1 + \kappa), \gamma = \ln(1 + \kappa) - \frac{1}{2} \sigma_J^2, r_n = r - \lambda \kappa + n \gamma / T, \sigma_n^2 = \sigma^2 + n \sigma_J^2 / T, and C_{BS} is the Black-Scholes call price with adjusted parameters.^[18] For more general cases or when closed forms are unavailable, Fourier transform methods, such as the fast Fourier transform (FFT) approach, invert the characteristic function to obtain option prices efficiently. These models extend to more complex derivatives, including barrier options, where jumps introduce the possibility of breaching barriers discontinuously, requiring integral equation solutions or series expansions for pricing.^[9] Jump diffusion also better captures empirical features of implied volatility surfaces, particularly the negative skew observed in equity options, as jumps induce asymmetry in the risk-neutral distribution, leading to higher implied volatilities for low strikes post-jump events. Calibration of model parameters to market option prices further allows extraction of implied jump risk measures, such as \lambda and \kappa, using least-squares optimization or entropy-based methods to fit the volatility smile and term structure.

Applications in Pattern Theory and Imaging

Bayesian Inference Frameworks

In Bayesian inference within pattern theory, jump diffusion processes serve as powerful samplers for exploring complex posterior distributions over configuration spaces, particularly those exhibiting multimodal structures. These processes integrate continuous diffusion dynamics for local exploration with discrete jumps for global transitions, enabling efficient sampling in high-dimensional, non-Euclidean spaces such as Lie groups or manifolds representing scene configurations.^[31]^[21] The mathematical setup typically involves a Langevin diffusion component for continuous refinement, governed by a stochastic differential equation (SDE) that drives the sampler toward regions of high posterior density. For a posterior distribution p(\mathbf{x}) on a configuration space \mathcal{X}, the Langevin diffusion evolves as

d\mathbf{X}_t = \nabla \log p(\mathbf{X}_t) \, dt + \sqrt{2} \, d\mathbf{W}_t,

where \mathbf{W}_t is a Wiener process, providing gradient-based local moves analogous to overdamped Langevin dynamics. This is augmented with Metropolis-Hastings jumps at discrete times, which propose global changes such as altering the dimensionality or topology of the configuration (e.g., adding or removing components in a scene representation). Jump proposals are accepted or rejected based on the posterior ratio, ensuring the overall process has the target posterior as its invariant distribution. On curved spaces like matrix Lie groups, the SDE is adapted to the manifold geometry using Stratonovich integration to preserve the posterior stationarity.^[31]^[32]^[33] In the context of jump-diffusion Markov chain Monte Carlo (MCMC), jumps facilitate global moves, such as abrupt configuration changes (e.g., reconfiguring object arrangements), while the diffusion handles fine-grained local refinements to escape local modes. This hybrid approach contrasts with pure MCMC methods like random-walk Metropolis, which struggle with slow mixing in multimodal posteriors due to high barriers between modes. Jump-diffusion MCMC achieves superior mixing by leveraging jumps to traverse these barriers, leading to more efficient empirical generation of posterior samples and estimators like conditional means. For instance, in spaces with variable model order (e.g., unions of Lie groups G = \bigcup_k G_k), jumps between subgroups G_k and G_{k \pm 1} allow seamless exploration of model uncertainty.^[31]^[32]^[33] A primary advantage over standard MCMC is enhanced ergodicity in multimodal posteriors typical of scene representations, where pure diffusion may trap in suboptimal local maxima, whereas jumps promote rapid traversal of the state space, reducing autocorrelation and improving effective sample size. This is particularly beneficial for Bayesian abduction tasks, where the goal is to infer latent configurations from observations. Empirical studies demonstrate convergence rates that scale favorably with dimensionality compared to dimension-matching MCMC alternatives.^[31]^[32] Grenander's framework in pattern theory formalizes this via jump-diffusion dynamics for empirical posterior generation, treating inference as a stochastic search over generative models of complex structures. Introduced as a synthesis of representations and inference, it uses the process to solve sequential jump-diffusion equations of generalized Langevin form, yielding posterior measures on countable unions of configuration spaces. This approach underpins abduction in pattern theory by enabling the construction of optimal estimators directly from simulated trajectories, without explicit likelihood maximization.^[34]^[32]^[33]

Computer Vision and Scene Understanding

Jump diffusion processes have been applied in computer vision to address challenges in image segmentation and scene understanding, particularly where traditional diffusion methods struggle with multimodal posteriors arising from complex, cluttered environments. These processes combine continuous diffusion for local refinement with discrete jumps to explore global configurations, enabling more robust inference in tasks like object recognition and boundary detection. Rooted in Bayesian frameworks, jump diffusion facilitates posterior sampling by modeling scene hypotheses as probabilistic configurations.^[23] In the 1990s and 2000s, seminal algorithms developed by Zhu and Mumford utilized jump-diffusion for range image segmentation, integrating geometric priors with data-driven likelihoods to delineate surfaces in 3D scans. This approach optimizes a Bayesian posterior by alternating between diffusive flows that smooth within segments and jumps that propose new segmentations, effectively handling occlusions and noise in range data. For instance, in forward-looking infrared (FLIR) scene understanding, jumps generate discrete target hypotheses such as object addition or removal, while diffusion refines continuous parameters like edges and poses, simulating scenes from CAD models to match observed imagery.^[23]^[35] Empirical evaluations demonstrate that jump diffusion enables effective posterior sampling in cluttered scenes, such as military surveillance imagery, where it outperforms pure diffusion methods by escaping local modes in multimodal distributions and achieving higher accuracy in target detection. In FLIR applications, this leads to improved conditional mean estimates for object positions and types, with reported success in resolving ambiguities from overlapping or partially occluded targets.^[35] Modern extensions integrate jump diffusion with deep learning for tasks like pose estimation on manifolds, such as orthogonal groups representing rotations. Classical formulations on Lie groups use jumps to switch between pose configurations and diffusion for fine-tuning, enhancing recognition in varying viewpoints. Recent adaptations, such as in visual tracking, embed jump-diffusion samplers within neural network pipelines to model visibility and motion uncertainties, improving robustness in dynamic scenes over purely data-driven deep methods.^[36]

Medical Imaging and Segmentation

Jump diffusion processes have been explored in medical imaging, drawing analogies to computer vision techniques for segmentation and tracking, with adaptations for physiological noise and constraints in clinical data. In diffusion MRI (dMRI), jump-diffusion models have been used to describe anomalous water dynamics in brain tissue, identifying fast and slow diffusion components that inform interpretations of dMRI signals for brain microstructure analysis.^[37] Foundational work from 1995 introduced conditional mean estimation via jump-diffusion processes for multiple target tracking and recognition, providing a basis for handling dynamic scenes with birth, death, and switching dynamics in noisy environments. These methods have been adapted for estimating moving features in medical sequences, such as in ultrasound or MRI, to improve precision in real-time segmentation tasks like monitoring organ deformation.^[38]^[39] These applications emphasize domain-specific modifications from broader pattern theory and vision frameworks to address tissue discontinuities and imaging artifacts.

References

[1]
[PDF] Jump-Diffusion Models for Asset Pricing in Financial Engineering
In Merton (1976) model, Y has a normal distribution, and in Kou (2002) it has a double exponential distribution. The double exponential distribution enables us ...
[2]
[PDF] Merton's Jump-Diffusion Model
The jump component is composed of lognormal jumps driven by a Poisson process. – It models the rare but large changes in the stock price because of the arrival ...
[3]
https://doi.org/10.1016/0304-405X(76)90022-2
[4]
https://link.springer.com/content/pdf/10.1007/978-981-97-6152-4_11.pdf
[5]
[PDF] An efficient jump-diffusion approximation of the Boltzmann equation
Dec 17, 2021 · A jump-diffusion process along with a particle scheme is devised as an accurate and efficient particle solution to the Boltzmann equation.
[6]
[PDF] Stochastic differential equations with jumps - arXiv
If Zt has a jump of size z, then Xt will have a jump of size a(Xt−)z. However, one might very well want Xt to have a jump whose size depends on Xt− and z, but ...
[7]
[PDF] Jump-diffusion models driven by Lévy processes
The basic principle is to replace the underlying Brownian motion of the Black-Scholes model with a type of jump-diffusion process. In this chapter, the basic ...
[8]
https://www.columbia.edu/~sk75/MagSci02.pdf
[9]
[PDF] A Jump-Diffusion Model for Option Pricing - Columbia University
Like the Black–Scholes model, the double exponential jump-diffusion model not only yields closed-form solutions for standard call and put options (see §5), but ...Missing: SDE | Show results with:SDE
[10]
http://www.dam.brown.edu/people/huiwang/research/kou_wang_jump.pdf
[11]
[PDF] © 2004 by CRC Press LLC - AltExploit
Cont, Rama. Financial modeling with jump processes / Rama Cont, Peter Tankov. p. cm. - (Chapman & HallICRC financial mathematics series).
[12]
[PDF] Exponential ergodicity of the jump-diffusion CIR process - EPFL
Exp. ergodicity results for JCIR processes ( will appear on. Proceeding of the conference at CAS in Oslo (2014), Springer. Verlag 2015).
[13]
https://www.columbia.edu/~mnb2/broadie/Assets/broadie_kaya_exact_sim_or_2006.pdf
[14]
Not Found | AIP Publishing
**Summary of Steele 1963 Paper on Rotational Brownian Motion:**
[15]
Molecular Reorientation in Liquids. I. Distribution Functions and ...
The Debye model of rotational Brownian motion and the rotational random jump model have been extended to allow for time fluctuations of the rotational ...
[16]
Stochastic Problems in Physics and Astronomy | Rev. Mod. Phys.
Stochastic Problems in Physics and Astronomy. S. Chandrasekhar. S. Chandrasekhar Yerkes Observatory, The University of Chicago, Williams Bay, Wisconsin.
[17]
[0907.0417] Microscopic origin of the jump diffusion model - arXiv
Jul 2, 2009 · The present paper is aimed at studying the microscopic origin of the jump diffusion. Starting from the N-body Liouville equation and making ...Missing: stellar | Show results with:stellar
[18]
Option pricing when underlying stock returns are discontinuous
January–March 1976, Pages 125-144. Journal of Financial Economics. Option pricing when underlying stock returns are discontinuous☆. Author links open overlay ...Missing: citation | Show results with:citation
[19]
The Pricing of Options and Corporate Liabilities
Anthony Neuberger The Black–Scholes paper: a personal perspective, Decisions in Economics and Finance 46, no.22 (Oct 2023): 713–730. https://doi.org/10.1007 ...
[20]
Hawkes jump-diffusions and finance: a brief history and review
A brief history of diffusions in Finance is presented, followed by an even briefer discussion of jump-diffusions that involve Poisson or Lévy jumps.
[21]
Representations of Knowledge in Complex Systems - 1994
This measure is synthesized by solving sequentially a jump–diffusion equation of generalized Langevin form. The jumps occur for the creation–annihilation of ...
[22]
[PDF] Pattern Theory: A Unifying Perspective - Applied Mathematics
The term "Pattern Theory" was introduced by Ulf Grenander in the 70s as a name for a field of applied mathematics which gave a theoretical setting.
[23]
[PDF] Range Image Segmentation by an Effective Jump-Diffusion Method
Abstract—This paper presents an effective jump-diffusion method for segmenting a range image and its associated reflectance image in the Bayesian framework.
[24]
https://doi.org/10.1016/S0378-3758(01)00195-1
[25]
Pattern Theory - Hardcover - Ulf Grenander; Michael Miller
Free delivery 25-day returnsPattern Theory provides a comprehensive and accessible overview of the modern challenges in signal, data, and pattern analysis in speech recognition, ...
[26]
An efficient jump-diffusion approximation of the Boltzmann equation
Dec 14, 2021 · The simulation results indicate that the Gamma-Boltzmann model yields a good approximation of the Boltzmann equation, provides a more accurate ...
[27]
An efficient jump-diffusion approximation of the Boltzmann equation
Oct 1, 2023 · A jump-diffusion process along with a particle scheme is devised as an accurate and efficient particle solution to the Boltzmann equation.
[28]
A Note on Estimating the Parameters of the Diffusion-Jump Model of ...
Merton [6] presents an option pricing model for the case where the stock price return dynamics are governed by a mixture of a diffusion and a jump pro? cess.
[29]
[PDF] parameter estimation in merton jump diffusion model
Jul 2, 2019 · This thesis models the US Dollar to Turkish Lira exchange rate using the Merton model, using Maximum Likelihood Estimation (MLE) for parameter ...
[30]
[PDF] Merton Jump-Diffusion Modeling of Stock Price Data
Sep 25, 2018 · Abstract. In this thesis, we investigate two stock price models, the Black-. Scholes (BS) model and the Merton Jump-Diffusion (MJD) model.<|control11|><|separator|>
[31]
None
### Summary of Mathematical Setup of Jump-Diffusion Processes for Bayesian Inference
[32]
[PDF] Representations of Knowledge in Complex Systems Ulf Grenander
Feb 20, 2007 · Ulf Grenander; Michael I. Miller. Journal of the Royal Statistical Society. Series B (Methodological), Vol. 56, No. 4. (1994), pp. 549-603 ...
[33]
[PDF] Pattern Theoretic Bayesian Inference for Multisensor Fusion. - DTIC
Feb 1, 1996 · As we describe below, a jump-diffusion process has the defining properties that it (i) executes jumps in the discrete state variables at.
[34]
Pattern Theory: From representation to inference | Oxford Academic
Pattern Theory provides a comprehensive and accessible overview of the modern challenges in signal, data, and pattern analysis in speech recognition.
[35]
Jump-diffusion processes for the automated understanding of FLIR ...
A jump-diffusion process empirically generates the posterior distribution. The jumps accommodate the discrete aspects of the estimation problem, such as adding ...
[36]
Jump–diffusion Markov processes on orthogonal groups for object ...
The two components of X(t), jumps and diffusions, are chosen in such a way that the resulting Markov process has the desired ergodic property: averages along ...
[37]
Anomalous water dynamics in brain: a combined diffusion magnetic ...
Aug 14, 2019 · The translational model used is the jump-diffusion model and the rotational model used is the continuous rotational diffusion on a circle.
[38]
[PDF] Brain Mapping Methods: Segmentation, Registration, and ...
classification performance, and use the Jump-Diffusion process for generating the model esti- mates in MCMC sampling (Srivastava et al., 2002). We can use ...
[39]
(PDF) Conditional-mean estimation via jump-diffusion processes in ...
Aug 6, 2025 · A new algorithm is presented for generating the conditional mean estimates of functions of target positions, orientations and type in ...
[40]
Jump-diffusion algorithm for multiple target recognition using laser ...
This creates the vanishing point effect in which objects that are further away from the sensor appear closer to the center of the detector. Objects will appear ...