We can now begin to value two-colour step barrier options in the following order:
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
defined by:
where the acronym “
” 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 is given by:where the are given by (3) and: - -
is the trivariate standard normal cumulative distribution function with correlation coefficients
Corollary 1. It suffices to multiply by all the first three arguments of each function and to substitute each operator by a operator in Proposition 1 to obtain an exact formula for defined as:where the acronym “” 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
and defined by: Corollary 3. It suffices to substitute each argument in each function of Proposition 1 by , , to obtain an exact formula for the early-ending variant defined by: Corollary 4. Let be the value of when the value of becomes “very high”, i.e., high enough for the probability to tend to zero; then, the difference provides the value of the following minor variant: 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
, is given by:
where
is the riskless interest rate assumed to be constant,
is the indicator function and
is the set constructed by the intersection of elements of the
algebra generated by the pair of processes
that characterises the probability
as given by the arguments of the probability operator in (5).
A simple application of the Cameron–Martin–Girsanov theorem yields:
where
is the measure under which
is a standard Brownian motion (the classical so-called risk-neutral measure), while
is the measure under which
and
are two independent standard Brownian motions.
To factor in a continuous dividend rate associated with each asset , simply replace by .
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
probability along with the pairs
and
.
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
and
, and different values of the knock-out barriers. All the initial values of the underlying assets
and the strike prices
are set at 100. Expiry is 1 year. The two time intervals
and
have equal length, i.e.,
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 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
and
, 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
is three times greater than when
. 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
and
, 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
and
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
and
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 denote the joint cumulative distribution function defined by:where the acronym “” stands for “Rainbow Up and Down”. Then, the exact value of
is given by:
where all the notations are identical, as in Proposition 1.
Corollary 1. It suffices to multiply by all the first three arguments of each function and substitute each operator by a operator as well as each operator by a operator in Proposition 2 to obtain an exact formula for defined as: 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 and defined by: 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
or
probability along with the pairs
and
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
and
, 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
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
and
.
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
) instead of the second and last one (asset
): 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
the following cumulative joint distribution at the core of reverse rainbow option valuation:
where the acronym “
” stands for “Reverse Rainbow Up and Up”.
Then, Proposition 3 provides the exact value of in the form of a triple integral.
Proposition 3. - -
is the univariate standard normal cumulative distribution function;
- -
the functions and 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 functions, according to the considered combination of events.
Remark 2. is the partial correlation between and conditional on , while is the partial correlation between and conditional on and , and is the conditional standard deviation of given , and .
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
, is given by:
where
- -
is the set constructed by the intersection of elements of the algebra generated by the pair of processes that characterises the probability as given by the arguments of the probability operator in (26);
- -
is the measure under which 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
is introduced, w.r.t. which the option’s moneyness is measured at expiry, while the assets
and
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
with the following differential:
The instantaneous pairwise correlations between the Brownian motions are denoted as .
The objective is to compute the probabilities
defined by:
Let be a vector of five coordinates where each , .
Let the function
,
, be defined by:
The following Proposition 4 combines all the probabilities defined in (32)–(35) into a single formula.
Proposition 4. The exact values of the probabilities , written in shorter notation as , are given by:where are as in Proposition 2, , and we have: Corollary 1. The corresponding knock-in probabilities can be inferred in the same way as in Proposition 1 and Proposition 2. Let the probabilities be defined by: Then, is given by (40).
Corollary 2. It suffices to substitute each argument in each function of Proposition 4 by , , to obtain an exact formula for the early-ending variant of .
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
, is given by:
where
is the set constructed by the intersection of elements of the
algebra generated by the pair of processes
that characterises the probability
as given by the arguments of the probability operator in (32), and the acronym “
” stands for “Outside Rainbow Up and Up”.
Using the following orthogonal decomposition of Brownian motion
:
where:
and
is a basis of three independent Brownian motions (
Guillaume 2018), the multidimensional Cameron-Martin-Girsanov theorem yields:
where:
The measure is the measure under which , and are three independent standard Brownian motions.
A simple and robust numerical evaluation of the function
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
,
, and
, and different values of the knock-out barriers. The parameters
,
,
,
, and
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
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
and
is negative and the correlation between
and both
and
is positive. The lowest option values are when the correlation between
and
is positive and the correlation between
and both
and
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
is high. This remains true under very different combinations of values for all the other parameters (volatilities of
and
, 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
, whereas the value of a rainbow step barrier option is not a monotonically increasing function of
, just like the value of a reverse rainbow step barrier option is not a monotonically increasing function of
. 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,
and
, are only concerned with barrier crossing during the option life, and the third one,
, 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
and
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
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
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
and
denote an upward and a downward barrier, respectively, on the time interval
. Similarly,
and
represent an upward and a downward barrier, respectively, on the time interval
. As in the previous sections, barrier crossing is monitored w.r.t. process
following Equation (1) on
and w.r.t. process
following Equation (2) on
. Our objective now is to find the value of the joint cumulative distribution function
defined by:
where the acronym “
” 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 is given by:where: 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
, is given by:
where the parameters
and
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
and
, and different values of the knock-out barriers. Expiry is 6 months and
is one quarter of a year. The parameters
,
, and
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
convergence is attained by stopping at 8, the number of iterations controlled by the absolute values of
and
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
function introduced in
Section 3 is slightly faster than the one using the trivariate standard normal cumulative distribution function
, 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
to
is followed by a relatively significant decrease in value from
to
, before a new increase in value from
to
. 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.