We propose a new accurate method for pricing European spread options by extending the lower bound approximation of Bjerksund and Stensland (2011) beyond the classical Black–Scholes framework. This is possible via a procedure requiring a univariate Fourier inversion. In addition, we are also able to obtain a new tight upper bound. Our method provides also an exact closed form solution via Fourier inversion of the exchange option price, generalizing the Margrabe (1978) formula. The method is applicable to models in which the joint characteristic function of the underlying assets forming the spread is known analytically. We test the performance of these new pricing algorithms performing numerical experiments on different stochastic dynamic models.

We present methodologies to price discretely monitored Asian options when the underlying evolves according to a generic Lévy process. For geometric Asian options we provide closed-form solutions in terms of the Fourier transform and we study in particular these formulas in the Lévy-stable case. For arithmetic Asian options we solve the valuation problem by recursive integration and derive a recursive theoretical formula for the moments to check the accuracy of the results. We compare the implementation of our method to Monte Carlo simulation implemented with control variates and using different parametric Lévy processes. We also discuss model risk issues.

We estimate a standard structural model of credit risk to draw insights about the premium demanded by investors for bearing default risk, using data on credit default swaps and market capitalization. We pin down the daily market value of assets for a set of non-financial firms and uncover cross-sectional heterogeneity in terms of the magnitude and time variation of the premium. By exploring the link between asset and default risk premia, we show that this heterogeneity closely depends on the relationship between the firm-specific market value of the assets and the business cycle.

We compute an analytical expression for the moment generating function of the joint random vector consisting of a spot price and its discretely monitored average for a large class of square-root price dynamics. This result, combined with the Fourier transform pricing method proposed by Carr and Madan [Carr, P., Madan D., 1999. Option valuation using the fast Fourier transform. Journal of Computational Finance 2(4), Summer, 61–73] allows us to derive a closed-form formula for the fair value of discretely monitored Asian-style options. Our analysis encompasses the case of commodity price dynamics displaying mean reversion and jointly fitting a quoted futures curve and the seasonal structure of spot price volatility. Four tests are conducted to assess the relative performance of the pricing procedure stemming from our formulae. Empirical results based on natural gas data from NYMEX and corn data from CBOT show a remarkable improvement over the main alternative techniques developed for pricing Asian-style options within the market standard framework of geometric Brownian motion.

We propose a consistent and computationally efficient 2-step methodology for the estimation of multidimensional non-Gaussian asset models built using Lévy processes. The proposed framework allows for dependence between assets and different tail behaviors and jump structures for each asset. Our procedure can be applied to portfolios with a large number of assets because it is immune to estimation dimensionality problems. Simulations show good finite sample properties and significant efficiency gains. This method is especially relevant for risk management purposes such as, for example, the computation of portfolio Value at Risk and intra-horizon Value at Risk, as we show in detail in an empirical illustration.