In financial markets, spread option is a derivative security with two underlying assets and the payoff of the spread option depends on the difference of these assets. We consider American style spread option which allows the owners to exercise it at any time before the maturity. The complexity of pricing American spread option is that the boundary of the corresponding partial differential equation which determines the option price is unknown and the model for the underlying assets is two-dimensional.In this dissertation, we incorporate the stochasticity to the interest rate and assume that it satisfies the Vasicek model or the CIR model. We derive the partial differential equations with terminal and boundary conditions which determine the American spread option with stochastic interest rate and formulate the associated free boundary problem. We convert the free boundary problem to the linear complimentarity conditions for the American spread option, so that we can go around the free boundary and compute the option price numerically. Alternatively, we approximate the option price using methods based on the Monte Carlo simulation, including the regression-based method, the Lonstaff and Schwartz method and the dual method. We make the comparisons among the option prices derived by the partial differential equation method and Monte Carlo methods to show the accuracy of the result.

Spread options are derivative securities, which are written on the difference between the values of two underlying market variables. They are very important tools to hedge the correlation risk. American style spread options allow the holder to exercise the option at any time up to and including maturity. Although they are widely used to hedge and speculate in financial market, the valuation of the American spread option is very challenging. Because even under the classic assumptions that the underlying assets follow the log-normal distribution, the resulting spread doesn't have a distribution with a simple closed formula. In this dissertation, we investigate the American spread option pricing problem. Several approaches for the geometric Brownian motion model and the stochastic volatility model are developed. We also implement the above models and the numerical results are compared among different approaches.