---
title: Time-Fractional Black–Scholes Equation
url: https://www.emergentmind.com/topics/time-fractional-black-scholes-equation-tfbse
type: topic
---

# Time-Fractional Black–Scholes Equation

Searching arXiv for recent and foundational papers on the Time-Fractional Black–Scholes Equation and related variants.
The **Time-Fractional Black–Scholes Equation (TFBSE)** is a fractional-in-time generalization of the Black–Scholes option-pricing PDE in which the first-order time derivative is replaced by a nonlocal fractional derivative, most commonly of Caputo type or an equivalent modified Riemann–Liouville form under suitable regularity. In the TFBSE, the spatial Black–Scholes operator is typically retained while temporal evolution acquires memory through a convolution kernel, so current option values depend on the full past history of the solution rather than only its instantaneous state. Across the literature, the TFBSE appears in several mathematically distinct but related forms: as an ad hoc time-fractional extension of the classical pricing PDE, as the pricing equation associated with subdiffusive or tempered-subdiffusive asset dynamics, and as the target of a wide range of numerical schemes including finite differences, compact methods, spline collocation, meshless RBF methods, neural solvers, and transform-based series representations [1907.00297] [1604.05178] [2207.09153] [2109.02028] [2303.10592] [2602.00201]. A persistent point of clarification is that not every “fractional Black–Scholes” model is a TFBSE: some papers introduce fractionality through stochastic volatility, long memory, or roughness in the volatility factor rather than through a time-fractional pricing operator [1509.01175].

## 1. Definition and model variants

In its standard one-asset form, the TFBSE replaces the classical Black–Scholes time derivative by a fractional derivative of order in \((0,1)\). A representative formulation is
\[
{}_{0}^{c}D_t^\alpha v(z,t) = -\frac{1}{2}\sigma^2 z^2 \frac{\partial^2 v(z,t)}{\partial z^2} -rz\frac{\partial v(z,t)}{\partial z} +r v(z,t),
\]
with European call terminal and boundary conditions [1907.00297]. After time reversal and log-price transformation, this becomes a constant-coefficient convection–diffusion–reaction equation,
\[
{}_{0}^{c}D_t^\alpha u(x,t) = \frac{1}{2}\sigma^2 u_{xx}(x,t)+\left(r-\frac{1}{2}\sigma^2\right)u_x(x,t)-ru(x,t),
\]
which is the form used by many numerical papers [1907.00297] [2109.02028] [2207.09153] [2602.00201].

A closely related form appears when the original problem is posed backward from maturity using a modified right Riemann–Liouville derivative and then transformed into a forward problem. In such settings, the transformed derivative is treated as equivalent to the Caputo derivative under \(C^{(1)}\) temporal regularity [2508.06780] [2602.00201] [2207.09153]. This equivalence is important because Caputo derivatives accommodate classical payoff-based initial conditions.

Several extensions of the TFBSE are represented in the literature. A two-asset time-fractional Black–Scholes PDE with Caputo derivative and mixed derivative term has been written as
\[
\prescript{C}{}{D}_t^\alpha c(S_1,S_2,t) +\frac{1}{2}\sigma_1^2 S_1^2 c_{S_1S_1} +\frac{1}{2}\sigma_2^2 S_2^2 c_{S_2S_2} +\rho \sigma_1 \sigma_2 S_1S_2 c_{S_1S_2} +rS_1 c_{S_1} +rS_2 c_{S_2} -r c=0,
\]
with \(0<\alpha\le 1\) [2010.13411]. Tempered variants replace the standard Caputo derivative by a tempered operator, typically to weaken the heaviest memory tails while retaining nonlocality [2303.10592] [2103.13679]. Time-space-fractional Black–Scholes-type models generalize further by using Caputo time derivatives together with radial Riemann–Liouville spatial derivatives, so the pure TFBSE is recovered only as a specialization [2011.07710].

Not all “fractional Black–Scholes” papers belong to this class. In particular, a model in which volatility is driven by a fractional Ornstein–Uhlenbeck process built from fractional Brownian motion yields fractional-power implied-volatility corrections to classical Black–Scholes, but it does not replace the time derivative in the pricing PDE by a fractional operator [1509.01175]. That distinction is foundational for the taxonomy of the field.

## 2. Fractional derivatives and probabilistic interpretation

The dominant operator in the TFBSE literature is the Caputo derivative. In the one-dimensional setting, it is written as
\[
{}_{0}^{c}D_t^\alpha g(t) = \frac{1}{\Gamma(1-\alpha)} \int_0^t g'(s)(t-s)^{-\alpha}\,ds, \qquad 0<\alpha<1,
\]
or in equivalent modified Riemann–Liouville form after appropriate time reversal [1907.00297] [2508.06780] [2207.09153]. The Caputo choice is standard because it preserves the interpretation of initial conditions in terms of option payoffs.

A more structural interpretation arises in subdiffusive and tempered-subdiffusive Black–Scholes models. In the subdiffusive case, the asset process is obtained by replacing the physical clock of geometric Brownian motion with an inverse \(\alpha\)-stable subordinator,
\[
Z_\alpha(t)=Z(S_\alpha(t)),
\]
where \(S_\alpha\) is independent of Brownian motion [1907.00297]. The resulting option price can be expressed as
\[
v(t)=\mathbb{E}[h(S_\alpha(t))],
\]
where \(h\) is the classical Black–Scholes price. Laplace-transform arguments then yield the Caputo TFBSE as the governing pricing equation [1907.00297]. In this reading, the fractional time derivative is not merely phenomenological; it is the PDE manifestation of random waiting times and market stagnation encoded by inverse-subordinator time change.

A more general construction uses the inverse of a strictly increasing subordinator with Laplace exponent
\[
\phi(\lambda)=\kappa\lambda+\int_0^\infty (1-e^{-\lambda x})\,\nu(dx),
\]
leading to a generalized time-fractional operator
\[
\partial_t^w f(t) := \frac{d}{dt}\int_0^t w(s)\big(f(t-s)-f(0)\big)\,ds,\qquad w(x)=\nu([x,\infty))
\]
[2511.10371]. In this framework, the option-pricing kernel solves
\[
(\kappa \partial_t+\partial_t^w)u(t,x)=\frac12 \sigma^2 x^2 u_{xx}(t,x),\qquad u(0,x)=(x-K)^+,
\]
and the TFBSE is derived from a rigorously specified arbitrage-free but generally incomplete sub-diffusive market model rather than postulated directly [2511.10371]. This suggests that the standard Caputo TFBSE is a special case within a broader inverse-subordinator theory.

Tempered subdiffusive models modify the inverse stable clock by using an inverse tempered stable subordinator, leading to the tempered operator
\[
\partial^{\alpha,\lambda}_{t} u(z,t)=e^{-\lambda t} D_t^\alpha\!\left(e^{\lambda t}(u(z,t)-u(z,0))\right)-\lambda^\alpha (u(z,t)-u(z,0)),
\]
which reduces to the ordinary TFBSE operator when \(\lambda=0\) [2103.13679]. Here tempering weakens the most extreme heavy-tail effects while preserving temporal nonlocality.

## 3. Transformations and canonical computational forms

The standard computational treatment of the TFBSE begins with time reversal and logarithmic transformation. Starting from a terminal-value problem in \((S,\tau)\), one sets
\[
t=T-\tau,\qquad x=\ln S,\qquad u(x,t)=v(e^x,T-t),
\]
or equivalent notation, thereby converting the problem into an initial-value equation on a log-price domain [1907.00297] [2508.06780] [2207.09153] [2109.02028] [2602.00201]. This removes the variable coefficients \(S^2\) and \(S\) from the Black–Scholes operator and yields a constant-coefficient PDE in \(x\).

For bounded-domain computations, the infinite log-price interval is truncated to \([x_{\min},x_{\max}]\) or \((B_x,B_y)\), and Dirichlet conditions are imposed based on asymptotic option behavior [1907.00297] [2508.06780] [2207.09153]. In some approaches, nonhomogeneous boundary values are eliminated by an affine lifting. For example, one may define
\[
z(x,\tau)=\frac{\Psi(\tau)-\Phi(\tau)}{x_r-x_l}(x-x_l)+\Phi(\tau),
\]
followed by a transformed unknown \(v\) satisfying homogeneous boundary conditions [2303.10592] [2109.02028]. In the tempered case, an additional exponential integrating factor is used:
\[
v(x,\tau):=k(x)[U(x,\tau)-z(x,\tau)],\qquad k(x)=\exp\!\left(\frac{c(x-x_l)}{\sigma^2}\right),
\]
which transforms the Black–Scholes operator into a diffusion–reaction equation with homogeneous boundaries [2303.10592].

An older high-order compact treatment uses a different route for a European put. After time reversal and boundary homogenization,
\[
W(s,t):=V(s,t)+\frac{K}{S}(s-S),
\]
and then
\[
U(s,t)=s^q W(s,t),\qquad q=\frac{r-d}{\sigma^2},
\]
to eliminate the first derivative term and obtain a diffusion-form TFBSE involving \(s^2U_{ss}\) and a reaction term [1604.05178]. This diffusion form is particularly convenient for nonuniform spatial meshes tailored to the degenerate left boundary and strike-region nonsmoothness.

These transformations are not mere algebraic conveniences. They determine the structure of the discrete operator, the appropriate boundary treatment, and the attainable order of accuracy, especially when combined with graded temporal meshes designed to handle startup singularities [2109.02028] [2303.10592].

## 4. Analytical solutions, transform methods, and formal series

Closed-form solutions to the TFBSE are generally unavailable except through subordinated semigroup representations, Laplace-domain formulas, or formal series. In the subdiffusive setting, the option price can be written as the classical Black–Scholes price averaged over the law of the inverse subordinator:
\[
v(t)=\mathbb{E}[h(S_\alpha(t))]=\int_0^\infty h(z)\,\varrho_\alpha(z,t)\,dz,
\]
which directly links the TFBSE to the classical model [1907.00297]. The more general sub-diffusive formulation yields
\[
u(t,x)=\mathbb E_x[(Y_{L_t}-K)^+],
\]
where \(Y\) is geometric Brownian motion in operational time and \(L_t\) is the inverse subordinator, and this function is proved to be the unique solution of the generalized TFBSE [2511.10371].

For a European call in the sub-diffusive model, the value at time zero can be written as
\[
V_0 = e^{-rT} \int_0^\infty \mathbb E\!\left[ \left( \mathcal S_0 e^{rT-\sigma^2 y/2+\sigma\sqrt y Z}-K \right)^+ \right] \mathbb P(L_{T-a}\in dy),
\]
which is an explicit mixture of Black–Scholes call values over the inverse-subordinator law [2511.10371]. A plausible implication is that generalized TFBSEs can often be solved more naturally by probabilistic subordination than by direct PDE methods.

Transform-based formal series also appear. A two-asset Caputo TFBSE has been treated using the Sumudu transform, leading to an integral equation and a recursive infinite-series solution
\[
c(u,v,t)=c(u,v,0)+\sum_{n=1}^\infty g_n(u,v)\frac{t^{n\alpha}}{\Gamma(1+n\alpha)},
\]
with
\[
g_{n+1}(u,v) = -\left[ \frac{1}{2}\sigma_1^2 g_{n,uu} +\frac{1}{2}\sigma_2^2 g_{n,vv} +\rho \sigma_1 \sigma_2 g_{n,uv} -r g_n \right]
\]
[2010.13411]. The same recurrence can be compactly interpreted as an operator series, although the paper itself does not formalize that viewpoint [2010.13411]. Its main limitation is that convergence is asserted but not rigorously proved.

ANN-based work also treats the TFBSE from a semidiscrete perspective. There the Caputo derivative is approximated by a history-weighted scheme, and each time step is solved as a spatial problem by a two-layer feed-forward network trained by Adam [2105.11240]. However, the paper’s explicit fractional Black–Scholes example is a manufactured problem with internal typographical inconsistencies, so its value is methodological rather than canonical [2105.11240].

## 5. Numerical methods

A large fraction of TFBSE research is numerical, and the field is organized around the tradeoff between temporal nonlocality, spatial accuracy, boundary treatment, and the weak initial singularity typical of fractional evolution equations.

The most classical family uses \(L1\)-type approximations for the Caputo derivative on uniform meshes. In the subdiffusive Black–Scholes model, a weighted finite difference method interpolates between implicit, explicit, and fractional Crank–Nicolson analogues through a parameter \(\theta\in[0,1]\), with the fully implicit case \(\theta=0\) unconditionally stable for all \(\alpha\in(0,1)\) and the fractional Crank–Nicolson case conditionally stable [1907.00297]. The same paper proves convergence with order \(2-\alpha\) in time and \(2\) in space under its stability conditions [1907.00297]. A tempered extension modifies the time operator and the convolution weights but retains the same general finite-difference architecture [2103.13679].

High-order compact finite differences on nonuniform meshes constitute another major line. A three-point compact approximation on smooth nonuniform meshes generated by Tavella–Randall or quadratic mappings has fourth-order local spatial accuracy, and when applied to the transformed TFBSE yields observed \(O(h^4+\tau)\) behavior, with fourth-order spatial convergence and roughly first-order temporal accuracy in experiments due to the startup singularity [1604.05178]. This work is important because it ties mesh design directly to two structural features of Black–Scholes operators: degeneracy near \(s=0\) and reduced regularity near the strike [1604.05178].

Later high-order methods explicitly address the weak initial singularity. A nonuniform Alikhanov formula combined with a fourth-order average approximation in space and sum-of-exponentials acceleration produces a fast variable-step solver with unconditional stability and global convergence
\[
\|e^n\| \le C\bigl(\tau^{\min\{\gamma\sigma,2\}}+h^4\bigr),
\]
so second-order temporal accuracy is recovered when the grading parameter satisfies \(\gamma\sigma\ge 2\) [2109.02028]. A tempered analogue combines a compact fourth-order spatial operator, graded time mesh \(\tau_n=T(n/M)^y\), nonuniform tempered \(L1\) approximation, and SOE history compression, achieving temporal order \(\min\{y\alpha,2-\alpha\}\), fourth-order space, and unconditional stability for both direct and fast schemes [2303.10592].

Spline-based collocation methods form a distinct class. Exponential B-spline collocation combined with an \(L1\)-type time discretization yields a tridiagonal system at each time step, unconditional stability by von Neumann analysis, and convergence of order \(O(h_t^{2-\mu}+h_\varkappa^2)\) [2207.09153]. The method has been applied not only to manufactured tests but also to European call, European put, and double barrier knock-out call options [2207.09153]. A more recent exponential B-spline method uses a Crank–Nicolson-type approximation for the fractional derivative, claims unconditional stability, and reports temporal order \(2-\mu\) and second-order spatial convergence on smooth test problems [2602.00201]. Modified cubic B-spline differential quadrature has also been used together with \(L1\) time discretization, with reported fourth-order spatial convergence and order \(2-\alpha\) in time under a matrix inverse bound derived via the Neumann series theorem [2508.06780].

Meshless methods are represented by a global RBF collocation approach for time-space-fractional Black–Scholes-type PDEs. Specializing to the pure time-fractional case by taking \(\beta=1\) yields a TFBSE solver using Caputo \(L1\)-type time stepping, cubic polyharmonic spline interpolation, and a preconditioned global collocation matrix [2011.07710]. The paper emphasizes flexibility in multidimensional settings and condition-number reduction, but it does not provide a full convergence theorem for the fully discrete RBF method [2011.07710].

ANN-based solvers, finally, replace the spatial linear solver by optimization over a neural network ansatz at each time step, with domain mapping for semi-infinite \(S\)-domains and fine-tuning across adjacent time levels [2105.11240]. This suggests a broadening of TFBSE computation beyond classical discretizations, though the validation remains limited [2105.11240].

## 6. Applications, interpretation, and conceptual boundaries

The TFBSE has been used to price a range of European-style derivatives. Standard applications include European calls and puts under time-fractional dynamics [1907.00297] [2207.09153]. More specialized applications include double barrier knock-out call options [2109.02028] [2207.09153] and two-asset basket-type payoffs in a Caputo two-dimensional Black–Scholes PDE [2010.13411]. In the generalized sub-diffusive setting, the pricing theory is extended to European-style contingent claims under a sub-diffusion equivalent martingale measure, with the call price represented by a shifted TFBSE and an explicit inverse-subordinator mixture formula [2511.10371].

A recurring interpretive theme is that the fractional order governs the strength of temporal memory or subdiffusion. In the subdiffusive Black–Scholes model, smaller \(\alpha\) corresponds to more pronounced waiting times or stagnation, interpreted as illiquidity or inactivity [1907.00297]. In the tempered subdiffusive extension, \(\lambda\) weakens the heaviest waiting-time tails, but the effect of a numerically small \(\lambda\) can remain significant when \(\alpha\) is small because \(\lambda^\alpha\) need not be small [2103.13679]. This matters for both modeling and numerical interpretation.

At the level of computed prices, several papers report that the influence of \(\alpha\) is strongest near the money or in barrier-sensitive regions. In the exponential B-spline collocation study, the impact of \(\mu\) on European call prices is mild far from the strike but significant near \(\xi\approx \widetilde K\), while for double barrier knock-out calls the price is reported to be inversely proportional to \(\mu\) when the stock price is near or above strike [2207.09153]. A plausible implication is that fractional memory effects manifest most clearly where payoff geometry and boundary interaction already make the pricing surface sensitive.

The main conceptual boundary in the literature concerns what counts as a TFBSE. A paper on fractional stochastic volatility driven by an fOU process with fBm input derives \(O(\delta)\) corrections to Black–Scholes prices and implied volatilities, including fractional-power maturity effects, but the pricing operator itself remains the classical Black–Scholes PDE plus perturbative corrections [1509.01175]. It belongs to the broader fractional-finance literature, but not to the TFBSE in the PDE sense. This distinction is necessary to avoid conflating memory in volatility dynamics with memory in the pricing equation.

Another internal boundary separates direct fractional-PDE modeling from probabilistically derived TFBSEs. The ad hoc Caputo replacement
\[
\partial_t \to {}_0^C D_t^\alpha
\]
is common and computationally tractable [2207.09153] [2602.00201], but the subdiffusive and generalized inverse-subordinator works provide a more rigorous foundation by deriving the pricing PDE from time-changed processes and martingale pricing [1907.00297] [2511.10371]. This suggests two parallel strands of the literature: phenomenological fractionalization and stochastic-clock derivation.

## 7. Open issues and research directions

The literature identifies several persistent technical issues. The first is the **weak initial singularity** of TFBSE solutions. Uniform-step schemes with formal order \(2-\alpha\) or higher often lose their advertised convergence unless sufficient regularity is assumed. This has driven the development of graded meshes, variable-step Alikhanov formulas, and nonuniform tempered \(L1\) methods [2109.02028] [2303.10592] [1604.05178]. The dominance of this issue suggests that temporal regularity, not merely fractional order, is central to the design of reliable solvers.

The second issue is the **history cost** of Caputo-type derivatives. Direct convolution implies \(O(N^2)\) temporal work in straightforward implementations. Some papers leave this unaddressed [2508.06780] [2207.09153], while others introduce sum-of-exponentials compression or recursive auxiliary variables to reduce cost to essentially linear growth in the number of time steps up to the SOE rank [2109.02028] [2303.10592]. Fast memory evaluation is therefore one of the main algorithmic frontiers of TFBSE computation.

A third issue is the relationship between **high spatial order** and **payoff nonsmoothness**. Fourth-order compact and average schemes perform strongly on manufactured smooth solutions [1604.05178] [2109.02028] [2303.10592], but practical option payoffs contain kinks at the strike. Some papers explicitly note that this can degrade the full formal order and motivate local mesh refinement or piecewise-uniform meshes [2303.10592]. This suggests that future high-order work must treat payoff regularization and adaptive resolution as first-class concerns.

A fourth issue is **financial derivation and arbitrage interpretation**. Many PDE papers simply replace the time derivative by a fractional operator, without deriving the modified pricing equation from self-financing arguments, anomalous diffusion, or time-changed martingale models [2010.13411] [2105.11240]. By contrast, the subdiffusive and generalized sub-diffusive works explicitly derive the TFBSE from time-changed stochastic models and construct appropriate martingale measures [1907.00297] [2511.10371]. A plausible implication is that future TFBSE research may increasingly distinguish between numerically useful fractional surrogates and fully specified no-arbitrage market models.

Finally, multidimensional and nonstandard-contract TFBSEs remain comparatively underdeveloped. Two-asset Caputo PDEs exist in formal series form [2010.13411], and time-space-fractional meshless frameworks can in principle handle higher-dimensional asset spaces [2011.07710], but rigorous, financially calibrated multidimensional TFBSE pricing remains sparse. The same is true for American-style products, where free boundaries would interact nontrivially with temporal memory.

Taken together, these works position the TFBSE as a mature but still evolving subject at the intersection of anomalous diffusion, inverse-subordinator stochastic modeling, and high-order nonlocal PDE numerics. The core idea is stable across the literature: the pricing operator retains Black–Scholes spatial structure while time evolution is made nonlocal. What varies is the derivative class, the probabilistic foundation, the contract geometry, and the numerical strategy used to resolve long memory efficiently and accurately [1907.00297] [2109.02028] [2303.10592] [2511.10371].

Source: https://www.emergentmind.com/topics/time-fractional-black-scholes-equation-tfbse