^{1}

^{2}

^{3}

We replicate and price European options on stocks modeled by time-changed geometric Brownian motion. The time change is obtained as the integrated intensity of random arrival times of price changes of the underlier over the life of the option. For European call options we obtain explicit hedging and pricing formulas. This approach is motivated by the need to connect option prices directly to the microstructure properties of the limit order book that determines tick-by-tick stock price changes. The continuous time model is obtained as an appropriate limit of discrete time random walks with random jump times, in the limit of infinitely many independent representative agents.

In this paper we derive unique European option prices and hedging strategies in a Black-Scholes-style market for which the underlying stock process is time-changed by a measure of market activity. The time change is the compensator of an arrival process, and the option prices are conditional on the integrated intensity of the arrival process over the life of the option. In particular, options, such as call options, that depend only on the terminal price of the underlier may be perfectly replicated as a function of the terminal value of the compensator. In the case of European call options, we obtain an explicit hedging and price formulas as functions of the integrated intensity.

We are motivated by efforts to model the dynamics of the limit order book that underlies stock price formation. From this perspective, stock price movements are the result of limit and market order arrivals collected in an electronic limit order book. Order arrival times are random and described by point processes, most notably by “self-exciting” Hawkes processes, e.g. [

As a model for the high frequency limit order book microstructure, we consider price changes as the limit of generalized random walks for which the i.i.d. space increments (e.g. log price changes) occur at random times given by a simple point process and associated counting process. Considering the effect on market price to be the summed contribution of many small such changes, we obtain as a scaled limiting price process a time-changed Brownian motion, where the time change is the compensator of the original counting process. With mild assumptions, the resulting time-changed Brownian motion is a continuous square integrable martingale with uncorrelated, but not necessarily independent, increments, as is consistent with stylized facts about equity prices.

The limiting process serves as the basis for a Black-Scholes style option pricing model that we call “activity-driven Black-Scholes”, to emphasize its connection to the market microstructure. It incorporates the characteristics of the chosen point process describing the “excitability” of order arrivals causing price changes in the underlying limit order book. This affords a perspective on a heteroscedastic stock price model for option pricing in which the fluctuating volatility can be intrinsically connected to a calibrated model of the underlying order book. When the point process is a homogeneous Poisson process with unit intensity, we recover the usual Brownian motion and Black-Scholes model.

When there is a second source of randomness—such as market activity—affecting the stock price, replication of an option in the ordinary Black-Scholes sense is not possible using only the underlying stock and a numeraire. But we can still obtain conditional pricing and replication. This is the content of our main theorem and corollaries in Section 4.

We mention two ways that this kind of conditional pricing and hedging can be useful. First, if there is an underlying calibrated limit order book model defining a market activity process, this will define a distribution of integrated intensities, and hence a calibrated distribution of options prices, distributed over the range of market activity outcomes. The width of this distribution is a way to measure the “volatility risk” of the option price. Second, the possible trajectories of market activity intensity can be considered a way to parametrize “scenarios” for use by banks in stress-testing options portfolios, as required by banking regulators. A single underlying market activity scenario Λ ( t ) can be used to underly the prices of a basket of different options that share the same scenario.

We next briefly summarize our results. We say that a simple counting process N is regular if the compensator Λ of N is continuous, strictly increasing, and Λ ( ∞ ) = ∞ . This is a mild assumption including most examples of interest.

For any regular counting process N with compensator Λ and for any integer n ≥ 1 there is a natural regular counting process N n with compensator Λ n = n Λ , which we use as follows. Given an i.i.d. sequence { ε i : i ≥ 1 } of random variables with E [ ε i ] = 0 , V a r ( ε i ) = σ 2 < ∞ , and independent of N n for all n, we define a sequence { S N n : n = 1 , 2 , 3 , ⋯ } of rescaled random walks over N by

S N n ( t ) = 1 n ∑ i = 1 N n ( t ) ε i .

These random walks over N can be compared to the standard sequence of rescaled random walks over deterministic (integer) times

W n ( t ) = 1 n ∑ i = 1 ⌊ n t ⌋ ε i ,

that are well-known (Donsker’s Theorem [

A similar fact is true for the random walks S N n over N: the sequence { S N n : n ≥ 1 } of stochastic processes defined above converges weakly to the scaled, time-changed Brownian motion σ B ∘ Λ .

Since the time change Λ need only be regular, we do not require it to be a Levy subordinator, or even Markov. This opens the door to “self-exciting” point processes such as those often used to model order arrivals in the high-frequency limit order book.

In general, B ( Λ ( t ) ) is a continuous square integrable martingale; its quadratic variation is Λ ; it is a standard Brownian motion if and only if N is Poisson with unit intensity; it has uncorrelated increments, but independent increments if and only if Λ is deterministic.

Our activity-driven Black-Scholes model for the stock price is

S ( t ) = S 0 exp ( σ B ( Λ ( t ) ) + μ Λ ( t ) ) , (1)

where B is a standard Brownian motion, along with the numeraire

B t = exp ( r t ) .

This is flexible enough to display most of the stylized facts of stock price returns. For certain contingent claims X paying off at the maturity time T, and conditional on the value of Λ ( T ) , we show X is attainable (may be replicated by a self-financing strategy), find an equivalent martingale measure Q for the discounted stock, and obtain a no-arbitrage pricing formula of the form

V t ( X ) = E Q [ B t B T X | F t ] .

This is stated in detail in Theorem 4, where we formalize conditioning on Λ ( T ) by the condition that Λ is “pinned” at T (Section 3.2).

As an application, we can give explicit conditional hedging and pricing formulas (Corollary 3) for a vanilla European call option, conditional on the integrated intensity Λ ( T ) = T ′ . The price formula at t = 0 turns out to be the classical Black-Scholes formula, but with the strike price K replaced by e r ( T ′ − T ) K , and the maturity T replaced by T ′ . The resulting call price is an increasing function of T ′ , which can be considered a “market activity” or “market clock” parameter similar in effect to the volatility.

Financial applications of time-changed Lévy processes have been extensively studied, most significantly in the work of Peter Carr and collaborators, for example in [

Often in the literature attention is directed toward particular Markov time change processes such as integrated affine processes ( [

The work of [

Some extensions of [

A related direction in the literature has been the study of option pricing for discrete time models motivated by the binomial tree model, e.g. [

The comprehensive book [

The remainder of the paper is organized as follows. In Section 2 we describe random walks over continuous time point processes and their rescalings and limits. Sections 3 and 4 describe our activity-driven Black-Scholes model, and the pricing and hedging of terminal-time-payoff options. It contains the main theorem, corollaries, and discussion. Section 5 summarizes, and the proof of the main theorem appears in Section 6.

In this section we define some terms and describe a class of random walks over counting processes suited to our purposes. For background, see for example [

We consider only simple, non-explosive point processes on [ 0, ∞ ) , i.e. sequences { T n } of [ 0, ∞ ) -valued random variables on a probability space ( Ω , F , P ) such that T 0 = 0 , T n < T n + 1 a.s. for all n, and lim n → ∞ T n = + ∞ . For such a point process { T n } , the corresponding counting process is

N ( t ) = ∑ n ≥ 1 1 ( T n ≤ t )

with N ( 0 ) = 0 . The natural filtration F t N = σ ( N s : s ≤ t ) of N is automatically right continuous ( [

In this paper we restrict attention to the (large) class of regular counting processes as defined in the introduction.

To set notation, we denote the classical one-dimensional random walk W ( t ) to be a piecewise constant cadlag stochastic process defined by

W ( t ) = ∑ i = 1 ⌊ t ⌋ ε i ,

where { ε i : i ≥ 1 } is an i.i.d. sequence of random variables, and we assume E [ ε i ] = 0 and V a r ( ε i ) = σ 2 < ∞ .

Motivated by the interpretation of W as a log-price process, we can think of W as the running sum of a sequence ε i of random contributions, or price shocks, delivered at times i = 1 , 2 , 3 , ⋯ by a representative agent at unit frequency as a result of trading.

Now instead of a single agent contributing at unit frequency, we can imagine n smaller agents also contributing at unit frequency, but contributing only ε i / n each independently, and where the deterministic contribution times are spread out uniformly, yielding

W n ( t ) = 1 n W ( n t ) = 1 n ∑ i = 1 ⌊ n t ⌋ ε i .

This process converges weakly as n → ∞ to the multiple σ B ( t ) of a standard Brownian motion B ( t ) (Donsker’s Theorem). It can be interpreted, in the limit, as the running total of the contributions of infinitely many infinitesimal agents all contributing independently with unit frequency, and forms the basis of the Black-Scholes model.

Now we would like to imagine a different representative agent for which the contribution times are the random jump times of a regular counting process N ( t ) , instead of the deterministic times t = 1 , 2 , 3 , ⋯ . The counting process may be quite general, including a Hawkes process or another self-exciting, non-Markov process. We assume it is independent of the i.i.d. sequence { ε i } . The corresponding random walk is denoted

S N ( t ) = W ( N ( t ) ) = ∑ i = 1 N ( t ) ε i .

If we denote by Λ ( t ) the compensator of N, then Λ may be interpreted as the time integral of an arrival rate intensity λ ( t ) , if Λ is absolutely continuous with respect to t.

Corresponding to the scaling of the unit frequency agent, we may also imagine we have n smaller agents, each contributing an independent amount ε i / n with the same arrival intensity. The integrated arrival intensity of n such agents will be n Λ ( t ) . It is straightforward to show that the corresponding counting process N n with compensator n Λ may be defined by

N n ( t ) = N ( Λ − 1 ( n Λ ( t ) ) ) .

(Here Λ − 1 ( t ) denotes the functional inverse of Λ ( t ) , well-defined since Λ is strictly increasing.)

We may then define the scaled random walk defined by a regular counting process N ( t ) to be

S N n ( t ) = 1 n W ( N n ( t ) ) = 1 n ∑ i = 1 N n ( t ) ε i .

As before, we can interpret this as the running total of the independent contributions of n small agents, each with rate described by the compensator Λ of N.

We note that, like the rescaled classical random walk, the random walk S N n is a martingale with respect to the natural filtration and there are explicit expressions for the variance and covariances in terms of N [

Just as the classical random walks W n ( t ) converge weakly to σ B ( t ) , so also the random walks S n converge to a time-changed Brownian motion. A general reference on weak convergence in this context is [

Theorem 1 Let ε i be an i.i.d. sequence of random variables with mean 0 and variance σ 2 . Given a regular counting process N with compensator Λ , the sequence { S N n : n ≥ 1 } of rescaled random walks over N converges weakly, in the Skorokhod metric, to the time-changed Brownian motion σ B ∘ Λ . Moreover, B and Λ are independent.

The strategy of proof is to establish that ( W n , 1 n N n ) converges weakly to ( σ B , Λ ) in the product space, and then apply the composition operator to both sides, using the Continuous Mapping Theorem and the continuity of composition.

Since we do not assume that Λ is a subordinator, the limit σ B ∘ Λ , while continuous, need not be a Lévy process, or even Markov. Theorem 1 is, in spirit, a generalization of Donsker’s theorem. The Donsker case is recovered when N is chosen to be a homogeneous Poisson process with unit intensity, in which case the compensator is Λ ( t ) = t and the weak limit becomes σ B ( t ) , a multiple of standard Brownian motion.

The process B ∘ Λ depends on N because Λ depends on N. To emphasize that dependence, we will use the notation

B N = B ∘ Λ , i .e . B N ( t ) = B ( Λ ( t ) ) .

To put this into words:

The limiting random walk B N is a time-changed Brownian motion, where the time change is the compensator of the point process describing the trading times of a representative infinitesimal agent.

Let { H t } be the filtration defined by

H t = σ ( B ( r ) , Λ − 1 ( u ) : r ≤ t , u ≤ t ) .

By the usual augmentation if necessary, we assume that these filtrations satisfy the usual conditions.

It is shown in [

Furthermore:

1) [ B N , B N ] t = Λ ( t ) ,

2) E [ B N ( t ) ] = 0 ,

3) c o v ( B N ( s ) , B N ( t ) ) = E [ N ( s ∧ t ) ] for all s , t ≥ 0 ,

4) B N has uncorrelated increments, and

5) B N has independent increments if and only if Λ is deterministic.

These properties make B N useful as a model for financial price processes. Financial returns time series can have close to zero autocorrelation, but squared returns show positive autocorrelation due to typically observed heteroscedasticity and the corresponding failure of independence of returns. Choice of a non-deterministic Λ naturally provides us with a stock price model consistent with these stylized facts of uncorrelated but dependent increments.

Here we introduce some brief preliminaries prior to stating our main results.

The Hawkes process (e.g. [

λ t = α + β ∫ 0 t μ ( t − s ) N ( d s ) > 0 ,

where α > 0 and the response function μ ( t ) is a positive function satisfying ∫ 0 ∞ μ ( t ) d t < ∞ . The (non-Markov) compensator is Λ ( t ) = ∫ 0 t λ t d t . Then the random walks S N n , n ≥ 1 , can be the basis of a model for the log stock price that will reflect the self-exciting nature of the Hawkes process.

Modelers who use this approach to describe the activity of the high-frequency limit order book then have access to the option price model described below if they have a parametrized distribution for a terminal value Λ ( T ) .

Definition 2 A continuous random process Λ ( t ) with Λ ( 0 ) = 0 is pinned at time T if Λ ( T ) = T ′ for some deterministic T ′ .

The Brownian bridge is a familiar example. Absolutely continuous, increasing, even self-exciting pinned processes are not hard to come by. For example, let r t be any positive adapted RCLL process, such as the intensity of a Hawkes process. Let R t = e − r t < 1 . For any T ′ , c > 0 with c < T ′ / T 2 , define

Λ ( t ) = T ′ T t − c ( T − t ) ∫ 0 t R s d s .

Then Λ ( 0 ) = 0 , Λ ( T ) = T ′ , Λ ( ∞ ) = ∞ , Λ is absolutely continuous, adapted to the natural filtration of r, and is increasing on [ 0, ∞ ) because

Λ ′ ( t ) = T ′ T + c ∫ 0 t R s d s − c ( T − t ) R t > T ′ T − c ( T − t ) > 0.

The following proposition demonstrates that regular counting processes with pinned compensators are easily constructed.

Proposition 3 Denote by P ( t ) the homogeneous Poisson process with rate 1, and let Λ be any absolutely continuous, strictly increasing process pinned at time T. For each t, denote by F t the sigma-algebra F t P ∨ F t Λ − 1 = σ ( P ( s ) , Λ − 1 ( s ) : s ≤ t ) . Then Λ ( t ) is an F -stopping time for each t, and

N ( t ) = P ( Λ ( t ) )

is a regular counting process on [ 0, T ] with compensator Λ with respect to F Λ .

Proof. Given s < t ≤ T , apply the optional stopping theorem (e.g. [

The familiar Black-Scholes option pricing framework for a bond and stock price process B t and S ( B S ) ( t ) , given by

B t = exp ( r t ) , S ( B S ) ( t ) = S 0 exp ( σ B ( t ) + μ t ) ,

is often motivated as a limit of the exponential of a binomial random walk of the form W n ( t ) discussed above. At short time scales, however, we tend to observe non-independent increments and some level of self-exciting behavior in the limit order book driving the price.

When N is any regular point process describing the arrival of orders to the order book, we can call the corresponding limiting market model an “activity-driven Black-Scholes model”, given by

B t = exp ( r t ) , S ( t ) = S ( B S ) ( Λ ( t ) ) = S 0 exp ( σ B N ( t ) + μ Λ ( t ) ) , (2)

where r , σ , μ are positive constants. Our interpretation of the time change Λ is the compensator of the arrival point process N of trades for a representative infinitesimal agent in the market.

This turns out to be a suitable framework for option pricing that can incorporate heteroscedasticity and other non-stationary stylized facts of stock price behavior. The goal of this section is to show that in certain cases, conditional on Λ ( T ) , we can also obtain a no-arbitrage pricing formula analogous to the standard Black-Scholes formula.

Recall the filtration H t = σ ( B ( r ) , Λ − 1 ( u ) : r , u ≤ t ) with respect to which Λ ( t ) is a stopping time for each t. If we define

G t = H Λ t ,

then both B N and Λ , hence S ( t ) , are G -adapted processes on ( Ω , F , P ) .

It’s convenient to define the additional filtrations, assumed complete:

F t B = σ ( B ( s ) : s ≤ t ) ⊂ H t ,

F t Λ − 1 = σ ( Λ − 1 ( s ) : s ≤ t ) ⊂ H t ,

F t Λ = σ ( Λ ( s ) : s ≤ t ) ,

F t S = σ ( S ( s ) : s ≤ t ) ⊂ G t .

To discuss the problem of option pricing, we review some standard terminology. We will say that a pair ( ϕ t , ψ t ) is a self-financing strategy if ϕ and ψ are G -predictable processes such that if V t = ϕ t S ( t ) + ψ t B t , then d V t = ϕ t d S ( t ) + ψ t d B t . Here we interpret ϕ t as the number of shares of stock held in a portfolio at time t, ψ t the same for the bond, and so V t is the time-t portfolio value.

If we fix a deterministic maturity time T, we call a random payoff X ∈ G T at time T a claim. A claim X is an attainable claim if there exists a self-financing strategy ( ϕ , ψ ) with

X = V T = ϕ T S ( T ) + ψ T B T ,

i.e. the self-financing strategy replicates the claim at the terminal time. When this is the case, ( ϕ , ψ ) is a replicating portfolio for X, and V t will be the no-arbitrage price of the claim X at any time t < T .

In the classical Black-Scholes model with a stock and a bond, every claim is attainable (the market is complete) due to the martingale representation property for Brownian motion. In the context of the market model (2), a suitable representation property for the time-changed Brownian motion B N is not available. Instead, we will show that certain classes of claims are attainable, and establish a familiar-looking conditional expectation formula for the claim price, with respect to a suitably defined measure. The remainder of this section spells out the details.

Market Assumptions:

1) ( Ω , F , P ) is a probability space supporting a Brownian motion B ( t ) and an independent regular point process N with absolutely continuous compensator Λ , with Λ ( t ) = ∫ 0 t λ u d u , where λ is left-continuous with right limits and for some ε > 0 , for all u, λ u > ε a.s.; and

2) there are two tradable assets, a stock S ( t ) and bond B t given by

B t = e r t , S ( t ) = S 0 exp ( σ B N ( t ) + μ Λ ( t ) ) , σ , μ > 0 , r ≥ 0. (3)

3) X ∈ G T ∩ L 1 ( Ω , F , P ) is a European option contract payoff at a fixed maturity T > 0 .

Theorem 4 (Pricing and Hedging a European option) Under the Market Assumptions above, suppose further that the compensator Λ of N is pinned at T. Let Q be the equivalent martingale measure defined by

d Q d P = exp ( − ∫ 0 Λ ( T ) γ t d B ( t ) − 1 2 ∫ 0 Λ ( T ) γ t 2 d t ) , (4)

where

γ t = 1 σ ( μ + σ 2 / 2 − r λ Λ − 1 ( t ) ) , 0 ≤ t ≤ Λ ( T ) . (5)

Then the discounted stock price Z t = B t − 1 S ( t ) (stopped at T) is a ( Q , G ) -martingale, and B ˜ ( t ) = B ( t ) + ∫ 0 t γ s d s is a standard Brownian motion with respect to ( Q , H ) .

Moreover, suppose in addition that the payoff X is Q-independent of F T Λ .

Then

1) X is attainable and its replicating portfolio is ( ϕ , ψ ) , with value

V t = ϕ t S ( t ) + ψ t B t ,

where

ϕ t = η Λ ( t ) , ψ t = ξ Λ ( t ) ,

and ( η , ξ ) is the replicating portfolio for the payoff X at time Λ ( T ) in the standard risk-neutral Black-Scholes model with stock Y ( t ) and bond B t given by

B t = e r t , Y ( t ) = S 0 exp ( σ B ˜ ( t ) − 1 2 σ 2 t ) , t ≤ Λ ( T ) . (6)

2) The unique no-arbitrage price V t of X at time t < T is

V t = ( B t / B T ) E Q [ X | G t ] . (7)

The proof is postponed to the Appendix. We remark that the Theorem provides explicit hedging and pricing, conditional on Λ ( T ) , of a general class of European option payoffs that might potentially depend on the whole history of S ( t ) . Moreover, the hedging strategy is attainable at any time t as long as Λ ( t ) is observable at time t. The hypothesis that Λ is pinned at T is equivalent to pricing that is conditional on the value Λ ( T ) . The additional assumption that the payoff X is Q-independent of F T Λ is made more concrete in the following Corollaries describing some interesting special cases.

Corollary 1 Under the Market Assumptions above, if N is an inhomogeneous Poisson process with positive intensity, then X is attainable and its unique no-arbitrage price at time t < T is

V t = ( B t / B T ) E Q [ X | G t ] ,

where the equivalent martingale measure is defined in Equation (4).

Proof. An inhomogeneous Poisson process has a deterministic intensity, hence deterministic compensator Λ . The conclusion is immediate from Theorem 4.

Corollary 2 Under the Market Assumptions above, suppose that the compensator Λ of the regular point process N is pinned at time T. Let X ∈ L 1 be the payoff of a European option of the form X = f ( S ( T ) ) for some function f.

Then X is attainable and the unique no-arbitrage price V t of X at time t < T is

V t = ( B t / B T ) E Q [ X | G t ] , (8)

where the measure Q is the equivalent martingale measure defined by Equation (4).

Moreover, the time 0 price of the option X = f ( S ( T ) ) is given by

V 0 = e − r T ∫ − ∞ ∞ f ( S 0 e − σ 2 T ′ / 2 + r T e σ y ) p ( y ; σ 2 T ′ ) d y , (9)

where T ′ = Λ ( T ) and p ( y ; σ 2 T ′ ) is the pdf of a normal distribution in y with mean zero and variance σ 2 T ′ .

Proof of Corollary 2. Let Λ ( T ) = T ′ . Since S ( T ) = exp ( σ B ( Λ ( T ) ) + m Λ ( T ) ) = exp ( σ B ( T ′ ) + m T ′ ) , it follows that S ( T ) , hence X, is independent of Λ and Λ − 1 . The first part then follows from Theorem 4.

The second part follows straightforwardly from (8) and the fact, as shown in the proof of Theorem 4, that

S ( T ) = S 0 exp ( σ B ˜ ( Λ ( T ) ) − ( 1 / 2 ) σ 2 Λ ( T ) + r T ) = S 0 exp ( σ B ˜ ( T ′ ) − ( 1 / 2 ) σ 2 T ′ + r T ) ,

where B ˜ ( T ′ ) is normally distributed with respect to Q with mean zero and variance T ′ .

Corollary 3 Under the same assumptions of Corollary 2, if X = ( S T − K ) + is a European call option with strike K and maturity T, then its time-0 price is

V 0 = C B S ( σ , r , S 0 , e r ( T ′ − T ) K , T ′ ) , (10)

where T ′ = Λ ( T ) and C B S ( σ , r , S 0 , K , T ) denotes the standard Black-Scholes call price at strike K and maturity T.

Furthermore, the hedging portfolio ( ϕ t , ψ t ) is given explicitly by

ϕ t = N ( d 1 Λ ( t ) ) and ψ t = − K e − r T ′ N ( d 2 Λ ( t ) ) , (11)

where N is the standard normal cumulative distribution function, and

d 1 Λ ( t ) = 1 σ T ′ − Λ ( t ) [ ln ( S t K + ( r + σ 2 / 2 ) ( T ′ − Λ ( t ) ) ) ] , (12)

d 2 Λ ( t ) = d 1 Λ ( t ) − σ T ′ − Λ ( t ) , (13)

and the time-t price of the option is

V t = ϕ t S t + ψ t e r t .

Corollary 3 follows from Corollary 2, part (a) of Theorem 4, and the standard Black-Scholes formulas for the replicating portfolio ( η , ξ ) of a European call in the market (6).

Notice that in the above results the quantity Λ ( t ) is observable at time t, as it is the quadratic variation of ( log S ) / σ up to time t. Therefore the price formulas and replicating portfolios depend only on the random variable T ′ = Λ ( T ) . We can view Equations (9) and (10) as the price conditional on Λ ( T ) , and therefore a scenario-based price, where the scenarios are parametrized by Λ ( T ) . These results would be relevant to stress-testing portfolio values under different possible future regimes of accumulated market activity over the life of the option.

Alternatively, if we have a separately calibrated limit order book model with arrival process N leading to a parametrized distribution for Λ ( T ) , we can take a further expectation of Equation (10) to obtain an unconditional option price.

It is straightforward to check that the expression in Equation (10) is increasing in T ′ = Λ ( T ) when the other parameters are held constant. Therefore one can think of the accumulated intensity (or realized quadratic variation) Λ ( T ) over the life of the option as playing a similar role as the volatility parameter σ in affecting the option price. The parameter Λ ( T ) is a “trading activity” parameter separate from volatility but having a similar qualitative effect on the option price. This model is therefore an alternative to a stochastic volatility model.

For a numerical illustration comparing Black-Scholes volatility σ as a parameter to the integrated arrival intensity Λ ( T ) , see

In

Motivated by the random jump times in the limit order book for a stock price and Theorem 1, we consider a generalized Black-Scholes model with a bond or cash account B t = e r t and a stock S ( t ) = S 0 exp ( σ B N ( t ) + μ Λ ( t ) ) . This market can be thought of as a continuous limit of a discrete stock price model where price changes are driven by the action of many small agents acting with integrated intensity Λ , reflecting underlying limit order book activity. In this sense the market model’s heteroskedastic features are derived from market-clock variations rather than imposed by an exogenously estimated stochastic volatility.

For certain classes of option payoffs X, and conditional on Λ ( T ) , we establish an option pricing formula in a familiar form

V t = ( B t / B T ) E Q [ X | H Λ ( t ) ] ,

where Q is the explicitly defined risk-neutral measure and H = F B ∨ F Λ − 1 . As an application, we can price a European call option with strike K and maturity T, conditional on the value Λ ( T ) = T ′ of the integrated intensity of the counting process N over the life of the option, as

C B S ( σ , r , S 0 , e r ( T ′ − T ) K , T ′ ) ,

where C B S is the usual Black-Scholes call option price formula. A separately calibrated risk-neutral distribution of Λ ( T ) could provide an unconditional call option price as a function of parameters of a limit order book model.

Given X, our goal is to construct a self-financing strategy in the stock and bond with time-t value given by V t of Equation (7).

In this proof we use the convention that when we call a process M ( t ) a martingale that is only defined on an interval [ 0, t 0 ] , we mean that the process stopped at t 0 is a martingale.

Since Λ is absolutely continuous with derivative λ bounded below by ε > 0 , the inverse function theorem tells us that

Λ − 1 ( t ) = ∫ 0 t λ ¯ u d u ,

where

λ ¯ u = 1 λ Λ − 1 ( u )

is positive and bounded above by 1 / ε . It follows that

γ t = 1 σ ( μ + σ 2 / 2 − r λ ¯ t )

is bounded and left continuous with right limits. By Girsanov’s Theorem ( [

B ˜ ( t ) = B ( t ) + ∫ 0 t γ s d s = B ( t ) + μ + σ 2 / 2 σ t − r σ ∫ 0 t λ ¯ u d u = B ( t ) + μ + σ 2 / 2 σ t − r σ Λ − 1 ( t )

is a standard Brownian motion with respect to ( Q , H ) . Following our convention, we write

B ˜ N ( t ) = B ˜ ( Λ ( t ) ) = B N ( t ) + μ + σ 2 / 2 σ Λ ( t ) − r σ t .

B ˜ N ( t ) is a square integrable ( Q , G ) -martingale, and since [ B ˜ N , B ˜ N ] = Λ , therefore so is the discounted stock price

Z t ≡ B t − 1 S ( t ) = S 0 exp ( σ B N ( t ) + μ Λ ( t ) − r t ) = S 0 exp ( σ B ˜ N ( t ) − 1 2 σ 2 Λ ( t ) ) .

We are assuming that the option payoff X is independent of F T Λ . This is equivalent to the independence of X from F Λ ( T ) Λ − 1 , since these two sigma-algebras are equal.

To construct a self-financing replicating portfolio, the difficulty is that the Predictable Representation Property (PRP, [

Our strategy is to apply the PRP to the Brownian motion model, and then change variables by means of the time change Λ and the optional stopping theorem.

Recall that B N , Λ , and S are adapted to the filtration G t = H Λ t , and B , Λ − 1 are adapted to H t .

Let

Y ( t ) = Z ( Λ − 1 ( t ) ) = S 0 exp ( σ B ( t ) − 1 2 σ 2 t ) .

Evidently Y is a ( Q , F B ˜ ) martingale.

Now define

E ( t ) = E Q [ B T − 1 X | H t ] .

Since

H = F B ∨ F Λ − 1 = F B ˜ ∨ F Λ − 1 ,

the independence of X from F Λ ( T ) Λ − 1 implies X ∈ F Λ ( T ) B ˜ and E ( t ) = E Q [ B T − 1 X | F t B ˜ ] for t ≤ Λ ( T ) .

Since E Q [ B T − 1 X | F t B ˜ ] is a ( Q , F B ˜ ) -martingale, by the Brownian martingale representation property (e.g. [

d E ( t ) = η t d Y ( t ) . (14)

Next we wish to compose this equation with Λ , which is justified by the following lemma.

Lemma 5 ( [

If H ∈ L ( X , H ) , then H C t − ∈ L ( X C , H C ) for all t. Moreover, with probability one, for all t ≥ 0 ,

∫ 0 C t H s d X s = ∫ 0 t H C s − d X C s .

Since Λ is a finite time-change with respect to H , and letting ϕ t = η Λ ( t ) , then ϕ t is H Λ -predictable and we obtain from Equation (14) that

d E ( Λ ( t ) ) = ϕ t d Z ( t ) . (15)

We may now consider a portfolio holding ϕ t shares of stock and ψ t = E ( Λ ( t ) ) − ϕ t Z ( t ) shares of the bond at time t.

The portfolio value process

V t = ϕ t S ( t ) + ψ t B t = B t E ( Λ ( t ) )

is self-financing by virtue of an easy computation using (15).

To complete the argument, the optional stopping theorem (e.g. [

E ( Λ ( t ) ) = E Q [ B T − 1 X | H Λ ( t ) ] . (16)

Our portfolio strategy is therefore a replicating strategy because

V T = B T E ( Λ ( T ) ) = B T E Q [ B T − 1 X | H Λ ( T ) ] = X .

Therefore the no-arbitrage price of X at any earlier time t must be the value of the replicating portfolio

V t = B t E ( Λ ( t ) ) = ( B t / B T ) E Q [ X | H Λ ( t ) ]

as desired. This completes the proof of Theorem 4.

The first author was partially supported by a Simons Collaboration Grant for Mathematicians.

The authors declare no conflicts of interest regarding the publication of this paper.

Kercheval, A., Salehy, N. and Salehy, N. (2022) Pricing and Hedging Options Conditional on Market Activity. Journal of Mathematical Finance, 12, 1-19. https://doi.org/10.4236/jmf.2022.121001