Mathematical and stochastic analysis of financial markets, derivative pricing and hedging.
Looking for a broader view? This category is part of:
Within a general semimartingale framework, we study the relationship between collective market efficiency and individual rationality. We derive a necessary and sufficient condition for the existence of (possibly zero-sum) exchanges among agents that strictly increase their indirect utilities and characterize this condition in terms of the compatibility between agents' preferences and collective pricing measures. The framework applies to both continuous- and discrete-time models and clarifies when cooperation leads to a strict improvement in each participating agent's indirect utility.
We study a speculative trading problem within the exploratory reinforcement learning (RL) framework of Wang et al. [2020]. The problem is formulated as a sequential optimal stopping problem over entry and exit times under general utility function and price process. We first consider a relaxed version of the problem in which the stopping times are modeled by the jump times of Cox processes driven by bounded, non-randomized intensity controls. Under the exploratory formulation, the agent's randomized control is characterized via the probability measure over the jump intensities, and their objective function is regularized by Shannon's differential entropy. This yields a system of the exploratory HJB equations and Gibbs distributions in closed-form as the optimal policy. Error estimates and convergence of the RL objective to the value function of the original problem are established. Finally, an RL algorithm is designed, and its implementation is showcased in a pairs-trading application.
Options with maturities below one week, hereafter "ultra-short-term" options, have seen a sharp increase in trading activity in recent years. Yet, these instruments are difficult to price jointly using classical pricing models due to the pronounced oscillations observed in the at-the-money implied-volatility term structure across ultra-short-term tenors. We propose Edgeworth++, a parsimonious jump-diffusion model featuring a nonparametric stochastic volatility component, which provides flexibility in capturing implied-volatility smiles for each tenor, combined with a deterministic shift extension, which allows the model to fit rich at-the-money implied-volatility shapes across tenors. We derive a local (in tenor) expansion of the process characteristic function suited to value ultra-short-term options. The expansion leads to fast and accurate option pricing in closed form via standard Fourier inversion. We discuss the benefits of the proposed approach relative to benchmarks.
2603.28948We study scaled trinomial models converging to the Black--Scholes model, and analyze exponential certainty-equivalent prices for path-dependent European options. As the number of trading dates $n$ tends to infinity and the risk aversion is scaled as $nl$ for a fixed constant $l>0$, we derive a nontrivial scaling limit. Our analysis is purely probabilistic. Using a duality argument for the certainty equivalent, together with martingale and weak-convergence techniques, we show that the limiting problem takes the form of a volatility control problem with a specific penalty. For European options with Markovian payoffs, we analyze the optimal control problem and show that the corresponding delta-hedging strategy is asymptotically optimal for the primal problem.
2603.28256We study the upper hedging price for contingent claims in market models with strong types of arbitrage: increasing profit, strong arbitrage, and arbitrage of the first kind. The existence of arbitrage may make the price smaller than if it did not exist. For example, when the asset price process has a reflecting boundary, which introduces increasing profit in the market model, the option prices are reduced to those of the corresponding options that knock-out at the boundary. Furthermore, we demonstrate that corporate stock price processes with increasing profit are obtained as a result of corporate stock issuance and repurchase plans.
Accurately characterizing the implied volatility curves is a central challenge in option pricing and risk management. The classical SABR model by Hagan et al. has been widely adopted in practice due to its well-defined stochastic volatility structure and its tractable closed-form approximation for Black implied volatility. However, under complex market conditions, its fitting accuracy for implied volatility curves remains limited. To address this issue, this paper proposes an extended model within the SABR framework, referred to as skew-SABR. Specifically, the proposed approach introduces an extension to the stochastic dynamics of the underlying asset price and its variance process, under which a corresponding Black implied volatility expression is derived. By further simplifying and reorganizing the resulting formula, the implied volatility can be expressed in a form that explicitly incorporates a skew parameter, thereby enabling a direct characterization of the asymmetry in the implied volatility curve. The resulting expression preserves the structural simplicity of the Hagan-SABR formula, while significantly enhancing the model's flexibility in capturing complex volatility smile patterns. From a theoretical perspective, the paper provides a systematic analysis of the model specification and the financial interpretation of its parameters. From an empirical perspective, a comprehensive comparison is conducted using data from the Chinese options market over the period 2018--2025. The skew-SABR model is evaluated against the classical Hagan-SABR model, the SVI parameterization, polynomial fitting, and spline-based methods. Numerical results show that, across different market regimes and a wide range of implied volatility curve shapes, the skew-SABR model consistently achieves high and stable fitting accuracy.
We develop a semi-static framework for the variance-optimal hedging of multi-asset derivatives exposed to correlation and covariance risk. The approach combines continuous-time dynamic trading in the underlying assets with a static portfolio of auxiliary contingent claims. Using a multivariate Galtchouk--Kunita--Watanabe decomposition, we show that the resulting global mean-variance problem decouples naturally into an inner continuous-time projection onto the space spanned by the underlying assets and an outer finite-dimensional quadratic optimization over the static hedging instruments. To systematically select suitable auxiliary claims, we leverage multidimensional functional spanning theory, establishing that otherwise unhedgeable cross-gamma exposures can be structurally mitigated through static strips of vanilla, product, and spread options. As a central application, we derive explicit semi-static replication formulas for covariance swaps and geometric dispersion trades. Our framework accommodates a broad class of asset dynamics, including quadratic and stochastic Volterra covariance models, as well as affine stochastic covariance models with jumps, yielding tractable semi-closed-form solutions via Fourier transform techniques. Extensive numerical experiments demonstrate that incorporating optimally weighted static strips of cross-asset instruments substantially reduces the mean-squared hedging error relative to purely dynamic benchmark strategies across various model classes.
We develop an averaging approach to robust risk measurement under payoff uncertainty. Instead of taking a worst-case value over an uncertainty neighborhood, we weight nearby payoffs more heavily under a chosen metric and average the baseline risk measure. We prove continuity in the neighborhood radius and provide a stable large-radius behavior. In Banach lattices, the approach leads to a convex risk measure and under separability of the space, a dual representation through a penalty term based on an inf-convolution taken over a Gelfand integral constraint. We also relate our veraging to aggregation at the distribution and quantile levels of payoffs, obtaining dominance and coincidence results. Numerical illustrations are conducted to verify calibration and sensitivity.
2603.22058This thesis develops equilibrium asset pricing models in incomplete markets with a large number of heterogeneous agents using mean field game theory. The market equilibrium is characterized by a novel form of mean field backward stochastic differential equations (BSDEs). First, we propose a theoretical model that endogenously derives the equilibrium risk premium. Agents with exponential preferences are heterogeneous in initial wealth, risk aversion, and unspanned stochastic terminal liability. We solve the optimal investment problem using the optimal martingale principle. The equilibrium is characterized by a mean field BSDE whose driver has quadratic growth in both the stochastic integrands and their conditional expectations. We prove the existence of solutions and show that the risk premium clears the market in the large population limit. Second, we extend the model to include consumption and habit formation, relaxing the time-separability assumption of utility functions. A similar mean field BSDE is derived, and its well-posedness and asymptotic behavior are examined. We also introduce an exponential quadratic Gaussian (EQG) reformulation to obtain equilibrium solutions in semi-analytic form. Finally, the model is extended to partially observable markets where agents must infer the risk premium from stock price observations to determine trading strategies. We provide semi-analytic expressions for the equilibrium via the EQG framework, and the equilibrium risk-premium process is constructed endogenously using Kalman-Bucy filtering theory. Numerical simulations are included to visualize the resulting market dynamics.
We investigate the data-driven discovery of parametric representations for implied volatility slices. Using symbolic regression, we search for simple analytic formulas that approximate the total implied variance as a function of log-moneyness and maturity. Our approach generates candidate parametrizations directly from market data without imposing a predefined functional form. We compare the resulting formulas with the widely used SVI parametrization in terms of accuracy and simplicity. Numerical experiments indicate that symbolic regression can identify compact parametrizations with competitive fitting performance.
Denoising diffusion probabilistic models (DDPMs) have emerged as powerful generative models for complex distributions, yet their use in arbitrage-free derivative pricing remains largely unexplored. Financial asset prices are naturally modeled by stochastic differential equations (SDEs), whose forward and reverse density evolution closely parallels the forward noising and reverse denoising structure of diffusion models. In this paper, we develop a framework for using DDPMs to generate risk-neutral asset price dynamics for derivative valuation. Starting from log-return dynamics under the physical measure, we analyze the associated forward diffusion and derive the reverse-time SDE. We show that the change of measure from the physical to the risk-neutral measure induces an additive shift in the score function, which translates into a closed-form risk-neutral epsilon shift in the DDPM reverse dynamics. This correction enforces the risk-neutral drift while preserving the learned variance and higher-order structure, yielding an explicit bridge between diffusion-based generative modeling and classical risk-neutral SDE-based pricing. We show that the resulting discounted price paths satisfy the martingale condition under the risk-neutral measure. Empirically, the method reproduces the risk-neutral terminal distribution and accurately prices both European and path-dependent derivatives, including arithmetic Asian options, under a GBM benchmark. These results demonstrate that diffusion-based generative models provide a flexible and principled approach to simulation-based derivative pricing.
2603.16108We develop a continuous-time general equilibrium framework for economies with a heterogeneous population -- modeled as a continuum -- that repeatedly optimizes over short horizons under relative-income (Duesenberry-type) criteria. The cross-section evolves through a Brownian flow on a type space, transporting an effective impatience field that captures time variation in preferences induced by demographic changes, aging, and broader social shifts. We establish three main results. First, we prove an optimal consumption--investment theorem for infinite heterogeneous populations in this Brownian-flow setting. Second, we define a \emph{short-horizon Duesenberry equilibrium} -- a sequential-trading (Radner-type) equilibrium in which agents repeatedly solve vanishing-horizon problems under a relative-income criterion -- and provide a complete characterization and existence proof under mild regularity conditions; notably, market completeness and absence of (state-tame) arbitrage emerge endogenously from the market clearing, rather than being imposed as hypotheses. Third, we derive sharp asset-pricing implications: in equilibrium, the market price of risk is pinned down by the volatility of aggregate \emph{total wealth} (financial plus human capital), implying that the equity premium is governed by the magnitudes and correlations of wealth and equity volatilities rather than by consumption volatility alone. This shifts the equity premium puzzle from an implausibly low consumption volatility to a question about the volatility of aggregate total wealth. The framework delivers explicit decompositions of the risk-free rate that are consistent with macro-finance stylized facts. All equilibrium quantities are expressed in terms of market primitives.
Short-horizon option book management relies on P&L expansions in a small set of risk factors. In practice, the quadratic term and common desk adjustments (smile corrections, execution cost add-ons) depend on the chosen factor coordinates, so predicted second-order P&L can change when moving between spot, forward, and log-forward parameterizations. We propose a local, model-agnostic framework that makes the quadratic term coordinate invariant. The usual Hessian is replaced by a covariant Hessian defined by an affine connection, yielding an invariant quadratic predictor. The connection is calibrated to match a desk target for quadratic P&L (Vanna-Volga for smile effects or, in principle, a local fit to realized P&L) while leaving first-order hedge Greeks unchanged. Execution frictions enter through a quadratic cost model for hedge trades. Combined with hedge ratios, this induces an equivalent quadratic penalty on factor moves, makes portfolio netting of costs explicit, and provides local liquidity-aware second-order sensitivities and rebalancing directions. Calibration reduces to small linear systems with clear identifiability conditions. Two FX barrier case studies (EURUSD, USDTRY) illustrate the workflow, and we briefly sketch extensions to other quadratic penalties (risk normalization, scenario/gap terms, and xVA/capital add-ons).
2603.14024Whenever dealing with horizons of different times scales, risk evaluation of losses may incur in both interest rate uncertainty and horizon risk as introduced in [11]. With the goal to capture both effects, we work with cash subadditive fully-dynamic risk measures. In this work we consider such measures obtained via the BSDE and the shortfall approaches. We stress that we consider BSDEs both with Lipschitz and quadratic drivers. We then introduce the hq-entropic risk measure on losses as an effective example of fully-dynamic risk measure serving the scope. Shortfall risk measures are extended to capture cash non-additivity. For our newly introduced h-generalized shortfall risk measures we provide a dual representation and we connect them to fully-dynamic certainty equivalent. To conclude, we can see that the hq-entropic risk measures on losses belong to the family h-generalized shortfall, but they are not of certainty equivalent type. We note that the classical entropic risk measure, besides being generated by a BSDE, is also both a shortfall and a certainty equivalent.
Martingale Optimal Transport (MOT) provides a framework for robust pricing and hedging of illiquid derivatives. Classical MOT enforces exact calibration of model marginals to the mid-prices of vanilla options. Motivated by the industry practice of fitting bid and ask marginals to vanilla prices, we introduce a relaxation of MOT in which model-implied volatilities are only required to lie within observed bid--ask spreads; equivalently, model marginals lie between the bid and ask marginals in convex order. The resulting Bid--Ask MOT (BAMOT) yields realistic price bounds for illiquid derivatives and, via strong duality, can be interpreted as the superhedging price when short and long positions in vanilla options are priced at the bid and ask, respectively. We further establish convergence of BAMOT to classical MOT as bid--ask spreads vanish, and quantify the convergence rate using a novel distance intrinsically linked to bid--ask spreads. Finally, we support our findings with several synthetic and real-data examples.
We investigate the performance of the Kelly rule in a setting in which the dynamics of the return is represented by a time change process. We find that in this general semi-martingale setting the Kelly rule does not maximize the average growth rate, unless the log-return is normally distributed. Namely, the investment position proposed by the Kelly rule is too large, and the investor could achieve a higher average growth rate by investing less aggressively. The higher the variance of the stochastic clock, the more material the failure of the Kelly rule. The ruin threshold proposed by Thorp (1969) is closer, even though examples based on stochastic clock variance estimates taken from the literature show that Kelly rule investment remains safely in the ruin-free region. Finally, the goal of keeping the investment below the ruin threshold for a family of stochastic clock distributions generates a long term investment problem that parallels the "acceptable investment" theory.
We establish a microstructural foundation of the rough Bergomi model. Specifically, we consider a sequence of order driven financial market models where orders to buy or sell an asset arrive according to a Poisson process and have a long lasting impact on volatility. Using a recently established C-tightness result for càdlàg processes we establish the weak convergence of the price-volatility process to a log-normal rough volatility model. Our weak convergence result is accompanied by weak error rates that employ a recently established Clark-Ocone formula for Poisson processes and turn our microstructure model into viable alternative to classical simulation schemes. The weak error rates strongly hinge on Poisson arrival dynamics and are novel to the rough microstructure literature.
We consider the pricing of derivatives written on accumulated marks, such as weather derivatives or aggregate loss claims, using a self-exciting marked point process. The jump intensity mean-reverts between events and increases at jump times by an amount proportional to the mark. The resulting state process, where the variable $U_t$ accumulates jump magnitudes, is a piecewise deterministic Markov process (PDMP). We derive the discounted pricing equation as a backward partial integro-differential equation (PIDE) in two spatial dimensions. To overcome the dimensionality, we propose an exponential (Laplace/Fourier) transform in the accumulated mark variable, which diagonalizes the translation operator and reduces the pricing problem to a family of one-dimensional PIDEs in the intensity variable along a Bromwich contour. For Gamma-mixture mark laws (under actuarial or Esscher-tilted measures), the nonlocal jump term is efficiently approximated by generalized Gauss--Laguerre quadrature. We solve the reduced PIDEs backward in time using a monotone IMEX finite difference scheme (implicit upwind drift and discounting, explicit jump operator) and recover option prices via numerical inversion. We provide a rigorous, term-by-term global error bound covering time and space discretization, quadrature, interpolation, and boundary effects, supported by numerical experiments and Monte Carlo benchmarks.
2603.12422We study the aggregate hazard rate of a heterogeneous population whose individual event intensities are modeled as Cox (doubly stochastic) processes. In the deterministic hazard setting, the observed pool hazard is the survival weighted mean of the individual hazards, and its time derivative equals the mean individual hazard drift minus a variance term. This yields a transparent structural explanation of burnout in mortgage pools. We extend this perspective to stochastic intensity models. The observed pool hazard remains a survival-weighted mean, but now evolves as an Ito process whose drift contains the mean drift of the individual hazards and a negative selection term driven by cross-sectional dispersion, together with a diffusion term inherited from the common factor. We formulate the general identity and discuss special cases relevant to mortgage prepayment modeling.
We find an approximate Nash equilibrium in a game between decentralized exchanges (DEXs) that compete for order flow by setting dynamic trading fees. We characterize the equilibrium via a coupled system of partial differential equations and derive tractable approximate closed-form expressions for the equilibrium fees. Our analysis shows that the two-regime structure found in monopoly models persists under competition: pools alternate between raising fees to deter arbitrage and lowering fees to attract noise trading and increase volatility. Under competition, however, the switching boundary shifts from the oracle price to a weighted average of the oracle and competitors' exchange rates. Our numerical experiments show that, holding total liquidity fixed, an increase in the number of competing DEXs reduces execution slippage for strategic liquidity takers and lowers fee revenue per DEX. Finally, the effect on noise traders' slippage depends on market activity: they are worse off in low-activity markets but better off in high-activity ones.