This thesis explores numerical methods for solving nonlinear partial differential equations (PDEs) that arise in option pricing problems. The goal is to develop or identify robust and efficient techniques that converge to the financially relevant solution for both one and two factor problems. To illustrate the underlying concepts, two nonlinear models are examined in detail: uncertain volatility and passport options. For any nonlinear model, implicit timestepping techniques lead to a set of discrete nonlinear equations which must be solved at each timestep. Several iterative methods for solving these equations are tested. In the cases of uncertain volatility and passport options, it is shown that the frozen coefficient method outperforms two different Newton-type methods. Further, it is proven that the frozen coefficient method is guaranteed to converge for a wide class of one factor problems. A major issue when solving nonlinear PDEs is the possibility of multiple solutions. In a financial context, convergence to the viscosity solution is desired. Conditions under which the one factor uncertain volatility equations are guaranteed to converge to the viscosity solution are derived. Unfortunately, the techniques used do not apply to passport options, primarily because a positive coefficient discretization is shown to not always be achievable. For both uncertain volatility and passport options, much work has already been done for one factor problems. In this thesis, extensions are made for two factor problems. The importance of treating derivative estimates consistently between the discretization and an optimization procedure is discussed. For option pricing problems in general, non-smooth data can cause convergence difficulties for classical timestepping techniques. In particular, quadratic convergence may not be achieved. Techniques for restoring quadratic convergence for linear problems are examined. Via numerical examples, these techniques are also shown to improve the stability of the nonlinear uncertain volatility and passport option problems. Finally, two applications are briefly explored. The first application involves static hedging to reduce the bid-ask spread implied by uncertain volatility pricing. While static hedging has been carried out previously for one factor models, examples for two factor models are provided. The second application uses passport option theory to examine trader compensation strategies. By changing the payoff, it is shown how the expected distribution of trading account balances can be modified to reflect trader or bank preferences.

An option, a type of a financial derivative, is a contract that creates an opportunity for a market player to avoid risks involved in investing, especially in equities. An investor desires to know the accurate value of an option before entering into a contract to buy/sell the underlying asset (stock). There are various techniques that try to simulate real market conditions in order to price or evaluate an option. However, most of them achieved limited success due to high uncertainty in price behavior of the underlying asset. In this study, I propose two new Firefly variant algorithms to compute accurate worth for European and American option contracts and compare them with popular option pricing models (such as Black-Scholes-Merton, binomial lattice, Monte-Carlo, etc.) and real market data. In my study, I have first modelled the option pricing as a multi-objective optimization problem, where I introduced the pay-off and probability of achieving that pay-off as the main optimization objectives. Then, I proposed to use a latest nature-inspired algorithm that uses the bioluminescence of Fireflies to simulate the market conditions, a first attempt in the literature. For my thesis, I have proposed adaptive weighted-sum based Firefly algorithm and non-dominant sorting Firefly algorithm to find Pareto optimal solutions for the option pricing problem. Using my algorithm(s), I have successfully computed complete Pareto front of option prices for a number of option contracts from the real market (Bloomberg data). Also, I have shown that one of the points on the Pareto front represents the option value within 1-2 % error of the real data (Bloomberg). Moreover, with my experiments, I have shown that any investor may utilize the results in the Pareto fronts for deciding to get into an option contract and can evaluate the worth of a contract tuned to their risk ability. This implies that my proposed multi-objective model and Firefly algorithm could be used in real markets for pricing options at different levels of accuracy. To the best of my knowledge, modelling option pricing problem as a multi-objective optimization problem and using newly developed Firefly algorithm for solving it is unique and novel. / October 2016

Following the financial crisis of 2008, the need for more robust techniques to quantify the capital charge for risk management has become a pressing problem. Under Basel II/III, banks are allowed to calculate the capital charge using internally developed models subject to regulatory approval. An interesting problem for the regulator is to compare the resulting figures against the required capital under worst case scenarios. The existing literature on the latter problem, which is based on the marginal problem, assumes that no a-priori information is known about the dependencies of contributing risks. These problems are linear optimisation problems over a constrained set of probability measures, discretisation of which leads to large scale LPs. But this approach is very conservative and cannot be implemented robustly in practice, due to the scarcity of historical data. In our approach, we take a less conservative strategy by incorporating dependence information contained in the data in a form that still leads to LPs, an important feature of such problems due to their high dimensionality. Conceptually, our model is the discretisation of an infinite dimensional linear optimisation problem over a set of probability measures. For some specific cases we can prove strong duality, opening up the approach of discretising the dual instead of the primal. This approach is preferable, as it yields better numerical results. In this work we also apply our model to model-free path-dependent option pricing. Use of delayed column generation techniques allows us to solve problems several orders of magnitude larger than via the standard simplex algorithm. For high-dimensional LPs we also implement Nesterov's smoothing technique to solve the problems.

Title: Option Pricing Author: Radek Moravec Department: Department of Probability and Mathematical Statistics Supervisor: doc. RNDr. Jan Hurt, CSc., Department of Probability and Mathematical Statistics In the present thesis we deal with European call option pricing using lattice approaches. We introduce a discrete market model and show a way how to find an arbitrage price of financial instruments on complete markets. It's equal to the discounted value of future expected cash flow. We present the binomial option pricing model and generalize it into multinomial model. We test the resulting formula on real market data obtained from NYSE and NASDAQ. We suggest a parameter estimate method which is based on time series of historical observations of daily close price. We compare calculated option prices with their real market value and try to explain the reasons of the differences. 1

[pt] A expansão do mercado de derivativos no mundo e principalmente no Brasil tem impulsionado seus usuários a aprimorar e desenvolver ferramentas de apreçamento mais eficientes. Com esse intuito, o presente trabalho tem por objetivo evidenciar qual janela de observações gera a curtose e a assimetria que mais suavize o sorriso da volatilidade utilizando-se do modelo Corrado-Su. Para tanto, as empresas escolhidas foram a Petrobrás PN e a Vale PNA, devido a suas ações e opções de compra serem as mais líquidas no mercado brasileiro. A análise dos dados apontou para uma maior suavização do sorriso da volatilidade por parte das janelas de dados de curto prazo sobre as longo prazo, e uma equivalência de desempenho das primeiras ao do modelo Black-Scholes. / [en] The expansion of the derivatives market in the world and especially in Brazil has driven its users to enhance and develop tools for more efficient pricing. With this purpose, this paper aims to point which window of observations generates the kurtosis and skewness that more soften the volatility smile using the Corrado-Su model. Therefore, the firms that were chosen were Petrobras PN and Vale PNA, because their stocks and options are the most liquid in Brazilian market. The data analysis indicated a greater smoothing volatility smile using the windows of observations of the short term instead of the long term, and a equivalent performance of the first ones to that of the Black-Scholes model.

Modelos bastante utilizados atualmente no apreçamento de derivativos de taxas de juros realizam, muitas vezes, premissas excessivamente restritivas com relação à volatilidade da série do ativo objeto. O método de Black and Scholes e o de Vasicek, por exemplo, consideram a variância da série como constante no tempo e entre as diferentes maturidades, suposição que pode não ser a mais adequada para todos os casos. Assim, entre as técnicas alternativas de modelagem da volatilidade que vêm sendo estudadas, destacam-se as regressões por kernel. Discutimos neste trabalho a modelagem não paramétrica por meio da referida técnica e posterior apreçamento das opções em um modelo HJM Gaussiano. Analisamos diferentes especificações possíveis para a estimação não paramétrica da função de volatilidade através de simulações de Monte Carlo para o apreçamento de opções sobre títulos zero cupom, e realizamos um estudo empírico utilizando a metodologia proposta para o apreçamento de opções sobre IDI no mercado brasileiro. Um dos principais resultados encontrados é o bom ajuste da metodologia proposta no apreçamento de opções sobre títulos zero cupom. / Many models which have been recently used for derivatives pricing make restrictive assumptions about the volatility of the underlying object. Black-Scholes and Vasicek models, for instance, consider the volatility of the series as constant throughout time and maturity, an assumption that might not be the most appropriate for all cases. In this context, kernel regressions are important technics which have been researched recently. We discuss in this framework nonparametric modeling using the aforementioned technic and posterior option pricing using a Gaussian HJM model. We analyze different specifications for the nonparametric estimation of the volatility function using Monte Carlo simulations for the pricing of options on zero coupon bonds and conduct an empirical study using the proposed methodology for the pricing of options on the Interbank Deposit Index (IDI) in the Brazilian market. One of our main results is the good adjustment of the proposed methodology on the pricing of options on zero coupon bonds.

The main contribution of this thesis is to derive the properties and present a closed from solution of the exotic options under some specific types of Levy processes, such as American put options, American call options, British put options, British call options and American knock-out put options under either double exponential jump-diffusion processes or one-sided exponential jump-diffusion processes. Compared to the geometric Brownian motion, exponential jump-diffusion processes can better incorporate the asymmetric leptokurtic features and the volatility smile observed from the market. Pricing the option with early exercise feature is the optimal stopping problem to determine the optimal stopping time to maximize the expected options payoff. Due to the Markovian structure of the underlying process, the optimal stopping problem is related to the free-boundary problem consisting of an integral differential equation and suitable boundary conditions. By the local time-space formula for semi-martingales, the closed form solution for the options value can be derived from the free-boundary problem and we characterize the optimal stopping boundary as the unique solution to a nonlinear integral equation arising from the early exercise premium (EEP) representation. Chapter 2 and Chapter 3 discuss American put options and American call options respectively. When pricing options with early exercise feature under the double exponential jump-diffusion processes, a non-local integral term will be found in the infinitesimal generator of the underlying process. By the local time-space formula for semi-martingales, we show that the value function and the optimal stopping boundary are the unique solution pair to the system of two integral equations. The significant contributions of these two chapters are to prove the uniqueness of the value function and the optimal stopping boundary under less restrictive assumptions compared to previous literatures. In the degenerate case with only one-sided jumps, we find that the results are in line with the geometric Brownian motion models, which extends the analytical tractability of the Black-Scholes analysis to alternative models with jumps. In Chapter 4 and Chapter 5, we examine the British payoff mechanism under one-sided exponential jump-diffusion processes, which is the first analysis of British options for process with jumps. We show that the optimal stopping boundaries of British put options with only negative jumps or British call options with only positive jumps can also be characterized as the unique solution to a nonlinear integral equation arising from the early exercise premium representation. Chapter 6 provides the study of American knock-out put options under negative exponential jump-diffusion processes. The conditional memoryless property of the exponential distribution enables us to obtain an analytical form of the arbitrage-free price for American knock-out put options, which is usually more difficult for many other jump-diffusion models.

This diploma thesis deals with problem of option pricing with stochastic volatility. At first, the Black-Scholes model is derived and then its biases are discussed. We explain shortly the concept of volatility. Further, we introduce three pricing models with stochastic volatility- Hull-White model, Heston model and Stein-Stein model. At the end, these models are reviewed.

Modelos bastante utilizados atualmente no apreçamento de derivativos de taxas de juros realizam, muitas vezes, premissas excessivamente restritivas com relação à volatilidade da série do ativo objeto. O método de Black and Scholes e o de Vasicek, por exemplo, consideram a variância da série como constante no tempo e entre as diferentes maturidades, suposição que pode não ser a mais adequada para todos os casos. Assim, entre as técnicas alternativas de modelagem da volatilidade que vêm sendo estudadas, destacam-se as regressões por kernel. Discutimos neste trabalho a modelagem não paramétrica por meio da referida técnica e posterior apreçamento das opções em um modelo HJM Gaussiano. Analisamos diferentes especificações possíveis para a estimação não paramétrica da função de volatilidade através de simulações de Monte Carlo para o apreçamento de opções sobre títulos zero cupom, e realizamos um estudo empírico utilizando a metodologia proposta para o apreçamento de opções sobre IDI no mercado brasileiro. Um dos principais resultados encontrados é o bom ajuste da metodologia proposta no apreçamento de opções sobre títulos zero cupom. / Many models which have been recently used for derivatives pricing make restrictive assumptions about the volatility of the underlying object. Black-Scholes and Vasicek models, for instance, consider the volatility of the series as constant throughout time and maturity, an assumption that might not be the most appropriate for all cases. In this context, kernel regressions are important technics which have been researched recently. We discuss in this framework nonparametric modeling using the aforementioned technic and posterior option pricing using a Gaussian HJM model. We analyze different specifications for the nonparametric estimation of the volatility function using Monte Carlo simulations for the pricing of options on zero coupon bonds and conduct an empirical study using the proposed methodology for the pricing of options on the Interbank Deposit Index (IDI) in the Brazilian market. One of our main results is the good adjustment of the proposed methodology on the pricing of options on zero coupon bonds.

Barrier options have become increasingly popular financial instruments due to the lower costs and the ability to more closely match speculating or hedging needs. In addition, barrier options play a significant role in modeling and managing risks in insurance and finance as well as in refining insurance products such as variable annuities and equity-indexed annuities. Motivated by these immediate applications arising from actuarial and financial contexts, the thesis studies the pricing of barrier options and some exotic variations, assuming that the underlying asset price follows the Black-Scholes model or jump-diffusion processes. Barrier options have already been well treated in the classical Black-Scholes framework. The first part of the thesis aims to develop a new valuation approach based on the technique of exponential stopping and/or path counting of Brownian motions. We allow the option's boundaries to vary exponentially in time with different rates, and manage to express our pricing formulas properly as combinations of the prices of certain binary options. These expressions are shown to be extremely convenient in further pricing some exotic variations including sequential barrier options, immediate rebate options, multi-asset barrier options and window barrier options. Many known results will be reproduced and new explicit formulas will also be derived, from which we can better understand the impact on option values of various sophisticated barrier structures. We also consider jump-diffusion models, where it becomes difficult, if not impossible, to obtain the barrier option value in analytical form for exponentially curved boundaries. Our model assumes that the logarithm of the underlying asset price is a Brownian motion plus an independent compound Poisson process. It is quite common to assign a particular distribution (such as normal or double exponential distribution) for the jump size if one wants to pursue closed-form solutions, whereas our method permits any distributions for the jump size as long as they belong to the exponential family. The formulas derived in the thesis are explicit in the sense that they can be efficiently implemented through Monte Carlo simulations, from which we achieve a good balance between solution tractability and model complexity.