Incorporating boundary conditions in a stochastic volatility model for the numerical approximation of bond prices

In this paper, we consider a two-factor interest rate model with stochastic volatility, and we assume that the instantaneous interest rate follows a jump-diffusion process. In this kind of problems, a two-dimensional partial integro-differential equation is derived for the values of zero-coupon bonds. To apply standard numerical methods to this equation, it is customary to consider a bounded domain and incorporate suitable boundary conditions. However, for these two-dimensional interest rate models, there are not well-known boundary conditions, in general. Here, in order to approximate bond prices, we propose new boundary conditions, which maintain the discount function property of the zero-coupon bond price. Then, we illustrate the numerical approximation of the corresponding boundary value problem by means of an alternative direction implicit method, which has been already applied for pricing options. We test these boundary conditions with several interest rate pricing models.     

Stochastic models for bond prices,function space integrals and immunization theory

This article presents several diffusion models for bond prices. By considering the spot interest rate as a state variable and invoking the no-arbitrage principle, the price of a default-free and non-callable pure discount bond is represented as a conditional expectation. The Ornstein—Uhlenbeck (O.U.) stochastic process is described, and used to model the spot rate. The O.U. process is then modified to exclude negative interest rates and the resulting bond-price partial differential equation is solved. By considering the yield rate as a state variable and using the Brownian bridge process, a simpler bond price model is obtained. Applications to immunization theory are presented.     

Repeated Red-Black ordering: a new approach

Hereafter, we describe and analyze, from both a theoretical and a numerical point of view, an iterative method for efficiently solving symmetric elliptic problems with possibly discontinuous coefficients. In the following, we use the Preconditioned Conjugate Gradient method to solve the symmetric positive definite linear systems which arise from the finite element discretization of the problems. We focus our interest on sparse and efficient preconditioners. In order to define the preconditioners, we perform two steps: first we reorder the unknowns and then we carry out a (modified) incomplete factorization of the original matrix. We study numerically and theoretically two preconditioners, the second preconditioner corresponding to the one investigated by Brand and Heinemann [2]. We prove convergence results about the Poisson equation with either Dirichlet or periodic boundary conditions. For a meshsizeh, Brand proved that the condition number of the preconditioned system is bounded byO(h ^–1/2) for Dirichlet boundary conditions. By slightly modifying the preconditioning process, we prove that the condition number is bounded byO(h ^–1/3).     

The transformation theorem for two-parameter pure jump martingales

Summary LetM be a martingale of pure jump type, i.e. the compensation of the process describing the total of the point jumps ofM in the plane.M can be uniformly approximated by martingales of bounded variation jumping only on finitely many axial parallel lines. Using this fact we prove a change of variables formula in which forC ⁴-functions f the processf(M) is described by integrals off ^(k) (M),k=1, 2, with respect to stochastic integrators of the types expected: a martingale, two processes behaving as martingales in one direction and as processes of bounded variation in the other, and one process of bounded variation. Hereby we are led to investigate two types of random measures not considered so far in this context. By combination with the integrators already known, they might complete the set needed for a general transformation formula.     

Լ��ģ���к��н��׶������÷��յĿ�ת��ծȯ�Ķ���

??This paper studies the price of convertible bonds with counterparty credit risk in a reduced-form model. We suppose that the default intensity process and the interest rate process follow the Vasicek model, and derive the price expression of convertible bonds using the method of measure changes. Moreover, we make some numerical analysis on the explicit formulae to demonstrate the sensitivity of a convertible bond price to changes in the parameters of the model.     

On a multiple credit rating migration model with stochastic interest rate

In this paper, we explore a pricing model for corporate bond accompanied with multiple credit rating migration risk and stochastic interest rate. The bond price volatility strongly depends on potentially multiple credit rating migration and stochastic change of interest rate. A free boundary problem of partial differential equation is presented, which is the equivalent transformation of the pricing model. The existence, uniqueness, and regularity for the free boundary problem are established to guarantee the rationality of the pricing model. Due to the stochastic change of interest rate, the discontinuous coefficient in the free boundary problem depends explicitly on the time variable but is convergent as time tends to infinity. Accordingly, an auxiliary free boundary problem is constructed, whose coefficient is the convergent limit of the coefficient in the original free boundary problem. With some constraint on the risk discount rate satisfied, we prove that a unique traveling wave exists in the auxiliary free boundary problem. The inductive method is adopted to fit the multiplicity of credit rating. Then we show that the solution of the original free boundary problem converges to the traveling wave in the auxiliary free boundary problem. Returning to the pricing model with multiple credit rating migration and stochastic interest rate, we conclude that the bond price profile can be captured by a traveling wave pattern coupling with a guaranteed bond price with face value equal to one at the maturity.     

The critical probability for random Voronoi percolation in the plane is 1/2

We study percolation in the following random environment: let Z be a Poisson process of constant intensity on ℝ², and form the Voronoi tessellation of ℝ² with respect to Z. Colour each Voronoi cell black with probability p, independently of the other cells. We show that the critical probability is 1/2. More precisely, if p>1/2 then the union of the black cells contains an infinite component with probability 1, while if p<1/2 then the distribution of the size of the component of black cells containing a given point decays exponentially. These results are analogous to Kesten's results for bond percolation in ℤ². The result corresponding to Harris' Theorem for bond percolation in ℤ² is known: Zvavitch noted that one of the many proofs of this result can easily be adapted to the random Voronoi setting. For Kesten's results, none of the existing proofs seems to adapt. The methods used here also give a new and very simple proof of Kesten's Theorem for ℤ²; we hope they will be applicable in other contexts as well. Research supported in part by NSF grant ITR 0225610 and DARPA grant F33615-01-C-1900 Research partially undertaken during a visit to the Forschungsinstitut für Mathematik, ETH Zürich, Switzerland     

Hyperuniform and rigid stable matchings

We study a stable partial matching τ of the d‐dimensional lattice with a stationary determinantal point process Ψ on R^d with intensity α>1. For instance, Ψ might be a Poisson process. The matched points from Ψ form a stationary and ergodic (under lattice shifts) point process Ψ^τ with intensity 1 that very much resembles Ψ for α close to 1. On the other hand Ψ^τ is hyperuniform and number rigid, quite in contrast to a Poisson process. We deduce these properties by proving more general results for a stationary point process Ψ, whose so‐called matching flower (a stopping set determining the matching partner of a lattice point) has a certain subexponential tail behavior. For hyperuniformity, we also additionally need to assume some mixing condition on Ψ. Furthermore, if Ψ is a Poisson process then Ψ^τ has an exponentially decreasing truncated pair correlation function.     

Numerical Evaluation of Dynamic Behavior of Ornstein–Uhlenbeck Processes Modified by Various Boundaries and its Application to Pricing Barrier Options

In financial engineering, one often encounters barrier options in which an action promised in the contract is taken if the underlying asset value becomes too high or too low. In order to compute the corresponding prices, it is necessary to capture the dynamic behavior of the associated stochastic process modified by boundaries. To the best knowledge of the authors, there is no algorithmic approach available to compute such prices repeatedly in a systematic manner. The purpose of this paper is to develop computational algorithms to capture the dynamic behavior of Ornstein-Uhlenbeck processes modified by various boundaries based on the Ehrenfest approximation approach established in Sumita et al. (J Oper Res Soc Jpn 49:256–278, 2006). As an application, we evaluate the prices of up-and-out call options maturing at time τ _M with strike price K _S written on a discount bond maturing at time T, demonstrating the usefulness, speed and accuracy of the proposed computational algorithms.     

Boundary Behavior of Harmonic Functions for Truncated Stable Processes

For any α∈(0,2), a truncated symmetric α-stable process in ℝ ^d is a symmetric Lévy process in ℝ ^d with no diffusion part and with a Lévy density given by c|x|^−d−α 1_{|x|<1} for some constant c. In (Kim and Song in Math. Z. 256(1): 139–173, [2007]) we have studied the potential theory of truncated symmetric stable processes. Among other things, we proved that the boundary Harnack principle is valid for the positive harmonic functions of this process in any bounded convex domain and showed that the Martin boundary of any bounded convex domain with respect to this process is the same as the Euclidean boundary. However, for truncated symmetric stable processes, the boundary Harnack principle is not valid in non-convex domains. In this paper, we show that, for a large class of not necessarily convex bounded open sets in ℝ ^d called bounded roughly connected κ-fat open sets (including bounded non-convex κ-fat domains), the Martin boundary with respect to any truncated symmetric stable process is still the same as the Euclidean boundary. We also show that, for truncated symmetric stable processes a relative Fatou type theorem is true in bounded roughly connected κ-fat open sets. The research of P. Kim is supported by the Korea Research Foundation Grant funded by the Korean Government (MOEHRD, Basic Research Promotion Fund) (KRF-2007-331-C00037). The research of R. Song is supported in part by a joint US-Croatia grant INT 0302167.     

Controlled diffusion processes on infinite horizon with the overtaking criterion

The optimal control of diffusion processes on the infinite time interval are studied. All the costs diverge to infinity and we employ the overtaking criterion, as well as considering minimal growth rate controls. The Bellman equation for the problem is considered and its solution is given an interpretation connected with the overtaking optimality notion. Control problems with a cost including a generalized discount factor (which is not integrable on [(, )) are also studied. Both cases, where the diffusion is inR ⁿ or where it is reflected from the boundary of a bounded set, are considered.This research was supported in part by the Institute for Mathematics and its Applications with funds provided by the National Science Foundation.     

Vasicek利率模型下带有零息票债券的投资-消费模型

应用随机最优控制理论研究Vasicek利率模型下的投资-消费问题,其中假设无风险利率是服从Vasicek利率模型的随机过程,且与股票价格过程存在一般相关性.假设金融市场由一种无风险资产、一种风险资产和一种零息票债券所构成,投资者的目标是最大化中期消费与终端财富的期望贴现效用.应用变量替换方法得到了幂效用下最优投资-消费策略的显示表达式,并分析了最优投资-消费策略对市场参数的灵敏度.     

Controllability and Hedgibility of Black-Scholes Equations with <Emphasis Type="Italic">N</Emphasis> Stocks

This paper is to investigate the controllability and observability properties of linear and certain nonlinear Black-Scholes (B-S) type equations consisting of N stocks in an appropriate bounded domain I of ℝ₊ ^N . In this model both the stock volatility and interest rate are influenced by the stock prices and the control which is related to the hedging ratio in option pricing of finance is distributed over a subdomain of I. The proof of the controllability result for the linear B-S equations relies on the suitable observability inequality for the associated adjoint problem, and for the nonlinear model, fixed point technique is applied. Our result leads to that the dynamic hedgibility in finance is proved in the context of controllability theory.     

Pricing Fixed-Income Securities in an Information-Based Framework

Abstract

The purpose of this article is to introduce a class of information-based models for the pricing of fixed-income securities. We consider a set of continuous-time processes that describe the flow of information concerning market factors in a monetary economy. The nominal pricing kernel is assumed to be given at any specified time by a function of the values of information processes at that time. Using a change-of-measure technique, we derive explicit expressions for the prices of nominal discount bonds and deduce the associated dynamics of the short rate of interest and the market price of risk. The interest rate positivity condition is expressed as a differential inequality. An example that shows how the model can be calibrated to an arbitrary initial yield curve is presented. We proceed to model the price level, which is also taken at any specified time to be given by a function of the values of the information processes at that time. A simple model for a stochastic monetary economy is introduced in which the prices of the nominal discount bonds and inflation-linked notes can be expressed in terms of aggregate consumption and the liquidity benefit generated by the money supply.     

Interest rate swap pricing with default risk under variance gamma process

Under the assumption that the dynamic assets price follows the variance gamma process, we establish a new bilateral pricing model of interest rate swap by integrating the reduced form model for swap pricing and the structural model for default risk measurement. Our pricing model preserves the simplicity of the reduced form model and also considers the dynamic evolution of the counterparty assets price by incorporating with the structural model for default risk measurement. We divide the swap pricing framework into two parts, simplifying the pricing model relatively. Simulation results show that, for a one year interest rate swap, a bond spread of one hundred basis points implies a swap credit spread about 0.1054 basis point.     

The Distribution of Condition Numbers of Rational Data of Bounded Bit Length

We prove that rational data of bounded input length are uniformly distributed (in the classical sense of H. Weyl, in [42]) with respect to the probability distribution of condition numbers of numerical analysis. We deal both with condition numbers of linear algebra and with condition numbers for systems of multivariate polynomial equations. For instance, we prove that for a randomly chosen n\times n rational matrix M of bit length O(n ⁴ log n) + log w , the condition number k(M) satisfies k(M) ≤ w n ^5/2 with probability at least 1-2w ^-1 . Similar estimates are established for the condition number μ_ norm of M. Shub and S. Smale when applied to systems of multivariate homogeneous polynomial equations of bounded input length. Finally, we apply these techniques to estimate the probability distribution of the precision (number of bits of the denominator) required to write approximate zeros of systems of multivariate polynomial equations of bounded input length. March 7, 2001. Final version received: June 7, 2001.     

Gibbsianness of fermion random point fields

We consider fermion (or determinantal) random point fields on Euclidean space ℝ^d. Given a bounded, translation invariant, and positive definite integral operator J on L²(ℝ^d), we introduce a determinantal interaction for a system of particles moving on ℝ^d as follows: the n points located at x₁,· · ·,x_n ∈ ℝ^d have the potential energy given by where j(x−y) is the integral kernel function of the operator J. We show that the Gibbsian specification for this interaction is well-defined. When J is of finite range in addition, and for d≥2 if the intensity is small enough, we show that the fermion random point field corresponding to the operator J(I+J)⁻¹ is a Gibbs measure admitted to the specification.     

Explicit results for a class of asset-selling problems

We consider the following optimal selection problem: There are n identical assets which are to be sold, one at a time, to coming bidders. The bids are i.i.d. where there are only two possible bid-values, with known probabilities. The stream of bidders constitutes a general renewal process, and rewards are continuously discounted at a constant rate. The objective is to maximize the total expected discounted revenue from the sale of the n assets. The optimal policy here is stationary, where the decision in question is only whether to accept a low bid or not; the answer is affirmative depending on a critical number n^* of remaining assets. In this paper we derive an explicit formula for n^*, being a function of the Laplace transform of the renewal distribution evaluated at the discount rate, the probability for a low bid, and the ratio between the two bid-values. We also specify the pertinent value functions. Applications of the model are discussed in detail, and extensions are made to include holding costs and to allow for optimal pricing.