BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Andrey Itkin

Andrey Itkin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
q-fin.CPq-fin.PRq-fin.MFq-fin.GNmath.APphysics.chem-phphysics.soc-phq-fin.RM

Catalog footprint

What is connected

20works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2026arXiv

American options valuation in time-dependent jump-diffusion models via integral equations and characteristic functions

Despite significant advancements in machine learning for derivative pricing, the efficient and accurate valuation of American options remains a persistent challenge due to complex exercise boundaries, near-expiry behavior, and intricate contractual features. This paper extends a semi-analytical approach for pricing American options in time-inhomogeneous models, including pure diffusions, jump-diffusions, and Levy processes. Building on prior work, we derive and solve Volterra integral equations of the second kind to determine the exercise boundary explicitly, offering a computationally superior alternative to traditional finite-difference and Monte Carlo methods. We address key open problems: (1) extending the decomposition method, i.e. splitting the American option price into its European counterpart and an early exercise premium, to general jump-diffusion and Levy models; (2) handling cases where closed-form transition densities are unavailable by leveraging characteristic functions via, e.g., the COS method; and (3) generalizing the framework to multidimensional diffusions. Numerical examples demonstrate the method's efficiency and robustness. Our results underscore the advantages of the integral equation approach for large-scale industrial applications, while resolving some limitations of existing techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Semi-analytical pricing of American options with hybrid dividends via integral equations and the GIT method

This paper introduces a semi-analytical method for pricing American options on assets (stocks, ETFs) that pay discrete and/or continuous dividends. The problem is notoriously complex because discrete dividends create abrupt price drops and affect the optimal exercise timing, making traditional continuous-dividend models unsuitable. Our approach utilizes the Generalized Integral Transform (GIT) method introduced by the author and his co-authors in a number of papers, which transforms the pricing problem from a complex partial differential equation with a free boundary into an integral Volterra equation of the second or first kind. In this paper we illustrate this approach by considering a popular GBM model that accounts for discrete cash and proportional dividends using Dirac delta functions. By reframing the problem as an integral equation, we can sequentially solve for the option price and the early exercise boundary, effectively handling the discontinuities caused by the dividends. Our methodology provides a powerful alternative to standard numerical techniques like binomial trees or finite difference methods, which can struggle with the jump conditions of discrete dividends by losing accuracy or performance. Several examples demonstrate that the GIT method is highly accurate and computationally efficient, bypassing the need for extensive computational grids or complex backward induction steps.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Multilayer heat equations and their solutions via oscillating integral transforms

By expanding the Dirac delta function in terms of the eigenfunctions of the corresponding Sturm-Liouville problem, we construct some new (oscillating) integral transforms. These transforms are then used to solve various finance, physics, and mathematics problems, which could be characterized by the existence of a multilayer spatial structure and moving (time-dependent) boundaries (internal interfaces) between the layers. Thus, constructed solutions are semi-analytical and extend the authors' previous work (Itkin, Lipton, Muravey, Multilayer heat equations: application to finance, FMF, 1, 2021). However, our new method doesn't duplicate the previous one but provides alternative representations of the solution which have different properties and serve other purposes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Semi-closed form prices of barrier options in the Hull-White model

In this paper we derive semi-closed form prices of barrier (perhaps, time-dependent) options for the Hull-White model, ie., where the underlying follows a time-dependent OU process with a mean-reverting drift. Our approach is similar to that in (Carr and Itkin, 2020) where the method of generalized integral transform is applied to pricing barrier options in the time-dependent OU model, but extends it to an infinite domain (which is an unsolved problem yet). Alternatively, we use the method of heat potentials for solving the same problems. By semi-closed solution we mean that first, we need to solve numerically a linear Volterra equation of the first kind, and then the option price is represented as a one-dimensional integral. Our analysis shows that computationally our method is more efficient than the backward and even forward finite difference methods (if one uses them to solve those problems), while providing better accuracy and stability.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Semi-closed form prices of barrier options in the time-dependent CEV and CIR models

We continue a series of papers where prices of the barrier options written on the underlying, which dynamics follows some one factor stochastic model with time-dependent coefficients and the barrier, are obtained in semi-closed form, see (Carr and Itkin, 2020, Itkin and Muravey, 2020). This paper extends this methodology to the CIR model for zero-coupon bonds, and to the CEV model for stocks which are used as the corresponding underlying for the barrier options. We describe two approaches. One is generalization of the method of heat potentials for the heat equation to the Bessel process, so we call it the method of Bessel potentials. We also propose a general scheme how to construct the potential method for any linear differential operator with time-independent coefficients. The second one is the method of generalized integral transform, which is also extended to the Bessel process. In all cases, a semi-closed solution means that first, we need to solve numerically a linear Volterra equation of the second kind, and then the option price is represented as a one-dimensional integral. We demonstrate that computationally our method is more efficient than both the backward and forward finite difference methods while providing better accuracy and stability. Also, it is shown that both method don't duplicate but rather compliment each other, as one provides very accurate results at small maturities, and the other one - at high maturities.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Semi-closed form solutions for barrier and American options written on a time-dependent Ornstein Uhlenbeck process

In this paper we develop a semi-closed form solutions for the barrier (perhaps, time-dependent) and American options written on the underlying stock which follows a time-dependent OU process with a log-normal drift. This model is equivalent to the familiar Hull-White model in FI, or a time dependent OU model in FX. Semi-closed form means that given the time-dependent interest rate, continuous dividend and volatility functions, one need to solve numerically a linear (for the barrier option) or nonlinear (for the American option) Fredholm equation of the first kind. After that the option prices in all cases are presented as one-dimensional integrals of combination of the above solutions and Jacobi theta functions. We also demonstrate that computationally our method is more efficient than the backward finite difference method used for solving these problems, and can also be as efficient as the forward finite difference solver while providing better accuracy and stability.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Filling the gaps smoothly

The calibration of a local volatility models to a given set of option prices is a classical problem of mathematical finance. It was considered in multiple papers where various solutions were proposed. In this paper an extension of the approach proposed in LiptonSepp2011 is developed by i) replacing a piecewise constant local variance construction with a piecewise linear one, and ii) allowing non-zero interest rates and dividend yields. Our approach remains analytically tractable; it combines the Laplace transform in time with an analytical solution of the resulting spatial equations in terms of Kummer's degenerate hypergeometric functions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

LSV models with stochastic interest rates and correlated jumps

Pricing and hedging exotic options using local stochastic volatility models drew a serious attention within the last decade, and nowadays became almost a standard approach to this problem. In this paper we show how this framework could be extended by adding to the model stochastic interest rates and correlated jumps in all three components. We also propose a new fully implicit modification of the popular Hundsdorfer and Verwer and Modified Craig-Sneyd finite-difference schemes which provides second order approximation in space and time, is unconditionally stable and preserves positivity of the solution, while still has a linear complexity in the number of grid nodes.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Nonlinear PDEs risen when solving some optimization problems in finance, and their solutions

We consider a specific type of nonlinear partial differential equations (PDE) that appear in mathematical finance as the result of solving some optimization problems. We review some existing in the literature examples of such problems, and discuss the properties of these PDEs. We also demonstrate how to solve them numerically in a general case, and analytically in some particular case.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Structural default model with mutual obligations

This paper considers mutual obligations in the interconnected bank system and analyzes their influence on joint and marginal survival probabilities as well as CDS and FTD prices for the individual banks. To make the role of mutual obligations more transparent, a simple structural default model with banks' assets driven by correlated multidimensional Brownian motion with drift is considered. This model enables a closed form representation for many quantities of interest, at least in a 2D case, to be obtained, and moreover, model calibration is provided. Finally, we demonstrate that mutual obligations have to be taken into account in order to get correct values for model parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Efficient Solution of Backward Jump-Diffusion PIDEs with Splitting and Matrix Exponentials

We propose a new, unified approach to solving jump-diffusion partial integro-differential equations (PIDEs) that often appear in mathematical finance. Our method consists of the following steps. First, a second-order operator splitting on financial processes (diffusion and jumps) is applied to these PIDEs. To solve the diffusion equation, we use standard finite-difference methods, which for multi-dimensional problems could also include splitting on various dimensions. For the jump part, we transform the jump integral into a pseudo-differential operator. Then for various jump models we show how to construct an appropriate first and second order approximation on a grid which supersets the grid that we used for the diffusion part. These approximations make the scheme to be unconditionally stable in time and preserve positivity of the solution which is computed either via a matrix exponential, or via P{á}de approximation of the matrix exponent. Various numerical experiments are provided to justify these results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Efficient solution of structural default models with correlated jumps and mutual obligations

The structural default model of Lipton and Sepp, 2009 is generalized for a set of banks with mutual interbank liabilities whose assets are driven by correlated Levy processes with idiosyncratic and common components. The multi-dimensional problem is made tractable via a novel computational method, which generalizes the one-dimensional fractional partial differential equation method of Itkin, 2014 to the two- and three-dimensional cases. This method is unconditionally stable and of the second order of approximation in space and time; in addition, for many popular Levy models it has linear complexity in each dimension. Marginal and joint survival probabilities for two and three banks with mutual liabilities are computed. The effects of mutual liabilities are discussed, and numerical examples are given to illustrate these effects.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High-Order Splitting Methods for Forward PDEs and PIDEs

This paper is dedicated to the construction of high-order (in both space and time) finite-difference schemes for both forward and backward PDEs and PIDEs, such that option prices obtained by solving both the forward and backward equations are consistent. This approach is partly inspired by Andreasen & Huge, 2011 who reported a pair of consistent finite-difference schemes of first-order approximation in time for an uncorrelated local stochastic volatility model. We extend their approach by constructing schemes that are second-order in both space and time and that apply to models with jumps and discrete dividends. Taking correlation into account in our approach is also not an issue.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Splitting and Matrix Exponential approach for jump-diffusion models with Inverse Normal Gaussian, Hyperbolic and Meixner jumps

This paper is a further extension of the method proposed in Itkin, 2014 as applied to another set of jump-diffusion models: Inverse Normal Gaussian, Hyperbolic and Meixner. To solve the corresponding PIDEs we accomplish few steps. First, a second-order operator splitting on financial processes (diffusion and jumps) is applied to these PIDEs. To solve the diffusion equation, we use standard finite-difference methods. For the jump part, we transform the jump integral into a pseudo-differential operator and construct its second order approximation on a grid which supersets the grid that we used for the diffusion part. The proposed schemes are unconditionally stable in time and preserve positivity of the solution which is computed either via a matrix exponential, or via P'ade approximation of the matrix exponent. Various numerical experiments are provided to justify these results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

To sigmoid-based functional description of the volatility smile

We propose a new static parameterization of the implied volatility surface which is constructed by using polynomials of sigmoid functions combined with some other terms. This parameterization is flexible enough to fit market implied volatilities which demonstrate smile or skew. An arbitrage-free calibration algorithm is considered that constructs the implied volatility surface as a grid in the strike-expiration space and guarantees a lack of arbitrage at every node of this grid. We also demonstrate how to construct an arbitrage-free interpolation and extrapolation in time, as well as build a local volatility and implied pdf surfaces. Asymptotic behavior of this parameterization is discussed, as well as results on stability of the calibrated parameters are presented. Numerical examples show robustness of the proposed approach in building all these surfaces as well as demonstrate a better quality of the fit as compared with some known models.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

USLV: Unspanned Stochastic Local Volatility Model

We propose a new framework for modeling stochastic local volatility, with potential applications to modeling derivatives on interest rates, commodities, credit, equity, FX etc., as well as hybrid derivatives. Our model extends the linearity-generating unspanned volatility term structure model by Carr et al. (2011) by adding a local volatility layer to it. We outline efficient numerical schemes for pricing derivatives in this framework for a particular four-factor specification (two "curve" factors plus two "volatility" factors). We show that the dynamics of such a system can be approximated by a Markov chain on a two-dimensional space (Z_t,Y_t), where coordinates Z_t and Y_t are given by direct (Kroneker) products of values of pairs of curve and volatility factors, respectively. The resulting Markov chain dynamics on such partly "folded" state space enables fast pricing by the standard backward induction. Using a nonparametric specification of the Markov chain generator, one can accurately match arbitrary sets of vanilla option quotes with different strikes and maturities. Furthermore, we consider an alternative formulation of the model in terms of an implied time change process. The latter is specified nonparametrically, again enabling accurate calibration to arbitrary sets of vanilla option quotes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

New solvable stochastic volatility models for pricing volatility derivatives

Classical solvable stochastic volatility models (SVM) use a CEV process for instantaneous variance where the CEV parameter $γ$ takes just few values: 0 - the Ornstein-Uhlenbeck process, 1/2 - the Heston (or square root) process, 1- GARCH, and 3/2 - the 3/2 model. Some other models were discovered in \cite{Labordere2009} by making connection between stochastic volatility and solvable diffusion processes in quantum mechanics. In particular, he used to build a bridge between solvable (super)potentials (the Natanzon (super)potentials, which allow reduction of a Schrödinger equation to a Gauss confluent hypergeometric equation) and existing SVM. In this paper we discuss another approach to extend the class of solvable SVM in terms of hypergeometric functions. Thus obtained new models could be useful for pricing volatility derivatives (variance and volatility swaps, moment swaps).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Pricing options on illiquid assets with liquid proxies using utility indifference and dynamic-static hedging

This work addresses the problem of optimal pricing and hedging of a European option on an illiquid asset Z using two proxies: a liquid asset S and a liquid European option on another liquid asset Y. We assume that the S-hedge is dynamic while the Y-hedge is static. Using the indifference pricing approach we derive a HJB equation for the value function, and solve it analytically (in quadratures) using an asymptotic expansion around the limit of the perfect correlation between assets Y and Z. While in this paper we apply our framework to an incomplete market version of the credit-equity Merton's model, the same approach can be used for other asset classes (equity, commodity, FX, etc.), e.g. for pricing and hedging options with illiquid strikes or illiquid exotic options.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Pricing options with VG model using FFT

We discuss various analytic and numerical methods that have been used to get option prices within a framework of the VG model. We show that some popular methods, for instance, Carr-Madan's FFT method could blow up for certain values of the model parameters even for an European vanilla option. Alternative methods - one originally proposed by Lewis, and Black-Scholes-wise method are considered that seem to work fine for any value of the VG parameters. Test examples are given to demonstrate efficiency of these methods. Convergency of all methods is also discussed.}

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Using pseudo-parabolic and fractional equations for option pricing in jump diffusion models

In mathematical finance a popular approach for pricing options under some Levy model is to consider underlying that follows a Poisson jump diffusion process. As it is well known this results in a partial integro-differential equation (PIDE) that usually does not allow an analytical solution while numerical solution brings some problems. In this paper we elaborate a new approach on how to transform the PIDE to some class of so-called pseudo-parabolic equations which are known in mathematics but are relatively new for mathematical finance. As an example we discuss several jump-diffusion models which Levy measure allows such a transformation.

0 citations0 reviewsNo code linkedOpen access