Exact simulation of normal tempered stable processes of OU type with applications

We study the Ornstein-Uhlenbeck process having a symmetric normal tempered stable stationary law and represent its transition distribution in terms of the sum of independent laws. In addition, we write the background driving Lévy process as the sum of two independent Lévy components. Accordingly, we can design two alternate algorithms for the simulation of the skeleton of the Ornstein-Uhlenbeck process. The solution based on the transition law turns out to be faster since it is based on a lower number of computational steps, as confirmed by extensive numerical experiments. We also calculate the characteristic function of the transition density which is instrumental for the application of the FFT-based method of Carr and Madan (J Comput Finance 2:61–73, 1999) to the pricing of a strip of call options written on markets whose price evolution is modeled by such an Ornstein-Uhlenbeck dynamics. This setting is indeed common for spot prices in the energy field. Finally, we show how to extend the range of applications to future markets.