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
at time
is
where
is the stock price of the company;
is the ESG score of the company;
is the ESG score of the market index
to which the company stock belongs; and the parameter
measures the ESG affinity of the market. As
will be negative for low ESG stocks, if the market ESG affinity is valued at
, the ESG-adjusted price
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.
- (i)
- (ii)
used arithmetic Brownian motion for the risky asset with a continuously compounded risk-free asset (
Choi et al. 2022); or
- (iii)
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
consists of a risky asset
, a riskless asset
, and a derivative
whose underlying asset is
. Unless otherwise noted, the derivative will be a European contingency claim (i.e., option), which we denote by ECC. The model assumes that
has the price dynamics of a continuous diffusion process determined by the stochastic differential equation (SDE)
where
,
, is a standard Brownian motion on a stochastic basis (filtered probability space)
on a complete probability space
. In (1),
,
, and
are
-adapted processes on
, satisfying the usual regularity conditions
2. We require
-almost surely (a.s.). The SDE (1) is a one-dimensional linear equation whose solution can be written
3
where
From (1), it is clear
4 that riskless asset
should have the price dynamics:
where
and
. Applying the Karatzas–Shreve solution to (3), we have
; therefore,
As
and
must have the same numèraire, the natural choice is to set
. We will use this choice throughout the article. A necessary condition for no arbitrage is the requirement
,
-a.s. for
. The no-arbitrage assumption leads to the requirement on the market price of risk,
As we require -a.s., (5) further requires that -a.s.
With (1) expressed as
and (3) expressed as
, we see that the BBSM model is properly viewed as a generalized Bachelier model that reduces to the BSM model when
. When
, the BBSM model reduces to the MB model of
Rachev et al. (
2024).
The ECC
has the price
, where
,
,
has continuous partial derivatives
and
on
;
is the terminal (maturity) time of
; and the option’s terminal payoff is
for some continuous function
. From It
’s formula,
3. Option Pricing with Riskless Bonds under the BBSM Model
The classical treatment of a replicating, self-financing portfolio
of
and
has the price process (
Duffie 2001, chap. 5E, p. 89):
where
and
are
-adapted processes. This portfolio treats the riskless asset as a financial product issued in units (i.e., a bond), with
being the unit price and
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
for
,
, with the boundary condition
.
Assuming that
and
satisfy regularity conditions sufficient to guarantee that (8) has a unique strong solution, the Feynman-Kac (FK) solution of the BBSM PDE (8) is
where
is a standard Brownian motion under the probability measure
.
We confirm that the probability measure
is the equivalent martingale measure (EMM) obtained from risk-neutral valuation. The FK solution (9a) suggests the deflator
We define the standard Brownian motion
, having the dynamics
on
. Under the EMM
, the deflated process
,
, is an
-martingale with
,
. Consequently,
As the equivalent martingale measure
is unique, the market
is complete. Therefore, for the option price process
,
, with terminal payoff
,
is an
-martingale under
. Thus,
Equations (11) and (12) obtained from the EMM are the same as the option price determined by the FK solution (9a), (9b) of the BBSM PDE. Thus, as claimed.
The valuation of the BBSM option price
with terminal payoff
is given by
where
is given by (2). In
Appendix A, we develop the explicit solution of the BBSM call option price for time-independent parameters.
Setting
(8) reduces to the BSM PDE (
Duffie 2001, sect. 5G)
for
,
5 , with boundary condition
Under this limit, (9a) and (9b) reduce to the usual FK BSM solution (
Duffie 2001, chap. 5H)
Similarly (13) reduces to the familiar valuation of the option price in the BSM model (
Duffie 2001, sect. 6H)
where
and the EMM
is determined by
having the dynamics
on
.
Setting
, (8) reduces to the MB PDE (
Rachev et al. 2024)
for
,
, with boundary condition
Under this limit, (9a) and (9b) reduces to the FK MB solution
Similarly, (13) produces the risk-neutral valuation of the option price in the MB model
with
and the EMM
is determined by
having the dynamics
on
. 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
. The replicating, self-financing portfolio can then be written:
We note that (16a) is equivalent to a self-financing portfolio consisting of a long position in one option minus
shares of the risky asset with per-share price
. As this portfolio should be riskless, it must equal the value
in the bank account (
Hull 2012, p. 255). The dynamics (16b) follow similarly.
Equating (16b) to (4), and solving for
to eliminate the Brownian motion term
, produces the coupled system of PDEs
for
,
, with boundary condition
. For a general functional form
, (17a) and (17b) must be solved simultaneously for
.
Under the BBSM model,
is linear in
and (17a) and (17b) reduce to the single risk-neutral PDE,
When
, (18) reduces to the PDE analyzed in
Rachev et al. (
2024). When
, (18) is identical to (8). We emphasize this result in the following lemma.
Lemma 1. BSM-like models with 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
We next confirm that
is the EMM obtained using the risk-neutral approach to option valuation using the deflator
. We define the standard Brownian motion
, as having the dynamics
on
. Under the measure
, the deflated process
,
, is an
-martingale with
,
. Consequently,
Taking the total derivative of (20), it is straightforward to show that under , is given by (19b). The EMM is unique and the market is complete.
From the solution (4), it is straightforward to show
Writing the price vector as
and the continuous rate vector as
, the results (20) and (21) can be summarized as
in agreement with the development of arbitrage pricing with dividends (
Duffie 2001, chap. 6L, pp. 123–25). Thus, the term
in the riskless bank account of (16a), (16b) has the effect of adding a continuous dividend rate
to the market
. 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
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
.
As the market is complete, the option price process
,
, with terminal payoff
,
must also be an
-martingale under
. From (19a),
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
. 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
obtained via multiplicative
and additive
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
when the asset
pays a dividend at a continuous rate with an instantaneous dividend yield
.
6 By the usual no-arbitrage argument, the price
of
obeys
Equation (24) describes the dynamics of the underlying asset from the point of view of the buyer of the ECC. The riskless bank account
has the price dynamics (3). The ECC
has the price
,
, where
,
,
has continuous partial derivatives
and
on
. The terminal time of
is
with payoff
Itô’s formula for the option price dynamics is
We consider the replicating, self-financing riskless portfolio of
and
having the price
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
The standard computation yields the risk-neutral BBSM PDE for the option price
for
,
, with boundary condition
The FK solution of (25) is
Setting
in (25) produces the familiar BSM PDE when the underlying asset pays out a dividend yield
(
Björk 2009, p. 162),
,
, with boundary condition
In this limit, (26a), (26b) results in the familiar FK solution (
Duffie 2001, chap. 5H) of the BSM PDE (27)
Setting
in (25) produces the MB PDE,
,
, with boundary condition
In this limit, (26a), (26b) produces the solution of the MB PDE (28),
where
We note that the risk-neutral measures 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,
the FK solution for the option price is
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
and riskless asset
whose price process are given, respectively, by (1) with
,
,
, and by (3) with
,
. For
,
,
, represents the price of a perpetual derivative,
, in the complete no-arbitrage market
. There are two major applications of
. The first is to use
when the trader does not have access to the riskless bank account
The trader can then form a self-financing portfolio of
and
to replicate
. The second is the use of
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
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
for the risk-neutral PDE. Using (29), (8) becomes
The right-hand side is deterministic, whereas the left-hand side contains the stochastic process
. Thus, both sides must separately vanish, giving the ordinary differential equations
Equation (30b) suggests solutions of the form
. Substituting into (30b) produces
Equation (32) admits a solution as long as
, which is possible as long as
, in which case
For the BSM model (
), the behavior of the function
for time-independent coefficients
,
was studied by
Lindquist and Rachev (
2024). In that case,
with
. This produced a one-parameter family of solutions to (32), one of which is the parameter value
, as previously used by
Shirvani et al. (
2020).
For the BBSM model, as long as
, (8) admits a perpetual derivative of the form (29), where
is given by (31) and
is given by
Under time-independent coefficients
,
, and
, an examination of
in (33) shows it has the form
. Thus, at fixed time value
, the analysis of the form of
in
Lindquist and Rachev (
2024) holds with the replacement
. For example, if
, then
changes monotonically from
to
. We summarize the principal result of this section in the following lemma.
Lemma 2. The requirement 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
. Under continuous-time option pricing models, the hedger (who takes a short position in
) 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
follows the binomial pricing tree
In (34):
- (i)
, , is the stock price at time , , , where is the fixed terminal time, ; and
- (ii)
for every , the , are independent, identically distributed Bernoulli random variables with determining the filtration and the stochastic basis on the complete probability space .
The riskless asset
has the discrete price dynamics
where
is an instantaneous simple riskless rate. We again invoke a common numèraire via
. Under the BBSM model, stock price changes (i.e., lag-1 differences) rather than normalized returns are of primary interest. The price changes
have the discrete dynamics
We would like to have weak convergence of the càdlàg process
7 generated by the binomial tree (34) in the Skorokhod space
to the continuous time process (1).
8 To achieve this, we must consider an instantaneous mean
and instantaneous variance
for the discrete price differences of the risky asset such that the conditional mean and conditional variance satisfy
As a consequence,
and
are determined by
The ECC
will have the discrete price dynamics
,
, with terminal value
. Consider a self-financing strategy,
, replicating the option price
:
Solving this system in the usual manner produces the recursion formula for the risk neutral valuation for the option price on the tree,
where
, and the risk-neutral probability is
The ratio appearing in (38) reinforces our choice of as the natural common numèraire. Enforcement of the “time value of money” requires for all . To ensure no opportunity for arbitrage, we require for all , ensuring that .
Using (36), the risk-neutral probability (38) can be written
The natural world probability satisfies ; however, as noted above, the requirement must be imposed.
The discrete version of (1) and (3) implies
,
and
, with the requirement that
and
. Thus, the BBSM risk-neutral probability (39) can be written
We note that
(and hence
) is positive for sufficiently small values of
. In general, as noted above, the condition
must be imposed. Arguments for the weak convergence of binomial-tree-generated càdlàg processes in the Skorokhod space
show (see
Hu et al. 2022)
that the discrete option price (37) converges weakly to (12).The choice
produces the risk-neutral valuation of the option (37) on a binomial tree in the BSM market model,
where the risk-neutral probability
is given by
9Here, is the discrete form of the BSM market price of risk , .
The choice
produces the risk-neutral valuation of the option (37) on a binomial tree in the MB market model,
with
given by
where
. For time-independent parameters, (42) becomes
Under time-independent parameter values,
, with
, and
For values
(i.e.,
),
which resembles the discounting in the BSM option pricing (40). However, for
(i.e.,
approaching the maturity time
), 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
and
. 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
and
. However, we raise a significant issue regarding these models, each of which requires knowledge of the risk-neutral probability measure
or, equivalently, the probability
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
(
Shreve 2004a, chap. 6). In our view, the measure
and the characterization of
must come from the appropriate market. As demonstrated in
Hu et al. (
2020b, Figure 3), the value of
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
. 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
where the dynamics of the price
, of the risky asset is given by (1) and the dynamics of the price
,
, of the riskless asset is given by (3). Here,
is a zero-coupon bond (see
Duffie (
2001, sect. 6M) and
Shreve (
2004b, sect. 5.6.1)), having maturity time
and price
. The payoff
for
is 1 if
and 0 for
. Under the BBSM model,
is given by (12)
In (45),
represents a fixed time horizon at, or prior to which, all bonds will mature. The family
constitutes the TSIR in the BBSM model.
10Setting in (45) produces the TSIR under the BSM model, . Setting in (45), we obtain the TSIR under the MB model, , where , , .
In (45) and its BSM and MB limit cases, the expected value is taken with respect to the risk-neutral measure
. Evaluation of (45) simplifies considerably under time-independent parameters to
where
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
and
from market option prices. Using the calibrated values for
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
in (45) is a martingale on
By the martingale representation theorem (
Duffie 2001, p. 336),
for some
-adapted
in
. Thus, as
, then
We consider the replicating self-financing portfolio having the price dynamics:
Solving for
and
in the usual manner,
produces the replicating, self-financing strategy for the zero-coupon bond in the BBSM market model. Setting
or
produces, respectively, the replicating, self-financing strategy for the zero-coupon bond under the BSM or MB market models. If
and
,
, are deterministic, then
is deterministic and
,
, 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 , two market agents ( and ) enter into a forward contract on the risky asset , with taking the long position and taking the short position. When the contract reaches maturity, ’s payoff is and ’s payoff is ), where F(t,T) is the -forward (delivery) price of asset evaluated under the BBSM model.
With
denoting the initial value for the forward contract, from (12),
Thus, the
-forward price of the asset
is
11The traditional definition of a forward contract under the BSM model assumes
. A more realistic scenario involves an initial forward contract price
. Then, for
the value (47) of the BBSM forward contract is
We conclude with the hedging strategy that should follow when . agrees to pay for the asset with value at . At initiation time , receives . As per (47), receives the precise amount to buy the asset and takes a short position in zero-coupon bonds. From (48), at terminal time , ’s position has the value ; delivers the asset and receives the forward price
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
,
,
,
be a sequence of trading instances (of the forward contracts). To maintain the futures contract, trader
should pay
. The classical setting of futures contracts assumes that
, allowing
to take on an unlimited number of futures contracts
12. In what follows, we follow the exposition in
Shreve (
2004b, sect. 5.6.2).
13For notational convenience, we write the accumulation factor
as
The instantaneous rate
,
, is an
-adapted process on
. We assume the regularity condition
To simplify the connection between the discrete- and continuous-time futures prices, we further assume that
,
, is continuous on
-a.s. Let
, and assume that
as
. We define the piecewise constant rate,
for
,
, with
. Using
, we define the discrete accumulation factor
Trader
enters a long futures contract at time
,
, when the market price of the futures contract is
. At time
, when the market price is
,
makes a profit or loss of
. From (12), the value at time
of the payment
is
As it costs nothing to enter the futures contract at time
,
Informally, passing to the limit as
(49) becomes
for some
-adapted process
,
, on the stochastic basis
. The same argument yields
As
, the final equality in each of (49), (50) and (51) also holds
. The
-futures price of asset
in the BBSM market is any
-adapted stochastic process
,
, satisfying (51) with
-a.s. (and thus
-a.s.). As shown in
Shreve (
2004b, sect. 5.6.2), the unique
-adapted process
, that satisfies (51) with
is
It is of interest to employ futures market data
14 to calibrate time-independent parameter values
and
By calibrating
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 (). We have concentrated on developing derivative pricing solutions for three cases: (i) all parameters are non-zero, (ii) the BSM limit in which , and (iii) the MB limit . Note, however, that there are parameter choices, such as , 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
to the historical price series of asset
. Let
be the price at the time
. Let
denote the mean value of
’s price over the time period
.
15 Repeating this process for times
will generate data for a robust regression model to fit the series
for the parameters
and
. A similar procedure can be used to calibrate
, where
is estimated from the variance of the data over the time period
. Again, robust regression should be used to fit the series
but now subject to the restriction
,
.
The calibration of
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
and
are estimated from spot prices as described above. The probability of upward price movement
can be directly estimated from price changes
,
. The parameters
and
can be estimated by minimizing
where for time
,
is the listed option price corresponding to the maturity time
and strike price
;
is the theoretical price computed from (37);
is the number of listed maturity times; and for each maturity time
,
is the number of listed strike prices. Note that this minimization must be performed subject to the restriction
.
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
such that the deflated process
is an
-martingale on
. In (52),
and
are constants. We assume
and
are time-dependent,
-adapted Itô processes. Further, we assume
is strictly positive (
Duffie 2001, chap. 6B) and invertible. Denoting the derivatives
and
,
Define the standard Brownian motion
, with dynamics given by
on
. Under the EMM
, the deflated process (52) satisfies
. Consequently,
where
. As there is a single central bank, the market is complete and an option having the risky asset as underlying with payoff
at maturity date
will have the price
From these results, we see that only the ratio
enters the pricing. Thus, (52) can be recast as
.
If we assume that the effect of the central bank’s actions on setting U.S. Treasury rates produces the riskless rate dynamics
, then in (53),
and the option price model (23) is identical to (53) under the change
.
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 ), 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
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) , ; (ii) , ; or (iii) , provides a better fit. Assuming time-independent parameters for simplicity, the BBSM model, with the addition of 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 .
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
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.
18While 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.