arxiv: v3 [q-fin.cp] 6 Jan 2011
|
|
- Gervais Francis
- 5 years ago
- Views:
Transcription
1 Phase transition in a log-normal Markov functional model Dan Pirjol Markit, 6 8 th Avenue, New York, NY 8 We derive the exact solution of a one-dimensional Markov functional model with log-normally distributed interest rates and constant volatility in discrete time. The model is shown to have two distinct limiting states, corresponding to small and asymptotically large volatilities, respectively. These volatility regimes are separated by a phase transition at some critical value of the volatility, at which certain expectation values display non-analytical behavior as a function of volatility. We investigate the conditions under which this phase transition occurs, and show that it is related to the position of the zeros of an appropriately defined generating function in the complex plane, in analogy with the Lee-Yang theory of the phase transitions in condensed matter physics. arxiv:7.69v3 [q-fin.cp] 6 Jan I. INTRODUCTION An important class of interest rate models, which includes many of the models currently used in practice, is the class of Markov-functional models [ 3]. The advantage of these models is that the value of discount bonds can be expressed as a functional of a low-dimensional Markov process. The specification of this functional dependence allows one to model the distribution of the forward rates with a prescribed probability distribution. This implies that these models can be calibrated exactly to a set of market instruments, such as caplets or swaptions. We consider a one-dimensional Markov-functional model with log-normally distributed forward Libors in discrete time. We show that the model can be solved exactly for a constant (time-independent) Libor volatility, and exact results can be found for the dependence of all discount bonds on the Markovian driver. The exact solution of the model is used to study its behaviorasafunction ofvolatility. We showthat themodel has two distinct regimes, in the low and large volatility limits, respectively. These regimes are separated by a sharp transition, occurring at some critical value of the volatility. We investigate the nature of this transition, and discuss the conditions under which it occurs. This volatility transition is similar to a first order phase transition in condensed matter physics, and is described by an analog of the Lee-Yang theory of phase transitions [8, 9]. The fundamental dynamical quantities are the zero coupon bonds P i,j P ti,t j. They are functions of a one-dimensional Markov process x(t), which will be assumed to be a simple Brownian motion with the usual properties E[x(t)] =,E[x (t)] = t. The model is defined by the probability distribution of the forward Libor rates L i (t i ) = τ i (P i,i+ ) for the (t i,t i+ ) period, with τ i t i+ t i. Wewillworkthroughoutinthet n forward measure, with numerairethe discount bond P t,tn. Specifically, in this measure the Libor rates L i will be assumed to be log-normally distributed L i = L i exp (ψx i ) ψ t i () For notational simplicity we denote the value of the Markov driver at time t i as x i x(t i ). We assume that the Libor volatility ψ is a constant, although a more general formulation of the model is possible, wherein ψ has term structure. For the purpose of illustrating the phenomenon considered here, it will be sufficient to consider a constant Libor volatility ψ. We denoted in Eq. () with L i the convexity-adjusted Libors; they are the expectation values of the Libor rates in the measure considered, and their determination is part of the solution of the model to be discussed below. The continuous time limit of this model follows the short rate process dr(t) r(t) = ψdx(t)+( d log r(t))dt (3) dt II. THE MODEL DEFINITION We consider a Markov functional model with discrete time evolution. The tenorstructure isafinite set ofdates = t < t < < t n () representing maturities equally spaced, e.g. by 3 or 6 months apart. where we introduced the convexity-adjusted forward short rate r(t) as the continuous time analog of L i. The model Eq.(3) describes a simple log-normally distributed short rate model without mean reversion [6, 7]. As discussedinref.[], meanreversioncanbeintroducedbyan appropriate choice of the time-dependence of ψ(t). Models with log-normally distributed rates in discrete time have been considered in [4, 5]. Denoting the numeraire-rebased zero coupon bond prices as pirjol@mac.com ˆP i,j = P i,j P i,n (4)
2 5 4 3 ψ=.4 ψ=.5 defined by appropriate interpolation. However, for timeindependent volatility ψ, it can be shown that the model can be solvedanalytically. In this case the solution of the recursion equation for the one-step rebased zero coupon bonds ˆP i,i+ has the general form ψ= ψ=. 9 ˆP i,i+ (x i ) = n i j= j ejψxi (jψ) t i (9) with j a set of constant coefficients. The convexityadjusted Libors are given by L i = ˆP,i ˆP,i+ N i τ i, () N i E[ˆP i,i+ f i (x i )] = n i j= j ejψ t i. ψ=.3 ψ= FIG.. The solution of the model for the convexity-adjusted Libors L i for a simulation with n = time steps (above), n = 4 (below), τ =.5 and constant forward short rate r = 5%, for several values of the volatility ψ. we note that the martingale condition for ˆP i,j can be expressed as ˆP i,j = E[ F i ] (5) P j,n Two particular cases of this relation are ˆP i,i+ = E[ˆP i+,i+ (+ L i+ τ i+ f i+ (x)) F i ] (6) ˆP,i = E[ˆP i,i+ (+ L i τ i f i (x))] (7) where we defined f i (x) = exp(ψx i ψ t i ). These relations can be solved recursively for ˆP i,i+ and L i, starting with the initial conditions L n τ n = ˆP,n, ˆPn,n = (8) and proceeding backwards in time. In the next section we present a method for solving these recursion relations in analytical form. III. ANALYTICAL SOLUTION In practice the expectation values in Eqs. (6) are computed by numerical integration over the functional dependence of ˆP i,i+ on the Markovian driver x i, which is The matrix of coefficients j has a triangular form (e.g. for n = 5) c (n ) c (n ) c (n ) ĉ = c (n 3) c (n 3) c (n 3) () c () c () c () c () 3 c () c () c () c () 3 c () 4 The coefficients j satisfy the recursion relation j = c (i+) j + L i+ τ i+ c (i+) t i+ j () e(j )ψ whichmustbesolvedsimultaneouslywitheq.()for L i. The initial condition is c (n ) =, L n τ n = ˆP,n. The recursion relation () can be solved backwards in time, for all i n, finding all coefficients in the matrix ĉ starting from the upper left corner and going downwards. Once the coefficients j and the convexity-adjusted Libors have been determined, all zero coupon bonds can be found as where P i,j (x i ) = ˆP i,j (x i ) ˆP i,i+ (x i )[+ L i τ i f i (x i )] (3) ˆP i,j (x i ) = E[ F i ] = E[ˆP j,j+ (+ P L e ψxj ψ t j ) F i ] j,n = n j k= n j + L c (j) k ekψxi (kψ) t i (4) k= c (j) k e(k+)ψxi (k +)ψ t i+kψ (t j t i). This completes the exact solution of the model.
3 3 The coefficients j satisfy certain general relations and sum rules, valid for arbitrary volatility. The first two coefficients, can be given in closed form = for i = n,n,, (5) n = L. (6) j=i+ A sum rule for the coefficients which will be useful in the following is n i j= j = ˆP,i+. (7) We illustrate the solution of the recursion in Figure, where we show results for the convexity adjusted Libors L i forseveralvaluesoftheliborvolatilityψ. Thetwonumerical examples considered assume a flat forward short rate r = 5% and a time discretization with quarterly time steps τ =.5. The terminal bond maturity is t n = 5 yr, and t n = yr, corresponding to n =,4 steps, respectively. From Fig. one observes that the convexity adjustment L fwd i L i is always positive, and increases with the volatility ψ. It is largest in the middle of the simulation interval, and with increasing volatility it becomes larger in a wider region expanding towards the beginning and the end of the simulation interval. This general behavior is expected on general grounds for the convexity adjusted rate L i provided that L i is positively correlated with the Libor associated with the payment delay L ti+,t n. The convexity adjusted rate appears to vanish in the middle of the simulation interval, for sufficiently large volatilities. The rate of vanishing has a sudden increase for volatilities larger than a certain value. This phenomenon is the main subject of this paper, and will be explained and quantified in Section IV below. A. Scaling The model is uniquely defined by the parameters {t i },{P,i },ψ. The zero coupon bonds P,i can be equivalently written in terms of the zero rates r i as P,i = exp( r i t i ). The model is invariant under a simultaneous rescaling of these parameters given by t i λt i (8) r i λ r i ψ λ / ψ. Under this transformation, the coefficients j and the expectation values N i are invariant, while the convexity adjusted Libors scale as L i λ Li. Due to this scaling invariance, the number of the relevant parameters of the model is reduced from three to two. B. Generating function In this section we present an efficient method for solvingthe recursionrelation() forthe coefficients j. We introduce the generating function at the time horizon t i (x) n i j= j xj (9) This function satisfies a recursion relation, expressing the generating function at time t i in terms of the generating function at the next time t i+ (x) = f (i+) (x)+ L i+ τxf (i+) (xe ψ t i+ ) () The initial condition for the recursion is f (n ) (x) =. The expectation value N i appearing in the expression for the convexity-adjusted Libor L i Eq. (6) is N i = (e ψ t i ) () The generating function (x) takes known values at x =, () =, () = ˆP,i+ () where the second relation follows from the sum rule Eq. (7). In the zero volatility limit ψ =, the generating function (x) can be found exactly (x) = Πn j=i+ (+Lfwd x). (3) where L fwd j = τ j (ˆP,j /ˆP,j+ ) is the forward Libor rate for the time period (t j,t j+ ). Expanding in powers of x this gives all the coefficients j in the zero volatility limit. As the volatility increases ψ >, the generating function (x) also changes, in such a way that the constraints Eqs. () are still satisfied. The recursion relation for (x) can be reformulated in such a way that it does not contain any reference to the convexity-adjusted Libors (x) = f (i+) (x)+(ˆp,i+ ˆP,i+ )x f(i+) (xe ψ t i+ ) f (i+) (e ψ t i+) (4) In the asymptotically large volatility limit ψ, this recursion relation can be again solved exactly, and thegeneratingfunction (x) is givenbythe asymptotic form (x) = +(ˆP,n )x+(ˆp,n ˆP,n )x + +(ˆP,i+ ˆP,i+ )x n i. (5) Note that the coefficients j have well-defined limiting values as ψ.
4 4 This showsthat the model considered has two verydifferent limiting regimes, corresponding to i) small volatility, and ii) large volatility. We will denote these regimes as the phases of the model. In each phase the generating function (x) has a well-defined expansion, given by Eqs. (3) and (5), respectively. We can use the results for the generating function in the zero and large volatility limits in order to obtain asymptotic expressions for the convexity adjusted Libors L i in the small and large volatility limits. This can be done using the relation between the expectation valuesn i definedin Eq.(), andthegeneratingfunction N i = (e ψ t i ). We start by considering first the asymptotics of L i in the small volatility limit ψ t i. The expansion of the generating function (x) around x = reads (x) = ˆP,i+ [+ n j=i+ +O((x ) )]. L fwd +L fwd (x ) (6) Using Eq. () this gives the small volatility asymptotics of the convexity-adjusted Libors L i, valid up to corrections of O((ψ t i ) ) L i L fwd i ( n j=i+ L fwd +L fwd (e ψ t i )) (7) ψ t i (Note that the useofthe zerovolatilitylimit ofthe generatingfunction (x) wassufficientin ordertoderivethis result. Thisisduetotheexactcondition () = ˆP,i+, which implies that the Taylor expansion of (x,ψ t i ) in powers of x and ψ t i contains only a linear term in x but not in ψ t i.) The low volatility approximation Eq. (7) has the familiar form of the convexity adjustment for a lognormally distributed rate. The growth of the convexity adjustment with the volatility has the familiar exponential form, proportional to exp(ψ t i ). Assuming a flat forward Libor curve, we have L fwd i = L fwd τ (erτ ). Then all terms in the sum overj are equal, and the result (7) simplifies as L i L fwd L fwd τ j i ( (n i ) +L fwd (e ψ t i ))(8) τ j This convexity adjustment is largest in the middle of the simulation interval, and vanishes near the boundaries. This agrees qualitatively with the main features of the convexity adjustment observed in Figure. Considernextthelargevolatilityasymptoticsψ t i of the convexity adjusted Libors L i. This follows from the large x asymptotics of the large volatility generating function f (i) (x) (x) (ˆP,i+ ˆP,i+ )x n i +O(x n i ).(9) Assuming again a flat forward Libor curve, this gives the large volatility asymptotics of the convexity adjusted Libors ( +L L i = L fwd fwd ) τ L fwd e (n i )ψ t i (3) τ ψ t i This shows that in the large volatility limit, the convexity adjusted rates L i drop off much faster with the volatility ψ. The decrease is still exponential, but it is much faster due to the additional factor n i (equal to the number of time steps to maturity) in the exponent L i exp( (n i )ψ t i ). The vanishing of the convexity adjusted Libors is faster in the middle of the simulation interval, just like in the small volatility case. As an aside, we note that all dynamical quantities of the model can be expressed formally in terms of the generating function (x). For example, the rebased bond prices are given by ˆP i,i+ (x) = exp( t i x) (e ψx ) (3) ˆP i,j (x) = exp( t i x)f (j) (e ψx ) (3) + L e ψx ψ t i f (j) (e ψx+ψ (t j t i) ), from which all zero coupon bond prices can be obtained using Eqs. (3). IV. LIBOR VOLATILITY TRANSITION As the volatility increases from zero to a large value, the coefficients j interpolate between the two limiting values, corresponding to low and asymptotically large volatilities, respectively. Equivalently, the generating function (x) changes between the two limiting expressions (x) and f(i) (x), in such a way that the two constraints Eq. () are still satisfied. As mentioned above, the convexity adjusted Libors L i appear to become vanishingly small at some value of the volatility, see Fig.. This phenomenon occurs first in the middle of the simulation time interval, and then it gradually extends also towards the boundaries. This is related to the expectation values N i which become very large as the volatility increases. To investigate this in more detail, we show in Fig. the plots of logn i as function of the volatility ψ for two simulations. We observe that the change is not gradual, but happens at a sharply defined value of the volatility, which will be called the critical volatility ψ cr. The transition becomes more sharp as n i increases. The critical volatility demarcates two regions of very different qualitative behaviour, in which the model has distinct limiting expressions for the functional dependence of discount bonds on the Markov driver.
5 ψ ψ FIG. 3. The zeros of the generating function (x) in the complex plane for several values of the volatility ψ =.5,.55,.6 around the critical value ψ cr =.53 (blue, violet, red). The black dots correspond to the large volatility asymptotic generating function f (i) (x), which is obtained for ψ =. The parameters are the same as in the upper plot of Fig., and the time slice t i considered is i =. The circle shown has radius exp(ψcrt i). FIG.. The log of N i as a function of the volatility ψ (solid line), and its large volatility approximation as N i = f (i) (exp(ψ t i)) (dashed line). Above: n =,i =. Below: n = 4,i = 3. Both cases correspond to τ =.5,r = 5%. The convexity-adjusted Libors L i become very small for volatilities above the critical value ψ > ψ cr. In practice, they can become so small that they are below machine precision. This phenomenon thus imposes a limit to the practical applicability of the model, and it is important to understand the conditions under which it occurs. In the following we investigate in some detail the nature of this transition, and formulate a criterion for finding the critical volatility ψ cr at each time horizon t i. We will show that the singular behaviour of the expectation values N i is related to the distribution of the zeros of the generating function (z) in the complex plane. The generating function (z) is a polynomial of z with positive coefficients, and thus does not have any zeros on the positive real axis z >. However, it is well known that the position of the zeros in the complex plane can influence the behaviour of the function along the real axis. The generating function (z) has n i zeros. At zero volatility, all zeros are at large real negative values z k = /(L fwd k τ k ), but they migrate in the complex plane as pairs of complex conjugate values as the volatility increases, and surround the origin, see Fig. 3. As the volatility increases to very large values, the zeros reach fixed positions, given by the zeros of the asymptotic generating function (x). Furthermore, as the polynomial order n i increases, the number of the zeros increases and they close in on the positive real axis, pinching it at some point z >. We will show that the function (z) is continuousat z, but its derivativehas ajump. Thusthe generating function has a cusp at this point. Recalling that N i is related to the generating function as N i = (e ψ t i ), see Eq. (), it follows that N i has a singular behaviour at the volatility ψ cr, given by the equation exp(ψ crt i ) = z. (33) This equation determines the critical volatility ψ cr at the time slice t i. Geometrically, this has the following meaning: the critical volatility at the time horizon t i is given by that value of ψ for which the zeros of the generating function (z) enter the circle of radius exp(ψ t i ). As ψ increases, the zeros move closer to the origin, while the circle of radius exp(ψ t i ) expands, such that at some intermediate value ψ cr, the zeros will cross the expanding circle. This picture is illustrated in Figure 3 on the example of the transition shown in the upper plot of Figure. The transition between the two volatility phases is similar to a phase transition in the Lee-Yang formalism [9], where critical points are associated with the values of the thermodynamic parameter (fugacity) at which the zeros of the grand canonical partition function pinch the real axis. For a general introduction to phase transitions see
6 6 [8]. As shown in [9], the partition function is continuous at the critical point, but its derivative has a jump, which is proportional to the density of zeros around this point. These results hold in the thermodynamical limit of an infinite volume; the analog of the thermodynamical limit in our case is n i, where n i is the number of time steps from maturity to the time slice considered t i. SuchabehaviourispreciselywhatisobservedinFig., where one can see that logn i is continuous everywhere, but its derivative has a jump at ψ cr. The analog of the partition function in our case is the generating function (z), and the jump occurs at z = exp(ψ crt i ) In the following we quantify this statement, and compute an explicit result for the discontinuity of the derivative of (z) at z. The generating function can be written explicitly in terms of its zeros z k as (z) = Π n i k= ( z/z k ) (34) As the polynomial order n f = n i increases, the roots arrange themselves on a closed curve around the origin, which can be parameterized in polar coordinates asz(θ) = ρ(θ)e iθ. Therootsappearincomplexconjugate pairs, which implies that the curve describing the zeros is symmetric under reflection on the real axis ρ(θ) = ρ( θ). The logarithm of the generating function (z) can be expressed in the limit n f as an integral log (z)= π dθg(θ)log( z ρ (θ) cosθ z ρ(θ) +) (35) where g(θ) is the density of roots at polar angle θ. It is normalized as π The derivative of log (z) is d dz logf(i) (z)= π dθg(θ) = n f (36) z ρ(θ)cosθ dθg(θ) z ρ(θ)cosθz +ρ (θ) (37) This is discontinuous across the boundary at z = ρ() with a jump given by the density of zeros at this point d dz logf(i) (z) z=ρ()+ǫ d dz logf(i) (z) z=ρ() ǫ = πg() ρ(). (38) These results are completely analogous to the expressions derived in [9] for the jump of the derivative of the grand canonical function at a critical point. A similar phenomenon occurs for any expectation value of the form similar to N i, with φ a real number E[ˆP i,i+ e φx φ t i ] = (e ψφti ) (39) The expectation value can be expressed in terms of the generatingfunction (x) asshown. The criticalvolatility corresponding to this expectation value is found in analogy to Eq. (33) and is given by exp(ψφt i ) = z x FIG. 4. The integrand in the expression (4) for the expectation value N i at time t i =.5 for several values of the volatility ψ: i).4 (solid line),.5 (dashed line),.5 (dotted line). The parameters of the simulation are n =,i =, r = 5%,τ =.5. Note the second maximum which appears for volatilities close to the critical volatility ψ cr. We close this section with a comment about the practical observability of the phase transition in usual implementations of the Markov functional model. The expectation value N i defined in () is written explicitly as N i = E[ˆP i,i+ f i (x)] (4) = dx i πti e x i t i ˆPi,i+ (x i )e ψxi ψ t i This integral is usually assumed to be dominated by contributions from small values of x i around the origin x i κ t i, where κ 3 5. The integral is then evaluated numerically either on a grid, or using Monte Carlo simulations. However, for volatilities ψ above the critical value, the integrand develops a second local maximum at large values of the Markov driver as seen in Figure 4, which will dominate the integral above the critical volatility. In the example of Figure 4, the secondary maximum appearsatx,whichisalmoststandarddeviations away from zero. Thus the integral (4) will receive significant contributions from a region in x which is sampled very inefficiently in Monte Carlo or grid methods. This implies that the model will not be simulated correctly in the large volatility phase, and the phase transition will be unobservable under these simulation methods. V. EXAMPLE: CONSTANT FORWARD SHORT RATE We illustrate the general results discussed in the previous section on the example of a forward yield curve with constant short rate r. The initial yield curve is P,i = exp( r t i ), and the numeraire rebased discount bonds are ˆP,i = exp(r (n i)τ).
7 7 In Figure 3 we show the zeros z k of the generating function (z) at the time slice i =, for a simulation with r = 5%, total simulation time t n = 5 yr, with time step τ =.5 and n = time steps. The colored dots show the exact zeros at three values of the volatility ψ =.5,.55,.6around the critical volatility ψ cr =.53, and the black dots are the zeros of the asymptotic generating function f (i) (z). We note that the asymptotic generating function f (i) (z) gives a reasonably good approximation for the position of the zeros in the large volatility phase. For this reason we discuss in some detail the position of the zeros of the asymptotic generating function f (i) (z) for which an analytical treatment is possible in the limit of a constant forward short rate. The asymptotic generating function for a constant forward short rate r is n i f (i) (x) = +[ e rτ ] (e rτ x) j (4) j= Its zeros are x k = e rτ z k, where z k are the zeros of the polynomial p nf (z) of degree n f = n i p nf (z) e rτ +z +z + +z n f (4) Under usual market conditions /( exp( r τ)) >, and by the Eneström-Kakeya theorem [], all zeros of this polynomial lie outside the open unit disk z k >. The structure of the zeros of p nf (z) can be studied by noting that this polynomial is the truncated Taylor series of the function F(z) = e rτ + z (43) The function F(z) has an exact zero at z = e rτ, and a pole at z =, which means that the convergence region is the circle z <. The theory of the zeros of a truncated Taylor expansion is a well studied subject in approximation theory. The main result is the Jentzsch-Szegö theorem [, ], according to which the zeros of the truncated Taylor series of a function F(z) either converge to the zeros of F(z) as n f, provided that they are inside its convergence region, or they accumulate on the boundary of the convergence region. For our case the latter situation applies, such that the zeros of p nf (z) will accumulate uniformly on the z = circle as the polynomialordern f increases. This implies that the zeros of f (i) (x) close in on the positive real axis at x = ρ e rτ as the polynomial order n f = n i becomes large. This phenomenon is visible already at moderate values of n f O(), as seen from Fig. 3. We can compute the properties of the generating function at the critical point x by applying the general results discussed above. Taking ρ(θ) = ρ e rτ and g(θ) = n f /(π), theintegralineq.(35)canbeperformed exactly with the result { logf (i) nf logρ, x ρ (x) = n f logx, x ρ (44) The jump of the derivative across the critical point x = ρ is d dx logf(i) (x = ρ+ǫ) d dx logf(i) (x = ρ ǫ)= n f ρ (45) In the large volatility phase the derivative of log (x) is very large, and is of the same order of magnitude as expected from the asymptotic form Eq. (9). Finally, we consider the case of practical interest of finite polynomial order n f = n i. It was observed in [3] that for the truncated Taylor series a + a z + +a n z n of a function F(z), a good approximation for the moduli of the zeros z k is obtained by neglecting all but the first and last terms (by considering the simpler polynomial p(z) a +a n z n ), provided that the function F(z) does not have zeros within the convergence region. In our case of the polynomial p nf (z), the corresponding polynomial p nf (z) = /( e rτ )+z n f has zeros z k = ( e rτ )/(n i ) e πi k n i, (46) k =,,,n i The critical volatility can be computed using equation Eq. (33) and is given by e rτ+ψ cr ti = ( e rτ )/(n i ), (47) ψcr i(n i )τ log( r τ ). The minimum value of ψ cr is reached at i = [n/], in the middle of the simulation interval, where i(n i ) is maximal. This agrees with the shape of the L i curves in Fig., where the critical volatility is first reached in the middle of the interval. Thus the practical applicability range of the model is restricted to volatilities smaller than the minimum critical volatility ψ < (ψ cr) min = [n/] τ log( ). (48) r τ This expression satisfies the general scaling properties of the model Eqs. (8). InTable weshowthe valuesofthe maximallyallowed volatility for several values of the total simulation tenor t n = nτ, the time step τ, and the short rate r. The maximum allowed volatility decreases with the size of the time step τ, with the tenor of the simulation t n, and with the short rate r.
8 8 t n = 5 t n = t n = t n = 3 r τ =.5 τ =.5 τ =.5 τ =.5 τ =.5 τ =.5 τ =.5 τ =.5 % 48.95% 65.% 4.48% 3.55%.4% 6.8% 8.6%.85% % 46.4% 6.7% 3.% 3.35%.5% 5.7% 7.67%.% 3% 44.4% 57.96%.% 8.98%.6% 4.49% 7.37% 9.66% 4% 4.9% 55.94%.46% 7.97%.73% 3.99% 7.5% 9.3% 5% 4.87% 54.3%.93% 7.6%.47% 3.58% 6.98% 9.5% TABLE I. The maximal Libor volatility ψ for which the model is everywhere below the critical volatility ψ cr, determined according to Eq.(47), for several choices of the total tenor t n, time step τ and the level of the interest rates r. VI. CONCLUSIONS We discussed in this paper the behaviour of a Markov functional model with discrete log-normally distributed Libors, as a function of the Libor volatility. The model can be solved exactly in the T forward measure in the limit of a time-independent volatility, by a backwards recursion relation. Analytical results can be obtained for the functional dependence of all discount bonds on the Markov driver. We showed that the model has two volatility regimes, corresponding to small and large volatility, with very different qualitative behaviour. As the Libor volatility increases, there is a transition between the two regimes, at an intermediate critical volatility. In the large volatility phase the convexity adjusted Libors are very small, and can be below machine precision. Thus the existence of the large volatility regime imposes a limit on the applicability of such a model, which is manifested as an upper bound on the allowed Libor volatility. We formulated the conditions under which this phase transition occurs, and showed that it is related to the position of the zeros of an appropriately defined generating function in the complex plane. The transition appears for long simulation times, and small time discretization steps. A similar phenomenon is expected to occur also for the practically relevant but analytically more complex case of time-dependent volatility ψ(t i ). Also, the discussion of this paper is limited to the T forward measure, and one expects that the details of the volatility dependencecouldchangeinadifferentmeasure,but nottheexistence of a phase transition. We hope to report progress along these directions in future work. Appendix A: The zeros of the generating function (x) We study here the distribution of the complex zeros of the asymptotic generating function f (i) (z) for a constant forward short rate r. This is related to the problem of finding the zeros of the polynomial p n (z) a+z +z + +z n (A) FIG. 5. The zeros of the polynomial p 9(z) for a = 8. The solid line denotes the curve ρ(θ) given in Eq. (A), and the dashed line shows the approximation for the moduli of the zeros z k = a /n obtained by retaining in p n(z) only the first and last terms. whereaisarealnumberlargerthan. Theyarethesame as the zeros of the equation z n+ +(a )z a = from which z = is excluded. Expressed in polar coordinates z = ρe iθ one finds that the zeros are on the curve ρ(θ) which is given by the solution to the equation ρ (n+) (θ)=(a ) ρ (θ) a(a )ρ(θ)cosθ +a (A) Figure 5 shows the exact roots of p 9 (z), along with the curve ρ(θ) (solid line), and the circle with radius a /n (dashed line), which was used in the main text. We note that the latter approximates the moduli of the zeros very well, as noted in [3]. The curve ρ(θ) intersects the real axis at some point ρ(). The solution of the equation (A) at θ = has always the solution ρ() =, and in addition it can have another solutions, or none, depending on the values of n,a. The solutions different from can be found as the solutions of the simpler equation h(ρ) = ρ n+ (a )ρ+a =. The function h(ρ) has a minimum at ρ =
9 9 [(a )/(n+)] /n. We distinguish the 3 cases, according to the value of h(ρ ).h(ρ ) < solutions for ρ.h(ρ ) = ρ = ρ 3.h(ρ ) > no solutions for ρ (A3) (A4) (A5) These cases are obtained for n < [n ],n = [n ],n > [n ] respectively, where n is the solution of the equation ( a n )n = ( a n+ )n+. (A6) For a = 8 one has [n ] =. A good approximation for the solutions in case is obtained by Taylor expanding h(ρ) around ρ to quadratic order. This gives ρ = ρ ± a n(n+)ρ n. (A7) nρn+ In Fig. 5 we show only the largest of the 3 solutions. The angular distribution of the roots can be obtained from the imaginary part of the equation for z, which reads ρ n sin[(n+)θ] = (a )sinθ. For the subset of even n, this has real solutions for ρ only if θ is in one of the regions π θ ((k ) n+,k π ), k =,, n (A8) n+ There is one solution in each of these intervals of equal angular opening, which means that in the large n limit, the angular distribution of the roots approaches a uniform distribution g(θ) = n/(π). In the same limit, the curveρ(θ)approachestheunitcircleρ(θ) =,asrequired by the Jentzsch-Szegö theorem. ACKNOWLEDGEMENTS I am grateful to Dyutiman Das and Adrian Ghinculov for useful discussions on this problem, and to Radu Constantinescu for comments on the manuscript. [] P. Hunt, J. Kennedy and A. Pellser, Markov-Functional Interest Rate Models, Finance and Stochastics, 4, (). [] J. B. Hunt and J. E. Kennedy, Financial Derivatives in Theory and Practice, Wiley Series in Probability and Statistics, 5. [3] P. Balland and L. P. Hughston, Markov Market Model Consistent with Cap Smile, Int. J. Th. Appl. Finance 3, 6-8 (). [4] P. Glasserman and X. Zhao, Arbitrage free discretization of log-normal forward Libor and swap rate models, Finance and Stochastics 4, () [5] K. Miltersen, L. Sandmann and D. Sondermann, Closed Form Solutions for Term Structure Derivatives with Lognormal Interest Rates, J. Finance 5, (997). [6] L. U. Dothan, On the Term Structure of Interest Rates, Journal of Financial Economics 6, (978). [7] D. Brigo and F. Mercurio, Interest Rate Models - Theory and Practice: With Smile, Inflation and Credit, Springer Verlag 6. [8] Eugene H. Stanley, Introduction to Phase Transitions and Critical Phenomena, Oxford University Press, 987. [9] T. D. Lee and C. N. Yang, Statistical Theory of Equations of State and Phase Transitions. II. Lattice Gas and Ising Model, Physical Review Letters 87, 4-49 (95); Phys. Rev. 87, 4 (95). [] P. B. Borwein and T. Erdelyi, Polynomials and Polynomial Inequalities, Graduate Texts in Mathematics 6, Springer Verlag, 995. [] R. Jentzsch, Untersuchungen zur Theorie der Folgen analytischer Funktionen, Acta Math. 4, 9-5 (98). [] G. Szegö, Über die Nullstellen von Polynomen, die in einem Kreis gleichmassig konvergieren, Sitzungsber. Ber. Math. Ges.,, (9). [3] S. Christiansen, Per A. Madsen, On Truncated Taylor series and the position of their spurious zeros, Applied Numerical Mathematics, 56, 9-4 (6).
Phase Transition in a Log-Normal Interest Rate Model
in a Log-normal Interest Rate Model 1 1 J. P. Morgan, New York 17 Oct. 2011 in a Log-Normal Interest Rate Model Outline Introduction to interest rate modeling Black-Derman-Toy model Generalization with
More informationMarket interest-rate models
Market interest-rate models Marco Marchioro www.marchioro.org November 24 th, 2012 Market interest-rate models 1 Lecture Summary No-arbitrage models Detailed example: Hull-White Monte Carlo simulations
More informationInterest Rate Volatility
Interest Rate Volatility III. Working with SABR Andrew Lesniewski Baruch College and Posnania Inc First Baruch Volatility Workshop New York June 16-18, 2015 Outline Arbitrage free SABR 1 Arbitrage free
More informationOptimal stopping problems for a Brownian motion with a disorder on a finite interval
Optimal stopping problems for a Brownian motion with a disorder on a finite interval A. N. Shiryaev M. V. Zhitlukhin arxiv:1212.379v1 [math.st] 15 Dec 212 December 18, 212 Abstract We consider optimal
More informationSemi-analytic Lattice Integration of a Markov Functional Term Structure Model
Semi-analytic Lattice Integration of a Markov Functional Term Structure Model... Christ Church College University of Oxford A thesis submitted for the degree of Master of Science in Mathematical Finance
More informationThe Yield Envelope: Price Ranges for Fixed Income Products
The Yield Envelope: Price Ranges for Fixed Income Products by David Epstein (LINK:www.maths.ox.ac.uk/users/epstein) Mathematical Institute (LINK:www.maths.ox.ac.uk) Oxford Paul Wilmott (LINK:www.oxfordfinancial.co.uk/pw)
More informationA THREE-FACTOR CONVERGENCE MODEL OF INTEREST RATES
Proceedings of ALGORITMY 01 pp. 95 104 A THREE-FACTOR CONVERGENCE MODEL OF INTEREST RATES BEÁTA STEHLÍKOVÁ AND ZUZANA ZÍKOVÁ Abstract. A convergence model of interest rates explains the evolution of the
More informationEFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS
Commun. Korean Math. Soc. 23 (2008), No. 2, pp. 285 294 EFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS Kyoung-Sook Moon Reprinted from the Communications of the Korean Mathematical Society
More informationRiccardo Rebonato Global Head of Quantitative Research, FM, RBS Global Head of Market Risk, CBFM, RBS
Why Neither Time Homogeneity nor Time Dependence Will Do: Evidence from the US$ Swaption Market Cambridge, May 2005 Riccardo Rebonato Global Head of Quantitative Research, FM, RBS Global Head of Market
More informationMarket Risk Analysis Volume I
Market Risk Analysis Volume I Quantitative Methods in Finance Carol Alexander John Wiley & Sons, Ltd List of Figures List of Tables List of Examples Foreword Preface to Volume I xiii xvi xvii xix xxiii
More informationHints on Some of the Exercises
Hints on Some of the Exercises of the book R. Seydel: Tools for Computational Finance. Springer, 00/004/006/009/01. Preparatory Remarks: Some of the hints suggest ideas that may simplify solving the exercises
More informationLibor Market Model Version 1.0
Libor Market Model Version.0 Introduction This plug-in implements the Libor Market Model (also know as BGM Model, from the authors Brace Gatarek Musiela). For a general reference on this model see [, [2
More informationMath 416/516: Stochastic Simulation
Math 416/516: Stochastic Simulation Haijun Li lih@math.wsu.edu Department of Mathematics Washington State University Week 13 Haijun Li Math 416/516: Stochastic Simulation Week 13 1 / 28 Outline 1 Simulation
More informationInterest rate volatility
Interest rate volatility II. SABR and its flavors Andrew Lesniewski Baruch College and Posnania Inc First Baruch Volatility Workshop New York June 16-18, 2015 Outline The SABR model 1 The SABR model 2
More informationSYLLABUS. IEOR E4728 Topics in Quantitative Finance: Inflation Derivatives
SYLLABUS IEOR E4728 Topics in Quantitative Finance: Inflation Derivatives Term: Summer 2007 Department: Industrial Engineering and Operations Research (IEOR) Instructor: Iraj Kani TA: Wayne Lu References:
More information16. Inflation-Indexed Swaps
6. Inflation-Indexed Swaps Given a set of dates T,...,T M, an Inflation-Indexed Swap (IIS) is a swap where, on each payment date, Party A pays Party B the inflation rate over a predefined period, while
More informationModel Validation for Interest Rate Models
Model Validation for Interest Rate Models Dan Pirjol 1 1 J. P. Morgan, New York 24 Oct. 2011 Outline Introduction Model Validation for Interest Rate Models Regulatory Mandates Types of Interest Rate Models
More informationTN 2 - Basic Calculus with Financial Applications
G.S. Questa, 016 TN Basic Calculus with Finance [016-09-03] Page 1 of 16 TN - Basic Calculus with Financial Applications 1 Functions and Limits Derivatives 3 Taylor Series 4 Maxima and Minima 5 The Logarithmic
More informationLecture 4. Finite difference and finite element methods
Finite difference and finite element methods Lecture 4 Outline Black-Scholes equation From expectation to PDE Goal: compute the value of European option with payoff g which is the conditional expectation
More informationMartingale Methods in Financial Modelling
Marek Musiela Marek Rutkowski Martingale Methods in Financial Modelling Second Edition \ 42 Springer - . Preface to the First Edition... V Preface to the Second Edition... VII I Part I. Spot and Futures
More informationEquity correlations implied by index options: estimation and model uncertainty analysis
1/18 : estimation and model analysis, EDHEC Business School (joint work with Rama COT) Modeling and managing financial risks Paris, 10 13 January 2011 2/18 Outline 1 2 of multi-asset models Solution to
More informationPDE Methods for the Maximum Drawdown
PDE Methods for the Maximum Drawdown Libor Pospisil, Jan Vecer Columbia University, Department of Statistics, New York, NY 127, USA April 1, 28 Abstract Maximum drawdown is a risk measure that plays an
More informationSaddlepoint Approximation Methods for Pricing. Financial Options on Discrete Realized Variance
Saddlepoint Approximation Methods for Pricing Financial Options on Discrete Realized Variance Yue Kuen KWOK Department of Mathematics Hong Kong University of Science and Technology Hong Kong * This is
More informationarxiv: v1 [q-fin.pm] 12 Jul 2012
The Long Neglected Critically Leveraged Portfolio M. Hossein Partovi epartment of Physics and Astronomy, California State University, Sacramento, California 95819-6041 (ated: October 8, 2018) We show that
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Other Miscellaneous Topics and Applications of Monte-Carlo Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationAN ANALYTICALLY TRACTABLE UNCERTAIN VOLATILITY MODEL
AN ANALYTICALLY TRACTABLE UNCERTAIN VOLATILITY MODEL FABIO MERCURIO BANCA IMI, MILAN http://www.fabiomercurio.it 1 Stylized facts Traders use the Black-Scholes formula to price plain-vanilla options. An
More informationA Hybrid Commodity and Interest Rate Market Model
A Hybrid Commodity and Interest Rate Market Model University of Technology, Sydney June 1 Literature A Hybrid Market Model Recall: The basic LIBOR Market Model The cross currency LIBOR Market Model LIBOR
More informationCALIBRATION OF THE HULL-WHITE TWO-FACTOR MODEL ISMAIL LAACHIR. Premia 14
CALIBRATION OF THE HULL-WHITE TWO-FACTOR MODEL ISMAIL LAACHIR Premia 14 Contents 1. Model Presentation 1 2. Model Calibration 2 2.1. First example : calibration to cap volatility 2 2.2. Second example
More informationIntroduction to Financial Mathematics
Department of Mathematics University of Michigan November 7, 2008 My Information E-mail address: marymorj (at) umich.edu Financial work experience includes 2 years in public finance investment banking
More informationL 2 -theoretical study of the relation between the LIBOR market model and the HJM model Takashi Yasuoka
Journal of Math-for-Industry, Vol. 5 (213A-2), pp. 11 16 L 2 -theoretical study of the relation between the LIBOR market model and the HJM model Takashi Yasuoka Received on November 2, 212 / Revised on
More informationFast and accurate pricing of discretely monitored barrier options by numerical path integration
Comput Econ (27 3:143 151 DOI 1.17/s1614-7-991-5 Fast and accurate pricing of discretely monitored barrier options by numerical path integration Christian Skaug Arvid Naess Received: 23 December 25 / Accepted:
More informationInterest Rate Modeling
Chapman & Hall/CRC FINANCIAL MATHEMATICS SERIES Interest Rate Modeling Theory and Practice Lixin Wu CRC Press Taylor & Francis Group Boca Raton London New York CRC Press is an imprint of the Taylor & Francis
More informationThe Lognormal Interest Rate Model and Eurodollar Futures
GLOBAL RESEARCH ANALYTICS The Lognormal Interest Rate Model and Eurodollar Futures CITICORP SECURITIES,INC. 399 Park Avenue New York, NY 143 Keith Weintraub Director, Analytics 1-559-97 Michael Hogan Ex
More informationRohini Kumar. Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque)
Small time asymptotics for fast mean-reverting stochastic volatility models Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque) March 11, 2011 Frontier Probability Days,
More informationProbability. An intro for calculus students P= Figure 1: A normal integral
Probability An intro for calculus students.8.6.4.2 P=.87 2 3 4 Figure : A normal integral Suppose we flip a coin 2 times; what is the probability that we get more than 2 heads? Suppose we roll a six-sided
More informationImplementing the HJM model by Monte Carlo Simulation
Implementing the HJM model by Monte Carlo Simulation A CQF Project - 2010 June Cohort Bob Flagg Email: bob@calcworks.net January 14, 2011 Abstract We discuss an implementation of the Heath-Jarrow-Morton
More informationCharacterization of the Optimum
ECO 317 Economics of Uncertainty Fall Term 2009 Notes for lectures 5. Portfolio Allocation with One Riskless, One Risky Asset Characterization of the Optimum Consider a risk-averse, expected-utility-maximizing
More informationAsymptotic methods in risk management. Advances in Financial Mathematics
Asymptotic methods in risk management Peter Tankov Based on joint work with A. Gulisashvili Advances in Financial Mathematics Paris, January 7 10, 2014 Peter Tankov (Université Paris Diderot) Asymptotic
More informationPractical example of an Economic Scenario Generator
Practical example of an Economic Scenario Generator Martin Schenk Actuarial & Insurance Solutions SAV 7 March 2014 Agenda Introduction Deterministic vs. stochastic approach Mathematical model Application
More informationThe ruin probabilities of a multidimensional perturbed risk model
MATHEMATICAL COMMUNICATIONS 231 Math. Commun. 18(2013, 231 239 The ruin probabilities of a multidimensional perturbed risk model Tatjana Slijepčević-Manger 1, 1 Faculty of Civil Engineering, University
More informationContinuous-Time Pension-Fund Modelling
. Continuous-Time Pension-Fund Modelling Andrew J.G. Cairns Department of Actuarial Mathematics and Statistics, Heriot-Watt University, Riccarton, Edinburgh, EH4 4AS, United Kingdom Abstract This paper
More informationMartingale Methods in Financial Modelling
Marek Musiela Marek Rutkowski Martingale Methods in Financial Modelling Second Edition Springer Table of Contents Preface to the First Edition Preface to the Second Edition V VII Part I. Spot and Futures
More informationValuation of performance-dependent options in a Black- Scholes framework
Valuation of performance-dependent options in a Black- Scholes framework Thomas Gerstner, Markus Holtz Institut für Numerische Simulation, Universität Bonn, Germany Ralf Korn Fachbereich Mathematik, TU
More informationM5MF6. Advanced Methods in Derivatives Pricing
Course: Setter: M5MF6 Dr Antoine Jacquier MSc EXAMINATIONS IN MATHEMATICS AND FINANCE DEPARTMENT OF MATHEMATICS April 2016 M5MF6 Advanced Methods in Derivatives Pricing Setter s signature...........................................
More informationELEMENTS OF MONTE CARLO SIMULATION
APPENDIX B ELEMENTS OF MONTE CARLO SIMULATION B. GENERAL CONCEPT The basic idea of Monte Carlo simulation is to create a series of experimental samples using a random number sequence. According to the
More informationVolatility Smiles and Yield Frowns
Volatility Smiles and Yield Frowns Peter Carr NYU CBOE Conference on Derivatives and Volatility, Chicago, Nov. 10, 2017 Peter Carr (NYU) Volatility Smiles and Yield Frowns 11/10/2017 1 / 33 Interest Rates
More informationSlides for Risk Management
Slides for Risk Management Introduction to the modeling of assets Groll Seminar für Finanzökonometrie Prof. Mittnik, PhD Groll (Seminar für Finanzökonometrie) Slides for Risk Management Prof. Mittnik,
More informationResearch Statement. Dapeng Zhan
Research Statement Dapeng Zhan The Schramm-Loewner evolution (SLE), first introduced by Oded Schramm ([12]), is a oneparameter (κ (0, )) family of random non-self-crossing curves, which has received a
More informationPricing Dynamic Solvency Insurance and Investment Fund Protection
Pricing Dynamic Solvency Insurance and Investment Fund Protection Hans U. Gerber and Gérard Pafumi Switzerland Abstract In the first part of the paper the surplus of a company is modelled by a Wiener process.
More informationTangent Lévy Models. Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford.
Tangent Lévy Models Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford June 24, 2010 6th World Congress of the Bachelier Finance Society Sergey
More informationLarge Deviations and Stochastic Volatility with Jumps: Asymptotic Implied Volatility for Affine Models
Large Deviations and Stochastic Volatility with Jumps: TU Berlin with A. Jaquier and A. Mijatović (Imperial College London) SIAM conference on Financial Mathematics, Minneapolis, MN July 10, 2012 Implied
More informationLIBOR models, multi-curve extensions, and the pricing of callable structured derivatives
Weierstrass Institute for Applied Analysis and Stochastics LIBOR models, multi-curve extensions, and the pricing of callable structured derivatives John Schoenmakers 9th Summer School in Mathematical Finance
More informationExtended Libor Models and Their Calibration
Extended Libor Models and Their Calibration Denis Belomestny Weierstraß Institute Berlin Vienna, 16 November 2007 Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Vienna, 16 November
More informationINTEREST RATES AND FX MODELS
INTEREST RATES AND FX MODELS 4. Convexity Andrew Lesniewski Courant Institute of Mathematics New York University New York February 24, 2011 2 Interest Rates & FX Models Contents 1 Convexity corrections
More informationDynamic Relative Valuation
Dynamic Relative Valuation Liuren Wu, Baruch College Joint work with Peter Carr from Morgan Stanley October 15, 2013 Liuren Wu (Baruch) Dynamic Relative Valuation 10/15/2013 1 / 20 The standard approach
More informationClassifying Solvency Capital Requirement Contribution of Collective Investments under Solvency II
Classifying Solvency Capital Requirement Contribution of Collective Investments under Solvency II Working Paper Series 2016-03 (01) SolvencyAnalytics.com March 2016 Classifying Solvency Capital Requirement
More informationVolatility Smiles and Yield Frowns
Volatility Smiles and Yield Frowns Peter Carr NYU IFS, Chengdu, China, July 30, 2018 Peter Carr (NYU) Volatility Smiles and Yield Frowns 7/30/2018 1 / 35 Interest Rates and Volatility Practitioners and
More informationModern Methods of Option Pricing
Modern Methods of Option Pricing Denis Belomestny Weierstraß Institute Berlin Motzen, 14 June 2007 Denis Belomestny (WIAS) Modern Methods of Option Pricing Motzen, 14 June 2007 1 / 30 Overview 1 Introduction
More informationWindow Width Selection for L 2 Adjusted Quantile Regression
Window Width Selection for L 2 Adjusted Quantile Regression Yoonsuh Jung, The Ohio State University Steven N. MacEachern, The Ohio State University Yoonkyung Lee, The Ohio State University Technical Report
More informationMODELLING 1-MONTH EURIBOR INTEREST RATE BY USING DIFFERENTIAL EQUATIONS WITH UNCERTAINTY
Applied Mathematical and Computational Sciences Volume 7, Issue 3, 015, Pages 37-50 015 Mili Publications MODELLING 1-MONTH EURIBOR INTEREST RATE BY USING DIFFERENTIAL EQUATIONS WITH UNCERTAINTY J. C.
More informationNEW AND ROBUST DRIFT APPROXIMATIONS FOR THE LIBOR MARKET MODEL
NEW AND ROBUT DRIFT APPROXIMATION FOR THE LIBOR MARKET MODEL MARK JOHI AND ALAN TACEY Abstract. We present four new methods for approximating the drift in the LIBOR market model. These are compared to
More informationOption Pricing for Discrete Hedging and Non-Gaussian Processes
Option Pricing for Discrete Hedging and Non-Gaussian Processes Kellogg College University of Oxford A thesis submitted in partial fulfillment of the requirements for the MSc in Mathematical Finance November
More informationProject 1: Double Pendulum
Final Projects Introduction to Numerical Analysis II http://www.math.ucsb.edu/ atzberg/winter2009numericalanalysis/index.html Professor: Paul J. Atzberger Due: Friday, March 20th Turn in to TA s Mailbox:
More information3.4 Copula approach for modeling default dependency. Two aspects of modeling the default times of several obligors
3.4 Copula approach for modeling default dependency Two aspects of modeling the default times of several obligors 1. Default dynamics of a single obligor. 2. Model the dependence structure of defaults
More informationEC316a: Advanced Scientific Computation, Fall Discrete time, continuous state dynamic models: solution methods
EC316a: Advanced Scientific Computation, Fall 2003 Notes Section 4 Discrete time, continuous state dynamic models: solution methods We consider now solution methods for discrete time models in which decisions
More informationApplication of an Interval Backward Finite Difference Method for Solving the One-Dimensional Heat Conduction Problem
Application of an Interval Backward Finite Difference Method for Solving the One-Dimensional Heat Conduction Problem Malgorzata A. Jankowska 1, Andrzej Marciniak 2 and Tomasz Hoffmann 2 1 Poznan University
More informationChapter 9 Dynamic Models of Investment
George Alogoskoufis, Dynamic Macroeconomic Theory, 2015 Chapter 9 Dynamic Models of Investment In this chapter we present the main neoclassical model of investment, under convex adjustment costs. This
More informationLOGNORMAL MIXTURE SMILE CONSISTENT OPTION PRICING
LOGNORMAL MIXTURE SMILE CONSISTENT OPTION PRICING FABIO MERCURIO BANCA IMI, MILAN http://www.fabiomercurio.it Daiwa International Workshop on Financial Engineering, Tokyo, 26-27 August 2004 1 Stylized
More informationComputer Exercise 2 Simulation
Lund University with Lund Institute of Technology Valuation of Derivative Assets Centre for Mathematical Sciences, Mathematical Statistics Fall 2017 Computer Exercise 2 Simulation This lab deals with pricing
More informationApproximation Methods in Derivatives Pricing
Approximation Methods in Derivatives Pricing Minqiang Li Bloomberg LP September 24, 2013 1 / 27 Outline of the talk A brief overview of approximation methods Timer option price approximation Perpetual
More informationSmoking Adjoints: fast evaluation of Greeks in Monte Carlo calculations
Report no. 05/15 Smoking Adjoints: fast evaluation of Greeks in Monte Carlo calculations Michael Giles Oxford University Computing Laboratory, Parks Road, Oxford, U.K. Paul Glasserman Columbia Business
More informationThe value of foresight
Philip Ernst Department of Statistics, Rice University Support from NSF-DMS-1811936 (co-pi F. Viens) and ONR-N00014-18-1-2192 gratefully acknowledged. IMA Financial and Economic Applications June 11, 2018
More informationLecture 3: Review of mathematical finance and derivative pricing models
Lecture 3: Review of mathematical finance and derivative pricing models Xiaoguang Wang STAT 598W January 21th, 2014 (STAT 598W) Lecture 3 1 / 51 Outline 1 Some model independent definitions and principals
More informationChapter 3: Black-Scholes Equation and Its Numerical Evaluation
Chapter 3: Black-Scholes Equation and Its Numerical Evaluation 3.1 Itô Integral 3.1.1 Convergence in the Mean and Stieltjes Integral Definition 3.1 (Convergence in the Mean) A sequence {X n } n ln of random
More informationPolicy iterated lower bounds and linear MC upper bounds for Bermudan style derivatives
Finance Winterschool 2007, Lunteren NL Policy iterated lower bounds and linear MC upper bounds for Bermudan style derivatives Pricing complex structured products Mohrenstr 39 10117 Berlin schoenma@wias-berlin.de
More informationHedging Under Jump Diffusions with Transaction Costs. Peter Forsyth, Shannon Kennedy, Ken Vetzal University of Waterloo
Hedging Under Jump Diffusions with Transaction Costs Peter Forsyth, Shannon Kennedy, Ken Vetzal University of Waterloo Computational Finance Workshop, Shanghai, July 4, 2008 Overview Overview Single factor
More informationThe Black Model and the Pricing of Options on Assets, Futures and Interest Rates. Richard Stapleton, Guenter Franke
The Black Model and the Pricing of Options on Assets, Futures and Interest Rates Richard Stapleton, Guenter Franke September 23, 2005 Abstract The Black Model and the Pricing of Options We establish a
More informationarxiv: v1 [q-fin.pr] 18 Feb 2010
CONVERGENCE OF HESTON TO SVI JIM GATHERAL AND ANTOINE JACQUIER arxiv:1002.3633v1 [q-fin.pr] 18 Feb 2010 Abstract. In this short note, we prove by an appropriate change of variables that the SVI implied
More informationEuropean call option with inflation-linked strike
Mathematical Statistics Stockholm University European call option with inflation-linked strike Ola Hammarlid Research Report 2010:2 ISSN 1650-0377 Postal address: Mathematical Statistics Dept. of Mathematics
More information1 Mathematics in a Pill 1.1 PROBABILITY SPACE AND RANDOM VARIABLES. A probability triple P consists of the following components:
1 Mathematics in a Pill The purpose of this chapter is to give a brief outline of the probability theory underlying the mathematics inside the book, and to introduce necessary notation and conventions
More informationA Continuity Correction under Jump-Diffusion Models with Applications in Finance
A Continuity Correction under Jump-Diffusion Models with Applications in Finance Cheng-Der Fuh 1, Sheng-Feng Luo 2 and Ju-Fang Yen 3 1 Institute of Statistical Science, Academia Sinica, and Graduate Institute
More informationANALYSIS OF THE BINOMIAL METHOD
ANALYSIS OF THE BINOMIAL METHOD School of Mathematics 2013 OUTLINE 1 CONVERGENCE AND ERRORS OUTLINE 1 CONVERGENCE AND ERRORS 2 EXOTIC OPTIONS American Options Computational Effort OUTLINE 1 CONVERGENCE
More informationComputational Finance Binomial Trees Analysis
Computational Finance Binomial Trees Analysis School of Mathematics 2018 Review - Binomial Trees Developed a multistep binomial lattice which will approximate the value of a European option Extended the
More informationCalibration of Interest Rates
WDS'12 Proceedings of Contributed Papers, Part I, 25 30, 2012. ISBN 978-80-7378-224-5 MATFYZPRESS Calibration of Interest Rates J. Černý Charles University, Faculty of Mathematics and Physics, Prague,
More informationDRAFT. 1 exercise in state (S, t), π(s, t) = 0 do not exercise in state (S, t) Review of the Risk Neutral Stock Dynamics
Chapter 12 American Put Option Recall that the American option has strike K and maturity T and gives the holder the right to exercise at any time in [0, T ]. The American option is not straightforward
More informationA SIMPLE DERIVATION OF AND IMPROVEMENTS TO JAMSHIDIAN S AND ROGERS UPPER BOUND METHODS FOR BERMUDAN OPTIONS
A SIMPLE DERIVATION OF AND IMPROVEMENTS TO JAMSHIDIAN S AND ROGERS UPPER BOUND METHODS FOR BERMUDAN OPTIONS MARK S. JOSHI Abstract. The additive method for upper bounds for Bermudan options is rephrased
More informationOptimal Security Liquidation Algorithms
Optimal Security Liquidation Algorithms Sergiy Butenko Department of Industrial Engineering, Texas A&M University, College Station, TX 77843-3131, USA Alexander Golodnikov Glushkov Institute of Cybernetics,
More informationDynamic Replication of Non-Maturing Assets and Liabilities
Dynamic Replication of Non-Maturing Assets and Liabilities Michael Schürle Institute for Operations Research and Computational Finance, University of St. Gallen, Bodanstr. 6, CH-9000 St. Gallen, Switzerland
More informationAll Investors are Risk-averse Expected Utility Maximizers
All Investors are Risk-averse Expected Utility Maximizers Carole Bernard (UW), Jit Seng Chen (GGY) and Steven Vanduffel (Vrije Universiteit Brussel) AFFI, Lyon, May 2013. Carole Bernard All Investors are
More informationCONSISTENCY AMONG TRADING DESKS
CONSISTENCY AMONG TRADING DESKS David Heath 1 and Hyejin Ku 2 1 Department of Mathematical Sciences, Carnegie Mellon University, Pittsburgh, PA, USA, email:heath@andrew.cmu.edu 2 Department of Mathematics
More informationContinuous-time Stochastic Control and Optimization with Financial Applications
Huyen Pham Continuous-time Stochastic Control and Optimization with Financial Applications 4y Springer Some elements of stochastic analysis 1 1.1 Stochastic processes 1 1.1.1 Filtration and processes 1
More informationI. Time Series and Stochastic Processes
I. Time Series and Stochastic Processes Purpose of this Module Introduce time series analysis as a method for understanding real-world dynamic phenomena Define different types of time series Explain the
More informationStrategies for Improving the Efficiency of Monte-Carlo Methods
Strategies for Improving the Efficiency of Monte-Carlo Methods Paul J. Atzberger General comments or corrections should be sent to: paulatz@cims.nyu.edu Introduction The Monte-Carlo method is a useful
More informationSmooth estimation of yield curves by Laguerre functions
Smooth estimation of yield curves by Laguerre functions A.S. Hurn 1, K.A. Lindsay 2 and V. Pavlov 1 1 School of Economics and Finance, Queensland University of Technology 2 Department of Mathematics, University
More informationsymmys.com 3.2 Projection of the invariants to the investment horizon
122 3 Modeling the market In the swaption world the underlying rate (3.57) has a bounded range and thus it does not display the explosive pattern typical of a stock price. Therefore the swaption prices
More informationValue at Risk and Self Similarity
Value at Risk and Self Similarity by Olaf Menkens School of Mathematical Sciences Dublin City University (DCU) St. Andrews, March 17 th, 2009 Value at Risk and Self Similarity 1 1 Introduction The concept
More informationMartingale Pricing Theory in Discrete-Time and Discrete-Space Models
IEOR E4707: Foundations of Financial Engineering c 206 by Martin Haugh Martingale Pricing Theory in Discrete-Time and Discrete-Space Models These notes develop the theory of martingale pricing in a discrete-time,
More informationLinear-Rational Term-Structure Models
Linear-Rational Term-Structure Models Anders Trolle (joint with Damir Filipović and Martin Larsson) Ecole Polytechnique Fédérale de Lausanne Swiss Finance Institute AMaMeF and Swissquote Conference, September
More informationEFFECT OF IMPLEMENTATION TIME ON REAL OPTIONS VALUATION. Mehmet Aktan
Proceedings of the 2002 Winter Simulation Conference E. Yücesan, C.-H. Chen, J. L. Snowdon, and J. M. Charnes, eds. EFFECT OF IMPLEMENTATION TIME ON REAL OPTIONS VALUATION Harriet Black Nembhard Leyuan
More informationQuantitative relations between risk, return and firm size
March 2009 EPL, 85 (2009) 50003 doi: 10.1209/0295-5075/85/50003 www.epljournal.org Quantitative relations between risk, return and firm size B. Podobnik 1,2,3(a),D.Horvatic 4,A.M.Petersen 1 and H. E. Stanley
More information