Next Article in Journal
The Effects of Monitoring Activities on Loan Defaults in Group-Based Lending Program: Evidence from Vietnam
Next Article in Special Issue
A Double Optimum New Solution Method Based on EVA and Knapsack
Previous Article in Journal
A Re-Appraisal of the Role of Monetary Policy: The Quantity Theory of Money through a Structural Vector Autoregressive Approach
Previous Article in Special Issue
Enhancing Model Selection by Obtaining Optimal Tuning Parameters in Elastic-Net Quantile Regression, Application to Crude Oil Prices
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Rainbow Step Barrier Options

by
Tristan Guillaume
Laboratoire Thema, CYU Cergy Paris Université, 33 Boulevard du Port, F-95011 Cergy-Pontoise Cedex, France
J. Risk Financial Manag. 2024, 17(8), 356; https://doi.org/10.3390/jrfm17080356
Submission received: 13 June 2024 / Revised: 30 July 2024 / Accepted: 5 August 2024 / Published: 13 August 2024
(This article belongs to the Special Issue Featured Papers in Mathematics and Finance)

Abstract

:
This article provides exact analytical formulae for various kinds of rainbow step barrier options. These are highly flexible and sophisticated multi-asset barrier options based on the following principle: the option life is divided into several time intervals on which different barriers are monitored w.r.t. different underlying assets. From a mathematical point of view, new results are provided for the first passage time of a multidimensional geometric Brownian motion to a boundary defined as a step function. The article shows how to implement the obtained option valuation formulae in a simple and very efficient manner. Numerical results highlight a strong sensitivity of rainbow step barrier options to the correlations between the underlying assets.

1. Introduction

Barrier options are characterised by the introduction in the option contract of a parameter called the “barrier”, which, in the most standard form, is a predefined reference value of the underlying asset S that may be located above the spot value of the underlying (upward barrier or “up-barrier”) or below it (downward barrier or “down-barrier”). The barrier may be of a “knock-out” type, i.e., the option expires worthless if the barrier is hit by S at any time during the option life (in which case the barrier is called “continuous”) or at one or a few predefined times (in which case the barrier is called “discrete”). The first passage time of an underlying to a knock-out barrier, before the option expiry, triggers the “deactivation” of the option. Alternatively, the barrier may be of a “knock-in” type, i.e., the option expires worthless unless the barrier has been hit at least once before expiry, an event called “activation”.
Barrier options are the oldest and the most widely traded non-vanilla options. They are embedded in a lot of popular structured derivatives in stock and interest rate markets (see, e.g., Bouzoubaa and Osseiran 2010). They are also extensively used as analytical tools in financial modelling, for instance, in the so-called “structural models” of default risk (see, e.g., Bielecki and Rutkowski 2004) or in the valuation of investments (theory of “real options”). Since their first appearance in the financial markets during the 1970s, there have been a huge number of variations in their original payoff, leading to an extraordinary variety of non-standard barrier options. Among the most well-known of them are the step barrier options, which divide the option lifetime into several time intervals on which the barrier takes on different values. In its standard form, a step barrier option features a piecewise constant barrier, i.e., a barrier defined as a step function. This allows modulation of the barrier level during the option life, thus offering increased flexibility and enhanced risk management capabilities, relative to a traditional barrier option. For instance, up-and-out barrier levels can be raised and down-and-out barrier levels can be lowered during time intervals when more protection is required, thus reducing the risk of deactivation. Likewise, up-and-in barrier levels can be lowered and down-and-in barrier levels can be raised in time intervals during which the implicit volatility of the underlying rises, thus, increasing the chances of activation.
The exact analytical valuation of step barrier functions was first achieved by Guillaume (2001). Further mathematical details on how closed-form valuation can be achieved, as well as exact results for more general deterministic step barriers, are provided in Guillaume (2015, 2016). In the last couple of years, there has been a renewed interest in step barrier options, as they stand out as an essential component of innovative forms of investment, such as autocallable structured products and other equity-linked products. This has led to a new series of academic contributions, such as Lee et al. (2019a, 2019b), Lee et al. (2021), Lee et al. (2022). Unlike previously cited references, these articles do not provide explicit formulae, except for a few simple and already known cases, nor do they handle outside step barriers as in Guillaume (2001), alternatively upward and downward steps as in Guillaume (2015) and exponentially moving step barriers as in Guillaume (2016). These recent contributions ignore previous results given in references they do not cite, such as Guillaume (2001), that actually solve the problems they discuss. They also claim to be able to analytically value a step barrier option with an arbitrary number of steps, but without explaining how they intend to solve the difficult problem known in numerical integration as the “curse of dimensionality”, nor even beginning to discuss the numerical implementation of their approach, which constitutes the main issue, though.
There are still a number of unsolved problems related to the valuation of step barrier options. In particular, multi-asset step barrier options are barely touched upon in the existing literature, apart from an isolated formula for an “outside” step barrier option given in Guillaume (2001), also called an “external” step barrier option, featuring one underlying asset w.r.t. which barrier crossing is monitored and another underlying asset w.r.t. which the moneyness of the option is measured at expiry (the reader may refer to Heynen and Kat 1994, or to Kwok et al. 1998, for background on outside barrier options in general). Yet, multi-asset contracts with step barriers are actively traded in today’s financial markets, as they allow investors to benefit from the advantages of diversification in terms of risk control and expansion of investment opportunities. A particularly important subset of these contracts is the so-called “rainbow” step barrier option. Broadly speaking, in the literature on options, the denomination “rainbow” applies to payoffs linked to the performances of two or more underlying assets (Chang et al. 2005; Gao and Wu 2022); metaphorically, each underlying represents a different color, so that the association of all of these factors makes up a rainbow. In the realm of barrier options, the rainbow step barrier option is characterised by the property that, at each time interval, the barrier is monitored w.r.t. a different underlying asset. A contract featuring a number n of underlying assets associated with n time intervals t 0 = 0 , t 1 , …, t n 1 , t n , on which n steps of a piecewise constant barrier are monitored, is called an n the colour rainbow barrier option. Each time interval is matched with a specific step of the barrier and a specific underlying asset. In the standard form of the contract, it is the n th asset associated with the n th last step of the barrier that is used to determine the moneyness of the contract at expiry. Rainbow step barrier options are typically priced by Monte Carlo simulation, even in a standard Black–Scholes model, because of the difficulties of the entailed analytical calculations and also because the dimension of the valuation problem quickly increases with the number of “colors”, leading to non-trivial issues of numerical evaluation of high-dimensional integrals. Due to these obstacles, the present article is restricted to two-colour rainbow step barrier options. Closed-form valuation is achieved not only for standard two-colour contracts but also for two-colour outside step barrier options involving a third correlated asset at expiry, and for two-colour contracts featuring a two-sided barrier (also known as a double barrier), i.e., both an upward and a downward barrier on each time interval. Numerical evaluation of the obtained analytical solutions is dealt with so that the valuation formulae derived in this paper can be immediately implemented and yield extremely accurate results in a few tenths of one second. Numerical results are provided, which reveal a strong and stable dependency of rainbow step barrier options on the correlations between the underlying assets, as well as the importance of the volatility of the asset used to measure moneyness when dealing with rainbow outside step barrier options. These findings suggest clear application ideas to traders and investors, whether for a hedging or speculative purposes. They highlight the specificity of rainbow step barrier options as instruments highly sensitive to correlation, in contrast to standard step barrier options, which are sensitive to volatility but, by constuction, cannot be sensitive to correlation.
This article is organised as follows: Section 2 states the main analytical results, provides numerical results, and discusses their implications; Section 3 gives the mathematical proofs of the analytical results presented in Section 2.

2. Formulae and Numerical Results

Let us begin with a few definitions. Let S 1 and S 2 be two GBMs (geometric Brownian motions) modelling two asset prices, whose differentials, under a given probability measure P , are given by:
d S 1 t = v 1 S 1 t d t + σ 1 S 1 t d B 1 t
d S 2 t = v 2 S 2 t d t + σ 2 S 2 t d B 2 t
where v 1 , v 2 , σ 1 , σ 2 + , and B 1 and B 2 are two standard Brownian motions whose correlation coefficient is denoted by ρ 1.2 .
The measure P is characterised by the pair v 1 , v 2 or, equivalently, by the pair:
μ 1 = v 1 σ 1 2 / 2 , μ 2 = v 2 σ 2 2 / 2
If we refer to the log-return processes X i t = ln S i t / S i 0 , i 1 , 2 , whose differentials under P are given by:
d X i t = μ i d t + σ i d B i t
Let H 1 , H 2 , K 1 , K 2 , K 3 be positive real numbers. The numbers H 1 , H 2 are the values of two knock-out continuous barriers. H 1 is monitored w.r.t. S 1 on a time interval t 0 = 0 , t 1 , while H 2 is monitored w.r.t. S 2 on a time interval t 1 , t 2 . The numbers K 1 , K 2 , K 3 are the values of three discrete knock-out barriers; K 1 is monitored w.r.t. S 1 at time t 1 , while K 2 and K 3 are monitored w.r.t. S 2 at times t 1 and t 2 , respectively.
We can now begin to value two-colour step barrier options in the following order:
-
both steps either upward or downward (Section 2.1);
-
one upward step and one downward step (Section 2.2);
-
reverse-type contract (Section 2.3);
-
outside or external two-colour step barrier (Section 2.4);
-
two-colour step double barrier (Section 2.5).

2.1. Valuation of Two-Colour Step Barrier Options When the Steps of the Barrier Are on the Same Side in Each Time Interval

Section 2.1 deals with the valuation of two-colour step barrier options when the steps of the barrier are either both upward or both downward. Our objective is to find the value of the joint cumulative distribution function P RUU μ 1 , μ 2 defined by:
P RUU μ 1 , μ 2 P sup 0 t t 1 S 1 t H 1 , S 1 t 1 K 1 , S 2 t 1 K 2 , sup t 1 t t 2 S 2 t H 2 , S 2 t 2 K 3
where the acronym “ RUU ” stands for “Rainbow Up and Up”.
The main result of Section 2.1 is given by the following Proposition 1.
Proposition 1. 
The exact value of  P RUU μ 1 , μ 2  is given by:
P RUU μ 1 , μ 2 = N 3 min k 1 , h 1 μ 1 t 1 σ 1 t 1 , min k 2 , h 2 μ 2 t 1 σ 2 t 1 , min k 3 , h 2 μ 2 t 2 σ 2 t 2 ; θ 1.2 , θ 1.3 , θ 2.3
exp 2 μ 1 h 1 σ 1 2 × N 3 min k 1 , h 1 2 h 1 μ 1 t 1 σ 1 t 1 , min k 2 , h 2 μ 2 t 1 σ 2 t 1 2 θ 1.2 h 1 σ 1 t 1 , min k 3 , h 2 μ 2 t 2 σ 2 t 2 2 θ 1.2 h 1 σ 1 t 2 ; θ 1.2 , θ 1.3 , θ 2.3
exp 2 μ 2 h 2 σ 2 2 × N 3 min k 1 , h 1 μ 1 t 1 σ 1 t 1 + 2 θ 1.2 μ 2 t 1 σ 2 , min k 2 , h 2 + μ 2 t 1 σ 2 t 1 , min k 3 , h 2 2 h 2 μ 2 t 2 σ 2 t 2 ; θ 1.2 , θ 1.3 , θ 2.3
+ exp 2 μ 1 σ 1 2 4 μ 2 θ 1.2 σ 1 σ 2 h 1 + 2 μ 2 h 2 σ 2 2 × N 3 min k 1 , h 1 2 h 1 μ 1 t 1 σ 1 t 1 + 2 θ 1.2 μ 2 t 1 σ 2 , min k 2 , h 2 + μ 2 t 1 σ 2 t 1 2 θ 1.2 h 1 σ 1 t 1 , min k 3 , h 2 2 h 2 μ 2 t 2 σ 2 t 2 + 2 θ 1.2 h 1 σ 1 t 2 ; θ 1.2 , θ 1.3 , θ 2.3
where the  μ i s  are given by (3) and:
-
N 3 b 1 , b 2 , b 3 ; c 12 , c 13 , c 23  is the trivariate standard normal cumulative distribution function with correlation coefficients  c 12 , c 13 , c 23
h 1 = ln H 1 S 1 0 , h 2 = ln H 2 S 2 0 , k 1 = ln K 1 S 1 0 , k 2 = ln K 2 S 2 0 , k 3 = ln K 3 S 2 0
θ 1.2 = ρ 1.2 , θ 1.3 = t 1 t 2 ρ 1.2 , θ 2.3 = t 1 t 2
Corollary 1. 
It suffices to multiply by  1  all the first three arguments of each  N 3 . , . , . ; . , . , .  function and to substitute each  min  operator by a  max  operator in Proposition 1 to obtain an exact formula for  P RDD μ 1 , μ 2  defined as:
P RDD μ 1 , μ 2 P inf 0 t t 1 S 1 t H 1 , S 1 t 1 K 1 , S 2 t 1 K 2 , inf t 1 t t 2 S 2 t H 2 , S 2 t 2 K 3
where the acronym “ RDD ” stands for “Rainbow Down and Down”.
Corollary 2. 
The term numbered (9) in Proposition 1 gives the value of the corrresponding knock-in probability denoted by P RUU I μ 1 , μ 2  and defined by:
P RUU I μ 1 , μ 2 P sup 0 t t 1 S 1 t H 1 , S 1 t 1 K 1 , S 2 t 1 K 2 , sup t 1 t t 2 S 2 t H 2 , S 2 t 2 K 3
Corollary 3. 
It suffices to substitute each argument  θ 2.3  in each  N 3 . , . , . ; . , . , .  function of Proposition 1 by  t 2 t 3 ,  t 3 t 2 , to obtain an exact formula for the early-ending variant  P EERUU μ 1 , μ 2  defined by:
P EERUU μ 1 , μ 2 P sup 0 t t 1 S 1 t H 1 , S 1 t 1 K 1 , S 2 t 1 K 2 , sup t 1 t t 2 S 2 t H 2 , S 2 t 3 K 3
Corollary 4. 
Let  p ^  be the value of  P RUU μ 1 , μ 2  when the value of  K 3  becomes “very high”, i.e., high enough for the probability  P S 2 t 2 K 3  to tend to zero; then, the difference  p ^ P RUU μ 1 , μ 2  provides the value of the following minor variant:
p ^ P RUU μ 1 , μ 2 = P sup 0 t t 1 S 1 t H 1 , S 1 t 1 K 1 , S 2 t 1 K 2 , sup t 1 t t 2 S 2 t H 2 , S 2 t 2 > K 3
End of Proposition 1.
Equipped with Proposition 1, one can value in closed form a two-colour step barrier option with two successive upward or two successive downward steps. Applying the theory of non-arbitrage pricing in a complete market (Harrison and Kreps 1979; Harrison and Pliska 1981), the value of a two-colour up-and-up knock-out put, denoted by V RUU , is given by:
V RUU = e r t 2 E Q K 3 1 A S 2 t 2 1 A
where r is the riskless interest rate assumed to be constant, 1 . is the indicator function and A is the set constructed by the intersection of elements of the σ algebra generated by the pair of processes S 1 t , S 2 t that characterises the probability P RUU μ 1 , μ 2 as given by the arguments of the probability operator in (5).
A simple application of the Cameron–Martin–Girsanov theorem yields:
V RUU = e r t 2 K 3 P RUU μ 1 Q , μ 2 Q S 2 0 P RUU μ 1 P 2 , μ 2 P 2
where
μ i Q = r σ i 2 2 ,   μ 1 P 2 = r σ 1 2 2 + σ 1 σ 2 ρ 1.2 ,   μ 2 P 2 = r + σ 2 2 2 ,
Q is the measure under which B i t + μ i r σ i t , t 0 is a standard Brownian motion (the classical so-called risk-neutral measure), while P 2 is the measure under which B 1 t σ 2 ρ 1.2 t , t 0 and B 2 t σ 2 1 ρ 1.2 2 t , t 0 are two independent standard Brownian motions.
To factor in a continuous dividend rate δ i associated with each asset S i , simply replace r by r δ i .
All the other two-colour rainbow barrier options subsequently mentioned in Section 2.1 and Section 2.2, whether they be knock-in or feature a mixture of a downward and an upward barrier, are identically valued, by taking the relevant P . probability along with the pairs μ 1 Q , μ 2 Q and μ 1 P 2 , μ 2 P 2 .
The numerical implementation of Proposition 1 is easy. Using Genz’s (2004) algorithm to evaluate the trivariate standard normal cumulative distribution function, the accuracy and efficiency required for all practical purposes can be achieved in computational times in the order of 0.1 s. Table 1 provides the prices of a few two-colour up-and-up knock-out put options, for various levels of the volatility and correlation parameters of the underlying assets S 1 and S 2 , and different values of the knock-out barriers. All the initial values of the underlying assets S i 0 and the strike prices K i are set at 100. Expiry is 1 year. The two time intervals t 0 , t 1 and t 1 , t 2 have equal length, i.e., t 1 = 6 months, but unequal time lengths can be handled just as well by the formulae. The riskless interest rate is assumed to be 2.5%.
In each cell, four prices are reported: the first one is the exact analytical value as obtained by implementing Proposition 1, while the prices in brackets are three successive approximations obtained by performing increasingly large Monte Carlo simulations. More specifically, these approximations rely on the conditional Monte Carlo method, which is well known for its accuracy and efficiency (Glasserman 2003). The number of simulations performed is 500,000 for the first approximation, 2,000,000 for the second approximation, and 10,000,000 for the third approximation. The pseudo-random numbers are drawn from the reliable Mersenne Twister generator (Matsumoto and Nishimura 1998).
In purely numerical terms, it can be clearly observed that the conditional Monte Carlo approximations gradually converge to the analytical values as more and more simulations are performed. A minimum of 10,000,000 simulations are necessary to guarantee a modest 10 3 convergence. This requires a computational time of approximately 35 s on a computer equipped with a Core i7 CPU. Much more accurate values can be obtained by means of Proposition 1 in only two-tenths of a second. This gap in accuracy and efficiency makes a particularly valuable difference when pricing large portfolios of options.
From a financial point of view, the most striking phenomenon observed in Table 1 is that the option price regularly and significantly increases with the value of the correlation coefficient between assets S 1 and S 2 , whatever the volatilities and the levels of the barriers. Roughly speaking, the price of an at-the-money two-colour up-and-up knock-out put option when ρ 1.2 = 0.6 is three times greater than when ρ 1.2 = 0.6 . This property can be exploited by traders who take positions on correlation, as the prices of these options will substantially increase if implicit correlation turns out to be underestimated by the markets. This property can also be harnessed by traders to construct hedges on sold derivatives that are sensitive to pairwise correlation. From an investor’s perspective, the observed phenomenon allows to define effective strategies to reduce the cost of hedging by tapping into negative correlation. Such a significant functional relation w.r.t. correlation is a major attraction of rainbow step barrier options relative to non-rainbow step barrier options, as the latter can only handle volatility effects.
Another noticeable fact in Table 1 is that lowering the up-and-out barriers seems much more effective in reducing the option’s price than lowering the volatilities of assets S 1 and S 2 , regardless of the sign and the magnitude of correlation. Indeed, looking at row 1 in Table 1, one can see that the options are relatively cheap, although the volatilities of both assets S 1 and S 2 are low, because the knock-out barriers are located quite near the spot prices of the underlying assets; and looking at row 2 in Table 1, one can see that the options are relatively expensive, although the volatilities of both assets S 1 and S 2 are high because the knock-out barriers are more distant. This shows that the barrier effect, which drives prices down as up-and-out barriers become lower and conversely drives prices up as the up-and-out barrier becomes higher, and prevails over the volatility effect, which exerts its influence in the opposite direction, i.e., a lower volatility pushes prices up by decreasing the probability of knocking out before expiry and a higher volatility pushes prices down by increasing the latter probability. This phenomenon can be explained by the ambivalent nature of volatility: on the one hand, less volatility means less risk of being deactivated before expiry, but on the other hand, it also means fewer chances of ending in-the-money at expiry; whichever of this positive and this negative effect weighs more on the option price depends on the relative values of barrier, strike, volatility and expiry parameters in a complex manner.

2.2. Valuation of Two-Colour Step Barrier Options Involving One Upward Step and One Downward Step

Section 2.2 deals with the case when the steps of the barrier are not on the same side in each time interval, i.e., either first downward, then upward, or first upward, then downward.
The main result of Section 2.2 is given by the following Proposition 2.
Proposition 2. 
Let  P RUD μ 1 , μ 2  denote the joint cumulative distribution function defined by:
P RUD μ 1 , μ 2 P sup 0 t t 1 S 1 t H 1 , S 1 t 1 K 1 , S 2 t 1 K 2 , inf t 1 t t 2 S 2 t H 2 , S 2 t 2 K 3
where the acronym “ RUD ” stands for “Rainbow Up and Down”.
Then, the exact value of P RUD μ 1 , μ 2 is given by:
P RUD μ 1 , μ 2 = N 3 min k 1 , h 1 μ 1 t 1 σ 1 t 1 , max k 2 , h 2 + μ 2 t 1 σ 2 t 1 , max k 3 , h 2 + μ 2 t 2 σ 2 t 2 ; θ 1.2 , θ 1.3 , θ 2.3
exp 2 μ 1 h 1 σ 1 2 × N 3 min k 1 , h 1 2 h 1 μ 1 t 1 σ 1 t 1 , max k 2 , h 2 + μ 2 t 1 σ 2 t 1 + 2 θ 1.2 h 1 σ 1 t 1 , max k 3 , h 2 + μ 2 t 2 σ 2 t 2 + 2 θ 1.2 h 1 σ 1 t 2 ; θ 1.2 , θ 1.3 , θ 2.3
exp 2 μ 2 h 2 σ 2 2 × N 3 min k 1 , h 1 μ 1 t 1 σ 1 t 1 + 2 θ 1.2 μ 2 t 1 σ 2 , max k 2 , h 2 μ 2 t 1 σ 2 t 1 , max k 3 , h 2 + 2 h 2 + μ 2 t 2 σ 2 t 2 ; θ 1.2 , θ 1.3 , θ 2.3
+ exp 2 μ 1 σ 1 2 4 μ 2 θ 1.2 σ 1 σ 2 h 1 + 2 μ 2 h 2 σ 2 2 × N 3 min k 1 , h 1 2 h 1 μ 1 t 1 σ 1 t 1 + 2 θ 1.2 μ 2 t 1 σ 2 , max k 2 , h 2 μ 2 t 1 σ 2 t 1 + 2 θ 1.2 h 1 σ 1 t 1 , max k 3 , h 2 + 2 h 2 + μ 2 t 2 σ 2 t 2 2 θ 1.2 h 1 σ 1 t 2 ; θ 1.2 , θ 1.3 , θ 2.3
where all the notations are identical, as in Proposition 1.
Corollary 1. 
It suffices to multiply by  1  all the first three arguments of each  N 3 . , . , . ; . , . , .  function and substitute each  min  operator by a  max  operator as well as each  max  operator by a  min  operator in Proposition 2 to obtain an exact formula for  P RDU μ 1 , μ 2  defined as:
P RDU μ 1 , μ 2 P inf 0 t t 1 S 1 t H 1 , S 1 t 1 K 1 , S 2 t 1 K 2 , sup t 1 t t 2 S 2 t H 2 , S 2 t 2 K 3
Corollary 2. 
The term numbered (23) in Proposition 2 provides the value of the corrresponding up-and-in, then down-and-in probability, denoted as  P RUD I μ 1 , μ 2  and defined by:
P RUD I μ 1 , μ 2 P sup 0 t t 1 S 1 t H 1 , S 1 t 1 K 1 , S 2 t 1 K 2 , inf t 1 t t 2 S 2 t H 2 , S 2 t 2 K 3
End of Proposition 2.
Equipped with Proposition 2, one can value in closed form a two-colour step barrier option with one upward step and one downward step, by taking the relevant P . or P . I probability along with the pairs μ 1 Q , μ 2 Q and μ 1 P 2 , μ 2 P 2 defined in (18), as explained in Section 2.1. Table 2 reports the prices of a few down-and-up two-colour knock-out put options by implementing Proposition 2 to obtain exact analytical values and by computing three successive conditional Monte Carlo approximations in the same way, as in Table 1.
In Table 2, the most salient feature is still the functional dependency of the option’s price on the correlation between assets S 1 and S 2 , but, this time, the direction is opposite to that in Table 1, i.e., the two-colour down-and-up knock-out put prices steadily decrease as ρ 1.2 goes from −60% to 60%. In a trader’s perspective, one could sum up the argument by saying that two-colour rainbow barrier options are a bet on a positive correlation when both barriers are on the same side (upward or downward), while they are a bet on a negative correlation when the barriers stand on opposite sides (up-and-down or down-and-up).
The barrier effect also prevails over the volatility effect in Table 2. Overall, two-colour down-and-up knock-out puts display maximum values that are a little higher, and minimum values that are a little lower than two-colour up-and-up knock-out puts, although up-and-out barriers and down-and-out barriers are designed with the exact same distance to the spot prices of S 1 and S 2 .

2.3. Valuation of Reverse Two-Colour Step Barrier Options

A two-colour rainbow barrier option is said to be reverse when the moneyness of the option is defined w.r.t. the first and former “colour” (i.e., asset S 1 ) instead of the second and last one (asset S 2 ): the option, so to speak, reverts back to asset one at expiry, hence the denomination. From a computational standpoint, this is not a trivial difference since it adds an additional dimension to the integral formulation of the problem. Let us define as P RRUU μ 1 , μ 2 the following cumulative joint distribution at the core of reverse rainbow option valuation:
P RRUU Rev μ 1 , μ 2 = P sup 0 t t 1 S 1 t H 1 , S 1 t 1 K 1 , S 2 t 1 K 2 , sup t 1 t t 2 S 2 t H 2 , S 1 t 2 K 3
where the acronym “ RRUU ” stands for “Reverse Rainbow Up and Up”.
Then, Proposition 3 provides the exact value of P RRUU μ 1 , μ 2 in the form of a triple integral.
Proposition 3. 
P RRUU μ 1 , μ 2 = 1 2 π 3 / 2 σ 2 1 σ 3 1.2 σ 1 2 σ 2 t 1 t 2 min k 1 , h 1 h 2 min k 2 , h 2 φ 2 x 1 φ 3 x 2 , x 3
e 1 2 x 1 μ 1 t 1 σ 1 t 1 2 1 2 σ 2 1 2 x 2 μ 2 t 1 σ 2 t 1 θ 1.2 x 1 μ 1 t 1 σ 1 t 1 2 1 2 σ 3 1.2 2 x 3 μ 2 t 2 σ 2 t 2 θ 1.3 x 1 μ 1 t 1 σ 1 t 1 θ 2.3 1 σ 2 1 x 2 μ 2 t 1 σ 2 t 1 θ 1.2 x 1 μ 1 t 1 σ 1 t 1 2
N 1 σ 4 1.2.3 k 3 μ 1 t 2 σ 1 t 2 θ 1.4 x 1 μ 1 t 1 σ 1 t 1 θ 2.4 1 σ 2 1 x 2 μ 2 t 1 σ 2 t 1 θ 1.2 x 1 μ 1 t 1 σ 1 t 1 θ 3.4 1.2 σ 3 1.2 x 3 μ 2 t 2 σ 2 t 2 θ 1.3 x 1 μ 1 t 1 σ 1 t 1 θ 2.3 1 σ 2 1 x 2 μ 2 t 1 σ 2 t 1 θ 1.2 x 1 μ 1 t 1 σ 1 t 1 d x 3 d x 2 d x 1
where: 
θ 2.4 1 = θ 2.4 θ 1.2 θ 1.4 1 θ 1.2 2 , θ 3.4 1.2 = θ 3.4 θ 1.3 θ 1.4 θ 2.3 1 θ 2.4 1 σ 3 1.2 , σ 4 1.2.3 = 1 θ 1.4 2 θ 2.4 1 2 θ 3.4 1.2 2
-
N .  is the univariate standard normal cumulative distribution function;
-
the functions  φ 2  and  φ 3  are defined by (80) and (81) in Section 3.
All the other notations in Proposition 3 have been previously defined.
Remark 1. 
Other types of reverse two-colour knock-out or knock-in barrier probability distributions are handled similarly by modifying the upper bounds of the integral and, possibly, the  φ i  functions, according to the considered combination of events.
Remark 2. 
θ 2.4 1  is the partial correlation between  X 2 t 1  and  X 1 t 2  conditional on  X 1 t 1 , while  θ 3.4 1.2  is the partial correlation between  X 2 t 2  and  X 1 t 2  conditional on  X 1 t 1  and  X 2 t 1 , and  σ 4 1.2.3  is the conditional standard deviation of  X 1 t 2  given  X 1 t 1 ,  X 2 t 1  and  X 2 t 2 .
End of Proposition 3.
The application of Proposition 3 to value a reverse two-colour step barrier option is now discussed. The no-arbitrage price of a reverse two-colour rainbow up-and-up knock-out put, denoted by V RRUU , is given by:
V RRUU R = e r t 2 E Q K 3 1 A S 1 t 2 1 A = e r t 2 K 3 P RRUU R μ 1 Q , μ 2 Q S 2 0 P RRUU R μ 1 P 1 , μ 2 P 1
where
μ 1 P 1 = r + σ 1 2 2 ,   μ 2 P 1 = r σ 2 2 2 + σ 1 σ 2 ρ 1.2
-
A is the set constructed by the intersection of elements of the σ algebra generated by the pair of processes S 1 t , S 2 t that characterises the probability P RRUU μ 1 , μ 2 as given by the arguments of the probability operator in (26);
-
P 1 is the measure under which B 1 t σ 1 t is a standard Brownian motion.
However, it is less easy to evaluate Proposition 3 than to evaluate Proposition 1 and Proposition 2. The problem at hand has two “nice” features from the standpoint of numerical integration: first, the dimension, equal to 3, is moderate; second, the integrand is continuous. The snag is the large number of parameters in each evaluation of the integrand in a quadrature process, especially the various conditional standard deviations at the denominators of the fractions, that may hinder fast convergence when they take on absolute values that become smaller and smaller. That is why it is recommended to use a subregion adaptive algorithm of numerical integration, as explained by Berntsen et al. (1991), that adapts the number of integrand evaluations in each subregion according to the rate of change of the integrand. Although more time-consuming than a fixed degree rule, it is more accurate to control the approximation error, as the subdivision of the integration domain stops only when the sum of the local error deterministic estimates becomes smaller than some prespecified requested accuracy. Adaptive integration can be enhanced by a Kronrod rule to reduce the number of required iterations (see, e.g., Davis and Rabinowitz 2007). These techniques are widely used in numerical integration, and it is easy to find available code or built-in functions in the usual scientific computing software.

2.4. Valuation of Two-Colour Outside Step Barrier Options

In this section, a third correlated asset S 3 is introduced, w.r.t. which the option’s moneyness is measured at expiry, while the assets S 1 and S 2 serve exclusively as the underlyings w.r.t. which barrier crossing is monitored. This is an important extension, as outside barrier options allow to manage volatility more consistently than standard (non-outside) barrier options, as explained, e.g., by Das (2006).
Let us consider a third asset S 3 with the following differential:
d S 3 t = v 3 S 3 t d t + σ 3 S 3 t d B 3 t
The instantaneous pairwise correlations between the Brownian motions B i s are denoted as ρ i . j .
The objective is to compute the probabilities p m μ 1 , μ 2 , μ 3 , m 1 , 2 , 3 , 4 defined by:
p 1 μ 1 , μ 2 , μ 3 = P sup 0 t t 1 S 1 t H 1 , S 1 t 1 K 1 , S 2 t 1 H 2 , sup t 1 t t 2 S 2 t H 2 , S 2 t 2 K 2 , S 3 t 2 K 3
p 2 μ 1 , μ 2 , μ 3 = P inf 0 t t 1 S 1 t H 1 , S 1 t 1 K 1 , S 2 t 1 H 2 , inf t 1 t t 2 S 2 t H 2 , S 2 t 2 K 2 , S 3 t 2 K 3
p 3 μ 1 , μ 2 , μ 3 = P sup 0 t t 1 S 1 t H 1 , S 1 t 1 K 1 , S 2 t 1 H 2 , inf t 1 t t 2 S 2 t H 2 , S 2 t 2 K 2 , S 3 t 2 K 3
p 4 μ 1 , μ 2 , μ 3 = P inf 0 t t 1 S 1 t H 1 , S 1 t 1 K 1 , S 2 t 1 H 2 , sup t 1 t t 2 S 2 t H 2 , S 2 t 2 K 2 , S 3 t 2 K 3
Let x = c 1 , c 2 , c 3 , c 4 , c 5 be a vector of five coordinates where each c i 1 , 1 , i 1 , , 5 .
Let the function Ψ 4 b 1 , b 2 , b 3 , b 4 ; x , b 1 , b 2 , b 3 , b 4 , be defined by:
Ψ 4 b 1 , b 2 , b 3 , b 4 ; x
= x 1 = b 1 x 2 = b 2 c 1 x 1 1 c 1 2 x 3 = b 3 c 4 x 2 1 c 1 2 c 4 c 1 x 1 1 c 2 2 1 2 π 3 / 2 exp x 1 2 2 x 2 2 2 x 3 2 2
N b 4 c 5 c 2 c 3 1 c 2 2 x 3 1 c 4 2 + x 2 c 4 1 c 1 2 x 1 c 3 + c 5 c 2 c 3 1 c 2 2 c 1 c 4 c 2 1 c 3 2 c 5 c 2 c 3 1 c 2 2 2 d x 3 d x 2 d x 1
The following Proposition 4 combines all the probabilities defined in (32)–(35) into a single formula.
Proposition 4. 
The exact values of the probabilities  p m μ 1 , μ 2 , μ 3 , m 1 , 2 , 3 , 4 , written in shorter notation as  p m , are given by:
p m = Ψ 4 δ 1 G 1 k 1 , h 1 μ 1 t 1 σ 1 t 1 , δ 2 h 2 μ 2 t 1 σ 2 t 1 , δ 2 G 2 k 2 , h 2 μ 2 t 2 σ 2 t 2 , δ 2 k 3 μ 3 t 2 σ 3 t 2 ; x 1
exp 2 μ 1 h 1 σ 1 2 × Ψ 4 δ 1 G 1 k 1 , h 1 2 h 1 μ 1 t 1 σ 1 t 1 , δ 2 h 2 μ 2 t 1 σ 2 t 1 θ 1.2 2 h 1 σ 1 t 1 , δ 2 G 2 k 2 , h 2 μ 2 t 2 σ 2 t 2 θ 1.2 2 h 1 σ 1 t 2 , δ 2 k 3 μ 3 t 2 σ 3 t 2 θ 1.4 2 h 1 σ 1 t 1 θ 3.4 1 θ 1.2 2 h 1 σ 1 t 2 θ 1.3 2 h 1 σ 1 t 1 ; x 1
exp 2 μ 2 h 2 σ 2 2 × Ψ 4 δ 1 G 1 k 1 , h 1 μ 1 t 1 σ 1 t 1 + θ 1.2 2 μ 2 t 1 σ 2 t 1 , δ 2 h 2 + μ 2 t 1 σ 2 t 1 , δ 2 G 2 k 2 , h 2 2 h 2 μ 2 t 2 σ 2 t 2 , δ 2 k 3 μ 3 t 2 σ 3 t 2 θ 3.4 1 2 h 2 σ 2 t 2 + θ 1.2 θ 1.3 2 μ 2 t 1 σ 2 t 1 + θ 1.2 θ 1.4 2 μ 2 t 1 σ 2 t 1 ; x 2
+ exp 2 μ 1 σ 1 2 4 μ 2 ρ 1.2 σ 1 σ 2 h 1 + 2 μ 2 h 2 σ 2 2 × Ψ 4 δ 1 G 1 k 1 , h 1 2 h 1 μ 1 t 1 σ 1 t 1 + θ 1.2 2 μ 2 t 1 σ 2 t 1 , δ 2 h 2 + μ 2 t 1 σ 2 t 1 θ 1.2 2 h 1 σ 1 t 1 , δ 2 G 2 k 2 , h 2 2 h 2 μ 2 t 2 σ 2 t 2 + θ 1.2 2 h 1 σ 1 t 2 , δ 2 k 3 μ 3 t 2 σ 3 t 2 + θ 1.4 θ 1.2 2 μ 2 t 1 σ 2 2 h 1 σ 1 t 1 θ 3.4 1 2 h 2 σ 2 t 2 θ 1.2 2 h 1 σ 1 t 2 + θ 1.3 θ 1.2 2 μ 2 t 1 σ 2 2 h 1 σ 1 t 1 ; x 2
where  k 1 , k 2 , h 1 , h 2  are as in Proposition 2,  k 3 = ln K 3 S 3 0 , and we have:
θ 1.2 = ρ 1.2 , θ 1.3 = t 1 t 2 ρ 1.2 , θ 1.4 = t 1 t 2 ρ 1.3 , θ 2.3 = t 1 t 2 , θ 3.4 = ρ 2.3 ,   θ 3.4 1 = θ 3.4 θ 1.3 θ 1.4 1 θ 1.3 2
δ 1 = 1   if   p m = p 1   or   p m = p 3 1   if   p m = p 2   or   p m = p 4 ,   δ 2 = 1   if   p m = p 1   or   p m = p 4 1   if   p m = p 2   or   p m = p 3
G 1 . , . = min . , .   if   p m = p 1   or   p m = p 3 max . , .   if   p m = p 2   or   p m = p 4 ,   G 2 . , . = max . , .   if   p m = p 1   or   p m = p 3 min . , .   if   p m = p 2   or   p m = p 4
x 1 = θ 1.2 , θ 1.3 , θ 1.4 , θ 2.3 , θ 3.4   if   p m = p 1   or   p m = p 2 θ 1.2 , θ 1.3 , θ 1.4 , θ 2.3 , θ 3.4   if   p m = p 3   or   p m = p 4
x 2 = θ 1.2 , θ 1.3 , θ 1.4 , θ 2.3 , θ 3.4   if   p m = p 1   or   p m = p 2 θ 1.2 , θ 1.3 , θ 1.4 , θ 2.3 , θ 3.4   if   p m = p 3   or   p m = p 4
Corollary 1. 
The corresponding knock-in probabilities can be inferred in the same way as in Proposition 1 and Proposition 2. Let the probabilities  p m I μ 1 , μ 2 , μ 3 , m 1 , 2 , 3 , 4  be defined by:
p 1 I μ 1 , μ 2 , μ 3 = P sup 0 t t 1 S 1 t H 1 , S 1 t 1 K 1 , S 2 t 1 H 2 , sup t 1 t t 2 S 2 t H 2 , S 2 t 2 K 2 , S 3 t 2 K 3
p 2 I μ 1 , μ 2 , μ 3 = P inf 0 t t 1 S 1 t H 1 , S 1 t 1 K 1 , S 2 t 1 H 2 , inf t 1 t t 2 S 2 t H 2 , S 2 t 2 K 2 , S 3 t 2 K 3
p 3 I μ 1 , μ 2 , μ 3 = P sup 0 t t 1 S 1 t H 1 , S 1 t 1 K 1 , S 2 t 1 H 2 , inf t 1 t t 2 S 2 t H 2 , S 2 t 2 K 2 , S 3 t 2 K 3
p 4 I μ 1 , μ 2 , μ 3 = P inf 0 t t 1 S 1 t H 1 , S 1 t 1 K 1 , S 2 t 1 H 2 , sup t 1 t t 2 S 2 t H 2 , S 2 t 2 K 2 , S 3 t 2 K 3
Then,  p m I μ 1 , μ 2 , μ 3  is given by (40).
Corollary 2. 
It suffices to substitute each argument  θ 3.4  in each  Ψ 4 . , . , . , . ; . , . , . , . , .  function of Proposition 4 by  ρ 2.3 t 2 t 3 ,  t 3 t 2 , to obtain an exact formula for the early-ending variant of  p m μ 1 , μ 2 , μ 3 .
End of Proposition 4.
Equipped with Proposition 4, one can value in closed form a two-colour outside step barrier option. More precisely, the value of a two-colour outside up-and-out put, denoted by V ORUU , is given by:
V ORUU = e r t 2 E Q K 3 1 A S 3 t 2 1 A
where A is the set constructed by the intersection of elements of the σ algebra generated by the pair of processes S 1 t , S 2 t that characterises the probability p 1 μ 1 , μ 2 as given by the arguments of the probability operator in (32), and the acronym “ ORUU ” stands for “Outside Rainbow Up and Up”.
Using the following orthogonal decomposition of Brownian motion B 3 t :
B 3 t = ρ 1.3 W 1 t + ρ 2.3 1 W 2 t + σ 3 1.2 W 3 t
where:
ρ 2.3 1 = ρ 2.3 ρ 1.2 ρ 1.3 1 ρ 1.2 2 ,   σ 3 1.2 = 1 ρ 1.3 2 ρ 2.3 1 2
and W 1 t , W 2 t , W 3 t is a basis of three independent Brownian motions (Guillaume 2018), the multidimensional Cameron-Martin-Girsanov theorem yields:
V ORUU = e r t 2 K 3 p 1 μ 1 Q , μ 2 Q , μ 3 Q S 3 0 p 1 μ 1 P 3 , μ 2 P 3 , μ 3 P 3
where:
μ 1 P 3 = r σ 1 2 2 + σ 1 σ 3 ρ 1.3 ,   μ 2 P 3 = r σ 2 2 2 + σ 2 σ 3 ρ 2.3 ,   μ 3 P 3 = r + σ 3 2 2
The measure P 3 is the measure under which B 1 t σ 3 ρ 1.3 t , B 2 t σ 3 ρ 2.3 1 t and B 3 t σ 3 σ 3 1.2 t are three independent standard Brownian motions.
A simple and robust numerical evaluation of the function Ψ 4 consists in selecting an appropriate cutoff value for the negative infinity lower bounds and then applying a fixed-degree quadrature rule. Given the smoothness of the integrand, even a low-degree rule will perform well. Table 3 provides the prices of a few two-colour outside up-and-down knock-out call options for various levels of the volatility and correlation parameters of the underlying assets S 1 , S 2 , and S 3 , and different values of the knock-out barriers. The parameters S i 0 , K i , t 1 , t 2 , and r are identical as those as in Table 1 and Table 2. In each cell, the first reported value is the exact analytical price, as obtained by implementing Proposition 4 by means of a classical 16-point Gauss–Legendre quadrature, while the numbers in the brackets are three successive Monte Carlo approximations, as explained in Section 2.1.
From a purely numerical standpoint, the pattern of convergence of conditional Monte Carlo approximations to the analytical values is as clear in Table 3 as in Table 1 and Table 2. This illustrates the robustness of our numerical integration scheme for the Ψ 4 function. The efficiency gap between Monte Carlo pricing and analytical pricing is even more pronounced than for non-outside rainbow step barrier options due to the presence of an additional stochastic process to simulate: the average computational time required by simulation is 42 s, whereas the evaluation of the analytical formula based on Proposition 4 only takes a few tenths of a second.
From a financial point of view, the prices in Table 3 display a very different pattern from those in Table 1 and Table 2. With regard to correlation, the highest option values attained are when the correlation between S 1 and S 2 is negative and the correlation between S 3 and both S 1 and S 2 is positive. The lowest option values are when the correlation between S 1 and S 2 is positive and the correlation between S 3 and both S 1 and S 2 is negative. On average across all volatilities and barrier levels in Table 3, options are approximately three times more expensive under the former correlation structure than under the latter one. In terms of volatility, the highest option values attained are when the volatility of asset S 3 is high. This remains true under very different combinations of values for all the other parameters (volatilities of S 1 and S 2 , barrier levels and correlation structure). Such an observation highlights the prominent role of the volatility of the asset chosen to determine the moneyness of the option at expiry. In particular, the value of a rainbow outside step barrier option is a monotonically increasing function of σ 3 , whereas the value of a rainbow step barrier option is not a monotonically increasing function of σ 2 , just like the value of a reverse rainbow step barrier option is not a monotonically increasing function of σ 1 . This is because a rainbow outside step barrier option allows to make a clear distinction between the functions of each underlying asset: two of them, S 1 and S 2 , are only concerned with barrier crossing during the option life, and the third one, S 3 , is only concerned with moneyness testing at the option expiry. That distinction is impossible to make when it comes to non-outside rainbow step barrier options, so that the impact of volatility becomes ambiguous and difficult to handle. It should be emphasised that, for the vast majority of parameters, the sensitivities of the rainbow outside step knock-out barrier options to σ 1 and σ 2 is negative, reflecting an increased risk of being deactivated before expiry. Only for quite specific correlation structures between the underlying assets and quite specific combinations of barrier values can these sensitivities be positive. A major advantage of closed form formulae such as those derived in this article is precisely to allow measurement of such sensitivities with high precision by mere differentiation of the formulae w.r.t. the relevant parameters.
One more noticeable difference in the reported numerical results between outside and non-outside rainbow step barrier options is that the volatility effect prevails over the barrier effect in Table 3, in contrast to Table 1 and Table 2. Indeed, in row 2 of Table 3, tight barrier levels do not preclude relatively high option prices thanks to the volatility of asset S 3 set at 60%. Likewise, in row 3 of Table 3, wider barrier levels do not preclude relatively low option prices due to the volatility of asset S 3 set at only 20%.

2.5. Valuation of Two-Sided, Two-Colour Step Barrier Options

In this section, a two-sided barrier is introduced in each time interval, i.e., the valuation of rainbow step double barrier options is handled.
Let H 1 and H 2 denote an upward and a downward barrier, respectively, on the time interval t 0 = 0 , t 1 . Similarly, H 3 and H 4 represent an upward and a downward barrier, respectively, on the time interval t 1 , t 2 . As in the previous sections, barrier crossing is monitored w.r.t. process S 1 following Equation (1) on t 0 = 0 , t 1 and w.r.t. process S 2 following Equation (2) on t 1 , t 2 . Our objective now is to find the value of the joint cumulative distribution function P RDKO μ 1 , μ 2 defined by:
P RDKO μ 1 , μ 2
= P sup 0 t t 1 S 1 t H 1 , inf 0 t t 1 S 1 t H 2 , S 1 t 1 min K 1 , H 1 , sup t 1 t t 2 S 2 t H 3 , inf t 1 t t 2 S 2 t H 4 , S 2 t 2 min H 3 , K 2
where the acronym “ RDKO ” stands for “Rainbow Double Knock Out”.
The main result of Section 2.5 is given by the following Proposition 5.
Proposition 5. 
The exact value of  P RDKO μ 1 , μ 2  is given by:
P RDKO μ 1 , μ 2 = n 1 = n 2 = exp 2 μ 1 σ 1 2 n 1 a 1 + 2 μ 2 σ 2 2 n 2 a 2 N 3 A 1 min k 1 , h 1 , A 2 h 3 , A 3 min h 3 , k 2 ; x 1 N 3 A 1 h 2 , A 2 h 3 , A 3 min h 3 , k 2 ; x 1 N 3 A 1 min k 1 , h 1 , A 2 h 4 , A 3 min h 3 , k 2 ; x 1 + N 3 A 1 h 2 , A 2 h 4 , A 3 min h 3 , k 2 ; x 1 N 3 A 1 min k 1 , h 1 , A 2 h 3 , A 3 h 4 ; x 1 + N 3 A 1 h 2 , A 2 h 3 , A 3 h 4 ; x 1 + N 3 A 1 min k 1 , h 1 , A 2 h 4 , A 3 h 4 ; x 1 N 3 A 1 h 2 , A 2 h 4 , A 3 h 4 ; x 1
n 1 = n 2 = exp n 1 a 1 2 μ 1 σ 1 2 4 ρ 1.2 μ 2 σ 1 σ 2 + 2 μ 2 σ 2 2 h 4 n 2 a 2 N 3 A 4 min k 1 , h 1 , A 5 h 3 , A 6 min h 3 , k 2 ; x 2 N 3 A 4 h 2 , A 5 h 3 , A 6 min h 3 , k 2 ; x 2 N 3 A 4 min k 1 , h 1 , A 5 h 4 , A 6 min h 3 , k 2 ; x 2 + N 3 A 4 h 2 , A 5 h 4 , A 6 min h 3 , k 2 ; x 2 N 3 A 4 min k 1 , h 1 , A 5 h 3 , A 6 h 4 ; x 2 + N 3 A 4 h 2 , A 5 h 3 , A 6 h 4 ; x 2 + N 3 A 4 min k 1 , h 1 , A 5 h 4 , A 6 h 4 ; x 2 N 3 A 4 h 2 , A 5 h 4 , A 6 h 4 ; x 2
n 1 = n 2 = exp 2 μ 1 σ 1 2 h 2 n 1 a 1 + 2 μ 2 σ 2 2 n 2 h 3 a 2 N 3 A 7 min k 1 , h 1 , A 8 h 3 , A 9 min h 3 , k 2 ; x 1 N 3 A 7 h 2 , A 8 h 3 , A 9 min h 3 , k 2 ; x 1 N 3 A 7 min k 1 , h 1 , A 8 h 4 , A 9 min h 3 , k 2 ; x 1 + N 3 A 7 h 2 , A 8 h 4 , A 9 min h 3 , k 2 ; x 1 N 3 A 7 min k 1 , h 1 , A 8 h 3 , A 9 h 4 ; x 1 + N 3 A 7 h 2 , A 8 h 3 , A 9 h 4 ; x 1 + N 3 A 7 min k 1 , h 1 , A 8 h 4 , A 9 d 2 ; x 1 N 3 A 7 h 2 , A 8 h 4 , A 9 d 2 ; x 1
+ n 1 = n 2 = exp h 2 n 1 a 1 2 μ 1 σ 1 2 4 ρ 1.2 μ 2 σ 1 σ 2 + 2 μ 2 σ 2 2 h 4 n 2 a 2 N 3 A 10 min k 1 , h 1 , A 11 h 3 , A 12 min h 3 , k 2 ; x 2 N 3 A 10 h 2 , A 11 h 3 , A 12 min h 3 , k 2 ; x 2 N 3 A 10 min k 1 , h 1 , A 11 h 4 , A 12 min h 3 , k 2 ; x 2 + N 3 A 10 h 2 , A 11 h 4 , A 12 min h 3 , k 2 ; x 2 N 3 A 10 min k 1 , h 1 , A 11 h 3 , A 12 h 4 ; x 2 + N 3 A 10 h 2 , A 11 h 3 , A 3 h 4 ; x 2 + N 3 A 10 min k 1 , h 1 , A 11 h 4 , A 12 h 4 ; x 2 N 3 A 10 h 2 , A 11 h 4 , A 12 h 4 ; x 2
where:
-   h 1 = ln H 1 S 1 0 > 0 , h 2 = ln H 2 S 1 0 < 0 , h 3 = ln H 3 S 2 0 , h 4 = ln H 4 S 2 0
-   a 1 = h 1 h 2 , a 2 = h 3 h 4
-   A 1 x = x 2 n 1 a 1 μ 1 t 1 σ 1 t 1 , A 2 x = x μ 2 t 1 σ 2 t 1 2 ρ 1.2 n 1 a 1 σ 1 t 1
-   A 3 x = x 2 n 2 a 2 μ 2 t 2 σ 2 t 2 2 ρ 1.2 n 1 a 1 σ 1 t 2
-   A 4 x = A 1 x 2 ρ 1.2 μ 2 t 1 σ 2 , A 5 x = x + μ 2 t 1 σ 2 t 1 2 ρ 1.2 n 1 a 1 σ 1 t 1
-   A 6 x = x 2 h 4 + 2 n 2 a 2 μ 2 t 2 σ 2 t 2 + 2 ρ 1.2 n 1 a 1 σ 1 t 2
-   A 7 x = x 2 h 2 + 2 n 1 a 1 μ 1 t 1 σ 1 t 1 , A 8 x = x μ 2 t 1 σ 2 t 1 2 ρ 1.2 h 2 n 1 a 1 σ 1 t 1
-   A 9 x = x 2 n 2 a 2 μ 2 t 2 σ 2 t 2 2 ρ 1.2 h 2 n 1 a 1 σ 1 t 2
-   A 10 x = A 7 x 2 ρ 1.2 μ 2 t 1 σ 2 , A 11 x = x + μ 2 t 1 σ 2 t 1 2 ρ 1.2 h 2 n 1 a 1 σ 1 t 1
-   A 12 x = x 2 h 4 + 2 n 2 a 2 μ 2 t 2 σ 2 t 2 + 2 ρ 1.2 h 2 n 1 a 1 σ 1 t 2
-   x 1 = ρ 1.2 , ρ 1.2 t 1 t 2 , t 1 t 2 , x 2 = ρ 1.2 , ρ 1.2 t 1 t 2 , t 1 t 2
All other notations have been defined in the previous sections.
End of Proposition 5.
Pricing two-colour double knock-out barrier options can be achieved through the same changes of probability measures as those applicable to two-colour single knock-out barrier options, i.e., the value of a two-colour double knock-out put, denoted as V RDKO , is given by:
V RDKO = e r t 2 K 3 P RDKO μ 1 Q , μ 2 Q S 2 0 P RDKO μ 1 P 2 , μ 2 P 2
where the parameters μ i Q and μ i P 2 are given by Equation (18).
Table 4 provides the prices of a few two-colour knock-out double barrier puts for various levels of the volatility and correlation parameters of the underlying assets S 1 and S 2 , and different values of the knock-out barriers. Expiry is 6 months and t 1 is one quarter of a year. The parameters S i 0 , K i , and r are identical to those in Table 1, Table 2 and Table 3. In each cell, the first number is the exact analytical value as derived from (71), while the numbers in the brackets are three successive Monte Carlo approximations, as explained in Section 2.1.
Thanks to the rapidly decaying exponential functions in the integrands, a level of 10 7 convergence is attained by stopping at 8, the number of iterations controlled by the absolute values of n 1 and n 2 in the double sum operators, which results in a total computational time of less than 1 s. For higher values of the volatility parameters than those in Table 4, however, the uniform convergence of the double sums in (56)–(59) may require a greater number of iterations and thus take more time. The implementation of Proposition 5 using the Φ 3 function introduced in Section 3 is slightly faster than the one using the trivariate standard normal cumulative distribution function N 3 , although the difference is relatively negligible for most practical purposes. Both methods of implementation yield prices equal to at least 4 decimals.
From a financial standpoint, a striking contrast between the numerical results in Table 4 and those of the previous sections is the much weaker dependency of the option value on the correlation structure, as illustrated by the smaller differences between the four option prices associated with each combination of volatility and barrier parameters. It seems that, the more volatility, the more dependency on the correlation structure, as suggested by the comparison between row 1 and row 2. Another noticeable difference is that the functional relation with the correlation structure is not monotonic. This is particularly clear in row 2 where a relatively significant increase in value from ρ 1.2 = 0.6 to ρ 1.2 = 0.2 is followed by a relatively significant decrease in value from ρ 1.2 = 0.2 to ρ 1.2 = 0.2 , before a new increase in value from ρ 1.2 = 0.2 to ρ 1.2 = 0.6 . This more complex and unstable dependency on correlation structure suggests that two-colour knock-out double barrier options are a less suitable instrument for correlation trading than two-colour knock-out single barrier options. However, one should remain wary of drawing hasty conclusions from the comparison of the results in Table 4 and those in the previous sections, as the option parameters are not identical, especially regarding volatility and expiry.

3. Proofs of Formulae

The proofs of Propositions 2 and 3 are only outlined as they essentially follow the same steps as the proof of Proposition 1.
Proof of Proposition 1. 
Since the log function is strictly increasing, we have:
P RUU μ 1 , μ 2 = P sup 0 t t 1 X 1 t h 1 , X 1 t 1 k 1 , X 2 t 1 k 2 , sup t 1 t t 2 X 2 t h 2 , X 2 t 2 k 3
Next, it can be noticed that, despite the non-zero correlation between X 1 t and X 2 t , the law of sup 0 t t 1 X 1 t conditional on X 1 t 1 and X 2 t 1 is equal to the law of sup 0 t t 1 X 1 t 1 conditional on X 1 t 1 .
Indeed, denoting the density function operator as f . and making use of the Markov property of X 2 t we have:
f sup 0 t t 1 X 1 t X 1 t 1 , X 2 t 1 = f X 2 t 1 sup 0 t t 1 X 1 t , X 1 t 1 f sup 0 t t 1 X 1 t , X 1 t 1 f X 1 t 1 , X 2 t 1
= f X 2 t 1 X 1 t 1 f sup 0 t t 1 X 1 t , X 1 t 1 f X 1 t 1 , X 2 t 1 = f X 1 t 1 , X 2 t 1 f X 1 t 1 f sup 0 t t 1 X 1 t , X 1 t 1 f X 1 t 1 , X 2 t 1
= f sup 0 t t 1 X 1 t , X 1 t 1 f X 1 t 1 = f sup 0 t t 1 X 1 t X 1 t 1
A translation from the time interval t 0 = 0 , t 1 to the time interval t 1 , t 2 , through the substitution of X 1 0 with X 2 t 1 , of X 1 t 1 with X 2 t 2 and of X 2 t 1 with X 1 t 2 , shows similarly that the law of sup t 1 t t 2 X 2 t conditional on X 2 t 1 , X 2 t 2 and X 1 t 2 is equal to the law of sup t 1 t t 2 X 2 t conditional on X 2 t 1 and X 2 t 2 .
Thus, by conditioning w.r.t. the absolutely continuous random variables X 1 t 1 , X 2 t 1 and X 2 t 2 , we can express the problem as the following integral:
P RUU μ 1 , μ 2 = min k 1 , h 1 min k 2 , h 2 min k 3 , h 2 φ 1 x 1 , x 2 , x 3 φ 2 x 1 φ 3 x 2 , x 3 d x 3 d x 2 d x 1
where
φ 1 x 1 , x 2 , x 3 = P X 1 t 1 d x 1 , X 2 t 1 d x 2 , X 2 t 2 d x 3 d x 3 d x 2 d x 1
φ 2 x 1 = P sup 0 t t 1 X 1 t h 1 X 1 t 1 d x 1 d x 1
φ 3 x 2 , x 3 = P sup t 1 t t 2 X 2 t h 2 X 2 t 1 d x 2 , X 2 t 2 d x 3 d x 2 d x 3
The functions φ 2 and φ 3 in (78) and (79) can be expanded by applying known formulae that can be found in Wang and Pötzelberger (1997):
φ 2 x 1 = 1 exp 2 h 1 x 1 h 1 σ 1 2 t 1
φ 3 x 2 , x 3 = 1 exp 2 h 2 x 2 x 3 h 2 σ 2 2 t 2 t 1
The function φ 1 derives from the trivariate normality of the triple X 1 t 1 , X 2 t 1 , X 2 t 2 . It is elementary to obtain the marginal distributions:
X 1 t 1 N μ 1 t 1 , σ 1 2 t 1 , X 2 t 1 N μ 2 t 1 , σ 2 2 t 1 , X 2 t 2 N μ 2 t 1 , σ 2 2 t 1
where N a , b 2 refers to the normal distribution with expectation a and variance b 2 .
Denoting by Z 1 , Z 2 , Z 3 three independent standard normal random variables, the pairwise covariances can be written as follows:
cov X 1 t 1 , X 2 t 1 = cov μ 1 t 1 + σ 1 t 1 Z 1 , μ 2 t 1 + σ 2 t 1 ρ 1.2 Z 1 + 1 ρ 1.2 2 Z 2 = σ 1 σ 2 ρ 1.2 t 1
cov X 2 t 1 , X 2 t 2 = cov μ 2 t 1 + σ 2 t 1 ρ 1.2 Z 1 + 1 ρ 1.2 2 Z 2 , μ 2 t 2 + σ 2 t 1 ρ 1.2 Z 1 + 1 ρ 1.2 2 Z 2 + σ 2 t 2 t 1 Z 3 = σ 2 2 t 1
cov X 1 t 1 , X 2 t 2 = cov μ 1 t 1 + σ 1 t 1 Z 1 , μ 2 t 2 + σ 2 t 1 ρ 1.2 Z 1 + 1 ρ 1.2 2 Z 2 + σ 2 t 2 t 1 Z 3 = σ 1 σ 2 ρ 1.2 t 1
where we have applied the bilinearity of the covariance operator, the independence of increments of Brownian motion, and the orthogonal decomposition of two-dimensional correlated Brownian motion. The correlation coefficients θ 1.2 , θ 1.3 , θ 2.3 in Proposition 1 ensue. Expanding the trivariate normal density function φ 1 x 1 , x 2 , x 3 as a product of normal conditional densities (Guillaume 2018), we obtain:
φ 1 x 1 , x 2 , x 3 = e 1 2 x 1 μ 1 t 1 σ 1 t 1 2 1 2 σ 2 1 2 x 2 μ 2 t 1 σ 2 t 1 θ 1.2 x 1 μ 1 t 1 σ 1 t 1 2 1 2 σ 3 1.2 2 x 3 μ 2 t 2 σ 2 t 2 θ 1.3 x 1 μ 1 t 1 σ 1 t 1 θ 2.3 1 σ 2 1 x 2 μ 2 t 1 σ 2 t 1 θ 1.2 x 1 μ 1 t 1 σ 1 t 1 2 2 π 3 / 2 σ 2 1 σ 3 1.2 σ 2 2 σ 1 t 1 t 2
where:
σ 2 1 = 1 θ 1.2 2 , θ 2.3 1 = θ 2.3 θ 1.2 θ 1.3 1 θ 1.2 2 , σ 3 1.2 = 1 θ 1.3 2 θ 2.3 1 2
The terms σ 2 1 , θ 2.3 1 and σ 3 1.2 have the following precise meanings:
-
σ 2 1 is the conditional standard deviation of X 2 t 1 given X 1 t 1 ;
-
θ 2.3 1 is the conditional correlation between X 2 t 1 and X 2 t 2 given X 1 t 1 ;
-
σ 3 1.2 is the conditional standard deviation of X 2 t 2 given X 1 t 1 and X 2 t 1 .
The rest of the proof, whose cumbersome details are omitted, then consists in solving the four integrals implied by (76). The final result takes the form of the linear combination of four N 3 functions written in Proposition 1.
Corollary 1 comes from the property of symmetry of Brownian paths.
Corollary 2 is a consequence of the fact that:
P RUU I μ 1 , μ 2 = min k 1 , h 1 min k 2 , h 2 min k 3 , h 2 φ 1 x 1 , x 2 , x 3 exp 2 h 1 x 1 h 1 σ 1 2 t 1 exp 2 h 2 x 1 x 2 h 2 σ 2 2 t 2 t 1 d x 3 d x 2 d x 1
Corollary 3 comes from the fact that the correlation coefficient between the random variables S 2 t 2 and S 2 t 3 is equal to t 2 t 3 .
Corollary 4 is a straightforward application of the law of total probability. □
Proof of Proposition 2. 
Using similar steps as in the proof of Proposition 1, one can express the problem at hand as the following integral:
P RUD μ 1 , μ 2 = min k 1 , h 1 max k 2 , h 2 max k 3 , h 2 φ 1 x 1 , x 2 , x 3 φ 2 x 1 φ 4 x 2 , x 3 d x 3 d x 2 d x 1
where the functions φ 1 and φ 2 are given by (86) and (80), respectively, and:
φ 4 x 2 , x 3 = P inf t 1 t t 2 X 2 t h 2 X 2 t 1 d x 2 , X 2 t 2 d x 3 d x 2 d x 3 = φ 3 x 2 , x 3
where the function φ 3 is given by (81).
Performing the necessary calculations, one can obtain the linear combination of four N 3 functions given in Proposition 2.
As in Proposition 1, Corollary 1 comes from the property of symmetry of Brownian paths.
Corollary 2 is a consequence of the fact that:
P RUD I μ 1 , μ 2 = min k 1 , h 1 max k 2 , h 2 max k 3 , h 2 φ 1 x 1 , x 2 , x 3 e 2 μ 1 σ 1 2 h 1 e 1 2 x 1 2 h 1 μ 1 t 1 σ 1 t 1 2 σ 1 2 π t 1 e 2 μ 2 σ 2 2 h 2 x 2 e 1 2 x 3 x 2 + 2 h 2 + μ 2 t 2 t 1 σ 2 t 2 t 1 2 σ 2 2 π t 2 t 1 d x 1 d x 2 d x 3
Proof of Proposition 3. 
One can express the problem at hand as the following integral:
P RRUU μ 1 , μ 2 = min k 1 , h 1 h 2 min k 2 , h 2 k 3 φ 2 x 1 φ 3 x 2 , x 3 φ 5 x 1 , x 2 , x 3 , x 4 d x 4 d x 3 d x 2 d x 1
where
φ 5 x 1 , x 2 , x 3 , x 4 = P X 1 t 1 d x 1 , X 2 t 1 d x 2 , X 2 t 2 d x 3 , X 1 t 2 d x 4 d x 4 d x 3 d x 2 d x 1
Plugging the quadrivariate normal joint density function of the set of random variables X 1 t 1 , X 2 t 1 , X 1 t 2 and X 2 t 2 , as a product of conditional density functions as explained in Guillaume (2018), and then factoring in the conditional cumulative distribution function of X 1 t 2 given the triple X 1 t 1 , X 2 t 1 , X 2 t 2 , Proposition 3 ensues. □
Proof of Proposition 4. 
Proof is given only for p 1 μ 1 , μ 2 , μ 3 and p 3 μ 1 , μ 2 , μ 3 , as p 2 μ 1 , μ 2 , μ 3 and p 4 μ 1 , μ 2 , μ 3 can then be deduced by the same symmetry argument as that already used in Corollary 1 of Proposition 1.
Following steps similar to the beginning of the proof of Proposition 1, one can express the problem at hand as the following two integrals:
p 1 μ 1 , μ 2 , μ 3 = x 1 = min k 1 , h 1 x 2 = min k 2 , h 2 x 3 = min k 3 , h 2 x 4 = k 4 φ 6 x 1 φ 7 x 1 , x 2 φ 8 x 2 , x 3 φ 9 x 1 , x 3 , x 4 d x 4 d x 3 d x 2 d x 1
p 3 μ 1 , μ 2 , μ 3 = x 1 = min k 1 , h 1 x 2 = max k 2 , h 2 x 3 = max k 3 , h 2 x 4 = k 4 φ 6 x 1 φ 7 x 1 , x 2 φ 10 x 2 , x 3 φ 9 x 1 , x 3 , x 4 d x 4 d x 3 d x 2 d x 1
where
φ 6 x 1 = P sup 0 t t 1 X 1 t h 1 , X 1 t 1 d x 1 d x 1
φ 7 x 1 , x 2 = P X 2 t 1 d x 2 X 1 t 1 d x 1 d x 1 d x 2
φ 8 x 2 , x 3 = P sup t 1 t t 2 X 2 t h 2 , X 2 t 2 d x 3 X 2 t 1 d x 2 d x 2 d x 3
φ 10 x 2 , x 3 = P inf t 1 t t 2 X 2 t h 2 , X 2 t 2 d x 3 X 2 t 1 d x 2 d x 2 d x 3
φ 9 x 1 , x 3 , x 4 = P X 3 t 2 d x 4 X 2 t 2 d x 3 , X 1 t 1 d x 1
The function φ 6 is obtained by differentiating the classical formula for the joint cumulative distribution of the maximum of a Brownian motion with drift and its endpoint over a closed time interval (see, e.g., Karatzas and Shreve 2000):
φ 6 x 1 = e 1 2 x 1 μ 1 t 1 σ 1 t 1 2 σ 1 2 π t 1 e 2 μ 1 σ 1 2 h 1 e 1 2 x 1 2 h 1 μ 1 t 1 σ 1 t 1 2 σ 1 2 π t 1 d x 1
The function φ 7 is easily derived from the bivariate normality of the pair X 1 t 1 , X 2 t 1 :
φ 7 x 1 , x 2 = e 1 2 1 ρ 1.2 2 x 2 μ 2 t 1 σ 2 t 1 ρ 1.2 x 1 μ 1 t 1 σ 1 t 1 2 σ 2 2 π t 1 1 ρ 1.2 2 d x 1 d x 2
To handle the function φ 8 , we notice that, by conditioning w.r.t. the filtration at time t 1 the same classical formula as the one used to derive φ 6 , we can obtain:
P sup t 1 t t 2 S 2 t H 2 , S 2 t 2 S 2 0 e x 3 S 2 t 1 = S 2 0 e x 2
= N ln S 2 0 e x 3 S 2 0 e x 2 μ 2 t 2 t 1 σ 2 t 2 t 1 H 2 S 2 0 e x 2 2 μ 2 σ 2 2 N ln S 2 0 e x 3 S 2 0 e x 2 2 ln H 2 S 2 0 e x 2 μ 2 t 2 t 1 σ 2 t 2 t 1
for any given x 2 , x 3 2 and H 2 > S 2 0 e x 3 .
Equation (103) can be rewritten as follows:
P sup t 1 t t 2 X 2 t h 2 , X 2 t 2 x 3 X 2 t 1 d x 2
= N x 3 x 2 μ 2 t 2 t 1 σ 2 t 2 t 1 exp 2 μ 2 σ 2 2 h 2 x 2 N x 3 x 2 2 h 2 x 2 μ 2 t 2 t 1 σ 2 t 2 t 1
Therefore, by differentiating (104) w.r.t. x 3 , we obtain:
φ 8 x 2 , x 3 = e 1 2 x 3 x 2 μ 2 t 2 t 1 σ 2 t 2 t 1 2 σ 2 2 π t 2 t 1 e 2 μ 2 σ 2 2 h 2 x 2 e 1 2 x 3 + x 2 2 h 2 μ 2 t 2 t 1 σ 2 t 2 t 1 2 σ 2 2 π t 2 t 1 d x 2 d x 3
By the symmetry of paths of Brownian motion, we have:
P inf t 1 t t 2 X 2 t h 2 , X 2 t 2 x 3 X 2 t 1 d x 2
= N x 3 + x 2 + μ 2 t 2 t 1 σ 2 t 2 t 1 exp 2 μ 2 σ 2 2 h 2 x 2 N x 3 + x 2 + 2 h 2 x 2 + μ 2 t 2 t 1 σ 2 t 2 t 1
Hence,
φ 10 x 2 , x 3 = φ 8 x 2 , x 3
The function φ 9 derives from the joint trivariate normality of the triple X 1 t 1 , X 2 t 2 , X 3 t 2 . The marginal distributions of the elements of this triple come from the known marginal distributions of S 1 t 1 , S 2 t 2 and S 3 t 2 . The pairwise correlations, as given by θ 1.3 , θ 1.4 and θ 3.4 in Proposition 4 can be easily determined using the same method as in the proof of Proposition 3. We obtain:
φ 9 x 1 , x 3 , x 4 = e 1 2 ϕ 4 1.3 2 x 4 μ 3 t 2 σ 3 t 2 θ 1.4 x 1 μ 1 t 1 σ 1 t 1 θ 3.4 1 x 3 μ 2 t 2 σ 2 t 2 θ 1.3 x 1 μ 1 t 1 σ 1 t 1 2 ϕ 4 1.3 σ 3 2 π t 2
where σ 3 1 = 1 θ 1.3 2 .
The rest of the proof, whose cumbersome details are omitted, then consists in solving the integrals implied by (94) and (95). The final result can be expressed as the linear combination of four Ψ 4 functions written in Proposition 4. The origin of the function Ψ 4 , which is a special form of quadrivariate normal cumulative distribution, lies in the FDD (Finite Dimensional Distribution) of the quadruple S 1 t 1 , S 2 t 1 , S 2 t 2 , S 3 t 2 . Indeed, a little algebra shows that, D 1 , D 2 , D 3 , D 4 + , we have:
P S 1 t 1 D 1 , S 2 t 1 D 2 , S 2 t 2 D 3 , S 3 t 2 D 4
= Ψ 4 ln D 1 S 1 0 μ 1 t 1 σ 1 t 1 , ln D 2 S 2 0 μ 2 t 1 σ 2 t 1 , ln D 3 S 2 0 μ 2 t 2 σ 2 t 2 , ln D 4 S 3 0 μ 3 t 2 σ 3 t 2 ; θ 1.2 , θ 1.3 , θ 1.4 , θ 2.3 , θ 3.4
Corollary 1 is a consequence of the fact that:
p 1 I μ 1 , μ 2 = x 1 = min k 1 , h 1 x 2 = min k 2 , h 2 x 3 = min k 3 , h 2 x 4 = k 4 e 2 μ 1 σ 1 2 h 1 e 1 2 x 1 2 h 1 μ 1 t 1 σ 1 t 1 2 σ 1 2 π t 1 φ 7 x 1 , x 2 e 2 μ 2 σ 2 2 h 2 x 2 e 1 2 x 3 + x 2 2 h 2 μ 2 t 2 t 1 σ 2 t 2 t 1 2 σ 2 2 π t 2 t 1 φ 9 x 1 , x 3 , x 4 d x 4 d x 3 d x 2 d x 1
and:
p 3 I μ 1 , μ 2 = x 1 = min k 1 , h 1 x 2 = max k 2 , h 2 x 3 = max k 3 , h 2 x 4 = k 4 e 2 μ 1 σ 1 2 h 1 e 1 2 x 1 2 h 1 μ 1 t 1 σ 1 t 1 2 σ 1 2 π t 1 φ 7 x 1 , x 2 e 2 μ 2 σ 2 2 h 2 x 2 e 1 2 x 3 x 2 + 2 h 2 + μ 2 t 2 t 1 σ 2 t 2 t 1 2 σ 2 2 π t 2 t 1 φ 9 x 1 , x 3 , x 4 d x 4 d x 3 d x 2 d x 1
Corollary 2 comes from the fact that the correlation coefficient between the random variables S 2 t 2 and S 3 t 3 is equal to ρ 2.3 t 2 t 3 . □
Proof of Proposition 5. 
Following steps similar to the beginnings of the previous proofs, one can express the problem at hand as the following integral:
P RDKO μ 1 , μ 2
= x 1 = h 2 min k 1 , h 1 x 2 = h 4 h 3 x 3 = h 4 min k 2 , h 3 φ 1 x 1 , x 2 , x 3 φ 11 x 1 , x 2 φ 12 x 2 , x 3 d x 3 d x 2 d x 1
where the function φ 1 is given by Equation (86) and:
φ 11 x 1 , x 2 = P sup 0 t t 1 X 1 t h 1 , inf 0 t t 1 X 1 t h 2 X 1 t 1 d x 1 d x 2 d x 1
φ 12 x 2 , x 3
= P sup t 1 t t 2 X 2 t h 3 , inf t 1 t t 2 X 2 t h 4 X 2 t 1 d x 2 , X 2 t 2 d x 3 d x 3 d x 2
From Pötzelberger and Wang (2001), one can plug:
φ 11 x 1 , x 2 = n = e 2 n a 1 x 1 n a 1 σ 1 2 t 1 e 2 h 1 n a 1 x 1 h 1 + n a 1 σ 1 2 t 1
φ 12 x 2 , x 3 = n = e 2 n a 2 x 3 x 2 n a 2 σ 2 2 t 2 t 1 e 2 h 3 x 2 n a 2 x 3 h 3 + n a 2 σ 2 2 t 2 t 1
The bulk of the proof, whose cumbersome details are omitted, then consists of solving the sixteen integrals implied by (112). The final result takes the form of the linear combinations of double sums of N 3 functions in Proposition 5.
An elementary adjustment identical to the one in Corollary 3 of Proposition 1 allows to value an early-ending variant of P RDKO μ 1 , μ 2 .
Alternatively, one can also expand the problem as the following integral:
P RDKO μ 1 , μ 2 = x 1 = h 2 min k 1 , h 1 x 2 = h 4 h 3 x 3 = h 4 min k 2 , h 3 φ 13 x 1 φ 7 x 1 , x 2 φ 14 x 2 , x 3 d x 3 d x 2 d x 1
where the function φ 7 is given by Equation (102) and:
φ 13 x 1 = P sup 0 t t 1 X 1 t h 1 , inf 0 t t 1 X 1 t h 2 , X 1 t 1 d x 1 d x 1
φ 14 x 2 , x 3
= P sup t 1 t t 2 X 2 t h 3 , inf t 1 t t 2 X 2 t h 4 , X 2 t 2 d x 3 X 2 t 1 d x 2 d x 3 d x 2
According to the classical formula for the distribution of the maximum, the minimum and the endpoint of a Brownian motion with drift over a closed time interval, which can be traced back to Anderson (1960), we have:
  P sup 0 t t 1 S 1 t < H 1 , inf 0 t t 1 S 1 t > H 2 , S 1 t 1 < S 1 0 e x 1 S 1 0
= n 1 = exp 2 μ 1 σ 1 2 n 1 a 1 N x 1 2 n 1 a 1 μ 1 t 1 σ 1 t 1 N h 2 2 n 1 a 1 μ 1 t 1 σ 1 t 1
n 1 = exp 2 μ 1 σ 1 2 h 2 n 1 a 1 N x 1 2 h 2 + 2 n 1 a 1 μ 1 t 1 σ 1 t 1 N h 2 + 2 n 1 a 1 μ 1 t 1 σ 1 t 1
for a given x 1 , and H 1 > S 1 0 e x 1 .
Mere differentiation of (120) w.r.t. x 1 yields the function φ 13 :
φ 13 x 1 = n 1 = e 2 μ 1 σ 1 2 n 1 a 1 1 2 σ 2 t 1 x 1 μ 1 t 1 2 n 1 a 1 2 σ 1 2 π t 1 n 1 = e 2 μ σ 1 2 h 2 n 1 a 1 1 2 σ 1 2 t 1 x 1 2 h 2 μ 1 t 1 + 2 n 1 a 1 2 σ 1 2 π t 1
To handle the function φ 14 , we notice that, by conditioning w.r.t. the filtration at time t 1 the same classical formula as the one used to derive φ 13 , we can obtain:
  P sup t 1 t t 2 S 2 t < H 3 , inf t 1 t t 2 S 2 t > H 4 , S 2 t 2 S 2 0 e x 3 S 2 t 1 = S 2 0 e x 2
= n 2 = exp 2 μ 2 σ 2 2 n 2 a 2 N ln S 2 0 e x 3 S 2 0 e x 2 2 n 2 a 2 μ 2 t 2 t 1 σ 2 t 2 t 1 N ln H 4 S 2 0 e x 2 2 n 2 a 2 μ 2 t 2 t 1 σ 2 t 2 t 1
n 2 = exp 2 μ 2 σ 2 2 ln H 4 S 2 0 e x 2 n 2 a 2 N ln S 2 0 e x 3 S 2 0 e x 2 2 ln H 4 S 2 0 e x 2 + 2 n 2 a 2 μ 2 t 2 t 1 σ 2 t 2 t 1 N ln H 4 S 2 0 e x 2 + 2 n 2 a 2 μ 2 t 2 t 1 σ 2 t 2 t 1
for any given x 2 , x 3 2 and H 3 > S 2 0 e x 3 .
Equations (122) and (123) can be rewritten as follows:
  P sup t 1 t t 2 X 2 t < h 3 , inf t 1 t t 2 X 2 t > h 4 , X 2 t 2 x 3 X 2 t 1 d x 2
= n 2 = exp 2 μ 2 σ 2 2 n 2 a 2 N x 3 x 2 2 n 2 a 2 μ 2 t 2 t 1 σ 2 t 2 t 1 N h 4 x 2 2 n 2 a 2 μ 2 t 2 t 1 σ 2 t 2 t 1
n 2 = exp 2 μ 2 σ 2 2 h 4 x 2 n 2 a 2 N x 3 x 2 2 h 4 x 2 + 2 n 2 a 2 μ 2 t 2 t 1 σ 2 t 2 t 1 N h 4 x 2 + 2 n 2 a 2 μ 2 t 2 t 1 σ 2 t 2 t 1
Therefore, by differentiating (124) and (125) w.r.t. x 3 , we obtain:
φ 14 x 2 , x 3 = n 2 = e 2 μ 2 σ 2 2 n 2 a 2 e 1 2 σ 2 2 t 2 t 1 x 3 x 2 2 n 2 a 2 μ 2 t 2 t 1 2 σ 2 2 π t 2 t 1
n 2 = e 2 μ 2 σ 2 2 h 4 x 2 n 2 a 2 e 1 2 σ 2 2 t 2 t 1 x 3 + x 2 2 h 4 + 2 n 2 a 2 μ 2 t 2 t 1 2 σ 2 2 π t 2 t 1
This second formulation leads to a formula identical to Proposition 5 except for the fact that the N 3 functions are replaced by Φ 3 functions defined as follows:
Φ 3 b 1 , b 2 , b 3 ; x = x = b 2 exp x 2 / 2 2 π N b 1 c 1 x 1 c 1 2   N b 3 c 2 x 1 c 2 2 d x
where b 1 , b 2 , b 3 3 and x is a vector with two real coordinates c 1 , c 2 1 , 1 .
The vectors of correlation coefficients x 1 and x 2 become:
x 1 = ρ 1.2 , t 1 t 2 , x 2 = ρ 1.2 , t 1 t 2
The origin of the function Φ 3 , which is a special form of trivariate normal cumulative distribution, lies in the FDD (finite dimensional distribution) of the triple S 1 t 1 , S 2 t 1 , S 2 t 2 . Indeed, a little algebra shows that, D 1 , D 2 , D 3 + , we have:
P S 1 t 1 D 1 , S 2 t 1 D 2 , S 2 t 2 D 3
= Φ 3 ln D 1 S 1 0 μ 1 t 1 σ 1 t 1 , ln D 2 S 2 0 μ 2 t 1 σ 2 t 1 , ln D 3 S 2 0 μ 2 t 2 σ 2 t 2 ; ρ 1.2 , t 1 t 2
Notice that the two-colour probability distributions of Section 2.1 and Section 2.2 can also be written as linear combinations of functions Φ 3 . □

4. Conclusions

This article has shown how to value in closed form an important kind of multi-asset step barrier option known as a rainbow step barrier option, under the condition that the number of “colours” is restricted to two, along with widespread variants such as a two-colour outside step barrier and a two-colour step double barrier. It may be feasible, albeit tedious, to find an analytical solution to an extended valuation problem with three or four colours, but the expected benefits, compared with a conditional Monte Carlo approximation method, would greatly depend on the degree of the quadrature required to numerically evaluate the resulting multidimensional integrals. It should be emphasised that, even if more sophisticated models (allowing, e.g., for stochastic volatility) or a greater number of underlying assets are needed, closed form solutions obtained in a low-dimensional Black–Scholes framework remain useful as fast and accurate benchmarks that can: (i) serve as control variates in a simulation; (ii) speed up the calibration process; (iii) facilitate the analysis and the understanding of the interactions between the variables, as well as of the sensitivities of the option value w.r.t. its main parameters, which is instrumental in devising appropriate hedging techniques or trading strategies.

Funding

This research received no external funding.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Anderson, Todd W. 1960. A modification of the sequential probability ratio test to reduce the sample size. Annals of Mathematical Statistics 31: 165–97. [Google Scholar] [CrossRef]
  2. Berntsen, Jarle, Terje O. Espelid, and Alan Genz. 1991. An adaptive algorithm for the approximate calculation of multiple integrals. ACM Transactions on Mathematical Software 17: 437–51. [Google Scholar] [CrossRef]
  3. Bielecki, Tomasz R., and Marek Rutkowski. 2004. Credit Risk: Modeling, Valuation and Hedging. Berlin/Heidelberg: Springer. [Google Scholar]
  4. Bouzoubaa, Mohamed, and Adel Osseiran. 2010. Exotic Options and Hybrids. Hoboken: Wiley Finance. [Google Scholar]
  5. Chang, Der-Chen, Eric C. Chang, Haitao Fan, and Duy-Minh Nhieu. 2005. Mathematical analysis of the two-color partial rainbow options. Applicable Analysis 84: 737–57. [Google Scholar] [CrossRef]
  6. Das, Satyajit. 2006. Structured Products: Volume 1. Hoboken: Wiley Finance, John Wiley and Sons (Asia). [Google Scholar]
  7. Davis, Philip J., and Philip Rabinowitz. 2007. Methods of Numerical Integration. New York: Dover Publications. [Google Scholar]
  8. Gao, Rong, and Xiaoli Wu. 2022. Pricing rainbow option for uncertain financial market. RAIRO Operations Research 56: 3973–89. [Google Scholar] [CrossRef]
  9. Genz, Alan. 2004. Numerical computation of rectangular bivariate and trivariate normal and t probabilities. Statistics and Computing 14: 151–60. [Google Scholar] [CrossRef]
  10. Glasserman, Paul. 2003. Monte Carlo Methods in Financial Engineering. New York: Springer. [Google Scholar]
  11. Guillaume, Tristan. 2001. Analytical valuation of options on joint minima and maxima. Applied Mathematical Finance 8: 209–35. [Google Scholar] [CrossRef]
  12. Guillaume, Tristan. 2015. On the computation of the survival probability of Brownian motion with drift in a closed time interval when the absorbing boundary is a step function. Journal of Probability and Statistics 2015: 391681. [Google Scholar] [CrossRef]
  13. Guillaume, Tristan. 2016. Computation of the survival probability of Brownian motion with drift in a closed time interval when the absorbing boundary is an affine or an exponential function of time. International Journal of Statistics and Probability 5: 119–38. [Google Scholar] [CrossRef]
  14. Guillaume, Tristan. 2018. Computation of the Quadrivariate and Pentavariate normal cumulative distribution functions. Communications in Statistics—Simulation and Computation 47: 839–51. [Google Scholar] [CrossRef]
  15. Harrison, J. Michael, and David M. Kreps. 1979. Martingales and arbitrage in multiperiod securities markets. Journal of Economic Theory 20: 381–408. [Google Scholar] [CrossRef]
  16. Harrison, J. Michael, and Stanley R. Pliska. 1981. Martingales and stochastic integrals in the theory of continuous trading. Stochastic Processes and Their Applications 11: 312–16. [Google Scholar] [CrossRef]
  17. Heynen, Ronald C., and Harry M. Kat. 1994. Crossing barriers. Risk Magazine 7: 46–50. [Google Scholar]
  18. Karatzas, Ioannis, and Steven E. Shreve. 2000. Brownian Motion and Stochastic Calculus. New York: Springer. [Google Scholar]
  19. Kwok, Yue-Kuen, Lixin Wu, and Hong Yu. 1998. Pricing Multi-Asset Options with an External Barrier. International Journal of Theoretical and Applied Finance 1: 523–41. [Google Scholar] [CrossRef]
  20. Lee, Hangsuck, Bangwon Ko, and Seongjoo Song. 2019a. Valuing step barrier options and their icicled variations. The North American Journal of Economics and Finance 49: 396–411. [Google Scholar] [CrossRef]
  21. Lee, Hangsuck, Gaeun Lee, and Seongjoo Song. 2022. Multi-step reflection principle and barrier options. The Journal of Futures Markets 42: 692–721. [Google Scholar] [CrossRef]
  22. Lee, Hangsuck, Hongjun Ha, and Minha Lee. 2021. Valuation of piecewise linear barrier options. The North American Journal of Economics and Finance 58: 101470. [Google Scholar] [CrossRef]
  23. Lee, Hangsuck, Soohan Ahn, and Bangwon Ko. 2019b. Generalizing the reflection principle of Brownian motion and closed form pricing of barrier options and autocallable investments. The North American Journal of Economics and Finance 50: 1–13. [Google Scholar] [CrossRef]
  24. Matsumoto, Makoto, and Takuji Nishimura. 1998. Mersenne Twister: A 623-dimensionally equidistributed uniform pseudorandom number generator. ACM Transactions on Modeling and Computer Simulation 8: 3–30. [Google Scholar] [CrossRef]
  25. Pötzelberger, Klaus, and Liqun Wang. 2001. Boundary crossing probability for Brownian motion. Journal of Applied Probability 38: 152–64. [Google Scholar] [CrossRef]
  26. Wang, Liqun, and Klaus Pötzelberger. 1997. Boundary crossing probability for Brownian motion and general boundaries. Journal of Applied Probability 34: 54–65. [Google Scholar] [CrossRef]
Table 1. Two-colour up-and-up knock-out put.
Table 1. Two-colour up-and-up knock-out put.
ρ 1.2 = 0.6 ρ 1.2 = 0.2 ρ 1.2 = 0.2 ρ 1.2 = 0.6
σ 1 = σ 2 = 20 %
H 1 = H 2 = 115
1.189
(1.123,1.168,
1.182)
2.163
(2.237, 2.141,
2.161)
3.134
(3.082, 3.116,
3.132)
4.141
(4.174, 4.127,
4.140)
σ 1 = σ 2 = 60 %
H 1 = H 2 = 125
4.065
(4.109, 4.072,
4.061)
6.440
(6.392, 6.449,
6.442)
8.769
(8.734, 8.755,
8.763)
11.213
(11.278, 11.196,
11.218)
σ 1 = 20 % , σ 2 = 60 %
H 1 = 115 , H 2 = 125
4.454
(4.411, 4.443,
4.454)
7.260
(7.355, 7.301,
7.264)
9.979
(9.912, 9.101,
9.975)
12.844
(12.957, 12.892,
12.848)
σ 1 = 60 % , σ 2 = 20 %
H 1 = 125 , H 2 = 115
1.096
(1.082, 1.108,
1.091)
1.921
(1.107, 1.953,
1.928)
2.753
(2.796, 2.745,
2.752)
3.617
(3.692, 3.599,
3.614)
Table 2. Two-colour down-and-up knock-out put.
Table 2. Two-colour down-and-up knock-out put.
ρ 1.2 = 0.6 ρ 1.2 = 0.2 ρ 1.2 = 0.2 ρ 1.2 = 0.6
σ 1 = σ 2 = 20 %
H 1 = 85 , H 2 = 115
4.299
(4.287, 4.305,
4.297)
3.293
(3.318, 3.286,
3.292)
2.307
(2.282, 2.298,
2.305)
1.303
(1.284, 1.309,
1.304)
σ 1 = σ 2 = 60 %
H 1 = 75 , H 2 = 125
10.221
(10.142, 10.228,
10.224)
7.610
(7.563, 7.597,
7.612)
5.212
(5.255, 5.204,
5.214)
2.862
(2.834, 2.854,
2.865)
σ 1 = 20 % , σ 2 = 60 %
H 1 = 85 , H 2 = 125
13.377
(13.385, 13.391,
13.376)
10.496
(10.472, 10.482,
10.948)
7.735
(7.783, 7.717,
7.731)
4.857
(4.894, 4.866,
4.857)
σ 1 = 60 % , σ 2 = 20 %
H 1 = 75 , H 2 = 115
3.344
(3.387, 3.358,
3.346)
2.401
(2.383, 2.413,
2.402)
1.545
(1.596, 1.530,
1.542)
0.748
(0.884, 0.787,
0.752)
Table 3. Outside two-colour up-and-down knock-out call.
Table 3. Outside two-colour up-and-down knock-out call.
ρ 1.2 = 0.6 ,
ρ 1.3 = ρ 2.3 = 0.4
ρ 1.2 = 0.6 ,
ρ 1.3 = ρ 2.3 = 0.4
ρ 1.2 = 0.6
ρ 1.3 = ρ 2.3 = 0.4
ρ 1.2 = 0.6
ρ 1.3 = ρ 2.3 = 0.4
σ 1 = σ 2 = 20 %
σ 3 = 20 %
H 1 = 115 , H 2 = 85
2.378
(2.452, 2.361,
2.375)
2.772
(2.914, 2.812,
2.779)
1.717
(1.585, 1.731,
1.720)
1.182
(1.193, 1.178,
1.180)
σ 1 = σ 2 = 20 %
σ 3 = 60 %
H 1 = 115 , H 2 = 85
5.769
(5.728, 5.781,
5.764)
7.053
(7.137, 7.036,
7.054)
4.522
(4.534, 4.541,
4.524)
2.897
(2.852, 2.923,
2.893)
σ 1 = σ 2 = 60 %
σ 3 = 20 %
H 1 = 125 , H 2 = 75
1.627
(1.592, 1.614,
1.628)
2.783
(2.848, 2.767,
2.788)
1.351
(1.320, 1.365,
1.353)
0.849
(0.915, 0.828,
0.842)
σ 1 = σ 2 = 60 %
σ 3 = 60 %
H 1 = 125 , H 2 = 75
3.864
(3.814, 3.872,
3.865)
7.554
(7.518, 7.535,
7.558)
3.823
(3.856, 3.829,
3.827)
2.076
(2.011, 2.091, 2.072)
σ 1 = 20 % , σ 2 = 60 %
σ 3 = 20 %
H 1 = 115 , H 2 = 75
1.534
(1.502, 1.526,
1.535)
2.573
(2.495, 2.556,
2.577)
1.188
(1.207, 1.179,
1.185)
0.621
(0.774, 0.684,
0.613)
σ 1 = 20 % , σ 2 = 60 %
σ 3 = 60 %
H 1 = 115 , H 2 = 75
3.629
(3.787, 3.662,
3.621)
6.697
(6.724, 6.684,
6.692)
3.217
(3.051, 3.252,
3.221)
1.518
(1.586, 1.476,
1.511)
σ 1 = 60 % , σ 2 = 20 %
σ 3 = 20 %
H 1 = 125 , H 2 = 85
2.572
(2.734, 2.548,
2.567)
2.989
(3.125, 3.016,
2.994)
1.931
(2.071, 1.965,
1.938)
1.496
(1.634, 1.454,
1.489)
σ 1 = 60 % , σ 2 = 20 %
σ 3 = 60 %
H 1 = 125 , H 2 = 85
6.248
(6.304, 6.237,
6.241)
7.953
(8.060, 7.984,
7.961)
5.323
(5.212, 5.348,
5.324)
3.707
(3.569, 3.726,
3.702)
Table 4. Two-colour rainbow double knock-out put.
Table 4. Two-colour rainbow double knock-out put.
ρ 1.2 = 0.6 ρ 1.2 = 0.2 ρ 1.2 = 0.2 ρ 1.2 = 0.6
σ 1 = σ 2 = 15 %
H 1 = H 3 = 120
H 2 = H 4 = 80
2.919
(2.985, 2.956,
2.916)
2.938
(2.792, 2.883,
2.932)
2.936
(2.974, 2.918,
2.931)
3.027
(3.191, 3.088,
3.032)
σ 1 = σ 2 = 30 %
H 1 = H 3 = 130
H 2 = H 4 = 70
4.370
(4.529, 4.281,
4.365)
5.791
(5.978, 5.697,
5.786)
4.427
(4.196, 4.443,
4.426)
5.175
(5.329, 5.092,
5.158)
σ 1 = 15 % , σ 2 = 30 %
H 1 = 120 , H 3 = 130
H 2 = 80 , H 4 = 70
4.726
(4.594, 4.752,
4.728)
4.744
(4.868, 4.785,
4.749)
4.724
(4.574, 4.771,
4.728)
5.042
(5.226, 5.018,
5.045)
σ 1 = 30 % , σ 2 = 15 %
H 1 = 130 , H 3 = 120
H 2 = 70 , H 4 = 80
2.614
(2.429, 2.576,
2.610)
2.683
(2.872, 2.612,
2.679)
2.771
(2.602, 2.742,
2.773)
2.942
(3.165, 3.036,
2.953)
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

Guillaume, T. Rainbow Step Barrier Options. J. Risk Financial Manag. 2024, 17, 356. https://doi.org/10.3390/jrfm17080356

AMA Style

Guillaume T. Rainbow Step Barrier Options. Journal of Risk and Financial Management. 2024; 17(8):356. https://doi.org/10.3390/jrfm17080356

Chicago/Turabian Style

Guillaume, Tristan. 2024. "Rainbow Step Barrier Options" Journal of Risk and Financial Management 17, no. 8: 356. https://doi.org/10.3390/jrfm17080356

APA Style

Guillaume, T. (2024). Rainbow Step Barrier Options. Journal of Risk and Financial Management, 17(8), 356. https://doi.org/10.3390/jrfm17080356

Article Metrics

Back to TopTop