Skip to main content

ORIGINAL RESEARCH article

Front. Appl. Math. Stat., 21 April 2023
Sec. Mathematical Finance

Three little arbitrage theorems

\nMauricio Contreras G.Mauricio Contreras G.1Roberto Ortiz H.,
Roberto Ortiz H.2,3*
  • 1Departamento de Física, Universidad Metropolitana de Ciencias de la Educación (UMCE), Santiago, Chile
  • 2Facultad de Ingeniería, Universidad Diego Portales, Santiago, Chile
  • 3Facultad de Ciencias Económicas y Administrativas (FACEA), Universidad Católica de la Santísima Concepción, Concepción, Chile

The authors proved three theorems about the exact solutions of a generalized or interacting Black–Scholes equation that explicitly includes arbitrage bubbles. These arbitrage bubbles can be characterized by an arbitrage number AN. The first theorem states that if AN = 0, then the solution at maturity of the interacting equation is identical to the solution of the free Black–Scholes equation with the same initial interest rate of r. The second theorem states that if AN ≠ 0, then the interacting solution can be expressed in terms of all higher derivatives of the solutions to the free Black–Scholes equation with an initial interest rate of r. The third theorem states that for a given arbitrage number, the interacting solution is a solution to the free Black–Scholes equation but with a variable interest rate of r(τ) = r + (1/τ)AN(τ), where τ = Tt.

1. Introduction

Since its introduction in 1973, the initial Black–Scholes model [1, 2] has undergone many changes, which have given rise to many different financial models, such as stochastic volatility models [3, 4] and the associated concept of the volatility smile [5, 6] and stochastic rate models [5, 7, 8], which account for the dynamics of the spot interest rate in the determination of option pricing, as well as, for example, the incorporation of jumps; this gives rise to integro differential equations for the price of the option [9]. All these generalizations are related to relaxing some of the assumptions of the initial Black–Scholes (BS) model.

One of the most important of these initial axioms is the non-arbitrage hypothesis, which is associated with equilibrium market dynamics. In analytic terms, the non-arbitrage hypothesis in the BS model can be characterized as follows: if B(t) and S(t) are the risk-free asset and underlying stock prices, the price dynamics of the bond and the stock in the usual BS model are given by the following equations:

dB(t)=rB(t)dt,    (1)
dS(t)=αS(t)dt+σS(t) dW(t),    (2)

where r is the constant risk-free interest rate, α and σ are the drift rate and the volatility of S, respectively, and W(t) is a Wiener process. To price the financial derivative π, it is assumed that it can be traded, so we can form a portfolio P based on the derivative π and the underlying stock S (no bonds are included). Considering only non-dividend-paying underlying assets and not considering consumption portfolios, the purchase of a new portfolio must be financed only by selling from the current portfolio.

Consider h(t)=(hS,hπ), where hS is the number of underlying assets S and hπ is the number of derivatives π present in the portfolio. Z(t)=(S,π) is the price vector of shares (underlying asset and option), and the value of the portfolio P at time t is given in equation (3),

P(t)=h(t)Z(t)=hS(t)S(t)+hπ(t)π(S(t),t),    (3)

and the dynamics of a self-financing portfolio with no consumption imply that [5, 7].

dP(t)=h(t)dZ(t)=h(t)dS(t)+hπ(t)dπ(t).    (4)

In other words, in a model without exogenous incomes or withdrawals, any change in value is due to changes in asset and derivative prices. Another important assumption for deriving the BS equation is that the market is efficient in the sense that it is free of arbitrage possibilities. This is equivalent to the fact that there exists a self-financed portfolio for which the portfolio return dP satisfies the standard equilibrium non-arbitrage condition [5, 7, 10].

dP(t)=rP(t)dt.    (5)

As standard textbooks show [5, 7, 10], equations (4), (5), and Itô's lemma for π(S, t) imply that the option price π satisfies the free or equilibrium Black–Scholes equation

πt+12σ2S22πS2+r(SπS-π)=0.    (6)

When deviations from this equilibrium state are considered (deviations can occur for several external reasons, including market imperfections such as transaction costs, asymmetric information issues, short-term volatility, and extreme discontinuities, among many others), the classical no-arbitrage assumption is violated. Thus, other types of models attempt to relax the no-arbitrage hypothesis to incorporate such non-equilibrium behavior. In fact, since the 1980s, economists have realized that future contracts are not always traded at a price predicted by the simple no-arbitrage relation in a real market. Substantial empirical evidence has supported this point many times and in different settings [1114]. However, economists have developed several alternative explanations for the variability of the arbitrage, such as differential tax treatments for spots and futures [15], and marking-to-market requirements for futures [16]. It has also been noted that there are certain factors that influence the arbitrage strategies and slow down the market's reaction to arbitrage. These factors include constrained capital requirements [17], position limits, and transaction costs [18].

2. Some arbitrage models

In this section, we analyze some models that explicitly incorporate the notion of arbitrage.

We start with the study of Ilinski [19] and quoted him in this first part, who states that arbitrage pricing theory (APT) [20] is a standard model for determining the expected rate of return on individual stocks and on a portfolio of stocks [21]. What the APT does, in simplified terms, is that the return on a risk-less portfolio should be equal to the risk-less rate of return, which can be taken equal to the rate of the return on a bank deposit (perfect capital market conditions are assumed). These ideas are summarized in equation (5).

As Ilinski [19] said, if any arbitrage possibility existed, then agents (arbitrageurs) would use the opportunity to make an abnormal risk-less profit, which itself will bring the system to the equilibrium and eliminate the arbitrage opportunity. Thus, the arbitrage cannot exist for long and its lifetime depends on the liquidity of the market. This fact, however, does not mean that arbitrage opportunities do not exist at all and cannot influence the asset pricing, violating the APT assumption. That is why Ilinski [19] tries to overcome the no-arbitrage assumption [equation 5] and suggests a model to account for the existence of virtual arbitrage opportunities and their influence on asset pricing in the framework of the APT. To this end, Ilinski [19] considers a risk-less portfolio P1 which is made up of (N + 1) assets. In the case of no-arbitrage, the portfolio P1 would satisfy the condition [equation 5]

dP1(t)=rP1(t)dt,    (7)

where r is the risk-less interest rate. Furthermore, Ilinski [19] introduces the idea of virtual arbitrage in derivative pricing. In this case the right side of equation (7) is changed by a factor R(t,P1)P1dt, so

dP1(t)=rP1(t)dt+R(t,P1)P1(t)dt,    (8)

where R(t,P1) represents the virtual arbitrage return. As Ilinski [19] said, to find an expression for R(t,P1), one can imagine that at some moment of time τ < t, a fluctuation of the return (an arbitrage opportunity) appeared in the market [this instantaneous arbitrage return is denoted by ν(τ, P1)]. Arbitrageurs would react to this circumstance and act in such a way that the arbitrage gradually disappears and the market returns to its equilibrium state, i.e., the absence of arbitrage. For small enough fluctuations, it is natural to assume that the arbitrage return R (in the absence of other fluctuations) evolves according to the following equation as follows [19]:

dRdt=-λR   (t>τ),   R(t=τ)=ν(τ,P1)

with a decay parameter λ, which is characteristic for the market. According to Ilinski [19], this parameter can be either estimated from a microscopic theory as in Ilinski [22] or can be found from the market using an analog of the fluctuation–dissipation theorem [19, 23].

The stochastic process ν(τ, P1) can be specified by assuming that its fluctuations at different times are independent and have the form of a white noise with a variance Σ2 which depends on the structure of the portfolio P as follows [19]:

ν(τ,P1)=0,   ν(τ,P1) ν(τ,P1)=Σ2δ(τ-τ).

A second way to introduce these same ideas, perhaps a more natural way, is to assume that there exist market short-life arbitrage statistical fluctuations x(t) which can be characterized, for example, as in previous studies Ilinski [13]. Ilinski [24] by an Ornstein–Uhlenbeck process with random noise η(t) of the form

dx(t)dt=-λx(t)+η(t),

where the decay rate λ is related to the characteristic life time τarb of the arbitrage fluctuation by λ=1τarb, and the white noise η(t) is defined by

η(t)=0,   η(t) η(t)=Σ2δ(t-t).

When these arbitrage fluctuations are present in the market, the portfolio returns dP1/P1 cannot be in the equilibrium given by equation (7) but would depend on both the spot rate r and the fluctuations x(t). Thus, one can assume that dP1/P1 is some function F = F(r, x(t)), with the condition F(r, 0) = r to recover the equilibrium case [equation (7)] for x = 0. In this way, on general grounds, one can generalize [equation 7] in the presence of arbitrage as follows:

dP1(t)/P1(t)=F(r,x(t))dt.    (9)

By expanding F in a Taylor series around x = 0 and keeping only the first terms, one has

dP1(t)/P1(t)=(F(r,0)+dFdx(r,0) x(t)+...)dt.    (10)

or

dP1(t)/P1(t)=(r+dFdx(r,0) x(t)+...)dt.    (11)

that is,

dP1(t)=rP1(t)dt+cx(t)P1(t)dt+...    (12)

where the constant c is c=dFdx(r,0). Now by considering the rescaled process x(t) → cx(t) instead of x(t), one can write finally

dP1(t)=rP1(t)dt+x(t)P1(t)dt+...    (13)

Thus, for fluctuations with x(t)/r < < 1, one can keep only the first two terms in the earlier expansion and the non-arbitrage hypothesis [equation 7] or [equation 5] becomes in this case an arbitrage hypothesis

dP1(t)=rP1(t)dt+x(t)P1(t)dt    (14)

with the same form as equation (8). Note that the no-arbitrage hypothesis (5) would be called really an equilibrium arbitrage hypothesis instead. After that, Ilinski and Stepanenko [25], using equation (14), proceed to derive (using a portfolio P made of one underlying asset S and one option π) the following PDE for the option price π [13].

πt+12σ2S22πS2+r(SπS-π)=x(t)(π-SπS),    (15)

which contains the effect of the arbitrage perturbation x(t) in the option dynamics.

Another effort in this line is Panayides [26] who [following Panayides and Fedotov [27]] considers a market that consists of a stock S, a bond B, and a European option π. The market is assumed to be affected by two sources of uncertainty, the random fluctuations of the return from the stock S, whose dynamics are given in equation (2), and a random arbitrage return from the bond B described by the equation

dBB=rdt+ξ(t)dt,    (16)

where the random process ξ(t) describes the fluctuations of the arbitrage return around the spot rate r. The random variations of arbitrage return ξ(t) are assumed to be on the scale of hours. This characteristic time is denoted by τarb, and it as an intermediate one between the time scale of a random stock return S and the lifetime T of the derivative (several months), so 0 < < τarb < < T.

To obtain the corresponding partial differential equation PDE for the option price π, one must consider an investor with a zero initial investment position by creating a portfolio P̄ consisting of one bond B, hS=-VS shares of the stock S, and one European option π. The value of this portfolio P̄ is given by

P̄=B-πSS+π.

The usual Black–Scholes dynamics of this portfolio can be obtained from the equations P̄/t=0, P̄=0, and ξ(t) = 0, which are equivalent of course to equation (5) in terms of the portfolio P (P̄=B+P).

When arbitrage is present, one can consider a generalization of P̄/t=0 by taking the simple non-equilibrium equation [26]

P̄t=-P̄τarb,

where τarb is the characteristic time during which the arbitrage opportunity ceases to exist. Using the self-financing condition dP̄=dB-VSdS+dV, Ito's lemma, and equations (2), (16), it is shown in Panayides [26] that the PDE equation for the option value π(S, t) is

πt+σ2S222πS2+rSπS-rπ+rP̄+ξ(t)P̄+ξ(t)(SπS-π)+P̄τarb=0,

which reduces to the classical Black–Scholes PDE when P̄=0 and ξ(t) = 0. By using the non-dimensional time

τ=T-tT,   0τ1,

and a small parameter

ε=τarbT<<1,

in the limit ε → 0, Panayides [26], following an approach suggested by Papanicolaou and Sircar [28], shows that the option price π finally obeys the following PDE

πτ=σ2S222πS2+rSπS-rπ+ξ(τε)(SπS-π),    (17)

with the initial condition π(S, 0) = max(SK, 0) for a call option with strike price K. Note that equation (17) has the same form as (15).

Contreras et al. [29], inspired by the ideas of Ilinski [13], proposed a generalization of the Black–Scholes model that incorporates market imperfections using arbitrage bubbles. In this case, by using portfolio (3), the corresponding portfolio return dP is assumed to follow a stochastic dynamics of the form [instead of (14)]

dP(t)=rP(t)dt+f(S,t)P(t)dW,    (18)

where the amplitude f(S, t) of the Wiener process dW is a given deterministic function and is called an arbitrage bubble. These dynamics mimic equation (2), which determines the dynamics of the underlying asset price S (essentially f is the portfolio's volatility). In addition, equation (18) is the minimal change that one can make to the no-arbitrage hypothesis without incorporating external structures (e.g., new independent Brownian motions).1 For a more recent review of the theoretical economic literature on asset price bubbles, see, for example, Heston et al. [32] and Jarrow [33]. The self-financing portfolio condition (4) and equation (18) give

hS(t)dS(t)+hπ(t)dπ(t)=rP(t)dt+f(S,t)P(t)dW,    (19)

The Itô law for the variation of the derivative price

dπ=πtdt+πSdS+122πS2dS2,    (20)

plus (2) implies that

dπ=(πt+αSπS+σ22S22πS2)dt+σSπSdW,    (21)

by replacing (21) and (2) in (19), one obtains

hS(αSdt+σSdW)+hπ[(πt+αSπS+σ22S22πS2)dt+σSπSdW]      =r(hSS+hππ)dt+f(hSS+hππ)dW.    (22)

By collecting dt and dW terms in the above equation, one gets the following system [29]:

hSS(α-r)+hπ(πt+αSπS+σ22S22πS2-rπ)=0hSS(σ-f)+hπ(σSπS-fπ)=0,    (23)

or in matrix form

(S(α-r)(πt+αSπS+σ22S22πS2-rπ)S(σ-f)(σSπS-fπ))(hShπ)=(00).    (24)

The condition for the existence of non-trivial portfolios, (hS, hπ) ≠ (0, 0) (that is, the determinant of the square 2 × 2 matrix in equation (24) must be equal to 0), finally gives the following interacting Black–Scholes equation in the presence of an arbitrage bubble:

πt+12σ2S22πS2+r(σ-αf(S,t)r)(σ-f(S,t))(SπS-π)=0.    (25)

The aforementioned equation can written as follows:

πt+12σ2S22πS2+r(SπS-π)+v(S,t)(SπS-π)=0,    (26)

with

v(S,t)=(r-α)f(S,t)σ-f(S,t),    (27)

which can be interpreted as the potential of an external time-dependent force generated by the arbitrage bubble f(S, t). Note that if f = 0, the potential v = 0, and thus equation (26) is reduced to the usual free Black–Scholes equation (6). Then, under market imperfections, the free Black–Scholes model becomes an interacting one. If one writes equation (26) in the form

πt+12σ2S22πS2+(r+v(S,t))(SπS-π)=0,    (28)

one can see that the arbitrage bubble f(S, t) also changes the constant interest r in the Black–Scholes equation to a time-dependent one r+v(S,t)=r+(r-α)f(S,t)σ-f(S,t). This effect would generate a back-reaction on the bond dynamics (1), so one would write

dB(t)=(r+v(S,t))B(t)dt    (29)

instead, in the same way as equation (16). Thus, from now on, we define the interacting Black–Scholes model in the presence of an arbitrage bubble f(S, t), by the set of equations (2), (18), (26), (27), and (29).

Note that equation (26) has the same form as equation (15) obtained by Ilinski [13, 24] and equation (17) obtained by Panayides [26] and includes these last models as special cases. It is interesting that three different models give the same PDE which incorporates the arbitrage effects. Equation (26) can be then thought of as a more general PDE model for arbitrage processes. Thus, a generic way to represent a non-risk-free portfolio is given by equation (18), where the function f(S, t) encapsulates all of the information about the market equilibrium's deviations regardless of their causes. Then, in principle, any non-equilibrium option behavior can be modeled endogenously in the framework of equation (26) by providing the appropriate bubble form f(S, t).

Contreras et al. [29], after analyzing empirical financial data, showed that the mispricing between the empirical and Black–Scholes prices can be described by a Heaviside-type function in time. This implies that the arbitrage bubble f = f(S, t) has the same time dependence and can be modeled by a pure time-dependent bubble of the form

f=f(t)={00tT1,f0T1<t<T2,0T2tT,    (30)

where f0 is the strength of the bubble and T is the option's maturity. Thus, arbitrage bubbles can be characterized by a finite timespan and a constant height that measures the deviation from the Black–Scholes model's prices.

Note that in general, a time-dependent arbitrage bubble f = f(t) can be determined approximately from the empirical financial data [34] by using semi-classical methods [35].

In the following, we consider a generic pure time-dependent arbitrage bubble f = f(t) and study its impact on the option price π(S, t). This formulation was chosen principally because the model for a pure time-dependent bubble f = f(t) is completely integrable. In fact, as we show in this article, the solution π(S, t) of the interacting equation (equation (26)) can be written in an exact closed form in terms of the free Black–Scholes solution. For the case of a general bubble f = f(S, t), one needs a full perturbative expansion in the (S, t) space, and it is not clear whether one can find such an exact simple integrable model.

Using this integrable model, three theorems on the behavior of the derivative asset price π(S, t) in the presence of the arbitrage bubble f(t) are inferred. A fundamental quantity in this model is the accumulated potential at time t, which defines what we call the arbitrage number AN(t).

The first arbitrage theorem states that the solution π(S, t = 0) of the interacting Black–Scholes equation evaluated at maturity has the same value as the solution C(S, t = 0, r)2 of the free Black–Scholes equation at maturity only if the accumulated potential at time t = 0, that is, AN(t = 0), is exactly zero.

The second arbitrage theorem implies that when AN(t = 0) ≠ 0, the solution of the interacting Black–Scholes equation differs from the solution of the free equation at maturity. In fact, it is a series of all the free Greeks. Here, the interacting solution's dynamical behavior is completely different from that given by the free dynamics.

The third theorem gives a closed-form solution π(S, t) of the interacting Black–Scholes equation as a composition of the free Black–Scholes solution C(S, t, r) and a time-dependent rate r(τ)=r+1τAN(τ), which depends on the accumulated potential; in fact, π(S, τ) = C(S, τ, r(τ)), where τ = Tt.

Thus, the arbitrage number permits us to classify arbitrage bubbles into two principal categories: bubbles with AN(t = 0) = 0, which are called innocuous at maturity (these arbitrage bubbles have no net effect on the derivative value at maturity), and bubbles with AN(t = 0) ≠ 0, which are called dangerous at maturity (in this case, the arbitrage bubbles change substantially the derivative value at maturity). Finally, to illustrate the effects of these theorems, several examples of different types of bubbles acting on the price of a call option are given; these examples clearly show the distinction between the two types of bubbles.

This method allows us to estimate the economic impacts that different arbitrage bubbles can generate in real time. Note that the duration of the bubble can range from a few seconds to the time remaining until the derivative asset becomes mature.

The rest of this article is organized as follows. Section 3 develops a quantum mechanical interpretation of the interacting Black–Scholes model. Section 4 presents the three arbitrage theorems. In Section 5, several examples are used to illustrate the theorems, and Section 6 contains the conclusion. All mathematical manipulations are included in the Appendices.

3. Quantum mechanical interpretation

Consider the Black–Scholes equation [equation (26)] in the presence of an arbitrage bubble f = f(t). This equation must be integrated with the final condition,

π(S,T)=Φ(S).    (31)

The function Φ is called the contract function and defines the type of option. Note that equation (26) must be integrated backward in time from the future time t = T to the present time t = 0. One can change the direction of time by using the change of variables given by

τ=T-t,    (32)

which implies that

τ=-t,    (33)

so equation (26) can be written as a forward τ-time Euclidean Schrödinger-like equation [36, 37]:

πτ=12σ2S22πS2+r(SπS-π)+v(τ)(SπS-π),    (34)

where

v(τ)=(r-α)f(τ)σ-f(τ),    (35)

with the Hamiltonian operator

Ȟ=12σ2Ť+rP^+v(τ)P^,    (36)

where

Ť=S22S2,    (37)

and

P^=(SS-Ǐ).    (38)

When the amplitude of the bubble f is zero, the potential function v(τ) is also zero, and the Hamiltonian reduces to

Ȟ0=12σ2Ť+rP^,    (39)

which gives the evolution of the usual free Black–Scholes model as follows:

πτ=12σ2S22πS2+r(SπS-π).    (40)

Note that equation (34) can be written as follows:

πτ=Ȟ0π+v(τ)P^π.    (41)

Due to the fact that P^ commutes with Ť and Ȟ0, equation (41) can be integrated to give

π(S,τ)=eȞ0τ+(0τv(τ)dτ)P^ π(S,0)=eȞ0τ+(0τv(τ)dτ)P^ Φ(S),    (42)

where π(S, 0) = Φ(S) is the contract function.

Now, we define the arbitrage number AN(τ) at time τ associated with the arbitrage bubble f(τ) using the integral

AN(τ)=0τv(τ)dτ=0τ(r-α)f(τ)σ-f(τ)dτ,    (43)

which represents the accumulated potential between 0 and τ, so

π(S,τ)=eȞ0τ+AN(τ)P^ Φ(S),    (44)

that is,

π(S,τ)=eAN(τ)P^eȞ0τ Φ(S).    (45)

If we denote by C(S, τ, r) the solution to the free (f = 0) Black–Scholes equation for a contract function Φ(S) (call, binary, etc.), which evolves with the free Hamiltonian (equation (39)) characterized by a constant interest rate r, we can write generically

C(S,τ,r)=eȞ0τ Φ(S).    (46)

Thus, the solution π(S, τ) to the interacting equation is as follows:

π(S,τ)=eAN(τ)P^C(S,τ,r).    (47)

By using the expansion given in Appendix for eAN(τ)P^ [see equations (136, 141)], the last expression can be written as follows:

π(S,τ)=n=0e-AN(τ)Qn(AN(τ))Sn n SnC(S,τ,r),    (48)

where (see Appendix)

Q0(x)=1,    (49)
Q1(x)=ex-1,    (50)
Qj(x)=m=jαm,jxmm!   j=2,3,4,,    (51)

and the coefficients αm,j are given by the recurrence relation

αn,1=1,      αn,n=1,      andαnm=mαn-1,m+αn-1,m-1.    (52)

Equation (48) makes it possible to write the solution π(S, t) of the interacting Black–Scholes equation in terms of all the derivatives of the free Black–Scholes solution C(S, τ, r), that is, in terms of all its Greeks.

4. Arbitrage theorems

The arbitrage theorems can be stated as follows:

First arbitrage theorem: Let f(τ) be a pure time-dependent arbitrage bubble that acts in the time interval 0 < τ < T. If the arbitrage number of f at time τ = T is zero, that is, if AN(T) = 0, then

π(S,T)=C(S,T,r).    (53)

This means that the solution π(S, T) of the interacting Black–Scholes equation (equation (34)) evaluated at maturity has the same value as the solution C(S, T, r) of the free Black–Scholes equation (equation (6)) at maturity only if the accumulated potential at time T is exactly zero. This fact is far from being a trivial conclusion, because the bubble f = f(τ) and its associated potential v(τ) in equation (35) are non-trivial in the time interval 0 < τ < T.

Proof: Consider the option price at time τ = T in equation (48), that is,

π(S,T)=n=0e-AN(T)Qn(AN(T))Sn n SnC(S,T,r).    (54)

It is because AN(T) = 0,

π(S,T)=n=0Qn(0)Sn n SnC(S,T,r),    (55)

and by noting that the functions Qn(x) satisfy Qn(0) = 0 for n = 1, 2, 3, ... and Q0(x) = 1, we have

π(S,T)=C(S,T,r).    (56)

Second arbitrage theorem: Let f(τ) be an arbitrage bubble that acts in the time interval 0 < τ < T. If the arbitrage number of f at time τ = T is non-zero, that is, if AN(T) ≠ 0, then the solution of the interacting Black–Scholes equation depends on all the Greeks of the free Black–Scholes solution and the bubble's behavior for 0 < τ < T, and

π(S,T)=n=0e-AN(T)Qn(AN(T))Sn n SnC(S,T,r).    (57)

Note that in this case, the solution of the interacting Black–Scholes equation [equation (34)] differs from the solution of the free equation [equation (40)] at maturity. In fact, it is a series of all the free Greeks. Here, the interacting solution's dynamical behavior is completely different from that given by the free dynamics.

Proof: Evaluate π(S, τ) in equation (48) at τ = T (see Appendix for details).

Third arbitrage theorem: Let f(τ) be an arbitrage bubble that acts in the time interval 0 < τ < T, and let AN(τ) be the arbitrage number of f at time τ. Then, the solution π(S, τ) of the interacting Black–Scholes equation [equation (34)] is just the solution to a free Black–Scholes equation with the variable interest rate

r(τ)=r+1τAN(τ),

that is,

π(S,τ)=C(S,τ,r(τ)),    (58)

and

C(S,τ,r+1τAN(τ))=n=0e-AN(τ)Qn(AN(τ))Sn n SnC(S,τ,r).    (59)

This theorem gives a closed-form solution of the interacting Black–Scholes equation (equation (34)) as a composition of the free Black–Scholes solution C(S, t, r) and a time-dependent rate r(τ)=r+1τAN(τ), which depends on the accumulated potential.

Proof: Consider equation (44), which is equivalent to

π(S,τ)=eȞ0τ+1τAN(τ)P^τ Φ(S)=e(12σ2Ť+rP^)τ+1τAN(τ)P^τ Φ(S),    (60)

that is,

π(S,τ)=e(12σ2Ť+(r+1τAN(τ))P^)τ Φ(S)=e(12σ2Ť+r(τ)P^)τ Φ(S)    (61)
=C(S,τ,r(τ)),    (62)

where r(τ)=r+1τAN(τ). By comparing equation (61) with equation (48), we obtain equation (59).

4.1. Risk-neutral valuation

Now, one can ask what happen with the risk-neutral valuation theorem in this case, that is, is the process

z(t)=π(S(t),t)B(t)    (63)

a martingale or not? In the usual arbitrage-free case, where π(S, t) = C(S, t, r) is the solution of the free Black–Scholes equation (6) for some contingent claim Φ(S), the process

z(t)=C(S(t),t,r)B(t),    (64)

is a martingale in the Q measure induced by the Feynman–Kac theorem

C(s,t,r)=e-r(T-t)Es,tQ[Φ(S(T))],    (65)

where the process S(u) for tuT satisfies the Q dynamics

dS=rSdu+σSdW,      S(t)=s,    (66)

and the bond B(t) satisfies equation (1). So, what would be the risk-neutral Q measure for the process in (63) for the interacting Black–Scholes model? A possible answer comes from equation (28), where one can see that the arbitrage bubble f(t) changes the constant interest r in the Black–Scholes equation to a time-dependent one r+v(t)=r+(r-α)f(t)σ-f(t). Heuristically, this would imply that the Q dynamics for S(t) in equation (66) would incorporate this change. In fact, the bond dynamics (29) in this case is

dB(t)=(r+v(t))B(t)dt.    (67)

Thus, one can state the following theorem:

Risk-neutral valuation theorem in presence of arbitrage bubbles

Define the Q dynamics for the underlying asset S in the presence of arbitrage bubble f(t) by the equation

dS(t)=(r+v(t))S(t)dt+σS(t) dW(t),    (68)

and consider the bond dynamics B(t) given by (67). Then, the process (63) is a martingale in this Q measure.

Proof: Consider the process

z(t)=G(π(t),B(t))=π(t)B(t),    (69)

with π(t) = π(S(t), t), then by the Itô lemma

dz(t)=Gπ dπ+GB dB+122Gπ2 dπ2+122GB2 dB2+2Gπ B dB dπ,    (70)

and using for in equation (20) the Q dynamics (68), one has that

dπ=(πt+(r+v(t))SπS+12σ2S22πS2) dt+SπSσ dW.    (71)

Now the interacting Black–Scholes equation (26) implies that

dπ=((r+v(t))π) dt+SπSσ dW.    (72)

Replacing (72) and (67) in equation (70) one gets, by noting that dB2 = 0, dB dπ = 0, Gπ = 1B, GB=-πB2, 2Gπ2=0, that

dz=Gπ dπ+GB dB=1B dπ-πB2 dB,    (73)

that is,

dz=1B((r+v(t))π dt+SπSσ dW)-πB2((r+v(t))B dt),    (74)

or

dz=πB(r+v(t)) dt+SBπSσ dW-πB(r+v(t)) dt=SBπSσ dW,    (75)

so dz has zero drift, which implies that E(z(t)) is time-independent, i.e., z(t) is a martingale.

5. Some examples

Example 1: To illustrate the consequences of the aforementioned theorems, we consider as the first example a simple, single square-time bubble [equation (30)] acting on a call option [29]. In this case, f(τ) is

f(τ)={00ττ1,f0τ1<τ<τ2,0τ2τT,    (76)

where τ1 = TT2 and τ2 = TT1. The potential function [equation (35)] for this bubble is

v(τ)={00ττ1,v0=(r-α)f0σ-f0,τ1<τ<τ2,0τ2τT,    (77)

and the non-zero arbitrage number for this square bubble is

AN(T)=0Tv(τ)dτ=τ1τ2(r-α)f0σ-f0 dτ=v0(τ2-τ1).    (78)

When r < α and f0 > σ, we speak of a positive bubble because the arbitrage number is positive. For f0 < σ, the arbitrage number becomes negative, so we speak of a negative bubble. For the case r > α, we have a positive bubble for f0 < σ and a negative bubble for f0 > σ. Note that a single positive or negative square bubble has an arbitrage number AN(T) ≠ 0. The interacting solution is given by [see equation (61)]

π(S,τ)={C(S,τ,r)0ττ1,C(S,τ,r+v0(τ-τ1)τ)τ1<τ<τ2,C(S,τ,r+v0(τ2-τ1)τ)τ2τT,    (79)

where, in this case, C(S, τ, r) is the corresponding call solution of the Black–Scholes equation.

Figure 1 shows the interacting solution (equation (79)) in terms of the time t for a call with K = 5, r = 0.1, σ = 0.3, T = 1, α = 0.2, τ1 = 0.4, and τ2 = 0.5 for f0 = 0.3045 (left figure) and f0 = 0.2955 (right figure). The curve on the left side of both figures at t = T = 1 or τ = 0 is the usual call contract function Φ(S) = max{0, SK}.

FIGURE 1
www.frontiersin.org

Figure 1. Interacting Black–Scholes solution π(S, t) for a call option with f0 = 0.3045 (left) and f0 = 0.2955 (right). Note that the interacting solution goes up when f0 > σ (left) and goes down when f0 < σ (right). In both cases, the option at maturity is very different from the free Black–Scholes solution C(S, t, r) due to the fact that the arbitrage number is non-zero.

Example 2: Consider the arbitrage process generated successively by a pair of positive and negative square bubbles given by

f(τ)={00ττ1,f0τ1<τ<τ2,f1τ2<τ<τ3,0τ3τT.    (80)

If the total arbitrage number of the bubble is zero, then we speak of a bubble–antibubble pair.

Figure 2 shows a bubble–antibubble pair f(t) for T = 1, τ1 = 0.4, τ2 = 0.5, and τ3 = 0.6. Here, f0 = 0.285 for the negative bubble, f1 = 0.3167 for the positive bubble, and σ = 0.3.

FIGURE 2
www.frontiersin.org

Figure 2. (Upper left) Bubble–antibubble pair as a function of time t = T − τ (blue curve) and the volatility σ (red curve). Note that the positive and negative bubbles are equidistant with respect to σ. (Upper right) Potential v in terms of t. (Lower figure) Arbitrage number in terms of t. Note that after the bubble–antibubble pair has finished acting, the arbitrage number is zero.

Figure 3 shows the effect of a bubble–antibubble pair on a call option with parameters K = 5, r = 0.1, and α = 0.2. The curve at τ = 0 (t = T = 1) is the usual contract function Φ(S) = max{0, SK}. For 0 < τ ≤ τ1, the dynamics are given by the free Black–Scholes solution. For τ1 ≤ τ ≤ τ2, the bubble acts and changes the dynamics to create an interacting solution. However, due to the fact that the arbitrage number of the bubble–antibubble pair is zero, for τ2 ≤ τ ≤ T, the interacting dynamics are again given by the free Black–Scholes model. Thus, when the arbitrage number is zero, the bubble has no effect on the option at maturity.

FIGURE 3
www.frontiersin.org

Figure 3. (Left) Effect of a bubble–antibubble pair on a call option. Note that for τ2 ≤ τ ≤ T, the option recovers its free initial dynamics due to the fact that the bubble arbitrage number is 0.

Example 3: Consider the case of three successive square bubbles of the form,

f(τ)={00ττ1,f0τ1<τ<τ2,f1τ2<τ<τ3,f2τ3<τ<τ4,0τ4τT.    (81)

Figure 4 shows the effect of two negative bubbles and one positive bubble (left figure), as well as that of one positive bubble and two negative bubbles (right figure), on a call option. In all cases, the total arbitrage number is zero. Here, K = 5, r = 0.1, σ = 0.3, T = 1, α = 0.2, τ1 = 0.2, τ2 = 0.4, τ3 = 0.5, and τ4 = 0.6. For the left figure, f0 = 0.24, f1 = 0.276, and f2 = 0.316, whereas for the right figure, f0 = 0.33, f1 = 0.24, and f2 = 0.28.

FIGURE 4
www.frontiersin.org

Figure 4. (Left) Effect of two negative bubbles and one positive bubble, respectively. (Right) Effect of one positive bubble and two negative bubbles. In all cases, the total arbitrage number is zero. Thus, the dynamics of the option at maturity are the usual free Black–Scholes dynamics.

Figure 5 shows the effect on a call option of one positive bubble, one negative bubble, and one positive bubble, respectively (left figure), and one negative bubble, one positive bubble, and one negative bubble, respectively (right figure). In all cases, the total arbitrage number is zero. Here, K = 5, r = 0.1, σ = 0.3, T = 1, α = 0.2, τ1 = 0.2, τ2 = 0.3, τ3 = 0.5, and τ4 = 0.6. For the left figure, f0 = 0.345, f1 = 0.273, and f2 = 0.326, whereas for the right figure, f0 = 0.27, f1 = 0.33, and f2 = 0.278.

FIGURE 5
www.frontiersin.org

Figure 5. (Left) One positive bubble, one negative bubble, and one positive bubble. (Right) One negative bubble, one positive bubble, and one negative bubble. The total arbitrage number is zero. Note again that the dynamics of the option at maturity are the usual free Black–Scholes dynamics.

Figure 6 shows the effect on a call option of one positive bubble and two negative bubbles (left figure) and two negative bubbles and one positive bubble (right figure), but in this case, the total arbitrage number AN ≠ 0. The parameters are the same as those used in Figure 5, but for the left figure, f0 = 0.345, f1 = 0.273, f2 = 0.294, and AN = −0.616, whereas for the right figure, f0 = 0.273, f1 = 0.285, f2 = 0.3015, and AN = 2.01. In this situation, the option dynamics at maturity are far away from the free dynamics. Thus, the bubble action has a net effect on the option price when AN ≠ 0.

FIGURE 6
www.frontiersin.org

Figure 6. (Left) Effect of one positive bubble and two negative bubbles for a total arbitrage number AN = −0.616. (Right) Effect of two negative bubbles and one positive bubble but with AN = 1.529. Note that after the bubble acts, the evolution of the option does not follow the dynamics of the free model; instead, its behavior changes radically.

6. Conclusion

Three theorems about the exact dynamical behavior of the option price when time-dependent arbitrage bubbles are incorporated explicitly in the Black–Scholes equation have been presented. These arbitrage bubbles (which can act at some instant t between 0 and the maturity T) can be characterized by an arbitrage number AN(τ) that corresponds to the accumulated external potential from time 0 to time τ. If AN(t = 0) = 0 or AN(τ = T) = 0, independent of the bubble's dynamic behavior, the option's values at maturity are given by the usual Black–Scholes dynamics with an interest rate r. In some sense, the option does not remember the past arbitrage process. If AN(t = 0) ≠ 0 or AN(τ = T) ≠ 0, the option price depends on the bubble and on all the higher derivatives (or Greeks) of the solution of the free Black–Scholes equation with an interest rate r, and it is equal to a solution of the free Black–Scholes equation with an interest rate r + (1/T)AN(τ = T). In this case, the option remembers the past arbitrage processes.

Thus, arbitrage bubbles with AN(t = 0) = 0 are innocuous because the system returns to the initial equilibrium state. However, bubbles with AN(t = 0) ≠ 0 are dangerous from a financial point of view, because these bubbles change the initial equilibrium trajectory of the option price.

We hope that these theorems will make it easier to understand an option's dynamic evolution when arbitrage processes are included. A generalization of these theorems and their consequences, for the case of a price-dependent bubble f = f(S, t), will be studied in future research.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.

Author contributions

MC and RO: conceptualization, investigation and formal analysis, visualization, writing—original draft, and editing. All authors contributed to the article and approved the submitted version.

Acknowledgments

The authors would like to thank the Industrial Engineering School of Diego Portales University for its economic support of the article published in this journal.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fams.2023.1138663/full#supplementary-material

Footnotes

1. ^Using a martingale approach, other authors have investigated the effect of bubbles on the valuation of the underlying asset and derivative assets. See, for example, Jarrow and Kwok [30] and Benth et al. [31].

2. ^Note that C(S, t = 0, r) is the solution of the equilibrium Black–Scholes equation for any derivative asset.

References

1. Black F, Scholes M. The pricing of options and corporate liabilities. In: Galai D, Crouhy M, Wiener Z, editors. World Scientific Reference on Contingent Claims Analysis in Corporate Finance: Volume 1: Foundations of CCA and Equity Valuation. Singapore: World Scientific (2019), p. 3–21. doi: 10.1142/9789814759588_0001

CrossRef Full Text | Google Scholar

2. Merton RC. Theory of rational option pricing. Bell J. Econ. Manage. Sci. (1973) 4:141–83. doi: 10.2307/3003143

CrossRef Full Text | Google Scholar

3. Heston SL. A closed-form solution for options with stochastic volatility with applications to bond and currency options. Rev Financ Stud. (1993) 6:327–43. doi: 10.1093/rfs/6.2.327

CrossRef Full Text | Google Scholar

4. Gatheral J. The Volatility Surface. Hoboken, NJ: Wiley (2006).

Google Scholar

5. Wilmott P. Paul Wilmott on Quantitative Finance. Hoboken, NJ: Wiley (2000).

Google Scholar

6. Dupire B. Pricing with a smile. Risk. (1994) 7:18–20.

Google Scholar

7. Björk T. Arbitrage Theory in Continuous Time. Oxford: Oxford University Press (2009).

Google Scholar

8. Ho TS, Lee SB. Term structure movements and pricing interest rate contingent claims. J Finance. (1986) 41:1011–29. doi: 10.1111/j.1540-6261.1986.tb02528.x

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Tankov P. Financial Modelling with Jump Processes. Boca Raton, FL: Chapman and Hall/CRC (2003). doi: 10.1201/9780203485217

CrossRef Full Text

10. Hull J. Options, Futures, and Other Derivative Securities, 2ed. Hoboken, NJ: Prentice Hall (1993).

11. Sofianos G. Index arbitrage profitability. J Deriv. (1993) 1:6–20. doi: 10.3905/jod.1993.407871

CrossRef Full Text | Google Scholar

12. Lo AW, MacKinlay AC. A Non-Random Walk Down Wall Street. Princeton, NJ: Princeton University Press (1999).

Google Scholar

13. Ilinski K. Physics of Finance: Gauge Modelling in Non-Equilibrium Pricing. Hoboken, NJ: John Wiley & Sons Inc (2001).

Google Scholar

14. Adrian T. Inference, arbitrage, and asset price volatility. J Financ Intermed. (2009) 18:49–64. doi: 10.1016/j.jfi.2008.06.001

CrossRef Full Text | Google Scholar

15. Cornell B, French KR. Taxes and the pricing of stock index futures. J Finance. (1983) 38:675–94. doi: 10.1111/j.1540-6261.1983.tb02496.x

CrossRef Full Text | Google Scholar

16. Figlewski S. Hedging performance and basis risk in stock index futures. J Finance. (1984) 39:657–69. doi: 10.1111/j.1540-6261.1984.tb03654.x

CrossRef Full Text | Google Scholar

17. Stoll HR, Whaley RE. Program trading and expiration-day effects. Financ Anal J. (1987) 43:16–28. doi: 10.2469/faj.v43.n2.16

CrossRef Full Text | Google Scholar

18. Brennan MJ, Schwartz ES. Arbitrage in stock index futures. J Bus. (1990) 63:S7–S31. doi: 10.1086/296491

CrossRef Full Text | Google Scholar

19. Ilinski K. Virtual arbitrage pricing theory. arXiv. (1999) [preprint]. doi: 10.48550/arXiv.cond-mat/9902045

CrossRef Full Text | Google Scholar

20. Ross SA. The arbitrage pricing theory of capital asset pricing. J Econ Theory. (1976) 13:334–60. doi: 10.1016/0022-0531(76)90046-6

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Elton EJ, Gruber MJ. Modern Portfolio Theory and Investment Analysis. Hoboken, NJ: John Wiley & Sons (1995).

22. Ilinski K. Physics of finance. In: Kertesz J, Kondor I, editors. Econophysics: An Emerging Science. Philadelphia, PA: Kluwer (1998).

Google Scholar

23. Callen HB. Thermodynamics. Hoboken, NJ: John Wiley & Sons (1960).

Google Scholar

24. Ilinski K. How to account for virtual arbitrage in the standard derivative pricing. arXiv. (1999) [preprint]. doi: 10.48550/arXiv.cond-mat/9902047

CrossRef Full Text | Google Scholar

25. Ilinski K, Stepanenko A. Derivative pricing with virtual arbitrage. arXiv. (1999) [preprint]. doi: 10.48550/arXiv.cond-mat/9902046

CrossRef Full Text | Google Scholar

26. Panayides S. Arbitrage opportunities and their implications to derivative hedging. Phys A: Stat Mech Appl. (2006) 361:289–96. doi: 10.1016/j.physa.2005.06.077

CrossRef Full Text | Google Scholar

27. Panayides S, Fedotov S. Stochastic arbitrage return and its implication for option pricing. Phys A: Stat Mech Appl. (2005) 345:207–17. doi: 10.1016/S0378-4371(04)00989-6

CrossRef Full Text | Google Scholar

28. Papanicolaou GC, Sircar KR. Stochastic volatility, smiles and asymptotics. Appl Math Finance. (1999) 6:107–45. doi: 10.1080/135048699334573

CrossRef Full Text | Google Scholar

29. Contreras M, Montalva R, Pellicer R, Villena M. Dynamic option pricing with endogenous stochastic arbitrage. Phys A: Stat Mech Appl. (2010) 389:3552–64. doi: 10.1016/j.physa.2010.04.019

CrossRef Full Text | Google Scholar

30. Jarrow RA, Kwok SS. Inferring Financial Bubbles from Option Data. School of Economics, University of Sydney (2020). Available online at: https://ideas.repec.org/p/syd/wpaper/2020-04.html

Google Scholar

31. Benth FE, Crisan D, Guasoni P, Manolarakis K, Muhle-Karbe J, Nee C, et al. A mathematical theory of financial bubbles. In: Henderson V, Sircar R, editors. Paris-Princeton Lectures on Mathematical Finance 2013. New York, NY: Springer (2013), p. 1–108. doi: 10.1007/978-3-319-00413-6

CrossRef Full Text | Google Scholar

32. Heston SL, Loewenstein M, Willard GA. Options and bubbles. Rev Financ Stud. (2007) 20:359–90. doi: 10.1093/rfs/hhl005

CrossRef Full Text | Google Scholar

33. Jarrow RA. Asset price bubbles. Ann Rev Financ Econ. (2015) 7:201–18. doi: 10.1146/annurev-financial-030215-035912

CrossRef Full Text | Google Scholar

34. Contreras M, Pellicer R, Santiagos D, Villena M. Calibration and simulation of arbitrage effects in a non-equilibrium quantum Black-Scholes model by using semiclassical methods. J Math Finance. (2016) 6:541–61. doi: 10.4236/jmf.2016.64042

CrossRef Full Text | Google Scholar

35. Contreras M, Pellicer R, Villena M, Ruiz A. A quantum model of option pricing: when Black-Scholes meets Schrödinger and its semi-classical limit. Phys A: Stat Mech Appl. (2010) 389:5447–59. doi: 10.1016/j.physa.2010.08.018

CrossRef Full Text | Google Scholar

36. Merzbacher E. Quantum Mechanics. Hoboken, NJ: Wiley. (1998).

Google Scholar

37. Sakurai J, Napolitano J. Modern Quantum Mechanics, 2nd ed. San Francisco, CA: Pearson. (2011).

Google Scholar

Keywords: Black–Scholes model, option pricing, derivatives, arbitrage, quantum mechanics

Citation: Contreras G. M and Ortiz H. R (2023) Three little arbitrage theorems. Front. Appl. Math. Stat. 9:1138663. doi: 10.3389/fams.2023.1138663

Received: 06 January 2023; Accepted: 17 March 2023;
Published: 21 April 2023.

Edited by:

Indranil SenGupta, North Dakota State University, United States

Reviewed by:

Emilio Barucci, Polytechnic University of Milan, Italy
Savin Treanta, Polytechnic University of Bucharest, Romania

Copyright © 2023 Contreras G. and Ortiz H. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Roberto Ortiz H., roberto.ortiz@udp.cl

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.