## Evaluating Microscopic and Macroscopic Models for Derivative Contracts on Commodity Indices ∗

Alberto Pedro Manzano-Herrero † Emanuele Nastasi ‡ Andrea Pallavicini ¶

§

Carlos Vázquez

First Version: June 20, 2024. This version: August 5, 2024

## Abstract

In this article, we analyze two modeling approaches for the pricing of derivative contracts on a commodity index. The first one is a microscopic approach, where the components of the index are modeled individually, and the index price is derived from their combination. The second one is a macroscopic approach, where the index is modeled directly. While the microscopic approach offers greater flexibility, its calibration results to be more challenging, thus leading practitioners to favor the macroscopic approach. However, in the macroscopic model, the lack of explicit futures curve dynamics raises questions about its ability to accurately capture the behavior of the index and its sensitivities. In order to investigate this, we calibrate both models using derivatives of the S&amp;P GSCI Crude Oil excess-return index and compare their pricing and sensitivities on pathdependent options, such as autocallable contracts. This research provides insights into the suitability of macroscopic models for pricing and hedging purposes in real scenarios.

JEL classification codes: C63, G13.

AMS classification codes: 65C05, 91G20, 91G60.

Keywords: Commodity Futures, Commodity Indices, Option Pricing, Stochastic Local Volatility.

∗ The opinions here expressed are solely those of the authors and do not represent in any way those of their employers.

† Department of Mathematics and CITIC, University of A Coruña, Spain alberto.manzano.herrero@udc.es ‡ swissQuant Group AG, Zurich, Switzerland nastasi@swissquant.com .

§ Intesa SanPaolo Milan, Italy, andrea.pallavicini@intesasanpaolo.com .

¶ Department of Mathematics and CITIC, University of A Coruña, Spain, carlos.vazquez.cendon@udc.es

## Contents

| 1 Introduction   | 1 Introduction             | 1 Introduction                                                                     |   3 |
|------------------|----------------------------|------------------------------------------------------------------------------------|-----|
| 2                | S&P GSCI indices           | S&P GSCI indices                                                                   |   4 |
|                  | 2.1                        | Evolution during a non-rolling period . . . . . . . .                              |   5 |
|                  | 2.2                        | Evolution during a rolling period . . . . . . . . . .                              |   6 |
| 3                | Modeling framework         | Modeling framework                                                                 |   7 |
|                  | 3.1                        | Micro model . . . . . . . . . . . . . . . . . . . . . .                            |   7 |
|                  | 3.2                        | Macro model . . . . . . . . . . . . . . . . . . . . .                              |   8 |
| 4                | Calibration and simulation | Calibration and simulation                                                         |   8 |
|                  | 4.1                        | Calibration of the micro model . . . . . . . . . . .                               |   9 |
|                  |                            | 4.1.1 The leverage functions . . . . . . . . . . . .                               |   9 |
|                  |                            | 4.1.2 Stochastic parameters and correlation structure . . . .                      |  11 |
|                  | 4.2                        | Calibration of the macro model . . . . . . .                                       |  13 |
|                  |                            | 4.2.1 The leverage function . . . . . . . . . . . .                                |  13 |
|                  |                            | 4.2.2 The stochastic volatility . . . . . . . . . . .                              |  15 |
|                  | 4.3                        | Simulation scheme for the micro and macro models                                   |  16 |
| 5                | Numerical investigations   | Numerical investigations                                                           |  17 |
|                  | 5.1 Calibration . . . .    | . . . . . . . . . . . . . . . . . .                                                |  18 |
|                  | 5.2                        | Pricing of some path dependent contracts . . . . .                                 |  18 |
|                  |                            | 5.2.1 Autocallable contracts . . . . . . . . . . . .                               |  19 |
|                  |                            | 5.2.2 Athena Jet on S&P GSCI Crude Oil Index Excess 5.2.3 .                        |  21 |
|                  |                            | Daily Knock In Option . . . . . . . . . .                                          |  22 |
|                  | 5.3                        | Sensitivities . . . . . . . . . . . . . . . . . . . . . .                          |  22 |
|                  |                            | . . . . . . . . . . . . .                                                          |  23 |
|                  |                            | 5.3.1 Delta . . . . . . . . 5.3.2 Vega . . . . . . . . . . . . . . . . . . . . . . |  24 |
|                  | 5.4                        | Computer implementation details . . . . . . . . . .                                |  25 |

6

Conclusions and further research

25

## 1 Introduction

Futures contracts are among the most liquid instruments in the commodities markets, making them the primary method for gaining exposure to underlying commodities. Since futures expire on specific dates, traders willing to maintain exposure to price changes in the underlying commodity must replace their expiring futures positions with new ones that have longer maturities. This process, known as a "rolling strategy" or "rolling hedge", involves shorting positions with the shortest maturities and taking long positions in futures with longer maturities. Specific implementations of rolling strategies can be found in indices like the S&amp;P GSCI indices . 1 Market participants seeking exposure in the commodities market can trade derivative contracts on these indices, thereby avoiding the complexities of directly trading futures, such as executing rolling strategies or managing margin requirements.

In the trading desks, the typical strategy to price derivative contracts on an index starts by making a model for the index. However, when building such models, the microscopic (micro) structure, i.e. the dynamics of the underlying components of the index, is often overlooked. This can be due to several reasons. From the numerical point of view, one reason comes from the vast increase in computational demands in order to perform a separate simulation for each of the individual components. From the modeling point of view, a first reason is motivated by the difficulty in specifying coherent relations between the components of the model. Another reason from the modeling point of view is that, although in our setting the underlying components and some derivatives contracts on them are very liquid, there are not enough liquid derivative contracts to calibrate the underlying components. The problem of not considering the micro structure is that, for certain products, the properties arising from the micro structure might have a direct impact on the prices that is not trivial to recover with a macroscopic (macro) model. Moreover, since in general the indices are not tradable assets, it is not clear how the macro models should be used for hedging purposes.

The main goal of this paper is building a model for futures prices that is able to calibrate plain-vanilla option prices on both commodity futures and indices and then compare its results with that of the more standard macro model on an index. Since contracts on indices are usually sensitive to smile effects and include path-dependency, the pricing model must be able to describe both curve and smile dynamics.

Local volatility (LV) models, as introduced by Dupire 1994; Derman et al. 1994, are well-known for their ability to perfectly match the market volatility surface for plain-vanilla options, thereby accurately reproducing the volatility smile. However, LV models can sometimes be unrealistic as they tend to flatten implied forward volatilities, a limitation highlighted in Rebonato 1999 and Hagan et al. 2002, making them potentially unsuitable for pricing contracts on indices. In contrast, Stochastic Volatility (SV) models, like the Heston model (see Heston 1993), where the variance of the asset follows its own stochastic differential equation (SDE), do not exhibit this flattening effect. Despite their advantages, SV models have a parametric form that prevents exact calibration to market prices of plain-vanilla options. In order to overcome these issues, Stochastic Local Volatility (SLV) models have been proposed,

aiming to combine exact calibration to plain-vanilla options with the stochastic nature of volatility. Emerging in the late 1990s and early 2000s (for example, see Said 1999; Lipton 2002; Ren et al. 2007), SLV models have since become a standard for pricing in many markets.

Examining our commodity settings more closely, we observe that accurately describing the dynamics of the futures term structure generally requires a distinct SV or LV model for the price of each futures contract, as noted in Pilz et al. 2011 and Chiminello 2015. While these models provide a precise depiction of futures term structure dynamics, they can be challenging to calibrate due to the typically limited number of quotes available in commodity markets. Consequently, we seek a more parsimonious approach. Following the methodology outlined by Nastasi et al. 2020 and Manzano et al. 2023, we employ a smaller family of SLV processes to govern the dynamics of the entire futures term structure.

In order to compare the micro and macro models, we will pay attention to the differences in prices and sensitivities of autocallable contracts, since they are path dependent. Autocallables are structured investment products which contain optionalities automatically triggered based on a specific event. These products are triggered on predetermined observation dates if an underlying asset or reference portfolio reaches or surpasses a barrier. More precisely, on the intermediate observation dates the autocallable pays a coupon if the underlying asset is above a certain coupon barrier level and it automatically redeems if it breaches an autocall barrier level. On the last observation date, the investor receives the principal amount, or a portion thereof, along with an optional payoff on the maturity date. The autocallability feature can be applied daily, monthly, yearly, or based on any schedule determined in accordance with the client. This feature is frequently offered in low-yield markets, thus providing the investor with the potential for an above-market yield, albeit with the risk of losing part of the initial capital. We refer the readers to Deng et al. 2011, Alm et al. 2013, Farkas et al. 2022 and references therein for a more detailed description of the autocallable product. As a result of the aforementioned structure, these products are clearly path-dependent. Moreover, they are actively traded in the market, what makes them an excellent choice to study the differences between the micro and macro models in a realistic setting.

The paper is organized as follows. In Section 2 we briefly describe the strategy to replicate a position on the index. Then, in Section 3 we illustrate the modeling framework for both the micro and macro models. Next, in Section 4 we explain the simulation and calibration procedures for both models. In Section 5 we calibrate the micro and macro models S&amp;P GSCI WTI Crude Oil index and investigate the differences between the micro and macro models by calculating the prices and sensitivities of different path-dependent contracts. Finally, we wrap up with the conclusions in Section 6.

## 2 S&amp;P GSCI indices

The following section describes S&amp;P indices and it is identical to Section 2 of our previous work Manzano et al. 2023.

S&amp;P GSCI indices are designed to replicate performances of actual commodity sectors. The idea behind the construction of the S&amp;P GSCI indices is to create an index that simulates a continuous investment on a basket of commodities (or a single commodity). The calculation of the S&amp;P GSCI indices takes into account the fact that a person holding positions in a contract near expiration would need to roll such positions forward as they approach settlement or

delivery dates. For this reason, the methodology for calculating S&amp;P GSCI indices includes a rolling procedure designed to replicate the rolling of actual futures positions. Since executing the rolling procedure on a single day could be difficult to implement or, if completed on a single day, could have an adverse impact on the market, such rolling takes place over a period of several days.

As the mechanics of the rolling depends particularly on which specific index we are dealing with, from now on we will focus on S&amp;P GSCI Excess-Return (ER) indices, which represent the daily return of a portfolio of commodity futures contracts. The rolling forward of the underlying futures contracts occurs once each month, on the fifth through ninth business day (the roll period) and the index is calculated as though these rolls occur at the end of each day during the roll period at the daily settlement prices. In the next subsections we explain the specific strategy of the non-rolling and rolling periods. As the strategy varies slightly for different indices and we will be working with the single-commodity index S&amp;P GSCI Crude Oil ER, we will describe how this index is calculated. This same procedure also covers how most of the single-commodity indices are built with a few exceptions (for example the S&amp;P GSCI Gold ER).

## 2.1 Evolution during a non-rolling period

On any business day during a non-rolling period the value of the S&amp;P GSCI Crude Oil ER index is equal to the product of the value of the index on the preceding business day times one plus the contract daily return on the business day on which the calculation is made.

For each futures maturity date T , i i = 1 , . . . , M F , we denote F t ( T i ) the futures price observed at time t ≤ T i . Moreover, let I t denote the value of the index at time t . The investment strategy implemented by the index consists in buying, at the beginning of each day t , a quantity Q t of futures contracts on the front month such that the nominal value of the investment is exactly I t . Therefore, the amount Q t of contracts we buy is given by

$$Q _ { t } \coloneqq \frac { I _ { t } } { F _ { t } ( T ^ { c } ) } \,,$$

where T c indicates the maturity of the front month. Due to the market movement of the futures price, at the end of the day our investment will have generated a profit-and-loss equal to

Such profit-and-loss is invested again in the strategy whose new value becomes

$$W _ { t + 1 } \coloneqq Q _ { t } \left ( F _ { t + 1 } ( T ^ { c } ) - F _ { t } ( T ^ { c } ) \right ) = I _ { t } \left ( \frac { F _ { t + 1 } ( T$$

$$I _ { t + 1 } = I _ { t } + W _ { t + 1 } = I _ { t } \, \frac { F _ { t + 1 } ( T ^ { c } ) } { F _ { t } ( T ^ { c } ) } \,.$$

If we repeat the same strategy the next day we get that the invested amount remains unchanged, so that

$$Q _ { t + 1 } = \frac { I _ { t + 1 } } { F _ { t + 1 } ( T ^ { c } ) } = \frac { I _ { t } } { F _ { t } ( T ^ { c } ) } \,. \quad \quad \quad$$

Proceeding recursively we have that the index value at n -th day after t can be calculated directly from the initial conditions as follows

$$I _ { t + n } = I _ { t } \, \frac { F _ { t + n } ( T ^ { c } ) } { F _ { t } ( T ^ { c } ) } \,.$$

## 2.2 Evolution during a rolling period

On a rolling period, we need to roll the nearest futures contract T c ≡ T 1 to the second nearest futures contract T f ≡ T 2 at a rate of 20% per day for the five days of the roll period. Notice that in other indices different from the S&amp;P GSCI Crude Oil ER index the rolling procedure might be done between contracts with other maturities (e.g. the GSCI Gold ER). We could also adapt our methodology to these indices in a straightforward way.

Until just before the end of the fifth business day, the entire S&amp;P GSCI Crude Oil ER index portfolio consists of the front futures contracts. At the end of the fifth business day, the portfolio is adjusted so that 20% of the held contracts are in the second futures contract, while 80% remain in the front one. The roll process continues on the sixth, seventh and eighth business days, with relative weights of front to second contracts of 60%/40%, 40%/60% and 20%/80%. At the end of the ninth business day, the last contract of the old front futures is exchanged, thus completing the roll and leaving the entire portfolio in what we have been calling the second futures contract. At this time, this former second futures becomes the new front contract and a new second futures is formed (with futures maturities further in the future) for use in the next month roll.

The last key point about the roll process is to specify exactly what the 80%/20% or other relative splits between front and second contracts mean. The roll percentages refer to contracts or quantities, not values. Taking the first day of the roll as an example, just before the roll takes place at the end of the day, the S&amp;P GSCI Crude Oil ER index consists of the front futures contract. That portfolio, constructed the night before and held throughout the fifth business day, has a dollar value. For the roll, that dollar value is distributed across the front and the second futures such that the number of contracts or the quantity of the front ones is 80% of the total and the quantity held of the second ones is 20% of the total.

We can illustrate the mechanics of the rolling period with an example to clarify the terms. Starting from the fifth business day of the month, due to the approaching expiry of the front month, the investment is gradually spread between the front and second futures. We call α t ( ) the investment percentage on the front futures contract F t ( T c ) and 1 -α t ( ) the investment percentage on the second futures contract F t ( T f ). Then, we define Q t within the rolling period as

$$Q _ { t } \coloneqq \frac { I _ { t } } { \alpha ( t ) F _ { t } ( T ^ { c } ) + ( 1 - \alpha ( t ) ) F _ { t } ( T ^ { f } ) } \,, \\ \text{$ \text{the numbers and unit} v of a fictitious contract made by the combination of}$$

which represents the purchased quantity of a fictitious contract made by the combination of the first and second futures. This investment generates the following profit-and-loss

$$\Pr _ { \ } L _ { t + 1 } \colon = I _ { t } \frac { \alpha ( t ) ( F _ { t + 1 } ( T ^ { c } ) - F _ { t } ( T ^ { c } ) ) + ( 1 - \alpha ( t ) ) ( F _ { t + 1 } ( T ^ { f } ) - F _ { t } ( T ^ { f } ) ) } { \alpha ( t ) F _ { t } ( T ^ { c } ) + ( 1 - \alpha ( t ) ) F _ { t } ($$

so that the new value of the strategy is given by

$$I _ { t + 1 } & = I _ { t } + \Pr \text{L} _ { t + 1 } = I _ { t } \frac { \alpha ( t ) F _ { t + 1 } ( T ^ { c } ) + ( 1 - \alpha ( t ) ) F _ { t + 1 } ( T ^ { f } ) } { \alpha ( t ) F _ { t } ( T ^ { c } ) + ( 1 - \alpha ( t ) ) F _ { t } ( T ^ { f } ) } \,. \quad \quad \quad$$

Contrary to the non-rolling period, if we now consider the evolution to the next day we obtain

that the quantities have changed due to the change in the investment weights

$$Q _ { t + 1 } & = \frac { I _ { t } } { \alpha ( t ) F _ { t } ( T ^ { c } ) + ( 1 - \alpha ( t ) ) F _ { t } ( T ^ { f } ) } \frac { \alpha ( t ) F _ { t + 1 } ( T ^ { c } ) + ( 1 - \alpha ( t ) ) F _ { t + 1 } ( T ^ { f } ) } { \alpha ( t + 1 ) F _ { t + 1 } ( T ^ { c } ) + ( 1 - \$$

so that the value of the index at time t + 2 cannot be evaluated only as a function of the initial conditions.

The description of the index strategy leads us to introduce a futures price model able to describe the complex path-dependent dynamics of rolling periods. In the next section we deal with this problem.

## 3 Modeling framework

The prices of a path depend options on a commodity index are expected to have a direct dependence not only on the paths followed by futures prices, but also on the correlation between futures. Hence, we are interested in the comparison of two different models. A macro model, i.e. , a model that tries to capture the dynamics of the index regardless of its underlying microstructure and a micro model, i.e. , a different model for each of the futures.

## 3.1 Micro model

In order to build the micro model we start by defining the dynamics of the futures prices. The index dynamics will be derived by implementing the index definition given in Section 2 in terms of the underlying futures. For each futures maturity date T , i i = 1 , . . . , M F , let F t ( T i ) be the futures price observed at time t ≤ T i . We describe the dynamics of futures prices under the risk-neutral measure as in Nastasi et al. 2020 and we consider deterministic interest rates. In line with the work in Manzano et al. 2023, we assume that each futures price follows a SLV model given by

$$d F _ { t } ( T _ { i } ) & = L ^ { F } ( t, T _ { i }, F _ { t } ( T _ { i } ) ) \sqrt { v _ { t } ^ { F, i } } \, d W _ { t } ^ { F, i }, \\ \mathfrak { m } _ { - } \text{ when } \, \mathfrak { m } \, f \, F / + \, T \, \mathfrak { m } \, \mathfrak { m } \, \quad.$$

$$\lambda \times \lambda \, \cdot \, d v _ { t } ^ { F, i } = \kappa ^ { F } ( \theta ^ { F } - v _ { t } ^ { F } ) + \chi ^ { F } \sqrt { v _ { t } ^ { v, i } } \, d W _ { t } ^ { v, i } \,, & & ( 1 1 ) \\ \intertext { e mean-reversion speed of the variance, } \theta ^ { F } \text{ is the long-term mean variance}$$

for i = 1 , . . . , M F , where { L F ( t, T , K i ) } i =1 ,...,M F are the so-called leverage functions and represent the local-volatility components of each SLV model. The common variance processes v F,i t satisfies the SDE

where κ F is the mean-reversion speed of the variance, θ F is the long-term mean variance and χ F is the volatility of the variance (also known as vol-of-vol). The initial value of the common variance v F,i 0 = v F 0 is not directly observable, so that we consider it as an additional parameter to be calibrated. For simplicity, all these parameters are assumed to be constant and the same across all futures. Moreover, { W F,i t } 1 ,...,M F and { W v,i t } 1 ,...,M F are standard Brownian motions under the risk-neutral measure with correlations

$$d \langle W ^ { F, i }, W ^ { F, j } \rangle _ { t } = \rho _ { i, j } ^ { F, F } \, d t \,, \quad d \langle W ^ { F, i }, W ^ { v, j } \rangle _ { t } = \delta _ { i, j } \rho ^ { F, v } \, d t$$

for i, j = 1 , . . . , M F . Notice that the variance processes v F,i t share the same correlation value between such process and the other futures Brownian motions since we aim at a parsimonious

description of futures dynamics when possible.

However, this model does not assume any specific form for the correlation matrix. Since the calibration of each entry in the correlation matrix would require many different contracts and a very complicated calibration strategy, we will assume that the correlation matrix follows the one parameter parametrization suggested in Rebonato 1999:

$$d \langle W ^ { F, i }, W ^ { F, j } \rangle _ { t } = \rho _ { i, j } ^ { F, F } \coloneqq e ^ { - \beta | T _ { i } - T _ { j } | }, \, i, j = 1, \dots, M ^ { F },$$

where i, j = 1 , . . . , M F and β &gt; 0 is called the de-correlation parameter. The assumption that the correlation matrix has the form of Equation (13) is somehow reasonable. On the one hand, because it always produces a valid correlation matrix in the sense that it produces a real, symmetric, positive-definite matrix. On the other hand, because in the market the correlation decreases for increasing maturity intervals and so does the proposed matrix. However, another feature of the correlation matrix that is typically observed in the markets is the fact that the correlation increases with large maturities and our model does not fulfill such restriction. In any case, for the shake of simplicity, we will restrict ourselves with correlation matrices of the form given in Equation (13).

## 3.2 Macro model

The dynamics of the index that we propose resembles the one of the micro-structure. Since the index is constructed as a combination of the different futures, we propose to model the index with the same dynamics of the underlying futures

$$d I _ { t } = L ^ { I } ( t, I _ { t } ) \sqrt { v _ { t } } I _ { t } \, d W ^ { I } _ { t },$$

where L I ( t, I t ) is the leverage function. The variance process v t satisfies the SDE:

$$d v _ { t } = \kappa ^ { I } ( \theta ^ { I } - v _ { t } ) + \chi ^ { I } \sqrt { v _ { t } } \, d W _ { t } ^ { v } \,,$$

Moreover, W I t is a standard Brownian motion under the risk-neutral measure and the instantaneous correlation between the W I t and W v t is given by

$$d \langle W ^ { I }, W ^ { v } \rangle _ { t } = \rho ^ { I, v } \, d t.$$

As we can see, the dynamics of both models are very similar, the most notable differences being the futures term and the correlation structure. In the case where the initial value of the futures is the same for all maturities and the correlation structure is just the identity matrix ρ F,F ij = δ i,j , the macro model collapses to the micro model. Note that, although they have the same structure, the parmaeters in the variance processes dv i t defined in Equation (15) are different from those in the variance process dv t defined in Equation (11).

## 4 Calibration and simulation

The next steps are the calibration and simulation of both models. The calibration will be different for the macro and micro models. In the micro model, we will calibrate the leverage functions to the prices of plain-vanilla options on the different futures and the correlation

structure plus the stochastic parameters on plain-vanilla options on the index. In the macro model, we will calibrate the leverage function and the stochastic parameters to the prices of plain-vanilla options on the index. In the next subsections we give more details on this calibration.

## 4.1 Calibration of the micro model

The calibration of the micro model is done in a iterative two-step process. On the one hand, the leverage function is calibrated to the prices of plain-vanilla options on the different futures by means of Markov projections via the application of the Gyöngy's Lemma. On the other hand, the correlation structure plus the stochastic parameters are calibrated to plain-vanilla options on the index using a global-local calibration strategy.

## 4.1.1 The leverage functions

If we look at the SLV literature, for instance in Guyon et al. 2012, we can find the description of a practical procedure to calibrate the leverage function by means of Markov projections via the application of the Gyöngy Lemma, see Gyöngy 1986. This Lemma states under which conditions the marginal densities of two semi-martingales are equivalent in law. Thus, if we are able to calibrate plain-vanilla options quoted by the market by means of a simpler model, e.g. a LV model, we can ensure that the SLV model does the same if such models satisfy the hypotheses of the Lemma.

In our case, the market quotes plain-vanilla options of futures prices, so that we introduce a LV model for each futures price, which is given by

$$d \hat { F } _ { t } ( T _ { i } ) = \hat { L } ^ { F } ( t, T _ { i }, \hat { F } _ { t } ( T _ { i } ) ) \, d \hat { W } _ { t } ^ { F, i } \,,$$

where the local-volatility functions { ˆ L F ( t, T , K i ) } i =1 ,...,M F are assumed to be Lipschitz, positive and at most of linear growth in price, so that there exists a solution for the SDE. Then, we can apply the Gyöngy's Lemma to perform the matching of the marginal distributions and we obtain

We can solve for the leverage functions to write the SLV model in terms of the local volatilities, thus getting

$$\text{own} \quad \hat { L } ^ { F } ( t, T _ { i }, K ) = L ^ { F } ( t, T _ { i }, K ) \, \sqrt { \mathbb { E } [ v _ { t } ^ { F, i } | F _ { t } ( T _ { i } ) = K ] } \,. \quad \text{(18)} \\ \text{live for the leverage functions to write the SLV model in terms of the local volatilities,}$$

$$d F _ { t } ( T _ { i } ) & = \hat { L } ^ { F } ( t, T _ { i }, F _ { t } ( T _ { i } ) ) \sqrt { \frac { v _ { t } ^ { F, i } } { \mathbb { E } [ v _ { t } ^ { F, i } | F _ { t } ( T _ { i } ) ] } } \, d W _ { t } ^ { i } \,. \\ 4. 1 \, \ W e n t o i c e \, \text{that the existence of solutions of } \text{McKean } S D E s \text{ such as the one in}$$

Remark 4.1. We notice that the existence of solutions of McKean SDEs such as the one in Equation (19) is a known open problem (see Guyon et al. 2012 for instance) and falls outside the scope of this paper.

The problem of calibrating the leverage functions has been transformed into the simpler problem of calibrating the local volatilities ˆ L F . This problem is considered and solved in a parsimonious way in Nastasi et al. 2020. Here, we adopt such solution, that we briefly describe in the following paragraphs.

First, we introduce a common driving factor ˆ s t for all future prices. We can understand it as a normalised 'spot' price. The dynamics of ˆ s t is given by

$$d \hat { s } _ { t } = a ( 1 - \hat { s } _ { t } ) \, d t + \hat { s } _ { t } \hat { L } ^ { s } ( t, \hat { s } _ { t } ) \, d \hat { W } _ { t } ^ { s } \,,$$

where a is the mean-reversion speed of the spot process, ˆ L s ( t, K ) is the local-volatility function of the spot price. This function is assumed to be Lipschitz, positive and bounded in price. Moreover, ˆ W ˆ s t is a standard Brownian motion under the risk-neutral measure. Then, we assume that futures prices in the LV model can be calculated starting from the normalised spot price as

$$& \Gamma \infty \infty \\ & \hat { F } _ { t } ( T _ { i } ) \coloneqq F _ { 0 } ( T _ { i } ) \, \mathbb { E$$

for i = 1 , . . . , M F , where F 0 ( T i ) is the term structure of futures prices observed in the market. Last equality in (21) can be derived by straightforward algebra due to the particular form of the spot price dynamics. Furthermore, we can apply the Itô Lemma and compare the result with equation (19), so that we get

$$\hat { L } ^ { F } ( t, T _ { i }, K ) \coloneqq \left ( K - F _ { 0 } ( T _ { i } ) \left ( 1 - e ^ { - a$$

where the effective strike k F can be defined as

$$k ^ { F } ( t, T _ { i }, K ) \coloneqq 1 - e ^ { a ( T _ { i } - t ) } \left ( 1 - \frac { K } { F _ { 0$$

Thus, the problem is now to calibrate the local volatility ˆ L s ( t, K ) of the spot price for a given mean-reversion speed a to the prices of plain-vanilla options on futures. Notice that the problem is now much simpler since instead of calibrating a different function for each futures we need to calibrate only one function. In order to do so, we write the price C F 0 ( T pv i , T , K i ) of plain-vanilla options with expiry T pv i on futures with maturity T i &gt; T pv i and strike K in terms of the spot price (in order to simplify notation, we assume that expiry and payment of the options are on the same date):

$$C _ { 0 } ^ { F } ( T _ { i } ^ { \text{pv} }, T _ { i }, K ) \coloneqq P _ { 0 } ( T _ { i } ^ { \$$

where P 0 ( T pv i ) is the zero-coupon bond with maturity T pv i and the normalized call prices c F ( t, k ) are defined as

$$c ^ { F } ( t, k ) \coloneqq \mathbb { E } _ { 0 } [ ( \hat { s } _ { t } - k ) ^ { + } ] \quad \quad \quad$$

and they satisfy the following extended version of the Dupire equation, see Nastasi et al. 2020:

$$\partial _ { t } c ^ { F } ( t, k ) = \left ( - a - a ( 1 - k ) \partial _ { k } + \frac { 1 } { 2 } k ^ { 2 } \hat { L } ^ { s } ( t, k ) ^ { 2 } \partial _ { k } ^ { 2 } \right ) c ^ { F } ( t, k ) \,.$$

In order to calibrate ˆ L s ( t, k ) we follow the same iterative strategy as in Nastasi et al. 2020. This strategy can be summarised in four main steps:

- 1. Solve Equation (26) for a fixed ˆ L s ( t, k ). We use a Cranck-Nicholson scheme in time and central finite differences in space (see Wilmott et al. 1995).
- 2. Compute the Black-Scholes volatilities ˆ σ F from the model prices c F .

- 3. Compare the model volatilities ˆ σ F with the market volatilities σ F and update the ˆ L s ( t, k ) accordingly.
- 4. Go back to step one until the target accuracy is met.

Once ˆ L s ( t, k ) is known we can calculate ˆ L F ( t, T , K i ) by means of equation (22) and then we get the leverage functions L F ( t, T , K i ) by solving equation (18).

## 4.1.2 Stochastic parameters and correlation structure

In Section 4.1.1 we already described how to calibrate the leverage functions L F ( t, T , K i ) to recover the prices of plain vanillas of the futures. The remaining degrees of freedom of the micro model:

$$p ^ { F } & = \{ a, \beta, \kappa ^ { F }, \theta ^ { F }, \chi ^ { F }, \rho ^ { F, v }, v _ { 0 } ^ { F } \} \,, \\. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &.. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &.$$

are adjusted to recover the implied volatilities of plain vanillas on the on GSCI ER index σ I market . For this purpose, we define the cost function:

$$D _ { \ell ^ { 1 } } ( \boldsymbol p ^ { F } ) = D _ { \ell ^ { 1 } } ( \boldsymbol \sigma _ { \text{market} } ^ { I }, \boldsymbol \sigma _ { \text{micro} } ^ { I } ( \boldsymbol p ^ { F } ) ) = \sum _ { j = 1 } ^ { J } \left | \sigma _ { \text{market} } ^ { I, j } - \sigma _ { \text{micro} } ^ { I, j } ( \boldsymbol p ^ { F } ) \right | \,, \quad \quad ( 2 8 ) \\ \text{where $\boldsymbol \sigma_{\text{micro} }^{\boldsymbol I}_{\text{micro} }$ are the implied volatilities given by the micro model when evaluated with pa-}.$$

where σ I micro are the implied volatilities given by the micro model when evaluated with parameters p F and J is the total number of vanillas on the index that have been considered.

In this setup, the calibration of the free parameters to fit the prices of plain-vanillas on the index can be formulated as the following unconstrained global optimization problem in a bounded domain:

$$\min _ { \mathbf p ^ { F } \in P \subseteq \mathbb { R } ^ { 7 } } D _ { \ell ^ { 1 } } ( \mathbf p ^ { F } )$$

where D d is the cost function defined on P F = Π 7 r =1 [ l F r , u F r ], with l F r and u F r being the lower and upper bounds in direction i , respectively. The solution vector p F ∗ contains the calibrated parameters and is defined as:

$$\mathbf p _ { * } ^ { F } = \underset { \mathbf p ^ { F } \in P ^ { F } \subseteq \mathbb { R } ^ { 7 } } { \arg \min } D _ { \ell ^ { 1 } } ( \mathbf p ^ { F } ),$$

Note that each evaluation of the cost function requires the numerical solution of the micro model.

In order to solve the unconstrained global optimization problem (29), we propose a simplified version of a two-phase calibration strategy (see Ferreiro-Ferreiro et al. 2020; Ferreiro et al. 2019 for a complete description). In particular, we start by defining an initial guess p F 0 of the optimal solution p F ∗ . We typically choose p F 0 randomly. Next, we run a global optimization algorithm starting with p F 0 . Once the global algorithm has finished, it provides an intermediate solution p F 1 . Then, we use p F 1 as the initial point for a local optimization algorithm. Finally, the local optimization algorithm gives us p F 2 . These steps are summarized in Algorithm 1.

## Algorithm 1 Overall optimization algorithm

| Input:                       |                                |
|------------------------------|--------------------------------|
| p F 0                        | // random seed for p F         |
| D d                          | // function to be minimized    |
| Output:                      |                                |
| p F 2                        | // approximated value of p F ∗ |
| Algorithm:                   |                                |
| ˇ p s 1 = ESCH( ˇ p s 0 )    | // global minimization.        |
| ˇ p s 2 = Subplex( ˇ p s 1 ) | // local minimization.         |

As the global optimization algorithm, we use the so-called ESCH (see Silva Santos et al. 2010). ESCH algorithm belongs to the broad category of Evolutionary Algorithms (EAs) (see Vikhar 2016) which is a class of heuristics inspired by natural selection in biological populations. We specifically use ESCH algorithm because it is available in the NLopt nonlinearoptimization package (see Johnson 2007 for more details). In Algorithm 2 we show the pseudocode for the ESCH method.

This evolutionary algorithm starts by dividing the population into two groups: parents and offspring. First, the parameters of the parents are randomly initialized from a uniform distribution. Then, some of the parents are selected and their information recombined to generate offspring. Some of the offspring would mutate one of their parameters. Next, we assign a score to the offspring based on their value in the function that we want to minimize. Finally, all individuals are ranked according to their fit score, the less fit individuals are removed from the population and only the best individuals are stored along the generations. The adopted recombination and the mutation operators are the so-called single point and the Cauchy distribution, respectively. The scheme for the evolution can be found in Algorithm 2. For a more detailed description we refer to the original article Silva Santos et al. 2010.

The advantage of starting by a global optimization algorithm is that it is able to escape from local minima. However, in general these algorithms in are computationally very expensive as they have a slow rate of convergence. Our purpose starting with a global optimization algorithm is to explore the space and end up with a good initial point for the local optimization algorithm. In the best case scenario the output of the algorithm should lay in the convex region defined by the global minima. Other possible choices for the global optimization routine are: Simulated Annealing (SA, see Vikhar 2016, Aarts et al. 1985), Differential Evolution (DE, see Storn et al. 1997) or Particle Swarm (PS, see Kennedy et al. 1995).

As the local optimization algorithm, we use a variation of the Nelder-Mead algorithm (see Nelder et al. 1965), called 'Subplex' developed in Rowan 1990. The Subplex algorithm starts by dividing the search space into subspaces. Then, the Nelder-Mead algorithm is used to minimize in each subspace. The subspaces in which the minimization has been larger are joined to form a new subspace. The process continues iteratively until the stopping criteria is met. An outline of the Subplex basic steps can be found in Algorithm 3. For a more detailed description we refer to the original thesis Rowan 1990.

Note that local search algorithms are not able to escape from local minima, although they

## Algorithm 2 ESCH pseudocode

Input:

R

// problem dimension

D d

// function to be minimized

x 0

// initial approximation to minimum

np

// number of parents

no

// number of offspring

Output:

x

// computed minimum

## Algorithm:

Initialize parents and offspring population

Parents fitness evaluation

while termination test not satisfied do

Crossover

Gaussian Mutation

Offspring fitness evaluation

Selection of the fittest as parents

end while

exhibit faster convergence rates than the global optimization ones. Therefore, the purpose of using the local optimization algorithm is to perform a fine grain optimization much faster than it would be possible with a global optimization method.

## 4.2 Calibration of the macro model

The calibration of the macro model will be different from the calibration of the micro model, since both the leverage function and the stochastic parameters need to be calibrated on the same surface of plain vanillas on the index. In a trading desk, these two are typically calibrated separately. First they build a local volatility model and calibrate it to the surface of plain vanillas on the index. Second, they calibrate a pure stochastic model on the volatility smile for a fixed maturity in the same surface. Finally, they calibrate the leverage function by means of Markov projections via the application of the Gyöngy's Lemma.

In our case there are two main differences. First, the calibration of the stochastic parameters is done on the whole surface of plain vanillas on the index and not on a fixed maturity. Second, the calibration is not performed on the market surface but on the surface produced by the micro model. In this way, we ensure that the marginal probabilities of both models are the same.

## 4.2.1 The leverage function

The procedure to calibrate the leverage function for the index model follows the same steps as in Section 4.1.1. First, we define the local-volatility model:

$$d \hat { I } _ { t } = \hat { L } ^ { I } ( t, \hat { I } _ { t } ) \hat { I } _ { t } \, d W _ { t } ^ { \hat { I } } \,.$$

| Algorithm 3 Subplex               | pseudocode                                                                                                                                                                                                                                                                                                                                    |
|-----------------------------------|-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Input:                            |                                                                                                                                                                                                                                                                                                                                               |
| R D d x 0 Scale α β γ δ ψ Ω nsmin | // problem dimension // function to be minimized // initial approximation to minimum // initial step size for the r coordinate directions // reflection coefficient // contraction coefficient // expansion coefficient // shrinkage coefficient // simplex reduction coefficient // step reduction coefficient // minimum subspace dimension |

where the local-volatility ˆ L I ( t, K ) is assumed to be Lipschitz, positive and bounded, so that there exists a solution for the SDE. Next, applying Gyöngy's Lemma we match the marginal distributions of the local-volatility model defined in Equation (31) with the SLV model defined in Equation (14):

$$L _ { I } ( t, K ) = \frac { \hat { L } ^ { I } ( t, K ) } { \sqrt { \mathbb { E } [ v _ { t } | I _ { t } = K ] } } \,. & & ( 3 2 ) \\ \text{Remark 4.1} \, w e \text{ stress that the existence of solutions of } \text{McKean SDEs} \\ \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Psi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi\Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phis \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Clone$$

Remark 4.2. As in Remark 4.1 we stress that the existence of solutions of McKean SDEs such as the one in Equation (32) is a known open problem (see Guyon et al. 2012 for instance) and falls outside the scope of this paper.

Now, we can solve for the leverage function to write the SLV model in terms of the local volatilities, thus getting:

$$d I _ { t } & = \hat { L } ^ { I } ( t, I _ { t } ) I _ { t } \sqrt { \frac { v _ { t } } { \mathbb { E } [ v _ { t } | I _ { t } ] } } d W _ { t } ^ { I } \,. \\ \intertext { a t i n d } \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \end{array}$$

With these manipulations, the problem of calibrating the leverage has been transformed in that of calibrating the local-volatility function η I ( t, K ). This local-volatility function satisfies a Dupire equation:

$$\partial _ { t } c ^ { I } ( t, k ) = \frac { 1 } { 2 } k ^ { 2 } \hat { L } ^ { I } ( t, k ) ^ { 2 } \partial _ { k k } c ^ { I } ( t, k ) \,,$$

where c I are normalised call prices defined by:

$$c ^ { I } ( t, k ) \coloneqq \frac { 1 } { I _ { 0 } } \mathbb { E } _ { t } [ ( I _ { t } - K ) ^ { + } ] \,,$$

In order to calibrate ˆ L I ( t, k ) we follow an analogous strategy as in Section 4.1.1. This strategy can be again summarised in four main steps:

- 1. Solve Equation (34) for a fixed ˆ L I ( t, k ). We use a Cranck-Nicholson scheme in time and central finite differences in space Wilmott et al. 1995.
- 2. Compute the Black-Scholes volatilities ˆ σ I macro from the model prices c I .
- 3. Compare the model volatilities ˆ σ I macro with the micro volatilities σ I micro and update the local-volatility function ˆ L I ( t, k ) accordingly.
- 4. Go back to step one until the target accuracy is met.

## 4.2.2 The stochastic volatility

In the micro model we calibrated the leverage and the stochastic parameters in a different set of options. However, since usually we don't have access to additional liquid options, for the macro model the stochastic parameters are typically calibrated on the same set of options as the leverage. For this purpose, we define the pure stochastic volatility model:

$$d \overline { I } _ { t } = \overline { I } _ { t } \sqrt { \overline { v } } _ { t } d \overline { W } _ { t } ^ { I } \,,$$

where the variance process v t satisfies the SDE:

$$d \overline { v } _ { t } = \overline { \kappa } ^ { I } ( \overline { \theta } ^ { I } - \overline { v } _ { t } ) + \overline { \chi } ^ { I } \sqrt { \overline { v } _ { t } } \, d \overline { W } _ { t } ^ { v, I } \,,$$

and the instantaneous correlation between the W I t and W v,I t is given by

$$d \langle \overline { W } ^ { I }, \overline { W } ^ { I, v } \rangle _ { t } = \overline { \rho } ^ { I, v } \, d t \,.$$

This model has five degrees of freedom:

$$\overline { \mathbf p } ^ { I } \coloneqq \{ \overline { \kappa } ^ { I }, \overline { \theta } ^ { I }, \overline { \chi } ^ { I }, \overline { \mathbf p } ^ { I, v }, \overline { v } ^ { I } _ { 0 } \} \,.$$

that will be calibrated on the implied volatilities of plain vanillas on the index given by the micro model σ I micro . For this purpose, we define the cost function:

$$D _ { \ell ^ { 1 } } ( \boldsymbol p ^ { I } ) = D _ { \ell ^ { 1 } } ( \boldsymbol \sigma ^ { I } _ { \text{micro} }, \overline { \boldsymbol \sigma } ^ { I } _ { \text{macro} } ( \overline { \boldsymbol p } ^ { I } ) ) = \sum _ { j = 1 } ^ { J } \left | \sigma ^ { I, j } _ { \text{micro} } - \overline { \boldsymbol \sigma } ^ { I, j } _ { \text{macro} } ( \overline { \boldsymbol p } ^ { I } ) \right | \,, \quad \quad ( 4 0 ) \\ \therefore \, \overline { \boldsymbol \sigma } ^ { I } _ { \text{macro} } \, \text{are the implied volatilities given by the pure stochastic macro model and } J \, \text{is}$$

where σ I macro are the implied volatilities given by the pure stochastic macro model and J is the total number of vanillas on the index that have been considered.

This setup closely resembles the one in Section 4.1.2 and the problem can be formulated as the following unconstrained global optimization problem in a bounded domain:

$$\min _ { \overline { p } ^ { I } \in \overline { P } ^ { I } \subseteq \mathbb { R } ^ { 5 } } D _ { \ell ^ { 1 } } ( \overline { p } ^ { I } ) \,,$$

where D d is the cost function defined on P I = Π 5 r =1 [ l I r , u I r ], with l I r and u I r being the lower and upper bounds in direction r , respectively. The solution vector p I ∗ contains the calibrated parameters and is defined as:

$$\Phi _ { ^ { * } } ^ { I } = \underset { \overline { p } ^ { I } \in \overline { P } \subseteq \mathbb { R } ^ { 5 } } { \arg \min } D _ { \ell ^ { 1 } } ( \overline { p } ^ { I } ),$$

Once again, in order to solve this unconstrained global optimization problem, we replicate the same hybrid global-local calibration as in Section 4.1.2, i.e., we start with a random seed p I 0 . Next, we first apply the ESCH global minimization algorithm to obtain an intermediate result p I 1 . Finally, we apply the Sublplex local minimization algorithm obtaining the solution p I 2 .

With the solution of the pure stochastic volatility model p I , we calibrate the macro model by simply fixing the values the stochastic parameters of I t with the values of I t :

$$\kappa ^ { I } = \overline { \kappa } ^ { I } \,, \theta ^ { I } = \overline { \theta } ^ { I } \,, \chi ^ { I } = \overline { \chi } ^ { I } \,, \rho ^ { I, v } = \overline { \rho } ^ { I, v } \,, v _ { 0 } ^ { I } = \overline { v } _ { 0 } ^ { I }.$$

## 4.3 Simulation scheme for the micro and macro models

For the simulation of both the micro and the macro model we resort to a Monte Carlo simulation since it is flexible, simple and efficient. We will focus on three main parts. First, the simulation of the conditional expectation in the diffusive term for both the micro and the macro model. Second the discretization of the variance process. Third, the simulation of the futures and index processes.

The simulation of the conditional expectation in the diffusive term for both the micro and the macro is achieved by means of the particle method described in Guyon et al. 2012.

In the case where we have N sample paths { ( F t ( T i ) n , v i,n t ) } n =1 ,...,N for the i -th future in the case where we have N sample paths { ( I n t , v n t ) } n =1 ,...,N for the macro model, the conditional expectations are approximated as:

and

$$\frac { 1 } { \sqrt { \mathbb { E } [ v _ { t } ^ { i } | F _ { t } ( T _ { i } ) ] } } = \sqrt { \frac { \sum _ { j = 1 } ^ { N } \delta ^ { \epsilon } ( F _ { t } ( T _ { i } ) ^ { n } - F _ { t } ( T _ { i } ) ^ { j } ) } { \sum _ { j = 1 } ^ { N } v _ { t } ^ { i, j } \delta ^ { \epsilon } ( F _ { t } ( T _ { i } ) ^ { n } - F _ { t } ( T _ { i } ) ^ { j } ) } } \,,$$

$$\frac { 1 } { \sqrt { \mathbb { E } [ v _ { t } | I _ { t } ] } } & = \sqrt { \frac { \sum _ { j = 1 } ^ { N } \delta ( I _ { t } ^ { n } - I _ { t } ^ { j } ) } { \sum _ { j = 1 } ^ { N } v _ { t } ^ { j } \delta ( I _ { t } ^ { n } - I _ { t } ^ { j } ) } }. \\ \intertext { n } \delta ^ { \epsilon } ( x ) \text{ is any mollifier of the Dirac delta.}$$

where the function δ /epsilon1 ( x ) is any mollifier of the Dirac delta.

The simulation of the variance processes in both models is achieved through the fulltruncated scheme as described in Lord et al. 2010. When we have a time grid { t m m } =1 ,...,M the full-truncated scheme applied to the micro model yields:

$$v _ { t _ { m + 1 } } ^ { i, n } = v _ { t _ { m } } ^ { i, n } + \kappa ^ { F } ( \theta ^ { F } - ( v _ { t _ { m } } ^ { i, n } ) ^ { + } ) \, \Delta t _ { m } + \chi ^ { F } \sqrt { ( v _ { t _ { m } } ^ { i, n } ) ^ { + } } \, \Delta W _ { t _ { m } } ^ { v, i, n } \,, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \,
d t o n t h o m o r o m o d o l.$$

and to the macro model:

$$v _ { t _ { m + 1 } } = v _ { t _ { m } } + \kappa ^ { I } ( \theta ^ { I } - ( v _ { t _ { m } } ) ^ { + } ) \, \Delta t _ { m } + \chi ^ { I } \sqrt { ( v _ { t _ { m } } ) ^ { + } } \, \Delta W ^ { v } _ { t _ { m } } \,. \quad \, \\ \text{full-truncated scheme avoids the negative values of the variance across at the discrete}$$

The full-truncated scheme avoids the negative values of the variance process at the discrete level while maintaining a good convergence to the continuous process.

Finally, for the discretization of the futures and the index we follow a standard EulerMaruyama scheme. For the micro model yields

$$& F _ { t _ { m + 1 } } ( T _ { i } ) ^ { n } = F _ { t _ { m } } ( T _ { i } ) ^ { n } + L ^ { F } ( t _ { m }, T _ { i }, F _ { t _ { m } } ( T _ { i } ) ^ { n } ) \, \sqrt { \frac { ( v _ { t _ { m } } ^ { i, n } ) ^ { + } \sum _ { j = 1 } ^ { N } \delta ^ { c } ( F _ { t _ { m } } ( T _ { i } ) ^ { n } - F _ { t _ { m } } ( T _ { i } ) ^ { j } ) } { \sum _ { j = 1 } ^ { N } ( v _ { t _ { m } } ^ { i, j } ) ^ { + } \, \delta ^ { c } ( F _ { t _ { m } } ( T _ { i } ) ^ { n } - F _ { t _ { m } } ( T _ { i } ) ^ { j } ) } } \, \Delta W _ { t _ { m } } ^ { i, n }, \\ & \text{and for the macro model}$$

and for the macro model

$$I _ { t _ { m + 1 } } ^ { n } = I _ { t _ { m } } ^ { n } + I _ { t _ { m } } ^ { n } \sqrt { v _ { t _ { m } } ^ { n } \frac { \sum _ { j = 1 } ^ { N } \delta ( I _ { t _ { m } } ^ { n } - I _ { t _ { m } } ^ { j } ) } { \sum _ { j = 1 } ^ { N } v _ { t _ { m } } ^ { j } \delta ( I _ { t _ { m } } ^ { n } - I _ { t _ { m } } ^ { j } ) } } d W _ { t _ { m } } ^ { I } \,. \\ \text{al volatility functions are obtained by making a piecewise-constant interpolation in}$$

The local volatility functions are obtained by making a piecewise-constant interpolation in time and a piecewise linear interpolation in price.

## 5 Numerical investigations

In this section we investigate the impact that the use of the two different models has on the prices and sensitivities of path-dependent options on an index. More precisely, we focus on derivatives contracts on the S&amp;P GSCI Crude Oil Index Excess Return composed by WTI Crude Oil futures.

## 5.1 Calibration

As explained in Section 4 the calibration of the micro model is performed against plain vanillas on the futures and the index quoted the 16th of December of 2019. Since not all parameters are relevant for the calibration we keep some of them fixed. For the micro model we fix the parameters in Table 1.

Table 1: Fixed values of p F for the calibration of the micro model.

| p F   | Value   |
|-------|---------|
| κ F   | 1 . 0   |
| θ F   | 1 . 0   |
| v F 0 | 1 . 0   |

The reason to fix these parameters is that the leverage function overruns their effect. Moreover, we choose κ F and θ F so that the Feller condition is not broken, thus avoiding numerical issues. From a more practical perspective, the impact of such parameters is already tested in Manzano et al. 2023.

For the macro model things are different. Although one could follow the same arguments and fix the parameters controlling the level, we have to recall that we are doing the calibration of the stochastic parameters by calibrating a pure stochastic model. For this reason, all parameters have a potential impact. However, we will also fix the value of κ I to avoid breaking the Feller condition, as indicated in Table 2.

Table 2: Fixed values of p I for the calibration of the macro model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00784v1-figure-9434da89-2886-4440-8c4a-53fa4c6f2ae1.jpg)

| p I   | Value   |
|-------|---------|
| κ I   | 1 . 0   |

After the previous considerations, the calibration of both models yields the results in Table

3.

| p F   | Seed   | Calibrated Value   |
|-------|--------|--------------------|
| a     | 0 . 3  | 0 . 338619         |
| β     | 0 . 1  | 0 . 172338         |
| χ     | 0 . 03 | 1 . 4              |
| ρ F,v | 0 . 2  | 0 . 40985          |

-

| p I     | Seed   | Calibrated Value   |
|---------|--------|--------------------|
| θ I     | 0 . 09 | 0 . 069918         |
| χ I     | 0 . 03 | 0 . 01277          |
| v I 0 I | 0 . 09 | 0 . 0637628        |
| ρ I,v   | 0 . 2  | 1 . 0              |

-

Table 3: Calibrated values of p F and p I against the mean of the plain vanillas on the S&amp;P GSCI Crude Oil ER index quoted on the 30 November 2019 and on 31 December 2019. For this calibration the hybrid global-local procedure described in Section 4 has been used.

## 5.2 Pricing of some path dependent contracts

So far we have performed a calibration of both the micro and the macro model. Now, we want to compare the results given by both models when pricing path-dependent options. However, it is difficult to directly compare the difference in prices between the micro and the

macro models. The standard procedure among traders to determine whether the difference is relevant consists in comparing the change in prices among models with the difference in prices that a change of market volatility of 1% would produce in the target contract. More specifically, for the micro model this means shifting the plain vanilla volatilities on the futures σ F micro by 1% and recomputing the leverage of the micro model accordingly:

$$\Delta V _ { \text{micro} } = V _ { \text{micro} } \left ( \sigma _ { \text{market} } ^ { F } + 1 \% \right ) - V _ { \text{micro} } ( \sigma _ { \text{market} } ^ { F } ) \. \quad \quad ( 5 0 ) \\ e \, \text{ macro model we proceed in a different way } \ \text{ Since we calibrated the macro model}$$

For the macro model we proceed in a different way. Since we calibrated the macro model to match the marginal probabilities of the micro one, instead of shifting the market volatilities σ I market by one percent, one should shift the micro volatilities σ I micro by one percent. Nevertheless, to keep things more consistent, we recalibrate the macro model on the implied volatilities of the index produced by the shifted micro model σ I micro ( σ F market +1%). In this way, both the micro and macro model retain the same marginal probabilities while, at the same time, a shift of 1% in the implied volatilities on the futures approximately produces a shift of 1% in the implied volatilities on the index:

$$\Delta V _ { \text{macro} } = V _ { \text{macro} } ( \sigma ^ { I } _ { \text{micro} } ( \sigma ^ { F } _ { \text{market} } + 1 \% ) ) - V _ { \text{macro} } ( \sigma ^ { I } _ { \text{micro} } ( \sigma ^ { F } _ { \text{market} } ) ) \. \quad \ \ ( 5 1 )$$

As a rule of thumbs, the traders consider that, when the model difference is 0 5 times greater . than the difference produced by a change of 1% in the surface of implied volatilities we can say that the difference between models is relevant, i.e.,

$$\frac { V _ { \text{micro} } - V _ { \text{macro} } } { \Delta V } \geq 0. 5 \longrightarrow V _ { \text{micro} } \neq V _ { \text{macro} }.$$

/negationslash

We divided the experiments in three sections. In Section 5.2.1 we price autocallable products with three different coupon structures. Next, in Section 5.2.2 we price a product that is actively traded in the market: the Athena Jet On S&amp;P GSCI Crude Oil Index Excess Return (see Paribas 2015). Finally, in Section 5.2.3 we price a daily barrier knock-in.

All experiments are performed with a confidence level of 95%.

## 5.2.1 Autocallable contracts

An autocallable is usually structured as a note paying coupons { γ , ..., γ 1 M A } on the observation dates { T , ..., T 1 M A } . The contract can be terminated before its expiry on T M A if a market event occurs. In case of early termination on T i a rebate β i is paid along with the capital redemption φ i . On the last date the whole capital or a part of it φ M A is returned along with an optional payoff β M A . If we consider for ease of exposition only autocallables on a single underlying asset, we can describe the market event as the first time the asset price S t is above a predefined level H i . Thus, the price of an autocallable is given by

$$V _ { t } \coloneqq \mathbb { E } _ { t } \left [ \sum _ { i = 1 } ^ { M _ { A } } \gamma _ { i } J _ { i } D ( t, T _ { i } ) + J _ { i - 1 } ( 1 - J _ { i } ) ( \beta _ { i } + \phi _ { i } ) D ( t, T _ { i } ) \right ], \\ \text{re the survival indicator } J _ { i } \text{ is defined for } i \in \{ 1. \dots, M _ { A } \} \text{ as}$$

where the survival indicator J i is defined for i ∈ { 1 , ..., M A } as

$$J _ { i } = 1 _ { \{ i < M _ { A } \} } \min _ { j \in \{ 1, \dots, i \} } 1 _ { \{ S _ { T _ { j } } < H _ { j } \} },$$

where we set J 0 := 1 since we assume that at inception the product is alive, while H j is the barrier level at time T j . We notice that on the last payment date the autocallable must be terminated, so that on the last date the indicator is always equal to zero. The rebate payments are usually equal to the corresponding coupon, namely

$$\beta _ { i } & = \gamma _ { i }.$$

On the other hand, the coupon payments γ i may be defined in different ways. Here, we briefly describe the most common ones:

- · Bullet coupons: In this case, the same quantity of money is paid on each date T i , namely

$$\gamma _ { i } ^ { \text{bullet} } \coloneqq \overline { \gamma }.$$

- · Digital coupons: In this case, a digital option is paid on each date T i , namely

$$\gamma _ { i } ^ { \text{digital} } \coloneqq \overline { \gamma } \mathbb { 1 } _ { \{ S _ { T _ { i } } > K _ { i } \} },$$

where K i is the digital option strike for the coupon paying at date T i .

- · Snowball coupons: In this case, a digital option with strike K is paid on each date T i . Moreover, if the option is triggered the amount paid is proportional to the number of coupons not paid since the last payment, namely

$$\gamma _ { i } ^ { \text{snowball} } \coloneqq N _ { i } \overline { \gamma } \mathbb { 1 } _ { \{ S _ { T _ { i } } > K _ { i } \} },$$

where the proportionality constant is defined as

$$N _ { i } \coloneqq i - \sup \{ 0, j \in \{ 1, \dots, i - 1 \} \, \colon S _ { T _ { j } } > K _ { j } \}.$$

The capital redemptions are only partial to allow to offer higher coupons. Indeed, they are usually defined as

$$\phi _ { i } \coloneqq 1 _ { \{ i < M _ { A } \} } + 1 _ { \{ i = M _ { A } \} } \left ( 1 _ { \{ S _ { T _ { M _ { A } } } \geq H _ { M _ { A } } \} } + 1 _ { \{ S _ { T _ { M _ { A } } } < H _ { M _ { A } } \} } \frac { S _ { T _ { M _ { A } } } } { H _ { M _ { A } } } \right ),$$

so that the capital is wholly redeemed if an early termination is triggered or if the contract is above the prespecified level H M A on the last date and it is partially redeemed if it is below H M A on the last date.

For the experiments we specifically choose an autocallable that expires nine months after the reference date and can be executed monthly. The barriers and strikes on each callable date are shown in Table 4. The barriers have been chosen in such a way that the probability of the contract being executed at any callabillity date is spread throughout the whole callability period. In this way, we try to emphasize the path-dependent properties of the contract, hence the differences among models.

Table 4: Barriers for the autocallable contracts.

|            |   H i |   K i |
|------------|-------|-------|
| 2020-01-17 | 1.1   | 1     |
| 2020-02-17 | 1.1   | 1     |
| 2020-03-17 | 1.075 | 0.975 |
| 2020-04-17 | 1.075 | 0.95  |
| 2020-05-19 | 1.075 | 0.925 |
| 2020-06-18 | 1.025 | 0.875 |
| 2020-07-17 | 0.95  | 0.775 |
| 2020-08-17 | 0.85  | 0.675 |
| 2020-09-16 | 0.7   | 0.5   |

With this structure we compute the prices of the product with the three aforementioned coupons for the micro and macro models. The prices and the normalised differences in prices for the micro and macro models are shown in Table 5. As we can see from the table, both models yield very similar results in terms of prices. Furthermore, the rule of thumbs from Equation (5.2) indicates that we cannot consider that both models are significantly different.

Table 5: Prices of autocallable contracts with bullet, snowball and digital coupons for the micro and macro models.

|                             | Bullet                    | Snowball                   | Digital                   |
|-----------------------------|---------------------------|----------------------------|---------------------------|
| Coupon                      | 0.005                     | 0.005                      | 0.005                     |
| V macro                     | 1 . 00996 ± 0 . 00018     | 1 . 00931 ± 0 . 00019      | 0 . 997769 ± 0 . 00020    |
| V micro                     | 1 . 00988 ± 0 . 00019     | 1 . 0091 ± 0 . 00021       | 0 . 99769 ± 0 . 00021     |
| ∆ V macro                   | - 0 . 002183 ± 0 . 000021 | - 0 . 0023147 ± 0 . 000023 | - 0 . 002329 ± 0 . 000022 |
| ∆ V micro                   | - 0 . 002246 ± 0 . 000024 | - 0 . 002399 ± 0 . 000027  | - 0 . 002495 ± 0 . 000025 |
| V micro - V macro ∆ V macro | 0 . 034 ± 0 . 121         | 0 . 099 ± 0 . 122          | 0 . 0020 ± 0 . 1244       |
| V micro - V macro ∆ V micro | 0 . 033 ± 0 . 117         | 0 . 095 ± 0 . 117          | 0 . 001915 ± 0 . 1210     |

## 5.2.2 Athena Jet on S&amp;P GSCI Crude Oil Index Excess Return Certificate

In this section we consider a real contract, the Athena Jet on S&amp;P GSCI Crude Oil Index Excess Return Certificate (see Paribas 2015). This contract has a maximum duration of one year but can expire after the first six months if the price of the S&amp;P GSCI Crude Oil Index Excess Return index is higher than or equal to its initial value. In that case it returns the nominal value of 1 euro and pays a 5% premium. Otherwise, the investor receives a premium of 2 5%. . At maturity (1 year) the certificate returns 1 euro and pays a premium equal to 1.5 times the performance of the underlying if the S&amp;P GSCI Crude Oil Index Excess Return index quotes at a value greater than or equal to its initial value. If the S&amp;P GSCI Crude Oil Index Excess Return index quotes between its initial value and a barrier, set at 70% of the initial value, the certificate protects the invested capital by returning 1 euro. If instead the S&amp;P GSCI Crude Oil Index Excess Return quotes below the Barrier, the investor receives an amount commensurate with the performance of the underlying index (with consequent loss

## Athena Jet on S&amp;P Certificate

Table 6: Prices of the Athena Jet Su on the S&amp;P GSCI Crude Oil Index Excess Return for the micro and macro models.

| V macro                     | 1 . 00207 ± 0 . 00048     |
|-----------------------------|---------------------------|
| V micro                     | 1 . 00348 ± 0 . 00047     |
| ∆ V macro                   | - 0 . 003718 ± 0 . 000091 |
| ∆ V micro                   | - 0 . 003529 ± 0 . 000047 |
| V micro - V macro ∆ V macro | - 0 . 38 ± 0 . 18         |
| V micro - V macro ∆ V micro | - 0 . 40 ± 0 . 19         |

Table 7: Prices of a daily knock-in option for the micro and macro models.

|                             | One year daily knock in   |
|-----------------------------|---------------------------|
| V macro                     | 0 . 03656 ± 0 . 00026     |
| V micro                     | 0 . 03828 ± 0 . 00029     |
| ∆ V macro                   | 0 . 003305 ± 0 . 000021   |
| ∆ V micro                   | 0 . 003425 ± 0 . 000027   |
| V micro - V macro ∆ V macro | 0 . 52 ± 0 . 12           |
| V micro - V macro ∆ V micro | 0 . 50 ± 0 . 11           |

on invested capital).

The results for this contract are shown in Table 6. As in the previous section, we observe that both models produce similar results in terms of prices and normalised differences. However, for this contract the differences become slightly more significant than in Table 5 to the point where we cannot completely discard the possibility that both models differ.

## 5.2.3 Daily Knock In Option

In order to wrap up this section, we check the difference in prices between micro and macro models on a contract that is heavily path dependent. For this purpose we have chosen a daily knock-in option. Table 7 shows a difference in prices between both models that is closer to being relevant.

## 5.3 Sensitivities

In this section we explore the effect that both models have in the sensitivities of the Athena Jet On S&amp;P Crude Oil Index Excess Return. This is interesting because the sensitivities for both the micro and the macro model refer to different quantities. More specifically, in Section 5.3.1 we compute the sensitivity of the micro and the macro model with respect to a change in the level of the future and the index, respectively. Next, in Section 5.3.2 we compute the

sensitivity of the micro and the macro model with respect to a change in the level of implied volatilities on the futures and the index, respectively.

## 5.3.1 Delta

The derivative of the Athena Jet price with respect to the underlying level represents a different thing for the micro and the macro model.

For the micro model, it represents the change in price with respect to a change in the value of the underlying futures curve. In the case of the first and second futures, its change also produces a change in the value of the index. In order to approximate the derivative we simply shift the values and redo the simulation, so that

$$\Delta _ { i } ^ { F } \coloneqq \frac { V _ { \text{micro} } ( ( 1 + 1 0 ^ { - 7 } ) \, F _ { t } ( T _ { i } ) ) - V _ { \text{micro} } ( F _ { t } ( T _ { i } ) ) } { 1 0 ^ { - 7 } F _ { t } ( T _ { i } ) }.$$

In the macro model, the delta refers to a change in the level of the index:

$$\Delta _ { \text{macro} } ^ { I } \coloneqq \frac { V _ { \text{macro} } ( ( 1 + 1 0 ^ { - 7 } ) \, I _ { t } ) - V _ { \text{macro} } ( I _ { t } ) } { 1 0 ^ { - 7 } I _ { t } }.$$

| ∆ F 1                  | ∆ I macro         |
|------------------------|-------------------|
| 0 . 0046046 0 . 000029 | 0 . 2848 0 . 0017 |

Table 8: Delta values for the micro and the macro computed according to Equations (57) and (58) respectively.

When we compute the deltas following Equations (57) and (58) we obtain the following the results in Table 8. We only show the derivative with respect to the first future because the deltas for the rest of the futures are negligible, actually they are of order O (10 -6 ) or less.

It is not surprising that the deltas from both models are different since we are calculating different things. Although this difference might not seem relevant it becomes very important in practice. The reason for this is that it is much easier to trade with futures than with the index, what potentially makes the macro model impractical. For this reason, we propose to define the following quantity in the micro model:

$$\Delta _ { \text{micro} } ^ { I } \coloneqq \Delta _ { 1 } ^ { F } \cdot \frac { 1 } { \frac { \partial I _ { t } } { \partial F _ { t } ( T _ { 1 } ) } } \,.$$

The new delta for the micro model defined in Equation (59) is the result of applying the chain rule to relate the impact that a change on the first future or a change of the index would have in the prices of our target contract. When the deltas are computed in that way the results are shown in Table 9. In this case, note that the deltas from both the micro and macro models are very similar. Since all the deltas with respect to futures are not meaningful except for the one on the first future and the fact that we have a deterministic relationship between the delta on the first future and the delta on the index, we can conclude that the macro model can be used to compute the delta.

| ∆ I micro         | ∆ I macro         |
|-------------------|-------------------|
| 0 . 2738 0 . 0017 | 0 . 2848 0 . 0017 |

±

±

Table 9: Delta values for the micro and the macro model computed according to Equations (59) and (58) respectively.

| T i        | ς F i               |
|------------|---------------------|
| 2020-01-21 | - 0 . 003 ± 0 . 034 |
| 2020-02-20 | - 0 . 027 ± 0 . 051 |
| 2020-03-20 | - 0 . 029 ± 0 . 063 |
| 2020-04-21 | - 0 . 024 ± 0 . 072 |
| 2020-05-19 | - 0 . 144 ± 0 . 068 |
| 2020-06-22 | - 0 . 315 ± 0 . 065 |
| 2020-07-21 | - 0 . 038 ± 0 . 051 |
| 2020-08-20 | - 0 . 053 ± 0 . 059 |
| 2020-09-22 | - 0 . 024 ± 0 . 052 |
| 2020-10-20 | - 0 . 003 ± 0 . 060 |
| 2020-11-20 | - 0 . 024 ± 0 . 041 |
| 2020-12-21 | 0 . 039 0 . 026     |

-

±

Table 10: Vega values for the micro model.

## 5.3.2 Vega

In order to complete the study we want to also compare the vegas in the micro and macro models. The vegas in the micro model are computed by shifting 1% the smile of the plain vanillas on the futures for a fixed maturity (here we denote by K the vector of strikes for that specific maturity):

$$\mathcal { V } _ { i } ^ { F } = \frac { V _ { \text{micro} } ( \sigma _ { \text{market} } ^ { F } ( T _ { i }, K ) + 0. 0 1 ) - V _ { \text{micro} } ( \sigma _ { \text{market} } ^ { F } ( T _ { i }, K ) ) } { 0. 0 1 }$$

More precisely, we select a futures contract and we only shift the plain vanillas on that specific future. This requires a recalibration of the leverage function. After the recalibration we obtain the results in Table 10. We can see that vega is small except for the futures that are close to the autocallability period. However, it may seem counter intuitive that the same effect does not happen when we approach the final payoff. The explanation is that when we change the probability of crossing the barrier in the middle of the contract, that change affects the number of paths that come to the end of the contract. This effect is much greater to that of simply changing the amount that you receive at the end.

At this point, it is difficult to find an equivalent for the macro model, since there is no direct mapping from the futures to the index. In any case, we suggest a similar procedure for the index, we fix a maturity and shift the market prices in terms of Black-Scholes volatilities by 1%, i.e. ,

$$\mathcal { V } _ { i } ^ { I } = \frac { V _ { \text{macro} } ( \sigma _ { \text{micro} } ^ { I } ( T _ { i } ^ { I }, K ^ { I } ) + 0. 0 1 ) - V _ { \text{macro} } ( \sigma _ { \text{micro} } ^ { I } ( T _ { i } ^ { I }, K ^ { I } ) ) } { 0. 0 1 }$$

| T i        | ς I i               |
|------------|---------------------|
| 2020-01-16 | 0 . 022 ± 0 . 036   |
| 2020-02-17 | - 0 . 009 ± 0 . 060 |
| 2020-03-16 | 0 . 013 ± 0 . 059   |
| 2020-06-16 | - 0 . 269 ± 0 . 036 |
| 2020-12-16 | 0 . 030 0 . 020     |

-

±

Table 11: Vega values for the macro model.

The 'vegas' V I i calculated with (5.3.2) are shown in Table 11. Both the results of Tables 11 and 10 are difficult to compare. The main point in common is that the vegas are greater when approaching the callability period and that both have the same sign. Apart from that, we didn't find any simple relationship between one and the other. We recall that options on the index are far less liquid than options on futures, thus making the micro model the best solution for hedging.

## 5.4 Computer implementation details

Concerning the hardware configuration, all tests have been performed in a Ubuntu server running over a virtualization layer (VMware) with 8 GB of RAM, 64 CPU cores (Intel(R) Xeon(R) CPU E5-2650 v4 at 2.2GHz for a total of 64 logical threads).

On the software side, we have done the implementation in the C ++ programming language. The GNU C ++ compiler has been used.

## 6 Conclusions and further research

In this paper we have presented two stochastic local volatility models for the pricing of derivative contracts on commodity indices:

- • A macro model that directly captures the dynamics of the index.
- • A micro model that captures the dynamics of the underlying futures curve.

Both models have been calibrated to market quotes by means of a hybrid global-local optimization algorithm.

The macro model produces results that are consistent with the more complex models that capture the behavior of the micro structure. In addition, it is very fast to execute and recalibrate since it requires the simulation of a single asset and the calibration of the stochastic parameters is based on a simple stochastic volatility model. However, it is not able to recover all the sensitivities that are needed for hedging and its calibration depends on data on the index, information that is scarce. Typically, the set of plain vanillas on the index is obtained from a consensus. The consensus is built from the data provided by different financial institutions on a monthly basis. Those financial institutions providing data that differs significantly from the consensus are expelled from the consensus data source. If we try to use it on a daily basis we will find that we either recalibrate the model once a month or we are forced to use another model that provides us a more stable stream of information.

The micro model is very robust. It produces prices for the derivatives on the index that are consistent with the models on the index. Its sensitivities refer to products that are actively

traded in the market what makes it better for hedging and trading purposes. Furthermore, it can be recalibrated on a daily basis, the stochastic parameters can be kept fixed and we can adjust the leverage function by recalibrating it to plain vanillas on the futures curve. Nevertheless, it is a very heavy model. The simulation of a contract depending on the index that expires one year into the future required the simulation of approximately 12 different SLV processes and the situation becomes worse as we go further into the future.

In view of the previous arguments, each model has its strengths and weaknesses and each of them can be effective on a different context.

Among the possible future research lines, we specifically mention two. On the one hand, we could price other indices, presumably multi-commodity ones. On the other hand, we could explore alternative models for the stochastic volatility. A possibility comes from extending the Heston model to admit powers of the stochastic volatility in the diffusion term of the stochastic volatility dynamics. This extension would be motivated by the fact that none of the parameters ( κ, θ, v 0 ) of the stochastic volatility affected the prices of the plain vanillas on the index. This suggests that a further modifications of the drift term will not have any impact on the prices of plain vanillas. However, the diffusion term had an impact. Therefore, by considering powers of the volatility in the diffusion term of the stochastic volatility dynamics we could enhance the calibration of the plain vanilla options on the index.

## Acknowledgements

A.P. Manzano-Herrero and C. Vázquez acknowledge the support from CITIC, as a center accredited for excellence within the Galician University System and a member of the CIGUS Network, receives subsidies from the Department of Education, Science, Universities, and Vocational Training of the Xunta de Galicia. Additionally, it is co-financed by the EU through the FEDER Galicia 2021-27 operational program (Ref. ED431G 2023/01). Also, both authors acknowledge the funding from Xunta de Galicia through the grant ED431C 2022/47, as well as the funding from Spanish Ministry of Science and Innovation with the grant PID2022141058OB-I00.

## References

Aarts, E. H. L. and P. J. M. Van Laarhoven (1985). 'Statistical cooling: A general approach to combinatorial optimization problems'. In: Philips Journal of Research 40.4, pp. 193-226. Alm, Thomas et al. (Sept. 2013). 'A Monte Carlo pricing algorithm for autocallables that allows for stable differentiation'. In: The Journal of Computational Finance 17, pp. 43-70. doi : 10.21314/JCF.2013.265 .

Chiminello, F. (2015). Oil goes local . Talk at Imperial College, London.

Deng, Geng, Craig Mccann, and Joshua Mallett (Aug. 2011). 'Modeling Autocallable Struc- tured Products'. In:

Journal of Derivatives &amp; Hedge Funds

17, pp. 326-340.

doi

:

10.1057/jdhf.2011.25

Derman, E. and I. Kani (1994). 'Riding on a smile'. In: Risk Magazine 7, pp. 32-39.

Dupire, Bruno (1994). 'Pricing with a Smile'. In: Risk Magazine 1, pp. 18-20.

and Urban Ulrych (2022). 'Pricing autocallables under

Farkas, Walter, Francesco Ferrari, local-stochastic volatility'. In: Swiss Finance Institute Research Paper 22-71.

.

Ferreiro, Ana M. et al. (2019). 'Parallel two-phase methods for global optimization on GPU'. In: Mathematics and Computers in Simulation 156, pp. 67-90. issn : 0378-4754. doi : https://doi.org/10.1016/j.matcom.2018.06.005 . url : https://www.sciencedirect.com/science/ar

Ferreiro-Ferreiro, Ana M. et al. (2020). 'A new calibration of the Heston Stochastic Local Volatility Model and its parallel implementation on GPUs'. In: Mathematics and Computers in Simulation 177, pp. 467-486.

Guyon, Julien and Pierre Henry-Labordère (2012). 'Being particular about calibration'. In: Risk Magazine 25.1, pp. 88-93.

Gyöngy, I. (1986). 'Mimicking the one-dimensional marginal distributions of processes having an Itô differential'. In: Probability Theory and Related Fields 71.4, pp. 501-516.

Hagan, Patrick S et al. (2002). 'Managing smile risk'. In: The Best of Wilmott 1, pp. 249-296. Heston, S. L. (1993). 'A closed-form solution for options with stochastic volatility with applications to bond and currency options'. In: The Review of Financial Studies 6.2, pp. 327343.

- Johnson, Steven G. (2007). The NLopt nonlinear-optimization package . https://github.com/stevengj/nlopt . Kennedy, J. and R. Eberhart (1995). 'Particle swarm optimization'. In: Proceedings of ICNN'95
- - International Conference on Neural Networks . Vol. 4, 1942-1948 vol.4. doi : 10.1109/ICNN.1995.488968 . Lipton, Alexander (2002). 'The vol smile problem'. In: Risk Magazine 15.2, pp. 61-66.

Lord, R., R. Koekkoek, and D. Van Dijk (2010). 'A comparison of biased simulation schemes for stochastic volatility models'. In:

Quantitative Finance

10.2, pp. 177-194.

Manzano, Alberto et al. (2023). 'Pricing commodity index options'. In:

Quantitative Finance

23.2, pp. 297-308. doi : 10.1080/14697688.2022.2138775 . eprint: https://doi.org/10.1080/14697688.2 url : https://doi.org/10.1080/14697688.2022.2138775 .

Nastasi, E., A. Pallavicini, and G. Sartorelli (2020). 'Smile modeling in commodity markets'.

In:

International Journal of Theoretical and Applied Finance

23.3, pp. 1-28.

Nelder, John A. and Roger Mead (1965). 'A Simplex Method for Function Minimization'.

In: The Computer Journal 7.4, pp. 308-313.

Paribas, BNP (2015). Certificate Athena Jet Su S&amp;P GSCI Crude Oil Index Excess Return .

url : https://www.borsaitaliana.it/cw-e-certificates/covered-warrant/brochureathenajet.pdf .

Pilz, K. and E. Schlögl (2011). 'A hybrid commodity and interest rate market model'. In: Quantitative Finance 4.13, pp. 543-560.

Rebonato, Ricardo (1999). Volatility and correlation: In the pricing of equity, FX and interestrate options . John Wiley &amp; Sons.

Ren, Yong., Dilip Madan, and Michael Qian Qian (2007). 'Calibrating and pricing with embedded local volatility models'. In: Risk Magazine 20.9, pp. 138-143.

Rowan, T. (1990). 'Functional Stability Analysis of Numerical Algorithms'. PhD thesis. Department of Computer Sciences, University of Texas at Austin.

Said, K. (1999). 'Pricing exotics under the smile'. In: Risk Magazine 12, pp. 72-75.

Silva Santos, Carlos Henrique da, Marcos Sergio Gonçalves, and Hugo Enrique HernándezFigueroa (2010). 'Designing Novel Photonic Devices by Bio-Inspired Computing'. In: IEEE Photonics Technology Letters 22.15, pp. 1177-1179. doi : 10.1109/LPT.2010.2051222 .

Storn, R. and K. Price (1997). 'Differential evolution - a simple and efficient heuristic for global optimization over continuous spaces'. In: Journal of Global Optimization 11.4, pp. 341-359.

Vikhar, Pradnya A. (2016). 'Evolutionary algorithms: A critical review and its future prospects'. In: 2016 International Conference on Global Trends in Signal Processing, Information Computing and Communication (ICGTSPICC) , pp. 261-265. doi : 10.1109/ICGTSPICC.2016.7955308 .

Wilmott, Paul, Sam Howison, and Jeff Dewynne (1995). The Mathematics of Financial Derivatives: A Student Introduction . Cambridge University Press.
Alberto Manzano Emanuele Nastasi Andrea Pallavicini Carlos Vázquez
Evaluating Microscopic and Macroscopic Models for Derivative Contracts on Commodity Indices
In this article, we analyze two modeling approaches for the pricing of
derivative contracts on a commodity index. The first one is a microscopic
approach, where the components of the index are modeled individually, and the
index price is derived from their combination. The second one is a macroscopic
approach, where the index is modeled directly. While the microscopic approach
offers greater flexibility, its calibration results to be more challenging,
thus leading practitioners to favor the macroscopic approach. However, in the
macroscopic model, the lack of explicit futures curve dynamics raises questions
about its ability to accurately capture the behavior of the index and its
sensitivities. In order to investigate this, we calibrate both models using
derivatives of the S\&P GSCI Crude Oil excess-return index and compare their
pricing and sensitivities on path-dependent options, such as autocallable
contracts. This research provides insights into the suitability of macroscopic
models for pricing and hedging purposes in real scenarios.