Next Article in Journal
Claim Prediction and Premium Pricing for Telematics Auto Insurance Data Using Poisson Regression with Lasso Regularisation
Next Article in Special Issue
Dynamics of Foreign Exchange Futures Trading Volumes in Thailand
Previous Article in Journal
Using the Fuzzy Version of the Pearl’s Algorithm for Environmental Risk Assessment Tasks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Asset Pricing in a Unified Bachelier–Black–Scholes–Merton Model

by
W. Brent Lindquist
1,
Svetlozar T. Rachev
1,
Jagdish Gnawali
1 and
Frank J. Fabozzi
2,*
1
Department of Mathematics and Statistics, Texas Tech University, Lubbock, TX 79409-4012, USA
2
Carey Business School, Johns Hopkins University, Baltimore, MD 21202, USA
*
Author to whom correspondence should be addressed.
Risks 2024, 12(9), 136; https://doi.org/10.3390/risks12090136
Submission received: 1 August 2024 / Revised: 19 August 2024 / Accepted: 21 August 2024 / Published: 27 August 2024
(This article belongs to the Special Issue Financial Derivatives: Market Risk, Pricing, and Hedging)

Abstract

:
We present a unified, market-complete model that integrates both Bachelier and Black–Scholes–Merton frameworks for asset pricing. The model allows for the study, within a unified framework, of asset pricing in a natural world that experiences the possibility of negative security prices or riskless rates. Unlike the classical Black–Scholes–Merton, we show that option pricing in the unified model differs depending on whether the replicating, self-financing portfolio uses riskless bonds or a single riskless bank account. We derive option price formulas and extend our analysis to the term structure of interest rates by deriving the pricing of zero-coupon bonds, forward contracts, and futures contracts. We identify a necessary condition for the unified model to support a perpetual derivative. Discrete binomial pricing under the unified model is also developed. In every scenario analyzed, we show that the unified model simplifies to the standard Black–Scholes–Merton pricing under specific limits and provides pricing in the Bachelier model limit. We note that the Bachelier limit within the unified model allows for positive riskless rates. The unified model prompts us to speculate on the possibility of a mixed multiplicative and additive deflator model for risk-neutral option pricing.

1. Introduction

This paper develops a unified Bachelier and Black–Scholes–Merton (BSM) market model for dynamic asset pricing. By incorporating a weighted mixture of arithmetic and geometric Brownian motion, the unified model enables the study of asset pricing in a natural world that experiences the possibility of negative security prices and riskless rates. The model provides a cohesive framework that leverages the strengths of both formulations while mitigating their respective weaknesses.
The seminal work of Bachelier (1900) introduced the precursor of what is now referred to as arithmetic Brownian motion. However, since the latter half of the 20th century, asset pricing theory has been dominated by that of the BSM model based on geometric Brownian motion (Black and Scholes 1973; Merton 1973). While there has been an extensive discussion of, and effort to, address the restrictive assumptions of the BSM model, its strengths (such as positive asset prices and positive riskless rates) have been challenged by the realities of 21st-century financial markets. What was perceived as a primary weakness of the arithmetic Brownian motion model (the possibility of negative asset prices and negative riskless rates) is now seen as having application to price processes for certain derivatives (such as commodity futures) as well as a dynamic theory of asset pricing that includes environmental, social and governance (ESG) scores (Rachev et al. 2024).
The ESG ratings industry continues to mature and its place in the U.S. financial system continues to be debated. Previous works (Pástor et al. 2021; Zerbib 2022; Pedersen et al. 2021) have explored how ESG ratings might affect asset valuation. As argued in Rachev et al. (2024), the ESG-adjusted stock price of company X at time t 0 is
A t = S t E S G = S t X 1 + γ E S G Z t X ; I R , Z t X ; I = Z t X Z t I Z t I   R ,
where S t X > 0 is the stock price of the company; Z t X [ 0,100 ] is the ESG score of the company; Z t I [ 0,100 ] is the ESG score of the market index I to which the company stock belongs; and the parameter γ E S G R measures the ESG affinity of the market. As Z t X ; I will be negative for low ESG stocks, if the market ESG affinity is valued at γ E S G < 1 / Z t X ; I , the ESG-adjusted price A t will be negative for such stocks. A Bachelier model, rather than the BSM model, is then required for modeling ESG-adjusted prices.
Negative asset valuation challenges economists and practitioners whose worldview is anchored in BSM theory. Under the Rachev et al. model, a company may have a negative ESG-adjusted stock price while retaining a positive financial stock price. Such a company has a negative valuation relative to its contribution to the long-term survival of the planet while retaining financial solvency. While stocks continue to be bought and sold under financial prices, the ESG-adjusted price adds sustainability as a third dimension in assessing the value of an asset.1 In this regard, the ESG incorporating model of Rachev et al. (2024) contrasts with models of Pástor et al. (2021), Zerbib (2022) and Pedersen et al. (2021), who incorporate ESG ratings strictly within the BSM world of positive prices. The unified model presented here widens the BSM worldview, incorporating BSM as one limit of a larger dynamic pricing model.
In contrast to the BSM model, which has driven continuous-time finance theory based upon stochastic differential equations, leading to important results such as the fundamental theorem of asset pricing (Delbaen and Schachermayer 1994; Delbaen and Schachermayer 1998) and the development of discrete pricing models (Cox et al. 1979; Jarrow and Rudd 1983), the academic literature on Bachelier’s model has been “scarce” and “scattered” (Choi et al. 2022). Prior to the work of Rachev et al. (2024), published Bachelier models have either:
(i)
used a risk-free rate of zero (e.g., Bachelier 1900; Shiryaev 1999; Schachermayer and Teichmann 2008);
(ii)
used arithmetic Brownian motion for the risky asset with a continuously compounded risk-free asset (Choi et al. 2022); or
(iii)
utilized an Ornstein–Uhlenbeck process for the risky asset (Brooks and Brooks 2017).
Of these three choices, only (i) results in a correct model. However, with no risk-free rate, model (i) has limited practical use. Choices (ii) and (iii) produce price solutions (either under the risk-neutral measure or under the natural measure) whose mean value will eventually diverge to infinity if a positive risk-free rate is chosen. Choosing a negative risk-free rate (for all time) is financially unrealistic, while we have noted that the choice of a zero risk-free rate is of limited use.
We refer to our unified model as the Bachelier–BSM (BBSM) model. In Section 2 we define the BBSM market pricing models. Appropriate parameter choices reduce the unified model to either the traditional BSM model based upon geometric Brownian motion or to a “modernized” Bachelier model based upon arithmetic Brownian motion (Rachev et al. 2024). In particular, the modernized Bachelier (MB) limit allows for either a negative or positive risk-free rate (with no solution divergence). We note that negative interest rates indeed materialized for significant periods as a result of the 2007–2009 Great Recession (Neufeld 2022) as central banks of several developed economies implemented negative interest rate policies to stimulate economic activity and combat deflationary pressures.
In contrast to the BSM model, the BBSM model produces different equivalent martingale measures (EMMs) and different option price solutions depending on whether the replicating, self-financing portfolio treats the riskless asset as a financial product issued in units (i.e., a bond) or as a single bank account. Pricing for a European contingency claim under these two replicating portfolios is presented in Section 3 and Section 4. In Appendix A and Appendix B, we develop closed-form call option price solutions for the BBSM and MB models assuming time-independent coefficients.
In Section 5, we examine option pricing in the case of a dividend-paying risky asset. The question of the existence of a perpetual derivative in the BBSM model is investigated in Section 6 where we show a necessary condition for its existence. In Section 7, we shift our attention to discrete option pricing in the BBSM model using binomial trees. In Section 8, we investigate the term structure of interest rates in this model by examining the pricing of zero-coupon bonds, forward contracts, and futures contracts. Section 9 explores the calibration of the model’s parameters to empirical data.
While the BSM model and the BBSM model under the replicating portfolio of Section 3 result in equivalent martingale measures obtained by multiplicative price deflators, the replicating portfolio of Section 4 produces deflated prices that can be described as a combination of multiplicative and additive deflation. In Section 10, we suggest an approach using real derivative (stock option and bond) prices to calibrate whether (and to what extent) the actions of the U.S. central bank result in such a deflator combination.
The final discussion of the model is presented in Section 11.

2. The Unified BBSM Market Model

Our unified BBSM market model A , B , C consists of a risky asset A , a riskless asset B , and a derivative C whose underlying asset is A . Unless otherwise noted, the derivative will be a European contingency claim (i.e., option), which we denote by ECC. The model assumes that A has the price dynamics of a continuous diffusion process determined by the stochastic differential equation (SDE)
d A t = a t + μ t A t d t + v t + σ t A t d B t   ϕ t d t + ψ t d B t , t 0 , A 0 > 0 ,
where B t , t [ 0 , ) , is a standard Brownian motion on a stochastic basis (filtered probability space) Ω , F , F = F t = σ B u , u t F , t 0 , P on a complete probability space Ω , F , P . In (1), a t R ,   μ t R , v t R 0 , and σ t R 0 are F -adapted processes on t 0 , satisfying the usual regularity conditions2. We require ψ t = v t + σ t A t > 0 ,   P -almost surely (a.s.). The SDE (1) is a one-dimensional linear equation whose solution can be written3
A t = A 0 η 0 , t + 0 t η s , t a s v s σ s d s + 0 t η s , t v s d B s   ,    t 0 ,
where
η s , t = exp s t μ u σ u 2 2 d u + s t σ u d B u .
From (1), it is clear4 that riskless asset B should have the price dynamics:
d β t = ρ t + r t β t d t χ t d t ,    β 0 > 0 ,
where ρ t R and r t R . Applying the Karatzas–Shreve solution to (3), we have η s , t = exp s t r u d u ; therefore,
β t = β 0 e 0 t r s d s + 0 t e s t r u d u ρ s d s .
As A t and β t must have the same numèraire, the natural choice is to set A 0 = β 0 . We will use this choice throughout the article. A necessary condition for no arbitrage is the requirement β t A t , P -a.s. for t > 0 . The no-arbitrage assumption leads to the requirement on the market price of risk,
Θ t = ϕ t χ t ψ t > 0 ,      P - a . s . ,       t   0 .
As we require ψ t > 0 , P -a.s., (5) further requires that ϕ t >   χ t , P -a.s.
With (1) expressed as d A t = ϕ ( t ,   A t ) d t + ψ ( t ,   A t ) d B t and (3) expressed as d β t = χ ( t ,   β t ) d t , we see that the BBSM model is properly viewed as a generalized Bachelier model that reduces to the BSM model when a t = v t = ρ t = 0 . When μ t = σ t = r t = 0 , the BBSM model reduces to the MB model of Rachev et al. (2024).
The ECC C has the price C t = f A t , t , where f x , t , x R , t 0 , T has continuous partial derivatives 2 f x , t / x 2 and f x , t / t on t 0 , T ; T is the terminal (maturity) time of C ; and the option’s terminal payoff is C T = g A T for some continuous function g : R R . From It ô ’s formula,
d f A t , t = f A t , t t + ϕ t f A t , t x + ψ t 2 2 2 f A t , t x 2 d t + f A t , t x ψ t d B t .

3. Option Pricing with Riskless Bonds under the BBSM Model

The classical treatment of a replicating, self-financing portfolio Π of A and B has the price process (Duffie 2001, chap. 5E, p. 89):
P t = a t A t + b t β t = C t = f A t , t ,
d P t = a t d A t + b t d β t = d C t = d f A t , t ,
where a t and b t are F -adapted processes. This portfolio treats the riskless asset as a financial product issued in units (i.e., a bond), with β t being the unit price and b t the number of units in the replicating portfolio. Solving (6), (7a) and (7b) in the usual manner produces the risk-neutral BBSM PDE for the option
f x , t t + ψ t 2 ( x , t ) 2 2 f x , t x 2 + χ t β t , t β t x f x , t x f x , t = 0 ,
for x = A t R , t 0 , T , with the boundary condition f x , T = g x .
Assuming that χ t / β t x : R × [ 0 , T ] R and ψ t : R × [ 0 , T ] R satisfy regularity conditions sufficient to guarantee that (8) has a unique strong solution, the Feynman-Kac (FK) solution of the BBSM PDE (8) is
f x , t = E t Q 1 exp t T χ u β u d u   g Z T Z t = x ,
d Z t = χ t β t   Z t   d t + ψ t   d B t Q 1 ,    t     0 ,
where B t Q 1 is a standard Brownian motion under the probability measure Q 1 .
We confirm that the probability measure Q 1 is the equivalent martingale measure (EMM) obtained from risk-neutral valuation. The FK solution (9a) suggests the deflator
D 1 t = exp 0 t χ u β u   d u D 2 t D 3 t ,     t 0 ,
D 2 t = exp 0 t r u   d u ,    D 3 t = exp 0 t ρ u β u   d u .
We define the standard Brownian motion B t Q ,   t 0 , having the dynamics d B t Q = d B t + ϕ t χ t β t A t / ψ t d t on P . Under the EMM Q ~ P , the deflated process Φ t = D 1 t   A t , t 0 , is an F -martingale with d Φ t = D 1 t ψ t d B t Q , A 0 > 0 . Consequently,
Φ t = E t Q Φ T A t = D 1 t 1 E t Q D 1 T A T .
From (10), under Q ,
d A t = χ t β t A t d t + ψ t d B t Q .
As the equivalent martingale measure Q is unique, the market ( S , B , C ) is complete. Therefore, for the option price process C t = f A t , t , t [ 0 , T ] , with terminal payoff C T = g A T , D t C t is an F -martingale under Q . Thus,
D 1 t C t = E t Q D 1 T C T      C t = D 1 t 1 E t Q D 1 T g A T .
Equations (11) and (12) obtained from the EMM Q are the same as the option price determined by the FK solution (9a), (9b) of the BBSM PDE. Thus, Q = Q 1 as claimed.
The valuation of the BBSM option price C 0 = f A 0 , 0 with terminal payoff C T = g A T is given by
C 0 = E Q 1 D 1 T g A T = E Q 1 exp 0 T χ u β u   d u g A T ,
where A T is given by (2). In Appendix A, we develop the explicit solution of the BBSM call option price for time-independent parameters.
Setting a t = v t =   ρ t = 0 (8) reduces to the BSM PDE (Duffie 2001, sect. 5G)
f x , t t + σ t 2 2 x 2 2 f x , t x 2 + r t x f x , t x f x , t = 0 ,
for   A t = x > 0 ,5  t 0 , T , with boundary condition f x , T = g x . Under this limit, (9a) and (9b) reduce to the usual FK BSM solution (Duffie 2001, chap. 5H)
f x , t = E t Q 2 e t T r u d u   g Z T Z t = x ,
d Z t = r t Z t d t + σ t Z t d B t Q 2 .
Similarly (13) reduces to the familiar valuation of the option price in the BSM model (Duffie 2001, sect. 6H)
C 0 = E Q 2 D 2 T g A T ,  
where
A t = A 0 exp 0 T r u σ u 2 2 d u + 0 T σ u d B u Q 2   .  
and the EMM Q 2 is determined by B t Q 2 , having the dynamics d B t Q 2 = d B t + μ t r t / σ t d t on P .
Setting μ t = σ t = r t = 0 , (8) reduces to the MB PDE (Rachev et al. 2024)
f x , t t + v t 2 2 2 f x , t x 2 + ρ t β t x f x , t x f x , t = 0 ,
for A t = x R , t 0 , T , with boundary condition f x , T = g x . Under this limit, (9a) and (9b) reduces to the FK MB solution
f x , t = E t Q 3 exp t T ρ u β u   d u   g Z T Z t = x ,
d Z t = ρ t β t   Z t   d t + v t   d B t Q 3 .
Similarly, (13) produces the risk-neutral valuation of the option price in the MB model
C 0 = E Q 3 D 3 T g A T ,  
with
A T =   A 0 η 0 , T + 0 T η s , T v s d B s Q 3 ,    η s , T = exp s T ρ u β u d u   ,
and the EMM Q 3 is determined by B t Q 3 , having the dynamics d B t Q 3 = d B t + a t ρ t β t A t / v t d t on P . In Appendix B, we develop the explicit solution of the FK MB call option price (15) assuming time-independent parameters.

4. Option Pricing with a Riskless Bank Account under the BBSM Model

We now consider the case in which the riskless asset is treated as a single bank account having total value β t . The replicating, self-financing portfolio can then be written:
P t = a t A t + β t = f A t , t ,
d P t = a t d A t + d β t = d f A t , t .
We note that (16a) is equivalent to a self-financing portfolio consisting of a long position in one option minus a t shares of the risky asset with per-share price A t . As this portfolio should be riskless, it must equal the value β t in the bank account (Hull 2012, p. 255). The dynamics (16b) follow similarly.
Equating (16b) to (4), and solving for a t to eliminate the Brownian motion term d W t , produces the coupled system of PDEs
f x , t t + ψ t 2 x , t 2 2 f x , t x 2 = χ t ( β t , t ) ,
  f x , t x f x , t x = β t ,
for x = A t , t 0 , T , with boundary condition f x , T = g x . For a general functional form χ t ( β t , t ) , (17a) and (17b) must be solved simultaneously for f x , t .
Under the BBSM model, χ t ( β t , t ) is linear in β t and (17a) and (17b) reduce to the single risk-neutral PDE,
f x , t t + 1 2 ψ t 2 2 f x , t x 2 + r t   x f x , t x f x , t ρ t = 0 .
When r t = 0 , (18) reduces to the PDE analyzed in Rachev et al. (2024). When ρ t = 0 , (18) is identical to (8). We emphasize this result in the following lemma.
Lemma 1.
BSM-like models with  χ t = r t β t  produce the same option prices, regardless of whether the riskless asset is packaged in units (7a), (7b) or as a bank account (16a), (16b).
The FK solution to (18) is
f x , t = E Q 2 exp t T r u d u g ( Z T ) t T exp t s r u d u ρ s d s Z t = x ,
d Z t = r t Z t d t + ψ t d B t Q 2 , t 0 .
We next confirm that Q 2 is the EMM obtained using the risk-neutral approach to option valuation using the deflator D 2 t . We define the standard Brownian motion B t Q 2 ,   t 0 , as having the dynamics d B t Q 2 = d B t + ϕ t r t A t / ψ t d t on P . Under the measure Q 2 ~ P , the deflated process Φ t = D 2 t   A t , t 0 , is an F -martingale with d Φ t = D 2 t ψ t d B t Q 2 , A 0 > 0 . Consequently,
Φ t = E t Q 2 Φ T A t = D 2 t 1 E t Q 2 D 2 T A T .
Taking the total derivative of (20), it is straightforward to show that under Q 2 , d A t is given by (19b). The EMM Q 2 is unique and the market A , B , C is complete.
From the solution (4), it is straightforward to show
β t = D 2 t 1 D 2 T β T t T D 2 s ρ s d s .
Writing the price vector as X t = A t , B t T and the continuous rate vector as D t = 0 , ρ t T , the results (20) and (21) can be summarized as
X t = D 2 t 1 E t   Q 2 D 2 T X T t T D 2 s d D s ,
in agreement with the development of arbitrage pricing with dividends (Duffie 2001, chap. 6L, pp. 123–25). Thus, the term ρ t in the riskless bank account of (16a), (16b) has the effect of adding a continuous dividend rate ρ t d t to the market A , B , C . The difference with the development in Duffie is that the dividend rate appears in the riskless asset, rather than in the risky asset. In contrast, from the solution (9a) (equivalently (12)) we see that the term ρ t appearing in the riskless bond price of the replicating portfolio (7a), (7b) does not result in the addition of a continuous dividend rate to the market A , B , C .
As the market is complete, the option price process C t = f A t , t , t [ 0 , T ] , with terminal payoff C T = g A T , D t C t must also be an F -martingale under Q 2 . From (19a),
C t = D 2 t 1 E t   Q 2 D 2 T g A T t T D 2 s ρ s d s .
Regardless of the asset generating the dividend rate, the option price (23) reflects that rate. While the development of arbitrage pricing with dividends in Duffie (2001, chap. 6L) is silent regarding the form of the self-financing trading strategy, we have identified an appropriate replicating, self-financing trading strategy in (16a), (16b) for the price processes (1), (3) and (6) with ρ t 0 . Duffie’s development is also silent on the dynamics of the price processes for the risky and riskless asset. Equations (16a) and (16b) constitute the only possible trading strategy with prices (1) and (3), producing the result (19a), (19b).
The solutions (22) and (23) can be interpreted as the results of a process Φ t = D m t   X t + D a t   obtained via multiplicative D m t   and additive D a t   deflators. We explore this idea further in Section 10.

5. Option Pricing When the Risky Asset Pays a Dividend Yield

We develop option pricing in the BBSM market ( A D , B , C D ) when the asset A D pays a dividend at a continuous rate with an instantaneous dividend yield D t = D ( A t , t ) R .6 By the usual no-arbitrage argument, the price A t D of A D obeys
d A t D = a t + μ t D t A t D d t + v t + σ t A t D d B t ϕ t D d t + ψ t D d B t ,   A 0 D > 0 .
Equation (24) describes the dynamics of the underlying asset from the point of view of the buyer of the ECC. The riskless bank account B has the price dynamics (3). The ECC C D has the price C t D = f A t D , t , t [ 0 , T ] , where f x , t , x R , t 0 , T has continuous partial derivatives 2 f x , t / x 2 and f x , t / t on t 0 , T . The terminal time of C D is T with payoff C T ( D ) = g ( A T D ) . Itô’s formula for the option price dynamics is
d C t D = f A t D , t t + ϕ t D f A t D , t x + 1 2 ψ t D 2 2 f A t D , t x 2 d t + ψ t D f A t D , t x d B t .
We consider the replicating, self-financing riskless portfolio of A D and B having the price
P t D = a t A t D + b t β t = C t D = f A t D , t .
In the case of dividend payouts, as the holder of the portfolio (option seller) receives the dividend, the change in value of the portfolio is
d P t D = a t d A t D + D t A t D d t + b t d β t = d C t D .
The standard computation yields the risk-neutral BBSM PDE for the option price
f x , t t + χ t β t D t x f x , t x + ψ t 2 ( x , t ) 2 2 f x , t x 2 χ t β t f x , t = 0 ,
for x = A t D R , t 0 , T , with boundary condition f x , T = g x . The FK solution of (25) is
f x , t = D 1 t 1 E Q 3 D 1 T g Z T Z t = x ,
d Z t = r t D t + ρ t β t Z t d t + v t + σ t Z t d B t Q 3 ,    t     0 .
Setting a t = v t = ρ t = 0 in (25) produces the familiar BSM PDE when the underlying asset pays out a dividend yield D t (Björk 2009, p. 162),
f x , t t + r t D t x f x , t x + σ t 2 x 2 2 2 f x , t x 2   r t f x , t = 0 ,
x > 0 , t 0 , T , with boundary condition f x , T = g x . In this limit, (26a), (26b) results in the familiar FK solution (Duffie 2001, chap. 5H) of the BSM PDE (27)
f ( B S M ) x , t = D 2 t 1 E Q 3 D 2 T g Z t     Z t = x ,
d Z t = r t D t Z t d t + σ t Z t d B t Q 3 ,    t     0 .
Setting μ t = σ t = r t = 0 in (25) produces the MB PDE,
f x , t t + ρ t β t D t x f x , t x + v t 2 2 2 f x , t x 2 ρ t β t f x , t = 0 ,
x R , t 0 , T , with boundary condition f x , T = g x . In this limit, (26a), (26b) produces the solution of the MB PDE (28),
f ( B S M ) x , t = D 3 t 1   E Q 3 D 3 T g Z t     Z t = x ,  
d Z t = ρ t β t D t Z t d t + v t d B t Q 3 ,    t     0 ,
where
D 3 t = e x p 0 t ρ u β u d u = D 1 t D 2 t .
We note that the risk-neutral measures Q 3 in the BBSM solution (26a) and (26b), and in the BSM and MB limits, will differ.
If we instead consider a replicating portfolio consisting of a bank account,
a t A t D + β t = f A t D , t ,
a t d A t D + D t A t D d t + d β t = d f A t D , t ,
the FK solution for the option price is
f x , t = D 2 t 1 E Q 4 D 2 T g Z T D t T D 2 s ρ s d s Z t D = x ,
d Z t D = r t D t Z t D d t + v t + σ t Z t D d B t Q 4 ,    t     0 .

6. A Perpetual Derivative Price Process

Shirvani et al. (2020) introduced a new class of hedging instruments, specifically perpetual derivatives (i.e., options with perpetual maturities). We consider a classical BSM market having a risky asset A B S M and riskless asset B B S M whose price process are given, respectively, by (1) with a t = v t = 0 , μ t = μ , σ t = σ > 0 , and by (3) with ρ t = 0 , r t = r . For δ = 2 r / σ 2 R , V t = A t δ , t 0 , represents the price of a perpetual derivative, V B S M , in the complete no-arbitrage market A B S M , B B S M , V B S M . There are two major applications of V B S M . The first is to use V B S M when the trader does not have access to the riskless bank account B B S M . The trader can then form a self-financing portfolio of A B S M and V B S M to replicate B B S M . The second is the use of V B S M in forming a hedging portfolio in trinomial option pricing. In trinomial option pricing, the hedger (the trader taking a short position in the option contract) requires three assets to form a self-financing portfolio that replicates the value of an option (having A B S M as underlying) on each node of the tree.
We investigate the development of a perpetual derivative in the more general case of the BBSM market model. The perpetual derivative for the BSM model was studied in Lindquist and Rachev (2024) by considering separable solutions of the form
f x , t = g x h t + w t
for the risk-neutral PDE. Using (29), (8) becomes
g x h t + χ t β t x g x g x + ψ 2 2 x , t g x h t = w t + χ t β t w t .
The right-hand side is deterministic, whereas the left-hand side contains the stochastic process x = A t . Thus, both sides must separately vanish, giving the ordinary differential equations
w t χ t β t w t = 0 ,
g x h t + χ t β t x g x g x + 1 2 ψ 2 x , t g x h t = 0 .
The solution to (30a) is
w t = D 1 t 1 w 0 .
Equation (30b) suggests solutions of the form g x = x γ . Substituting into (30b) produces
h t = η x , t h t ,    η x , t = 1 γ χ t β t + 1 2 ψ 2 x , t γ x 2 .
Equation (32) admits a solution as long as η x , t = η t , which is possible as long as v t = 0 , in which case
η t = 1 γ 2 χ t β t σ t 2 + γ σ t 2 2 ξ γ σ t 2 / 2 .
For the BSM model ( a t = v t = ρ t = 0 ), the behavior of the function ξ γ for time-independent coefficients r t = r , σ t = σ was studied by Lindquist and Rachev (2024). In that case, ξ γ = 1 γ δ + γ with δ 2 r / σ 2 . This produced a one-parameter family of solutions to (32), one of which is the parameter value γ = δ = 2 r / σ 2 , as previously used by Shirvani et al. (2020).
For the BBSM model, as long as v t = 0 , (8) admits a perpetual derivative of the form (29), where w ( t ) is given by (31) and h t is given by
h t = h 0 e x p 1 δ 0 t χ s β s + δ σ s 2 2 d s .
Under time-independent coefficients r t = r , σ t = σ , and ρ t = ρ , an examination of ξ γ in (33) shows it has the form ξ γ = 1 γ d ( t ) + γ . Thus, at fixed time value t , the analysis of the form of ξ γ in Lindquist and Rachev (2024) holds with the replacement δ d ( t ) . For example, if r > 0 , then d ( t ) changes monotonically from d 0 = δ + 2 ρ / β 0 σ 2 to lim t d t = δ . We summarize the principal result of this section in the following lemma.
Lemma 2.
The requirement  v t = 0  is a necessary condition for the BBSM to support a perpetual derivative having a solution of the form (29). As a corollary to this, the MB model cannot admit such a perpetual derivative.

7. Discrete Binomial Option Pricing

We consider binomial option pricing in the BBSM market model A , B , C . Under continuous-time option pricing models, the hedger (who takes a short position in C ) is assumed to be able to trade continuously. This trading fiction leads to the option price puzzle (Hu et al. 2020a) under which, as observed in (13) and (14), the option price has no dependence on the natural-world drift terms of the risky asset. Under a discrete-time binomial model, trading must occur at discrete time steps.
We start with the simplest form of binomial option pricing. The price of the risky asset S follows the binomial pricing tree
A k + 1 , n = A k + 1 , n u = A k , n + u k , n ,   i f   ξ k + 1 , n = 1 , A k + 1 , n d = A k , n + d k , n ,   i f   ξ k + 1 , n = 0 .
In (34):
(i)
A k , n , k = 0,1 , , n , n N = 1 , 2 , , is the stock price at time k , k = 0,1 , , n , = n = T / n , where T is the fixed terminal time, A 0 > 0 ; and
(ii)
for every n N , the ξ k , n   k = 1,2 , , n , are independent, identically distributed Bernoulli random variables with P ξ k , n = 1 = 1 P ξ k , n = 1 = p n determining the filtration F n = F k n = σ ξ j , n ,   j = 1 , , k ,   k = 1 , , n ,   F 0 n = , Ω ,   ξ 0 , n = 0 , and the stochastic basis Ω , F , F n , P on the complete probability space Ω , F , P .
The riskless asset B has the discrete price dynamics
β k + 1 , n = β k , n + χ k , n ,    k = 0,1 , , n 1 ,     β 0 , n > 0 ,
where χ k , n R is an instantaneous simple riskless rate. We again invoke a common numèraire via β 0 , n = A 0 , n . Under the BBSM model, stock price changes (i.e., lag-1 differences) rather than normalized returns are of primary interest. The price changes
c k + 1 , n = A k + 1 , n A k , n ,    k = 0 , , n 1 ,      c 0 , n = 0 ,
have the discrete dynamics
c k + 1 , n = c k + 1 , n u = u k , n , w . p . p n , c k + 1 , n d = d k , n , w . p . 1 p n .
We would like to have weak convergence of the càdlàg process7 generated by the binomial tree (34) in the Skorokhod space D 0 , T to the continuous time process (1).8 To achieve this, we must consider an instantaneous mean ϕ k , n and instantaneous variance ψ k , n 2 for the discrete price differences of the risky asset such that the conditional mean and conditional variance satisfy
E k , n c k + 1 , n E c k + 1 , n F k n = ϕ k , n ,
V a r k , n c k + 1 , n V a r c k + 1 , n F k n = ψ k , n 2 .
As a consequence, u k , n and d k , n are determined by
u k , n = ϕ k , n + 1 p n p n   ψ k , n ,       d k , n = ϕ k , n p n 1 p n   ψ k , n .
The ECC C will have the discrete price dynamics C k , n = f A k , n , k Δ , k = 0 , , n 1 , with terminal value C n , n = g A n , n . Consider a self-financing strategy, P k , n = a k , n A k , n + b k , n β k , n , replicating the option price C k , n :
a k , n A k , n + b k , n β k , n = C k , n ,
a k , n A k + 1 , n u + b k , n β ( k + 1 ) , n = C k + 1 , n u ,
a k , n A k + 1 , n d + b k , n β ( k + 1 ) , n = C k + 1 , n d .
Solving this system in the usual manner produces the recursion formula for the risk neutral valuation for the option price on the tree,
C k , n = 1 D   q k , n C k + 1 , n u + ( 1 q k , n ) C k + 1 , n d ,
where D = β ( k + 1 ) , n / β k , n = 1 + χ k , n / β k , n , and the risk-neutral probability is
q k , n = A k , n β k , n χ k , n d k , n u k , n d k , n .
The ratio A k , n / β k , n appearing in (38) reinforces our choice of β 0 , n = A 0 , n as the natural common numèraire. Enforcement of the “time value of money” requires D > 0 for all k = 0 ,   ,   N . To ensure no opportunity for arbitrage, we require u k , n > A k , n β k , n χ k , n >   d k , n for all k = 0 ,   ,   N , ensuring that q k , n ( 0,1 ) .
Using (36), the risk-neutral probability (38) can be written
q k , n = p n φ k , n χ k , n A k , n β k , n ψ k , n p n 1 p n   .
The natural world probability satisfies p n 0,1 ; however, as noted above, the requirement q k , n 0,1 must be imposed.
The discrete version of (1) and (3) implies φ k , n = a k , n + μ k , n A k , n , ψ k , n = v k , n + σ k , n A k , n and χ k , n = ρ k , n + r k , n β k , n , with the requirement that ψ k , n > 0 and χ k , n < φ k , n . Thus, the BBSM risk-neutral probability (39) can be written
q k , n = p n a k , n + μ k , n A k , n ρ k , n + r k , n β k , n A k , n β k , n v k , n + σ k , n A k , n p n 1 p n   .
From (35) and (3),
D = β ( k + 1 ) , n β k , n = 1 + r k , n + ρ k , n β k , n ,      β 0 , n > 0 .
As 0 ,
1 D ~ e x p r k , n + ρ k , n β k , n + o .
We note that 1 / D (and hence D ) is positive for sufficiently small values of . In general, as noted above, the condition D > 0 must be imposed. Arguments for the weak convergence of binomial-tree-generated càdlàg processes in the Skorokhod space D 0 , T  show (see Hu et al. 2022) that the discrete option price (37) converges weakly to (12).
The choice a k , n = v k , n = ρ k , n = 0 produces the risk-neutral valuation of the option (37) on a binomial tree in the BSM market model,
C k , n ( B S M ) = 1 1 + r k , n q k , n C k + 1 , n u + ( 1 q k , n ) C k + 1 , n d ,
where the risk-neutral probability q k , n is given by9
q k , n ( B S M ) = p n θ k , n p n 1 p n   ,      θ k , n       μ k , n r k , n σ k , n .
Here, θ k , n is the discrete form of the BSM market price of risk θ t = μ t r t / σ t , t 0 .
The choice μ k , n = σ k , n = r k , n = 0 produces the risk-neutral valuation of the option (37) on a binomial tree in the MB market model,
C k , n ( M B ) = 1 1 + ρ k , n / β k , n q k , n C k + 1 , n u + ( 1 q k , n ) C k + 1 , n d ,
with q k , n given by
q k , n ( M B ) = p n a k , n ρ k , n A k , n β k , n v k , n p n 1 p n   ,
where β k , n = β k 1 , n + ρ k 1 , n = β 0 + j = 1 k ρ j , n . For time-independent parameters, (42) becomes
q k , n ( M B ) = p n a ρ   A k , n β 0 + ρ k ν p n 1 p n   .
Under time-independent parameter values, β k , n = β 0 + ρ k , with β 0 = A 0 , and
D ( M B ) = 1 + ρ k , n β k , n = 1 + β 0 ρ + k = 1 + ρ β 0 1 1 + ρ k β 0 .
For values ρ k / β 0 1 (i.e., k β 0 / ρ ),
D M B ~   1 + ρ β 0 + O k 2 ,
which resembles the discounting in the BSM option pricing (40). However, for ρ k / β 0 1 (i.e., k approaching the maturity time T ), the approximation (44) is no longer valid. Thus, in MB binomial option pricing, one must use (41) or its time-independent parameter version (43).

8. The Term Structure of Interest Rates

We consider the price structure of a zero-coupon bond, a forward contract, and a futures contract under the BBSM market model. The prices require knowledge of the dynamics of the interest rates r t and ρ t . A complete description of a model for the term structure of interest rates (TSIR) requires a relationship between interest rates and maturity terms. There are numerous models that have been developed for interest rates (see, for example, Shreve (2004a, chap. 6) or Duffie (2001, chap. 7)); any of these models can be considered for r t and ρ t . However, we raise a significant issue regarding these models, each of which requires knowledge of the risk-neutral probability measure Q or, equivalently, the probability p t ( Q ) of an upward risk-neutral price movement. Almost always, the implementation of these models involves an ad hoc choice. As examples, both Ho–Lee and Black–Derman–Toy models use p t ( Q ) = 1 / 2 (Shreve 2004a, chap. 6). In our view, the measure Q and the characterization of p t ( Q ) must come from the appropriate market. As demonstrated in Hu et al. (2020b, Figure 3), the value of p t ( Q ) is often far from 1/2. We suggest that fitting an interest rate model to option prices based on a broad-based ETF as underlying, for example, the SPDR S&P 500 ETF Trust (SPY), is the appropriate way to develop values for p t ( Q ) . We leave this line of investigation for future work and concentrate here on developing the indicated derivative pricing under the BBSM model.

8.1. Pricing a Zero-Coupon Bond

We consider the BBSM market A , B , Z T , where the dynamics of the price A t ,   t 0 , of the risky asset is given by (1) and the dynamics of the price β t , t 0 , of the riskless asset is given by (3). Here, Z T is a zero-coupon bond (see Duffie (2001, sect. 6M) and Shreve (2004b, sect. 5.6.1)), having maturity time T > 0 and price B ( t ,   T ) . The payoff g B t ,   T for Z T is 1 if t = T and 0 for t < T . Under the BBSM model, B ( t ,   T ) is given by (12)
B t ,   T = D 1 t 1 E t   Q 1 D 1 T   D 1 t 1 M t , 0 t T T .
In (45), T represents a fixed time horizon at, or prior to which, all bonds will mature. The family B t ,   T , 0 t T T constitutes the TSIR in the BBSM model.10
Setting ρ t = 0 in (45) produces the TSIR under the BSM model, B ( B S M ) t ,   T = D 2 t 1 E t   Q 2 D 2 T . Setting r t = 0 in (45), we obtain the TSIR under the MB model, B ( M B ) t ,   T = D 3 t 1 E t   Q 3 D 3 T , where β t = β 0 + 0 t ρ u d u , t 0 , β 0 > 0 .
In (45) and its BSM and MB limit cases, the expected value is taken with respect to the risk-neutral measure Q 1 . Evaluation of (45) simplifies considerably under time-independent parameters to
B t ,   T = E t   Q 1 exp t T r + ρ β u d u , 0 t T T ,
where β t is given by (A3) in Appendix A. Using the solution to (12) under time-independent parameters developed in Appendix A, we can calibrate the parameters r and ρ from market option prices. Using the calibrated values for r , we can compute the theoretical value of zero-coupon bonds given by (46). Those values indicate the views of the option traders on the values of the zero-coupon bonds and thus on the TSIR. A large difference between the resultant theoretical values of the zero-coupon bonds and their market value may indicate a dislocation between the option and bond markets.
The process M t in (45) is a martingale on Ω , F , F , Q 1 . By the martingale representation theorem (Duffie 2001, p. 336),
M t = M 0 + 0 t θ B u d B u Q 1 = E Q D 1 T + 0 t θ B u d B u Q 1 ,   0 t T ,
for some F -adapted θ B t in L 2 . Thus, as B t , T = D 1 t 1 M t , then
d B t ,   T = χ t β t D 1 t 1 E Q 1 D 1 T + 0 t θ B u d B u Q 1   d t + D 1 t 1 θ B t d B t Q 1
= χ t β t B t ,   T d t + D 1 t 1 θ B t d B t Q 1 .
We consider the replicating self-financing portfolio having the price dynamics:
P t = a t A t + b t β t = B t ,   T ,
d P t = a t d A t + b t d β t = d B t ,   T , 0 t T .
Solving for a t and b t in the usual manner,
a t = D 1 t 1 θ B t ψ t , b t = 1 β t B t ,   T D 1 t 1 θ B t ψ t A t ,
produces the replicating, self-financing strategy for the zero-coupon bond in the BBSM market model. Setting ρ t = a t = v t = 0 or r t = μ t = σ t = 0 produces, respectively, the replicating, self-financing strategy for the zero-coupon bond under the BSM or MB market models. If r t and ρ t , t 0 , are deterministic, then D 1 t 1 is deterministic and γ B t = 0 , t 0 , in the BBSM framework (and its BSM and MB limits).

8.2. Pricing a Forward Contract

We define a forward contract in the BBSM market in a classical sense. At time t , two market agents ( and ) enter into a forward contract on the risky asset A , with taking the long position and taking the short position. When the contract reaches maturity, ’s payoff is V l T ,   T = A T F ( t , T ) and ’s payoff is V s T ,   T = V l T ,   T ), where F(t,T) is the T -forward (delivery) price of asset A evaluated under the BBSM model.
With V t ,   T denoting the initial value for the forward contract, from (12),
V t ,   T = D 1 t 1   E t Q D 1 T A T F t , T = A t F t , T B t ,   T , t 0 , T .
Thus, the T -forward price of the asset A is11
F t , T ; V t ,   T = A t V t ,   T B t ,   T , 0 t T .
The traditional definition of a forward contract under the BSM model assumes V t ,   T = 0 . A more realistic scenario involves an initial forward contract price V ( t ,   T ) V t   0 . Then, for t s T the value (47) of the BBSM forward contract is
V s ,   T ; V t = A s F t , T ; V t B s ,   T .
We conclude with the hedging strategy that should follow when V t   0 . agrees to pay F ( t , T ;   V t ) for the asset with value A T at T . At initiation time t , receives V t . As per (47), receives the precise amount to buy the asset and takes a short position in F t , T ; V t zero-coupon bonds. From (48), at terminal time T , ’s position has the value A T F t , T ; V t B T ,   T = A T F t , T ; V t ; delivers the asset and receives the forward price F t , T ; V t .

8.3. Pricing a Futures Contract

We view a futures contract in the BBSM market model as a “mark to market” sequence of forward contracts. Let t i , n , i = 0 , ,   n , 0 = t 0 , n < t 1 , n < < t n , n = T , n N = 1,2 ,     be a sequence of trading instances (of the forward contracts). To maintain the futures contract, trader should pay φ t i , n t i + 1 , n t i , n 0 . The classical setting of futures contracts assumes that φ t i , n n = 0 , allowing to take on an unlimited number of futures contracts12. In what follows, we follow the exposition in Shreve (2004b, sect. 5.6.2).13
For notational convenience, we write the accumulation factor φ 0 , t   D 1 t 1 as
φ 0 , t = exp 0 t r u d u ,   r t = χ t β t ,   t 0 .
The instantaneous rate r t R , t 0 , is an F -adapted process on Ω , F , F , P . We assume the regularity condition
s u p r t + 1 r t ;   t 0 , T 0 , , P - a . s .
To simplify the connection between the discrete- and continuous-time futures prices, we further assume that r t , t 0 , is continuous on 0 , T , P -a.s. Let Δ n = max i = 0 ,   , n 1 t i + 1 , n t i , n , and assume that Δ n 0 as n . We define the piecewise constant rate, r t n = r t i , n for t t i , n , t i + 1 , n , i = 0 , ,   n 1 , with r T n = r T . Using r t n , we define the discrete accumulation factor
φ 0 , t n = exp 0 t r u n d u , t 0 .
Then
φ 0 , t k + 1 , n n = exp 0 t k + 1 , n r u n d u = exp   i = 0 k r t i , n t i + 1 , n t i , n   .
Trader enters a long futures contract at time t k , n , k = 0 , , n 1 , when the market price of the futures contract is Φ t k , n n . At time t k + 1 , n , when the market price is Φ t k + 1 , n n , makes a profit or loss of Φ t k + 1 , n n Φ t k , n n R . From (12), the value at time t = 0 of the payment Φ t k + 1 , n n Φ t k , n n is
E Q 1 φ 0 , t k + 1 , n n Φ t k + 1 , n n Φ t k , n n = E Q exp i = 0 k r t i , n t i + 1 , n t i , n Φ t k + 1 , n n Φ t k , n n .
As it costs nothing to enter the futures contract at time t = 0 ,
E Q i = 0 n 1 1 φ 0 , t k + 1 , n n Φ t k + 1 , n n Φ t k , n n = E Q i = 0 n 1 exp i = 0 k r t i , n t i + 1 , n t i , n Φ t k + 1 , n n Φ t k , n n = 0 ,   Q - a . s .
Informally, passing to the limit as n , (49) becomes
E Q 0 T 1 φ 0 , u d Φ u = E Q 0 T e 0 u r s d s d Φ u = 0 ,   Q - a . s .
for some F -adapted process Φ t , t [ 0 , T ] , on the stochastic basis Ω , F , F , Q . The same argument yields
E t Q t T 1 φ 0 , u d Φ u = E t Q t T e 0 u r s d s d Φ u = 0 , Q - a . s . f o r   a l l   t 0 , T .
As Q ~ P , the final equality in each of (49), (50) and (51) also holds P - a . s . The T -futures price of asset A in the BBSM market is any F -adapted stochastic process Φ t , t 0 , T , satisfying (51) with Φ T = A T , Q   -a.s. (and thus P -a.s.). As shown in Shreve (2004b, sect. 5.6.2), the unique F -adapted process Φ t , t 0 , T , that satisfies (51) with Φ T = A T is
Φ t = E t Q A T .
It is of interest to employ futures market data14 to calibrate time-independent parameter values r t = r and ρ t = ρ . By calibrating r and ρ from option contracts (European puts and calls), we can investigate whether there is market dislocation between the futures and option markets.

9. Calibration of Model Parameters

The unified BBSM model has a rich parameter space ( a t ,   μ t , v t ,   σ t ,   ρ t ,   r t ). We have concentrated on developing derivative pricing solutions for three cases: (i) all parameters are non-zero, (ii) the BSM limit in which a t = v t = ρ t = 0 , and (iii) the MB limit μ t = σ t = r t = 0 . Note, however, that there are parameter choices, such as μ t = σ t =   ρ t = 0 , which result in a risky asset whose price follows arithmetic Brownian motion while the riskless asset price is continuously compounded (geometric). This parameter choice produces a model in which the dynamics of the risky and riskless assets are incompatible with each other.
We address how the model parameters can be calibrated to empirical data. Consider, for example, the calibration of ϕ t = a + μ A t to the historical price series of asset A . Let A k be the price at the time k . Let ϕ k denote the mean value of A ’s price over the time period [ k τ + 1 ,   k ] .15 Repeating this process for times k , k 1 ,   k 2 ,   . ,   k n + 1 will generate data for a robust regression model to fit the series
ϕ t k = a + μ A t k ,    k = 0 ,   ,   n 1 ,  
for the parameters a and μ . A similar procedure can be used to calibrate ψ k 2 = v + σ A k 2 , where ψ k 2 is estimated from the variance of the data over the time period [ k τ + 1 ,   k ] . Again, robust regression should be used to fit the series
ψ t k = v + σ A t k ,    k = 0 ,   ,   n 1 ,
but now subject to the restriction ψ t k 0 , k = 0 ,   ,   n 1 .
The calibration of χ t = ρ + r β t must be conducted in a risk-neutral setting, i.e., from option prices. Consider (under time-independent parameters) the recursion formula for the risk-neutral valuation for the option price on the binomial tree (37), which employs risk-neutral probability (39). The parameters a ,   μ ,   v and σ are estimated from spot prices as described above. The probability of upward price movement p n can be directly estimated from price changes A t k A t k 1 , k = 0 ,   ,   n 2 . The parameters ρ and r can be estimated by minimizing
min ρ , r C t T i , K j ; ρ , r C t A , e m p ( T i , K j ) C t A , e m p ( T i , K j ) , i = 1 ,   ,   n T , j = 1 ,   ,   m T i ,
where for time t , C t A , e m p ( T i , K j ) is the listed option price corresponding to the maturity time T i and strike price K j ; C t T i , K j ; ρ , r is the theoretical price computed from (37); n T is the number of listed maturity times; and for each maturity time T i , m T i is the number of listed strike prices. Note that this minimization must be performed subject to the restriction χ t = ρ + r B t < ϕ t = a + μ A t .

10. Option Pricing under Combined Multiplicative and Additive Deflators

Consider the position that “the U.S. central bank’s view of the market is the closest to the risk-neutral world of pricing securities”.16 We assume that the risky asset price dynamics are given by (1) and that the dynamics of the central bank’s riskless assets (U.S. treasuries) are unknown outside of the central bank. Under these assumptions, we propose to find an EMM Q 4 ~ P such that the deflated process
Φ t = a D m t A t   + b D a t ,    t 0 ,
is an F -martingale on Q 4 . In (52), a R and b R are constants. We assume D m t and D a t are time-dependent, F -adapted Itô processes. Further, we assume D m t is strictly positive (Duffie 2001, chap. 6B) and invertible. Denoting the derivatives D m t = d D m t / d t and D a t = d D a t / d t ,
d Φ t = a D m t A t + a D m t ϕ t + b D a t d t + a D m t ψ t d B t μ ϕ d t + σ ϕ d B t .
Define the standard Brownian motion B t Q 4 ,   t 0 , with dynamics given by d B t Q 4 = d B t + μ ϕ / σ ϕ d t on P . Under the EMM Q 4 , the deflated process (52) satisfies d Φ t = a D m t ψ t d B t Q 4 . Consequently,
Φ t = E t   Q 4 Φ T A t = D m t 1 E t   Q 4 D m T A T + b a D a T D a t ,
d A t = D m t 1 D m t A t b a D m t 1 D a t d t + ψ t d B t Q 4 ,
where D m t 1 d D m t 1 / d t . As there is a single central bank, the market is complete and an option having the risky asset as underlying with payoff g ( A T ) at maturity date T will have the price
C t = D m t 1 E t   Q 4 D m T g ( A T ) + b a D a T D a t .
From these results, we see that only the ratio c = b / a enters the pricing. Thus, (52) can be recast as Φ t = D m t A t   + c D a t .
If we assume that the effect of the central bank’s actions on setting U.S. Treasury rates produces the riskless rate dynamics d β t = c ρ t + r t β t d t , then in (53),
D m T = exp 0 t r u   d u ,    D a t = 0 t D m s ρ s d s   ,
and the option price model (23) is identical to (53) under the change ρ t c ρ t .
Support for a multiplicative deflator, under which the transition from the natural world to the risk-neutral world involves a multiplicative change of numèraire, has strong theoretical underpinnings. The existence of a multiplicative deflator is equivalent to the no-arbitrage condition (Duffie 2001, p. 4). Under the fundamental theorem of asset pricing, no-arbitrage is equivalent to the existence of an equivalent martingale measure (Delbaen and Schachermayer 1994). As noted by Duffie (2001, chap. 2C), no arbitrage requires positive asset prices.
The existence of additive deflation has not been subjected to theoretical scrutiny. We make the following informal observations. Under the MB limit of the unified BBSM model (and possibly under the unified model as long as ρ t 0 ), there is no limit on how negative prices can go, and consequently, there is no relationship between the existence of a risk-neutral measure and no arbitrage. In fact, this observation suggests the argument that the arbitrage opportunity should play no role in the existence of a risk-neutral measure. It is generally accepted in real financial markets that the risk-neutral world is the world of the central bank (R. Roll, private communication). Decisions made by the U.S. Federal Open Market Committee effectively determine the deflator model for the U.S. economy. It would appear that committee members pay no attention to the possibility of arbitrage in making their decisions. Rather, the enforcement of no arbitrage in the U.S. financial system appears to be a by-product of the efficiency of its markets.
Is it possible, therefore, that the connection between no-arbitrage and risk-neutral measures is a spurious consequence of a strict focus on the BSM model? The forms (1), (3) and (6) of the asset price dynamics under the unified BBSM model suggest a continuous diffusion process with time-varying parameters.17 If we accept the idea that asset pricing should be considered fair, under this general model, there exists the possibility of an equivalent martingale measure under a linear transformation c D a t + D m t A t of prices, containing both drift (additive) and scale (multiplicative) deflation terms.
Testing such an assumption on a large ensemble of empirical option and bond prices would determine which of the models (i) D a t = 0 , D m t 0 ; (ii) D a t 0 , D m t = 0 ; or (iii) D a t 0 , D m t 0 provides a better fit. Assuming time-independent parameters for simplicity, the BBSM model, with the addition of c a / b in (52), provides seven fitting parameters under which to determine the admixture of BSM and MB models and additive and multiplicative deflators. Under the ESG-adjusted valuation discussed in the Introduction, the parameter set is further augmented by γ E S G .

11. Discussion

We have developed the unified BBSM model in response to 21st-century finance which, under new financial crises, has experienced negative commodity prices and interest rates. As noted in the Introduction, we are also motivated by a desire to embed ESG ratings (a quantifier of corporate sustainability action) as a third dimension to the usual financial risk, financial reward worldview of finance. This contrasts with proposed models that treat ESG ratings as yet another financial risk factor. A consequence of our approach is negative ESG-adjusted valuations, which require a Bachelier-type model. We have effected this through a larger unified BBSM model, which contains the BSM model as well as a “modernized” Bachelier model as particular parameter limits. As discussed in the Introduction, our MB model corrects significant deficiencies of existing Bachelier models.
Current models for the TSIR seek the EMM only within the market of riskless bonds and their derivatives. (In particular, as demonstrated by Hu et al. (2020b, Figure 3), the use of q = 1 / 2 for the EMM in the binomial Ho–Lee and Black–Derman–Toy TSIR models is not acceptable). As shown in Section 8.1, the appropriate EMM used in determining the value of the zero-coupon bond is fixed by the stock option market. Thus, our BBSM model correctly notes that EMMs for bonds and stock options must be the same; otherwise, there is an effective assumption of two central banks, one governing risk-neutral valuation of equities and the other governing risk-neutral valuation of default-free assets.
Our study of the BBSM shows that certain results of the traditional BSM model are highly constrained.
  • It is well known that the BSM limit restricts prices and interest rates to be positive. We have not fully explored the BBSM model to determine in what parameter neighborhood of the BSM model this positivity holds; it does not hold for the full parameter space.18
  • While the BSM model produces the same option pricing results under the two common variations of the replicating, self-financing portfolio (Section 3 and Section 4), in general, the two portfolios produce different option prices. In particular, option pricing under the portfolio of Section 4 can be interpreted in terms of a combination of multiplicative and additive deflators. In Section 10, we explored this possibility in greater generality.
  • While the BSM model supports a perpetual derivative solution (of separable form), a large subset of the BBSM model parameter space, which includes the MB model, does not.
It is critical to note that our unified BBSM market model is complete, resulting in a unique EMM.19 While it is well intentioned to address the stylized facts of asset prices using Lévy processes to model the underlying price dynamics, such attempts result in uncountably many equivalent martingale measures, making the choice amongst them arbitrary. While the semimartingale approach to asset pricing based on the Strasbourg school of stochastic processes is beautiful, it has a fundamental drawback: all interesting non-Gaussian models exhibit infinitely many EMMs. Consequently, we submit that
Any meaningful market model should be complete, and arbitrage free, guaranteed by the existence of a unique EMM.
In our view, this completeness is cemented by fact that the U.S. Federal Reserve System (the “Fed”) is dominant in the global financial system. (i.e., there is effectively only one global central bank). In support of this statement, we note the following points that speak to the global impact of the Fed’s policies and actions: (i) the status of the U.S. dollar as the global reserve currency; (ii) the size and influence of the U.S. economy; (iii) the Fed’s role as a lender of last resort; (iv) the Fed’s reputation for independence and expertise; and (v) the network effect of the dominance of the U.S. dollar.

Author Contributions

Conceptualization, S.T.R.; methodology, W.B.L. and S.T.R.; validation, W.B.L., J.G. and S.T.R.; formal analysis, S.T.R. and W.B.L.; writing—original draft preparation, S.T.R. and W.B.L.; writing—review and editing, W.B.L. and F.J.F.; supervision, S.T.R. and F.J.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A. Explicit Solution of the BBSM Call Option Price for Time-Independent Parameters

We examine the explicit option solution to the full BBSM model for time-independent parameters:
a t = a R , μ t = μ R ,   v t = v R 0 ,   σ t = σ R 0 ,   ρ t = ρ R ,   r t = r R , with   v + σ A t > 0   P - a . s .   and   a + μ A t ρ + r β t > 0   P - a . s .
From (12), the option price at t = 0 is
C 0 = e r T E Q 1 g A T exp ρ 0 T 1 β u d u ,
where, under the riskless measure Q 1 , A T is given by (14). With time-independent coefficients,
A T = A 0 η 0 , T v σ 0 T η s , T d s + v 0 T η s , T d B s Q 1   ,
with
η s , T = exp r σ 2 2 T s + ρ s T 1 β u d u + σ B T Q 1 B s Q 1   .
From (4), the riskless asset price is
β t = β 0 + ρ r e r t ρ r i f r 0 , β 0 + ρ t i f r = 0 .  
Again, we assume the natural numéraire β 0 = A 0 .
When r 0 , then
ρ s T 1 β u d u = ln A 0 + ρ / r e r T ρ / r A 0 + ρ / r e r s ρ / r r T s ,
as long as ρ 0 and A 0 + ρ / r e r u > ρ / r for u [ s , T ] .20 Then,
η s , T = A 0 + ρ / r e rT ρ / r A 0 + ρ / r e rs ρ / r exp σ 2 2 T s + σ B T Q 1 B s Q 1 .
Thus, from (A1), the r 0 BBSM risk-neutral price of an option is given by
C 0 = A 0 A 0 + ρ / r e r T ρ / r   E Q 1 g A T ,
where A T is given by (A2), with η s , T given by (A4).
When a = v = ρ = 0 , (A5) reduces to the standard BSM risk-neutral option price
C 0 = e r T E Q 2 g A T ,
with
A T = A 0 exp r σ 2 2 T + σ B T Q 2 .
When μ = σ = r = 0 (and using the normalization β 0 = A 0 ),
ρ s T 1 β u d u = ρ s T 1 β 0 + ρ u d u = ln A 0 + ρ T A 0 + ρ s .
Then,
C 0 = A 0 A 0 + ρ T E Q 3 g A T ,
with
A T =   A 0 + ρ T 1 + v 0 T 1 A 0 + ρ s d B s Q 3 ,
in agreement with the MB FK solution (A7), (A8).

Appendix B. Explicit Solution of the MB Call Option Price for Time-Independent Parameters

We consider the explicit solution of the MB call option price (15) in the case of time-independent parameters a t = a R , v t = v R + , ρ t = ρ R , with ρ < a (to ensure no arbitrage). From (2) and (4), the risky and riskless asset prices are
A t = A 0 + a t + v B t , t 0 , A 0 > 0 ,
β t = β 0 + ρ t , t 0 , β 0 > 0 .
The FK solution (9a) and (9b) becomes
f x , t = β 0 + ρ t β 0 + ρ T   E Q 3   g Z T     Z t = x , t 0 , T ,
d Z t = ρ β t   Z t   d t + v   d B t Q 3 .
With the choice β 0 = A 0 , the t = 0 price of the option is
C 0 = f A 0 , 0 = A 0 A 0 + ρ T   E Q 3   g Z T Z 0 = A 0 , t 0 , T .
Similar to (2), the solution to (A6) is
Z T = A 0 + ρ T 1 + v 0 T 1 A 0 + ρ u   d B u Q 3 .
This agrees with Equation (4) in Shiryaev (1999, p. 736), which was developed for the case ρ = 0 .21
If we consider pricing a call option, g Z T = m a x Z T K , 0 , its price (A7) is
C 0 = A 0   E Q 3 m a x 0 ,   1 + v 0 T 1 A 0 + ρ u   d B u Q 3 K A 0 + ρ T .
when ρ = 0 , and this agrees with Equation (6) in Shiryaev (1999, p. 736).
We note that the process (see Karatzas and Shreve (1988, Theorem 1.69, chap. 29))
X t = 0 t 1 A 0 + ρ u   d B u Q 3 ,
is Gaussian with mean m t = E Q 3 X t = 0 and covariance
Σ 2 s , t = E Q 3 B s Q 3 B t Q 3 = 0 s t 1 A 0 + ρ u 2   d u
= 1 ρ 1 A 0 1 A 0 + ρ   s t , s t m i n s , t .
Denoting the variance by Σ 2 s = Σ 2 s , s , (A9) can be written
C 0 = E m a x 0 ,   A 0 A 0 K A 0 + ρ T + v A 0 Σ T N 0 ,
where
A 0 A 0 K A 0 + ρ T + v A 0 Σ T N 0
is a normal random variable with mean value A 0 A 0 K / A 0 + ρ T and variance v 2 A 0 2 Σ 2 ( T ) . As (A12) is of the form a + b N 0 where a and b are constants, the computation of the expectation (A11) is straightforward, giving
C 0 = a Φ a b + b ϕ a b ,
where ϕ z and Φ z are, respectively, the standard normal probability density and the cumulative density functions. Thus, the explicit solution to the price of a call option in the MB model is (A13), where
a = A 0 A 0 K / A 0 + ρ T , b = v A 0 1 ρ 1 A 0 1 A 0 + ρ   T .
An analysis of the MB call option price in the limit ρ = 0 requires revisiting (A10) and noting that X t = B t Q 3 / A 0 , which is a mean-zero Gaussian process with variance t / A 0 2 . Consequently, (A11) becomes
C 0 = E Q 3 m a x 0 ,   A 0 K + v T N 0 ,
having explicit solution
C 0 = A 0 K Φ A 0 K v T + v T ϕ A 0 K v T ,
in agreement with Equation (12) in Shiryaev (1999, p. 737).

Notes

1
It becomes an academic exercise to consider whether ESG-valuation and financial valuation are correlated, and to what degree. In our view, that exercise is irrelevant to the critical importance of sustainability.
2
See Duffie (2001, sect. 5G and Appendix E). The regularity conditions are satisfied if ρ x , t and v x , t , x R , t 0 , are measurable and satisfy the Lipschitz and growth conditions in x R and 0 t ρ s d s < , 0 t v s 2 d s < for all t 0 . To simplify the exposition, we will assume that ρ t and v t , t 0 have trajectories that are continuous and uniformly bounded on [ 0 , ) .
3
See Equations (6.30), (6.31), and (6.32) in Karatzas and Shreve (1988, pp. 360–61).
4
For model consistency, the riskless asset must have the same dynamics as the risky asset, but with no stochastic component.
5
The restriction A t = x > 0 is required for the risky asset price under the BSM model.
6
Over the time period ( t , t + d t ] , the asset-holder receives a dividend D t A t D d t .
7
A stochastic càdlàg process has paths that are right-continuous with left limits.
8
See Hu et al. (2020a) for a similar approach to binomial option pricing.
9
This is consistent with the binomial option pricing formula given in Hu et al. (2020a).
10
We follow the exposition of the term structure presented in Shreve (2004b, chap. 10) and Chalasani and Jha (1997).
11
The determination of the T -forward price is based on no-arbitrage assumptions (Whaley 2012). Assuming V t ,   T 0 is not merely a simple transformation of the formula for the delivery price. If we assume V t ,   T = 0 , each trader entering the contract can multiply the value of the contract by any number, increasing the risk of their position without “paying for it”. In other words, assuming t ,   T = 0 is unrealistic in business practice.
12
In reality, will be responsible for all margin requirements (initial, variation, and maintenance) to the exchange where the futures contract is traded. Engaging in an unbounded number of futures contracts is impossible. In general, the assumption φ ( t i , n ) ( t i + 1 , n t i , n ) > 0 can be treated as a transaction cost, and our analysis can follow the exposition in Janecek and Shreve (2010, Definition 5.6.4, p. 244).
13
In our development the risk-neutral measure Q is determined by option pricing under the BBSM market model, while the exposition in Shreve assumes that Q is determined by option pricing under the BSM market model.
14
See https://www.marketwatch.com/market-data/futures (accessed on 1 August 2024).
15
The period τ should be relatively short, for example, 10 trading days. An exact value will need to be determined in practice.
16
“The central bank is the closest thing to an investor/issuer with zero risk aversion”. (R. Roll, private communication).
17
This is reinforced by the development of asset pricing in Duffie (2001, pp. 86, 106).
18
As a simple example, evaluating (4) assuming ρ t = ρ R and r t = r R are time independent shows that the riskless asset price β t remains positive only if r β 0 + ρ 1 e r t > 0 for all t 0 .
19
Whose form is “a function of” the parameters of the model.
20
Analysis of the condition A 0 + ρ / r e r u > ρ / r for u [ s , T ] shows that it holds for all 0 s u T < , except when r > 0 , ρ < 0 , and A 0   < ρ / r . In this case, T must be restricted to
T 1 r ln ρ / r ρ / r A 0 .
21
There is a slight error in Equation (4) of Shiryaev (1999, p. 736) where the placement of the real-world measure P T and the risk-neutral measure P ~ T should be reversed.

References

  1. Bachelier, Louis. 1900. Théorie de la Spéculation. Annales Scientifiques de l’École Normale Supérieure 3: 21–86. [Google Scholar] [CrossRef]
  2. Björk, Thomas. 2009. Arbitrage Theory in Continuous Time, 3rd ed. Oxford: Oxford University Press. [Google Scholar]
  3. Black, Fischer, and Myron Scholes. 1973. The Pricing of Options and Corporate Liabilities. Journal of Political Economy 81: 637–54. [Google Scholar] [CrossRef]
  4. Brooks, Robert, and Joshua A. Brooks. 2017. An Option Valuation Framework Based on Arithmetic Brownian Motion: Justification and Implementation Issues. Journal of Financial Research 40: 401–27. [Google Scholar] [CrossRef]
  5. Chalasani, Prasad, and Somesh Jha. 1997. Steven Shreve: Stochastic Calculus and Finance. Available online: http://efinance.org.cn/cn/FEshuo/stochastic.pdf (accessed on 29 June 2023).
  6. Choi, Jaehyuk, Minsuk Kwak, Chyng Wen Tee, and Yumeng Wang. 2022. A Black–Scholes User’s Guide to the Bachelier Model. Journal of Futures Markets 42: 959–80. [Google Scholar] [CrossRef]
  7. Cox, John. C., Stephen A. Ross, and Mark Rubinstein. 1979. Options Pricing: A Simplified Approach. Journal of Financial Economics 7: 229–63. [Google Scholar] [CrossRef]
  8. Delbaen, Freddy, and Walter Schachermayer. 1994. A General Version of the Fundamental Theorem of Asset Pricing. Mathematische Annalen 300: 463–520. [Google Scholar] [CrossRef]
  9. Delbaen, Freddy, and Walter Schachermayer. 1998. The Fundamental Theorem of Asset Pricing for Unbounded Stochastic Processes. Mathematische Annalen 312: 215–50. [Google Scholar] [CrossRef]
  10. Duffie, Darrell. 2001. Dynamic Asset Pricing Theory, 3rd ed. Princeton: Princeton University Press. [Google Scholar]
  11. Hu, Yuan, Abootaleb Shirvani, Stoyan Stoyanov, Young Shin Kim, Frank J. Fabozzi, and Svetlozar T. Rachev. 2020a. Option Pricing in Markets with Informed Traders. International Journal of Theoretical and Applied Finance 23: 2050037. [Google Scholar] [CrossRef]
  12. Hu, Yuan, Abootaleb Shirvani, W. Brent Lindquist, Frank J. Fabozzi, and Svetlozar T. Rachev. 2020b. Option Pricing Incorporating Factor Dynamics in Complete Markets. Journal of Risk and Financial Management 13: 321. [Google Scholar] [CrossRef]
  13. Hu, Yuan, W. Brent Lindquist, Svetlozar T. Rachev, Abootaleb Shirvani, and Frank J. Fabozzi. 2022. Market Complete Option Valuation using a Jarrow-Rudd Pricing Tree with Skewness and Kurtosis. Journal of Economic Dynamics and Control 137: 104345. [Google Scholar] [CrossRef]
  14. Hull, John C. 2012. Options, Futures and Other Derivatives, 8th ed. Essex: Pearson Education. [Google Scholar]
  15. Janecek, Karel, and Steven E. Shreve. 2010. Futures Trading with Transaction Costs. Illinois Journal of Mathematics 54: 1239–84. [Google Scholar] [CrossRef]
  16. Jarrow, Robert A., and Andrew Rudd. 1983. Option Pricing. Homewood: Dow Jones-Irwin. [Google Scholar]
  17. Karatzas, Ioannis, and Steven E. Shreve. 1988. Brownian Motion and Stochastic Calculus, 2nd ed. Berlin: Springer. [Google Scholar]
  18. Lindquist, W. Brent, and Svetlozar T. Rachev. 2024. Alternatives to Classical Option Pricing. arXiv arXiv:2403.17187. [Google Scholar]
  19. Merton, Robert C. 1973. Theory of Rational Option Pricing. Bell Journal of Economics and Management Science 4: 141–83. [Google Scholar] [CrossRef]
  20. Neufeld, Dorothy. 2022. Mapped: Interest Rates by Country in 2022. World Population Review. Available online: https://advisor.visualcapitalist.com/mapped-interest-rates-by-country-in-2022/ (accessed on 29 June 2023).
  21. Pástor, Ľuboš, Robert F. Stambaugh, and Lucian A. Taylor. 2021. Sustainable Investing in Equilibrium. Journal of Financial Economics 142: 550–71. [Google Scholar] [CrossRef]
  22. Pedersen, Lasse Heje, Shaun Fitzgibbons, and Lukasz Pomorski. 2021. Responsible Investing: The ESG-Efficient Frontier. Journal of Financial Economics 142: 572–97. [Google Scholar] [CrossRef]
  23. Rachev, Svetlozar T., Nancy Asare Nyarko, Blessing Omotade, and Peter Yegon. 2024. Bachelier’s Market Model for ESG Asset Pricing. In Directional and Multivariate Statistics: A Volume in Honour of Professor Ashis Sen-Gupta. Edited by S. Kumar. Berlin: Springer. [Google Scholar] [CrossRef]
  24. Schachermayer, Walter, and Josef Teichmann. 2008. How Close are the Option Pricing Formulas of Bachelier and Black–Merton–Scholes? Mathematical Finance 18: 155–70. [Google Scholar] [CrossRef]
  25. Shirvani, Abootaleb, Stoyan V. Stoyanov, Svetlozar T. Rachev, and Frank J. Fabozzi. 2020. A New Set of Financial Instruments. Frontiers in Applied Mathematics and Statistics 6: 606812. [Google Scholar] [CrossRef]
  26. Shiryaev, Albert N. 1999. Essentials of Stochastic Finance; Facts, Models, Theory. Singapore: World Scientific. [Google Scholar]
  27. Shreve, Steven E. 2004a. Stochastic Calculus for Finance I: The Binomial Asset Pricing Model. Berlin: Springer. [Google Scholar]
  28. Shreve, Steven E. 2004b. Stochastic Calculus for Finance II. Continuous-Time Finance. Berlin: Springer. [Google Scholar]
  29. Whaley, Robert E. 2012. No-Arbitrage Price Relations for Forwards, Futures, and Swaps, Derivatives Valuation. Encyclopedia of Financial Models, Wiley Online Library. Available online: https://onlinelibrary.wiley.com/doi/epdf/10.1002/9781118182635.efm0031 (accessed on 29 June 2023).
  30. Zerbib, Oliver David. 2022. A Sustainable Capital Asset Pricing Model (S-CAPM): Evidence from Environmental Integration and Sin Stock Exclusion. Review of Finance 26: 1345–88. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lindquist, W.B.; Rachev, S.T.; Gnawali, J.; Fabozzi, F.J. Dynamic Asset Pricing in a Unified Bachelier–Black–Scholes–Merton Model. Risks 2024, 12, 136. https://doi.org/10.3390/risks12090136

AMA Style

Lindquist WB, Rachev ST, Gnawali J, Fabozzi FJ. Dynamic Asset Pricing in a Unified Bachelier–Black–Scholes–Merton Model. Risks. 2024; 12(9):136. https://doi.org/10.3390/risks12090136

Chicago/Turabian Style

Lindquist, W. Brent, Svetlozar T. Rachev, Jagdish Gnawali, and Frank J. Fabozzi. 2024. "Dynamic Asset Pricing in a Unified Bachelier–Black–Scholes–Merton Model" Risks 12, no. 9: 136. https://doi.org/10.3390/risks12090136

APA Style

Lindquist, W. B., Rachev, S. T., Gnawali, J., & Fabozzi, F. J. (2024). Dynamic Asset Pricing in a Unified Bachelier–Black–Scholes–Merton Model. Risks, 12(9), 136. https://doi.org/10.3390/risks12090136

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop