arxiv: v1 [math-ph] 4 Feb 2008

Size: px
Start display at page:

Download "arxiv: v1 [math-ph] 4 Feb 2008"

Transcription

1 On the local energy decay of solutions of the Dirac equation in the non-extreme Kerr-Newman metric, I arxiv:8.4v1 [math-ph] 4 Feb 8 M. Winklmeier, O. Yamada November 8, 18 Abstract 375.ptWe investigate the local energy decay of solutions of the Dirac equation in the non-extreme Kerr-Newman metric. First, we write the Dirac equation as a Cauchy problem and define the Dirac operator. It is shown that the Dirac operator is selfadjoint in a suitable Hilbert space. With the RAGE theorem, we show that for each particle its energy located in any compact region outside of the event horizon of the Kerr-Newman black hole decays in the time mean. 1 Introduction The Kerr-Newman metric is the most general stationary solution of Einstein s field equations and has the physical interpretation of a massive charged rotating black hole. The physical parameters M, Q and a = J/M have the interpretation as mass, electric charge and angular momentum parameter. As in flat spacetime, a spin- 1 particle is described by a Dirac equation. In the Kerr- Newman metric, the Dirac equation is given by the coupled system of partial differential equations see, e.g., [Pag76], [Cha98] where R := R+Â Ψ = 1 imr R t,ϕ + imr R t,ϕ R t,ϕ imr, R t,ϕ + imr  := D L t,ϕ + D L t,ϕ L t,ϕ + D D L t,ϕ winklmeier@math.unibe.ch yamadaos@se.ritsumei.ac.jp 1

2 and D := amcosθ, 3 R t,ϕ ± := r ± i [ r +a i t +ai ] ϕ +eqr on r +,, 4 L t,ϕ ± := θ + cotθ [ asinθ i t + 1 sinθ i ] on,π. 5 ϕ The physical parameters m and e are the mass and the electric charge of the spin- 1 - particle which is described by the wave function Ψ. The functions and Σ are given by r = r Mr+a +Q = r M +a +Q M, Σr,θ = r +a cos θ. According to the black hole parameters M, Q and a, three cases can arise: If a + Q M <, then the function has exactly two distinct zeros; this case is called the non-extreme Kerr-Newman metric. If a + Q M =, then the function has exactly one zero; this case is referred to as the extreme Kerr-Newman metric. In both cases the metric is interpreted physically as the spacetime generated by a charged rotating massive black hole. If a +Q M >, then the function has no zeros; in this case the Kerr-Newman metric has no interpretation as the metric of a black hole. It is a remarkable fact that the Dirac equation 1 can be separated into ordinary differential equations. This was first shown by Chandrasekhar [Cha76].The first step is to separate the dependence on the coordinates t and ϕ by the ansatz Ψr,θ,ϕ,t = e iωt Ψr,θ,ϕ, 6 ω R. We call ω an energy eigenvalue of equation 1 if it has a solution of the form 6 which is square integrable in the sense explained below in section. The separation of the Dirac equation into ordinary differential equations is achieved by the ansatz X rs + θ Ψr,θ,ϕ,t = e iωt e iκϕ X + rs θ X + rs + θ X rs θ with κ Z+ 1. One obtains two differential equations, the so-calledradial equation for the radial coordinate r and the so-called angular equation for the angular coordinate θ. Both equations can be written as eigenvalue problems for selfadjoint operators in appropriate Hilbert spaces. The spectrum of the corresponding angular operator consists of simple and discrete eigenvalues which are unbounded from below and above see, e.g., [Win6]. Belgiorno and Martellini [BM99] showed that the spectrum of the radial equation comprises all of the real axis. Recently, it was shown that in the non-extreme Kerr-Newman metric no embedded eigenvalues exist. In the extreme Kerr-Newman case, however, embedded eigenvalues can existsee[sch4] and[wy6]. In this paper we investigate the local energy decay of solutions of the Dirac equation1. We assume that the non-extreme case holds, that is, there are no energy eigenvalues.

3 We do not employ the separation ansatz described above, but we rewrite the Dirac equation as a Cauchy problem. After some transformations, we obtain the equation i t Ψx,θ,ϕ,t = Sx,θ 1 HΨx,θ,ϕ,t 7 where Sx,θ is a bounded and boundedly invertible 4 4-matrix and H is the formal Dirac operator associated with the Dirac equation 1. We show that the expression H on the right hand side has a selfadjoint realisation H in an appropriate Hilbert space Theorem.1. The operator H can be written as the sum over partial wave operators H κ, κ Z. In order to apply the RAGE theorem to the Dirac operator, we show in Theorem 3.3 that it enjoys a Rellich type property. One of our main results in this paper is Theorem 4. which shows the local energy decay of the wave functions in the time mean: Let U κ be the group associated with the skew-selfadjoint operator is 1 H κ. Then, any function f which is subjected to the time evolution given by 7 satisfies lim T [ 1 T R T T R ] U κ tfx,θ dxdθ dt =. C 4 for every R >. Since the expression in the square brackets is related to the energy of the particle described by f at time t in the region Ω := R,R,π, this result is related to the decay of the particle s energy in Ω. The above result remains valid eqr if the bounded electric potential r +a given by the background metric is substituted by a possibly unbounded potential q satisfying q C 1 qr r +,, lim r r exists, q r = Or, as r. The local energy decay has also been investigated by Finster, Kamran, Smoller and Yau [FKSY3]. They restrict the Dirac operator to an annulus outside the event horizon and impose Dirichlet type boundary conditions so that the restricted operator has purely discrete spectrum. Hence the propagator associated with it is a sum of projections. Extending the radii of the annulus to the event horizon and infinity, respectively, the sum of the propagator representation becomes an integral. Our approach, however, deals directly with the Dirac operator in the exterior of the black hole. Using the fact that the partial wave operators H κ have no point spectrum, the local energy decay of wave functions with initially compact support follows by methods from scattering theory. In particular, we do not impose any additional boundary conditions. Transformation of the Dirac equation to a Cauchy problem To investigate the time evolution of solutions of the time-dependent Dirac equation we rewrite it as a Cauchy problem. In the following, we use the Pauli matrices to abbreviate notation: 1 σ 1 =, σ 1 = i i, σ 3 =

4 In addition, let 1 I σj I =, β =, Σ 1 I j =, j = 1,,3. σ j For the ϕ-coordinate we always assume the boundary condition Ψr,θ,ϕ = Ψr,θ,ϕ+ π. In the rest of the paper we use the notation Ω 3 :=,,π [,π, Ω :=,,π. Theorem.1. Assume that the function has at least one real zero that is, either the extreme or the non-extreme Kerr-Newman case holds and denote the largest zero by r +. Define the new radial coordinate x, by dx dr = r +a r. 8 To keep notation simple, we often write r instead rx. A function Ψ is a solution of the Dirac equation 1 if and only if i t Ψ = S 1 HΨ 9 where Ψx,θ,ϕ,t = sinθ 1/ Ψrx,θ,ϕ,t, 1 S = I + a sinθ r +a σ, σ H = H 1 +H, h1 H 1 =, h h1 1 = σ 3 i r r +a σ 1 i θ +σ 1 sinθ i, ϕ H = 1 I r +a ai I ϕ +eqr mr I r +a + am cosθ ii I r +a. ii Moreover, the operator H defined by H Ψ = HΨ, DH = C Ω is essentially selfadjoint in the Hilbert space L Ω 3 4 and the operator S 1 H is essentially selfadjoint in the Hilbert space L S Ω 3 4 := L Ω 3 4 given by Ψ, Φ S = Ψ, SΦ = R with scalar product, S accordingly, the norm on L S Ω 3 will be denoted by S. Ψx,θ,ϕ, Sx,θΦx,θ,ϕ C 4dxdθdϕ; 4

5 Let H be the closure of H. Then S 1 H is the closure of S 1 H, and we call S 1 H the time-independent Dirac operator in the Kerr-Newman metric. Remark.. Observe that for each x,θ,,π the matrix Sx,θ is bounded and boundedly invertible since a sinθ rx +a a rx rx +a 1. Therefore, the norms on L Ω 3 4 and L S Ω 3 4 are equivalent. Proof of Theorem.1. First we rearrange equation 1 such that all time derivatives are on the left hand side and all other terms are on the right hand side, thus we obtain [ iσ 3 + a sinθ σ ] 1 iσ 3 r +a i = σ 1 t Ψx,θ,ϕ,t [ I + I x + r +a iσ iσ θ σ 1 r +a i sinθ σ 1 ϕ + 1 r +a iσ 3 ai iσ 3 ϕ +eqr + imr σ3 r +a + am ] cosθ σ3 Ψx,θ,ϕ,t. σ3 r +a σ3 Multiplication from the left by iσ 3 iσ 3 yields [ I + a ] sinθ σ r +a i σ t Ψx,θ,ϕ,t [ σ3 = i σ 3 x + σ1 r +a i σ 1 θ σ + r +a i sinθ σ ϕ 1 I r +a ai I ϕ +eqr mr I r +a + am cosθ I r +a = H 1 +H Ψx,θ,ϕ,t, ii ] ii Ψx,θ,ϕ,t thus the first assertion of the theorem is proved. By Remark. the norms on L Ω 3 4 and LS Ω 3 4 are equivalent. Hence, H is essentially selfadjoint in L Ω 3 4 if and only if S 1 H is so in LS Ω 3 4. Therefore, it suffices to show that H is essentially selfadjoint. Since ϕ with the boundary condition stated at the beginning of this section has the complete system of eigenfunctions e iκϕ κ Z+1/ and all the differential expressions 5

6 in the equation above commute with ϕ, each Ψ L Ω 3 R has a representation Ψx,θ,ϕ,t = κ Z+1/ Ψ satisfies 9 if and only if each Ψ κ satisfies e iκϕ Ψ κ x,θ,t. i t Ψκ = S 1 H κ Ψ κ 1 where H κ, H κ 1 and H κ are obtained from H, H 1 and H by substituting i ϕ with κ. Moreover, H is essentially selfadjoint if and only if each partial wave operator H κ the restriction of H to the partial wave space represented by e iκϕ is essentially selfadjoint. In Proposition.4 below, it is shown that h κ 1,, defined by h κ 1, Ψκ = h κ 1 Ψκ, Dh κ 1, = C Ω, 13 is essentially selfadjoint in L Ω with the usual scalar product. Since both H 1 and S 1 have diagonal structure and H κ defines a bounded selfadjoint operator on LS Ω 4, the assertion on the essential selfadjointness of H κ, κ Z, and therefore of H follows. Lemma.3. For κ Z+1/ let 1 A κ := 1 θ + κ sinθ 1 1 and define the unitary matrix W := 1. i The operator A κ is the angular part of the Dirac equation 1 arising from Chandrasekhar s separation process in the case a =. i For each κ Z + 1/, the operator A κ with domain C,π is essentially selfadjoint in L,π, dθ. We denote its closure again by A κ. It is compactly invertible and its spectrum consists of simple eigenvalues only, given by λ κ m := signm κ 1 + m, m = ±1, ±,. We denote the corresponding normalized eigenfunctions by g κ m. ii For every m Z\{} we have λ κ m = λ κ m and σ 3 g κ m = g κ m. iii WA κ W 1 = W 1 A κ W = iσ 3 A κ. iv W 1 h κ 1 W = W 1 σ 3 W i x r +a A κ = σ 3 i x r +a A κ. 6

7 Proof. For i, we refer to [Win6]. Statements iii and iv are easily verified by direct computation. To prove ii we observe that all eigenvalues of A κ are simple and that 1 A κ σ 3 g κ m = σ 3σ 3 A κ σ 3 g κ m = 1 1 = 1 1 A κ g 1 m κ 1 A 1 κ gm κ 1 = λκ m g 1 m κ = λκ m σ 3gm κ. = λκ m σ 3gm κ. Proposition.4. The minimal operator h κ 1, defined in 13 by h κ 1, ψ = hκ 1 ψ = σ 3 i r r +a σ 1 i θ +σ 1 sinθ i ψ, ϕ Dh κ 1, = C,,π is essentially selfadjoint in the space L Ω 4. Proof. Obviously, h κ 1 is hermitian and h κ 1, is symmetric. Hence, the lemma is proved if we have shown that kerh κ 1, ±i = {}. The equality hκ 1, ψ = ±iψ is equivalent to [ W 1 h κ 1, W ]W 1 ψ = ±iw 1 ψ with the matrix W as defined in Lemma.3. If we set ϕ := W 1 ψ, then, by Lemma.3, the equation above is equivalent to [σ 3 i ] + x r +a A κ ϕ = ±iϕx,θ. 14 Since the set of eigenfunctions of the operator A κ, denoted by g κ n as in Lemma.3, is complete in L,π, dθ, we can expand the function ϕ in 14 as ϕx,θ = m Z\{} ξ m xg κ mθ with functions ξ m L,,dx. Thus 14 yields ±i m Z\{} ξ m xg m θ = m Z\{} = m Z\{} i x σ 3ξ m xg κ mθ+ m Z\{} i x ξ mxg m θ+ λκ m r +a ξ mxgmθ. κ r +a ξ mxa κ g κ mθ 7

8 Taking the scalar product with g µ θ in C and integrating with respect to θ yields on the left hand side ±i m Z\{} and on the right hand side m Z\{} = Hence we obtain ξ m g κ m θg µ θ C dθ = ±i = ±iξ µ m Z\{} ξ m g κ m θg µ θ dθ C i x ξ mxg m θ+ λκ m r +a ξ mxgmθ, κ g µ θ dθ C m Z\{} + m Z\{} i x ξ mx λ κ m r +a ξ mx g m θ, g µ θ C dθ = i x ξ mx+ λκ m r +a ξ mx. g m θ, g µ θ C dθ g m θ, g µ θ C dθ ±iξ µ x = i x ξ µx+ λκ µ r +a ξ µx. 15 An analogous procedure with g µ gives ±iξ µ x = i x ξ µx+ λκ µ r +a ξ µx. 16 Combining equations 15 and 16 we obtain ±i ξ µ = ξ µ Themultiplicationoperator λκ µ that the operator i d dx i d dx λ κ µ r +a i d dx i d dx λκ µ r +a ξ µ. 17 ξ µ r +a isboundedon L,,dxanditiswellknown with domain C, is essentially selfadjoint in L,,dx, see, for instance, [Wei87, Ch.6.B]. Hence the same is true for λ κ m r +a d i dx and equation 17 has only the trivial solution in L,,dx. d i dx λκ µ r +a This, however, implies that also equation 14 has only the trivial solution and the lemma is proved. Let H κ be the unique selfadjoint extension of H κ see 1. The following result is proved in [Sch4]. Proposition.5. For each half-integer κ, the spectrum σh κ coincides with R, and H κ has no eigenvalues. 8

9 3 Rellich property Definition 3.1. Let Ω be an open set in R n and T a selfadjoint operator on L Ω. We say that T has the Rellich property, if the following is true: Let F DT be such that there exists a positive constant K such that f + Tf K, f F. Then F is precompact in L loc Ω for every bounded open set Ω Ω, that is, every sequence f n n F has a subsequence f nk k which convergesstronglyin L Ω. As usual, we identify functions f F with the corresponding functions f Ω L Ω. The classical result that differential operators in L R d have the Rellich property is cited in the appendix. The aim of this section is to show that H κ has the Rellich property. First we proof the following technical lemma needed for Theorem 3.3. Lemma 3.. Let W and A κ, κ Z+1/ be defined as in Lemma.3 and let H κ be the Dirac operator in the κ-th partial wave space see Theorem.1. Then there exists a positive constant K such that R x Ψ jx,θ j=1 C + r +a A κψ j x,θ C dθdx K S 1 H κ Ψ S + Ψ S for all Ψ = t Ψ 1,Ψ DS 1 H κ where Ψ j L S Ω, j = 1,. Proof. First assume Ψ C Ω 4. Observe that 1 I + asinθ r +a βσ 1, r,θ,,π. With the notation of Section we write S 1 H κ as the sum of an unbounded operator S 1 H κ 1 and a bounded S 1 H κ and obtain Recall that 18 S 1 H κ Ψ S 1 S 1 H κ 1 Ψ S S 1 H κ Ψ S = 1 S 1/ H κ 1 Ψx,θ C 4dxdθ S 1 H κ Ψ S R 1 H κ 1 Ψx,θ C 4 4dxdθ S 1 H κ Ψ S. 19 H κ 1 = = R σ3 i σ3 + x [ σ3 i + σ3 x W 1 A κ W r +a r +a iaκ i A κ W 1 A κ W by Lemma.3. Since the first matrix in the above representation is unitary, we obtain ]. 9

10 H κ R 1 Ψx,θ C dxdθ = R i x Ψx,θ+ iaκ r +a Ψx, θ dθdx ia κ C 4 = j=1 R i x Ψ jx,θ+ r +a ia κψ j x,θ dθdx C [ π ] = j=1 R x Ψ jx,θ + C r +a A κψ j x,θ C dθdx π + Ψ j x,θ, R xr +a A κ Ψ j x,θ dθdx C where in the last step we have used that i x is symmetric and A κ commutes with i x and performed integration by parts. There are no boundary terms since Ψ has compact support. There exists a positive constant C such that d dxr +a = dr d dxdr r +a = d r +a dr r +a = r +a r r +a + r M r +a C r +a 1 since the term in is continuous and tends to zero for x ±. Moreover, for every ε > Ψ j, A xr +a κ Ψ j L Ω Ψ j L Ω xr +a A κ Ψ j L Ω ε Ψ j L Ω 1 ε A xr +a κ Ψ j. L Ω Hence, from 19,, 1 and we obtain S 1 H κ Ψ 1 4 j=1 R [ x Ψ jx,θ + 1 ] C ε C r +a A κψ j x,θ C dθdx ε 4 Ψ L Ω 4 S 1 H κ Ψ S. 1

11 C If we choose ε small enough such that ε < 1 and observe that Ψ L Ω 4 S 1/ Ψ S then we can choose an K > large enough such that the estimate in the assertion holds. For any Ψ DH κ there exists a sequence Ψ n n such that Ψ n Ψ, H κ Ψ n H κ Ψ in LS Ω. Noticethatforthis sequencealsotheleft hand sideof 18convergesbythe dominated convergence theorem. Therefore, assertion holds for all Ψ DH κ. Note that K does not depend on Ψ but only on the L -bound C of d dr r +a. Theorem 3.3 Rellich property. Let κ be any half integer, R > and set Let Ψ n n DH κ such that Ω,R := R,R π,π, 3 Ψ n S + H κ Ψ n S K, n N, for some constant K >. Then there exists a subsequence Ψ nl l and a Φ L Ω,R 4 such that in L S Ω,R 4. Ψ nl Φ as l Remark 3.4. As usual, we identify elements Ψ L Ω 4 with elements of Ψ L Ω,R 4 by restriction. The assertion of the Theorem implies that there exists a subsequence Ψ nl l and a Φ L Ω 4 such that χ R Ψ nl χ R Φ as l in L S Ω 4. Proof of Theorem 3.3. Since the norms on LS Ω 4 and L Ω 4 are equivalent, it suffices to consider convergence in the latter space. Let Ψ n =: ψn, 1 ψn with ψn j L Ω, j = 1,. Lemma 3. and the existence of a positive constant δ > such that r +a δ, x R,R imply that there is an K > such that R π R x ψj n x,θ + A κ ψn j x,θ C dxdθ C j=1 K Ψ n S + Hκ Ψ n S KK, 4 for n N. Since L,π, dθ has a basis of orthonormal eigenfunctions g κ m m of A κ with corresponding eigenvalues λ κ m, m Z\{}, see Lemma.3. Hence ψ 1 n can be expanded in the double Fourier series ν,m u n = ν,mα ν,m,n e iνπx/r g κ mθ, α ν,m,n = 1 R R α ν,m,n = 1 R R R R ψ 1 n x,θ, e iνπx/r g κ m θ C dxdθ ψ 1 n dxdθ 1 R Ψ S. 11

12 Moreover, inequality 4 yields [ νπ/r +λ κ m ] α ν,m,n KK. ν,m For fixed ν,m the sequence α ν,m,n n is a bounded sequence, hence it contains a subsequence such that α ν,m,nl nl is a Cauchy sequence for any ν,m by a diagonal series argument. Thus for arbitray L > it follows that 1 R R R ψ 1 n l ψ 1 n j dxdθ = 1 R R R = νπ/r +λ κ m L + νπ/r +λ κ m L νπ/r +λ κ m L ψ 1 n l ψ 1 n j dxdθ = ν,m α ν,m,nl α ν,m,nj α ν,m,nl α ν,m,nj [ νπ/r +λ κ m ] α ν,m,nl α ν,m,nj νπ/r +λ κ m α ν,m,nl α ν,m,nj + KK 1+L, which shows that ψ nl l is a Cauchy sequence in L Ω. With the same argument we find that also the sequence ψ n l l, contains a convergent subsequence in L Ω. 4 Weak local energy decay Since H κ, κ Z\{} is selfadjoint in L S Ω 4, it is generator of a unitary group U κ t := exp its 1 H κ. 5 Theorem 4.1. Let κ Z\{}, R > and Ψ DH κ. For any sequence t n n R there exist a subsequence t nl l and Φ L Ω 4 such that in L Ω 4. U κ t nl Ψ Φ as l Proof. Let Ψ n := U κ t n Ψ. Then, for all n N, Ψ n S + S 1 H κ Ψ n S = U κ t n Ψ S + S 1 H κ U κ t n Ψ S = U κ t n Ψ S + U κ t n S 1 H κ Ψ S = Ψ S + S 1 H κ Ψ S. Thus the assertion follows by Theorem 3.3 with K = Ψ + H κ Ψ. Theorem 4.. For every R > and Φ LS Ω the time mean of the localization of Φ in R,R,π is to zero: [ 1 T R ] π lim U κ tφx,θ dxdθ dt =. T T C 4 T R 1

13 Proof. It is well-known that in the non-extreme Kerr-Newman case the operator H κ has no eigenvalues [Sch4], [WY6], see also Proposition.5. Let χ R x be the characteristic function of R, R. Then, the Rellich property Theorem 3.3 implies that χ R H κ +i 1 is a compact operator. Therefore, the assertion follows from the RAGE theorem, see [RS79]. The decay of the partial waves implies the decay of the solutions of the original Cauchy problem 9: Theorem 4.3. Let Ψ L S Ω 3 4. Its Fourier expansion Ψx,θ,ϕ = converges strongly in L S Ω 3 4 and κ Z+1/ κ Z+1/ e iκϕ Ψ κ x,θ Ψ κ x,θ S <. Then Ψx,θ,ϕ,t = e iκϕ U κ tψ κ x,θ κ Z+1/ is the unique weak solution of i t Ψ = S 1 H Ψ, Ψ = Ψ cf. 9 in L S Ω 3 4 and satisfies lim T for any R >. [ 1 T R T T R ] S 1/ Ψx,θ,ϕ,t dxdθdϕ dt = C 4 Remark 4.4. The above results remain still valid if the electric potential eqr in 4 is substituted by any real-valued C 1 qr -function qr defined on [r +, such that lim r r exists and q r = O1 as r. A Equivalent representations of the time-independent Dirac operator Remark A.1. Ifwedonotapplythetransformations8and1tothewavefunction Ψ, the time-independent Dirac equation takes the form i t Ψ = S 1 Ĥ Ψ 6 with Ĥ = Ĥ1 +H with H as in Theorem.1 and 13

14 ĥ1 Ĥ 1 =, ĥ1 ĥ 1 = σ 3 r +a i r r +a σ 1 i Ĥ is formally symmetric in the Hilbert space with the scalar product Ψ, Φ S = r + θ + cotθ L S := L r +,,π [,π 1 +σ sinθ i. ϕ Ψr,θ,ϕ, Sr,θ,ϕ Φr,θ,ϕ C 4 r +adr sinθdθdϕ. Remark A.. In Theorem.1 we have equipped the Hilbert space L,,π,π with a scalar product, S such that H is formally symmetric. Alternatively, we can work with the usual scalar product, but then we have to define the Dirac operator in the following way: H Ψ = S 1/ HS 1/ Ψ, D H = C,,π [,π4 on L,,π [,π 4 equipped with the usual scalar product Ψ, Φ = R Ψx,θ,ϕ, Φx,θ,ϕ C 4dxdθdϕ. Note that S leaves the space C,,π [,π 4 invariant. B Rellich s theorem Theorem B.1 Rellich s theorem. Let d N, K R d compact and Ω a bounded open neighbourhood of K. Further let s > t N. Then for each c > the set F := {u H s R d : u H s R d c, suppu K} is precompact in H t Ω, i.e., for each sequence f n n F there exists an f H t R d and a subsequence f nk k such that f nk f in H t Ω. Proof. See, e.g., [Miz73]. Acknowledgements M. W. is grateful for the hospitality at Ritsumeikan University, Kusatsu. The research work on this paper started during her visit at Ritsumeikan University, Kusatsu, supported by Open Research Center Project for Private Universities: matching fund subsidy from MEXT,

15 References [BM99] [BS6] F. Belgiorno and M. Martellini. Quantum properties of the electron field in Kerr-Newman black hole manifolds. Phys. Lett. B, 4531-:17, D. Batic and H. Schmid. The Dirac propagator in the Kerr-Newman metric. Progr. Theoret. Phys., 1163: , 6. [Cha76] S. Chandrasekhar. The solution of Dirac s equation in Kerr geometry. Proc. Roy. Soc. London Ser. A, : , [Cha98] S. Chandrasekhar. The mathematical theory of black holes. Oxford Classic Texts in the Physical Sciences. The Clarendon Press Oxford University Press, New York, Reprint of the 199 edition. [FKSY3] F. Finster, N. Kamran, J. Smoller, and S.-T. Yau. The long-time dynamics of Dirac particles in the Kerr-Newman black hole geometry. Adv. Theor. Math. Phys., 71:5 5, 3. [Miz73] [Pag76] [RS79] [Sch4] S. Mizohata. The theory of partial differential equations. Cambridge University Press, New York, Translated from the Japanese by Katsumi Miyahara. D. N. Page. Dirac equation around a charged, rotating black hole. Phys. Rev., D14: , Michael Reed and Barry Simon. Methods of modern mathematical physics. III. Academic Press [Harcourt Brace Jovanovich Publishers], New York, Scattering theory. H. Schmid. Bound state solutions of the Dirac equation in the extreme Kerr geometry. Math. Nachr., 74/75:117 19, 4. [Wei87] J. Weidmann. Spectral theory of ordinary differential operators, volume 158 of Lecture Notes in Mathematics. Springer-Verlag, Berlin Heidelberg New York London Paris Tokyo, [Win6] M. Winklmeier. The angular part of the Dirac equation in the Kerr-Newman metric: Estimates for the eigenvalues. PhD thesis, Universität Bremen, 6. [WY6] M. Winklmeier and O. Yamada. Spectral analysis of radial Dirac operators in the Kerr-Newman metric and its applications to time-periodic solutions. J. Math. Phys., 471:153, 17, 6. 15

M5MF6. Advanced Methods in Derivatives Pricing

M5MF6. 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 information

Conditional Square Functions and Dyadic Perturbations of the Sine-Cosine decomposition for Hardy Martingales

Conditional Square Functions and Dyadic Perturbations of the Sine-Cosine decomposition for Hardy Martingales Conditional Square Functions and Dyadic Perturbations of the Sine-Cosine decomposition for Hardy Martingales arxiv:1611.02653v1 [math.fa] 8 Nov 2016 Paul F. X. Müller October 16, 2016 Abstract We prove

More information

GUESSING MODELS IMPLY THE SINGULAR CARDINAL HYPOTHESIS arxiv: v1 [math.lo] 25 Mar 2019

GUESSING MODELS IMPLY THE SINGULAR CARDINAL HYPOTHESIS arxiv: v1 [math.lo] 25 Mar 2019 GUESSING MODELS IMPLY THE SINGULAR CARDINAL HYPOTHESIS arxiv:1903.10476v1 [math.lo] 25 Mar 2019 Abstract. In this article we prove three main theorems: (1) guessing models are internally unbounded, (2)

More information

"Pricing Exotic Options using Strong Convergence Properties

Pricing Exotic Options using Strong Convergence Properties Fourth Oxford / Princeton Workshop on Financial Mathematics "Pricing Exotic Options using Strong Convergence Properties Klaus E. Schmitz Abe schmitz@maths.ox.ac.uk www.maths.ox.ac.uk/~schmitz Prof. Mike

More information

In Discrete Time a Local Martingale is a Martingale under an Equivalent Probability Measure

In Discrete Time a Local Martingale is a Martingale under an Equivalent Probability Measure In Discrete Time a Local Martingale is a Martingale under an Equivalent Probability Measure Yuri Kabanov 1,2 1 Laboratoire de Mathématiques, Université de Franche-Comté, 16 Route de Gray, 253 Besançon,

More information

Stability in geometric & functional inequalities

Stability in geometric & functional inequalities Stability in geometric & functional inequalities A. Figalli The University of Texas at Austin www.ma.utexas.edu/users/figalli/ Alessio Figalli (UT Austin) Stability in geom. & funct. ineq. Krakow, July

More information

arxiv: v2 [math.lo] 13 Feb 2014

arxiv: v2 [math.lo] 13 Feb 2014 A LOWER BOUND FOR GENERALIZED DOMINATING NUMBERS arxiv:1401.7948v2 [math.lo] 13 Feb 2014 DAN HATHAWAY Abstract. We show that when κ and λ are infinite cardinals satisfying λ κ = λ, the cofinality of the

More information

A class of coherent risk measures based on one-sided moments

A class of coherent risk measures based on one-sided moments A class of coherent risk measures based on one-sided moments T. Fischer Darmstadt University of Technology November 11, 2003 Abstract This brief paper explains how to obtain upper boundaries of shortfall

More information

Rohini Kumar. Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque)

Rohini 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 information

American Option Pricing Formula for Uncertain Financial Market

American Option Pricing Formula for Uncertain Financial Market American Option Pricing Formula for Uncertain Financial Market Xiaowei Chen Uncertainty Theory Laboratory, Department of Mathematical Sciences Tsinghua University, Beijing 184, China chenxw7@mailstsinghuaeducn

More information

SELF-ADJOINT BOUNDARY-VALUE PROBLEMS ON TIME-SCALES

SELF-ADJOINT BOUNDARY-VALUE PROBLEMS ON TIME-SCALES Electronic Journal of Differential Equations, Vol. 2007(2007), No. 175, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) SELF-ADJOINT

More information

4: SINGLE-PERIOD MARKET MODELS

4: SINGLE-PERIOD MARKET MODELS 4: SINGLE-PERIOD MARKET MODELS Marek Rutkowski School of Mathematics and Statistics University of Sydney Semester 2, 2016 M. Rutkowski (USydney) Slides 4: Single-Period Market Models 1 / 87 General Single-Period

More information

THE NUMBER OF UNARY CLONES CONTAINING THE PERMUTATIONS ON AN INFINITE SET

THE NUMBER OF UNARY CLONES CONTAINING THE PERMUTATIONS ON AN INFINITE SET THE NUMBER OF UNARY CLONES CONTAINING THE PERMUTATIONS ON AN INFINITE SET MICHAEL PINSKER Abstract. We calculate the number of unary clones (submonoids of the full transformation monoid) containing the

More information

CONSISTENCY AMONG TRADING DESKS

CONSISTENCY 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 information

Risk minimizing strategies for tracking a stochastic target

Risk minimizing strategies for tracking a stochastic target Risk minimizing strategies for tracking a stochastic target Andrzej Palczewski Abstract We consider a stochastic control problem of beating a stochastic benchmark. The problem is considered in an incomplete

More information

Game Theory: Normal Form Games

Game Theory: Normal Form Games Game Theory: Normal Form Games Michael Levet June 23, 2016 1 Introduction Game Theory is a mathematical field that studies how rational agents make decisions in both competitive and cooperative situations.

More information

arxiv: v1 [math.lo] 27 Mar 2009

arxiv: v1 [math.lo] 27 Mar 2009 arxiv:0903.4691v1 [math.lo] 27 Mar 2009 COMBINATORIAL AND MODEL-THEORETICAL PRINCIPLES RELATED TO REGULARITY OF ULTRAFILTERS AND COMPACTNESS OF TOPOLOGICAL SPACES. V. PAOLO LIPPARINI Abstract. We generalize

More information

Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility model

Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility model Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility model 1(23) Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility

More information

MITCHELL S THEOREM REVISITED. Contents

MITCHELL S THEOREM REVISITED. Contents MITCHELL S THEOREM REVISITED THOMAS GILTON AND JOHN KRUEGER Abstract. Mitchell s theorem on the approachability ideal states that it is consistent relative to a greatly Mahlo cardinal that there is no

More information

The ruin probabilities of a multidimensional perturbed risk model

The 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 information

MATH3075/3975 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS

MATH3075/3975 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS MATH307/37 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS School of Mathematics and Statistics Semester, 04 Tutorial problems should be used to test your mathematical skills and understanding of the lecture material.

More information

Martingales. by D. Cox December 2, 2009

Martingales. by D. Cox December 2, 2009 Martingales by D. Cox December 2, 2009 1 Stochastic Processes. Definition 1.1 Let T be an arbitrary index set. A stochastic process indexed by T is a family of random variables (X t : t T) defined on a

More information

On Existence of Equilibria. Bayesian Allocation-Mechanisms

On Existence of Equilibria. Bayesian Allocation-Mechanisms On Existence of Equilibria in Bayesian Allocation Mechanisms Northwestern University April 23, 2014 Bayesian Allocation Mechanisms In allocation mechanisms, agents choose messages. The messages determine

More information

arxiv: v1 [q-fin.pm] 13 Mar 2014

arxiv: v1 [q-fin.pm] 13 Mar 2014 MERTON PORTFOLIO PROBLEM WITH ONE INDIVISIBLE ASSET JAKUB TRYBU LA arxiv:143.3223v1 [q-fin.pm] 13 Mar 214 Abstract. In this paper we consider a modification of the classical Merton portfolio optimization

More information

Square-Root Measurement for Ternary Coherent State Signal

Square-Root Measurement for Ternary Coherent State Signal ISSN 86-657 Square-Root Measurement for Ternary Coherent State Signal Kentaro Kato Quantum ICT Research Institute, Tamagawa University 6-- Tamagawa-gakuen, Machida, Tokyo 9-86, Japan Tamagawa University

More information

Part 3: Trust-region methods for unconstrained optimization. Nick Gould (RAL)

Part 3: Trust-region methods for unconstrained optimization. Nick Gould (RAL) Part 3: Trust-region methods for unconstrained optimization Nick Gould (RAL) minimize x IR n f(x) MSc course on nonlinear optimization UNCONSTRAINED MINIMIZATION minimize x IR n f(x) where the objective

More information

A Translation of Intersection and Union Types

A Translation of Intersection and Union Types A Translation of Intersection and Union Types for the λ µ-calculus Kentaro Kikuchi RIEC, Tohoku University kentaro@nue.riec.tohoku.ac.jp Takafumi Sakurai Department of Mathematics and Informatics, Chiba

More information

L 2 -theoretical study of the relation between the LIBOR market model and the HJM model Takashi Yasuoka

L 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 information

Martingale Pricing Theory in Discrete-Time and Discrete-Space Models

Martingale 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 information

No-arbitrage theorem for multi-factor uncertain stock model with floating interest rate

No-arbitrage theorem for multi-factor uncertain stock model with floating interest rate Fuzzy Optim Decis Making 217 16:221 234 DOI 117/s17-16-9246-8 No-arbitrage theorem for multi-factor uncertain stock model with floating interest rate Xiaoyu Ji 1 Hua Ke 2 Published online: 17 May 216 Springer

More information

Drunken Birds, Brownian Motion, and Other Random Fun

Drunken Birds, Brownian Motion, and Other Random Fun Drunken Birds, Brownian Motion, and Other Random Fun Michael Perlmutter Department of Mathematics Purdue University 1 M. Perlmutter(Purdue) Brownian Motion and Martingales Outline Review of Basic Probability

More information

Pricing Dynamic Solvency Insurance and Investment Fund Protection

Pricing 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 information

symmys.com 3.2 Projection of the invariants to the investment horizon

symmys.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 information

An overview of some financial models using BSDE with enlarged filtrations

An overview of some financial models using BSDE with enlarged filtrations An overview of some financial models using BSDE with enlarged filtrations Anne EYRAUD-LOISEL Workshop : Enlargement of Filtrations and Applications to Finance and Insurance May 31st - June 4th, 2010, Jena

More information

GPD-POT and GEV block maxima

GPD-POT and GEV block maxima Chapter 3 GPD-POT and GEV block maxima This chapter is devoted to the relation between POT models and Block Maxima (BM). We only consider the classical frameworks where POT excesses are assumed to be GPD,

More information

arxiv: v1 [math.st] 18 Sep 2018

arxiv: v1 [math.st] 18 Sep 2018 Gram Charlier and Edgeworth expansion for sample variance arxiv:809.06668v [math.st] 8 Sep 08 Eric Benhamou,* A.I. SQUARE CONNECT, 35 Boulevard d Inkermann 900 Neuilly sur Seine, France and LAMSADE, Universit

More information

Equivalence between Semimartingales and Itô Processes

Equivalence between Semimartingales and Itô Processes International Journal of Mathematical Analysis Vol. 9, 215, no. 16, 787-791 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ijma.215.411358 Equivalence between Semimartingales and Itô Processes

More information

COMBINATORICS OF REDUCTIONS BETWEEN EQUIVALENCE RELATIONS

COMBINATORICS OF REDUCTIONS BETWEEN EQUIVALENCE RELATIONS COMBINATORICS OF REDUCTIONS BETWEEN EQUIVALENCE RELATIONS DAN HATHAWAY AND SCOTT SCHNEIDER Abstract. We discuss combinatorial conditions for the existence of various types of reductions between equivalence

More information

Chapter 3: Black-Scholes Equation and Its Numerical Evaluation

Chapter 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 information

A Property Equivalent to n-permutability for Infinite Groups

A Property Equivalent to n-permutability for Infinite Groups Journal of Algebra 221, 570 578 (1999) Article ID jabr.1999.7996, available online at http://www.idealibrary.com on A Property Equivalent to n-permutability for Infinite Groups Alireza Abdollahi* and Aliakbar

More information

Lecture 4. Finite difference and finite element methods

Lecture 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 information

Hints on Some of the Exercises

Hints 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 information

Shape of the Yield Curve Under CIR Single Factor Model: A Note

Shape of the Yield Curve Under CIR Single Factor Model: A Note Shape of the Yield Curve Under CIR Single Factor Model: A Note Raymond Kan University of Chicago June, 1992 Abstract This note derives the shapes of the yield curve as a function of the current spot rate

More information

Computation of one-sided probability density functions from their cumulants

Computation of one-sided probability density functions from their cumulants Journal of Mathematical Chemistry, Vol. 41, No. 1, January 27 26) DOI: 1.17/s191-6-969-x Computation of one-sided probability density functions from their cumulants Mário N. Berberan-Santos Centro de Química-Física

More information

Characterization of the Optimum

Characterization 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 information

Stochastic calculus Introduction I. Stochastic Finance. C. Azizieh VUB 1/91. C. Azizieh VUB Stochastic Finance

Stochastic calculus Introduction I. Stochastic Finance. C. Azizieh VUB 1/91. C. Azizieh VUB Stochastic Finance Stochastic Finance C. Azizieh VUB C. Azizieh VUB Stochastic Finance 1/91 Agenda of the course Stochastic calculus : introduction Black-Scholes model Interest rates models C. Azizieh VUB Stochastic Finance

More information

On the Lower Arbitrage Bound of American Contingent Claims

On the Lower Arbitrage Bound of American Contingent Claims On the Lower Arbitrage Bound of American Contingent Claims Beatrice Acciaio Gregor Svindland December 2011 Abstract We prove that in a discrete-time market model the lower arbitrage bound of an American

More information

Minimal Variance Hedging in Large Financial Markets: random fields approach

Minimal Variance Hedging in Large Financial Markets: random fields approach Minimal Variance Hedging in Large Financial Markets: random fields approach Giulia Di Nunno Third AMaMeF Conference: Advances in Mathematical Finance Pitesti, May 5-1 28 based on a work in progress with

More information

CHARACTERIZATION OF CLOSED CONVEX SUBSETS OF R n

CHARACTERIZATION OF CLOSED CONVEX SUBSETS OF R n CHARACTERIZATION OF CLOSED CONVEX SUBSETS OF R n Chebyshev Sets A subset S of a metric space X is said to be a Chebyshev set if, for every x 2 X; there is a unique point in S that is closest to x: Put

More information

DENSITY OF PERIODIC GEODESICS IN THE UNIT TANGENT BUNDLE OF A COMPACT HYPERBOLIC SURFACE

DENSITY OF PERIODIC GEODESICS IN THE UNIT TANGENT BUNDLE OF A COMPACT HYPERBOLIC SURFACE DENSITY OF PERIODIC GEODESICS IN THE UNIT TANGENT BUNDLE OF A COMPACT HYPERBOLIC SURFACE Marcos Salvai FaMAF, Ciudad Universitaria, 5000 Córdoba, Argentina. e-mail: salvai@mate.uncor.edu Abstract Let S

More information

Received May 27, 2009; accepted January 14, 2011

Received May 27, 2009; accepted January 14, 2011 MATHEMATICAL COMMUNICATIONS 53 Math. Coun. 6(20), 53 538. I σ -Convergence Fatih Nuray,, Hafize Gök and Uǧur Ulusu Departent of Matheatics, Afyon Kocatepe University, 03200 Afyonkarahisar, Turkey Received

More information

On the h-vector of a Lattice Path Matroid

On the h-vector of a Lattice Path Matroid On the h-vector of a Lattice Path Matroid Jay Schweig Department of Mathematics University of Kansas Lawrence, KS 66044 jschweig@math.ku.edu Submitted: Sep 16, 2009; Accepted: Dec 18, 2009; Published:

More information

THE MARTINGALE METHOD DEMYSTIFIED

THE MARTINGALE METHOD DEMYSTIFIED THE MARTINGALE METHOD DEMYSTIFIED SIMON ELLERSGAARD NIELSEN Abstract. We consider the nitty gritty of the martingale approach to option pricing. These notes are largely based upon Björk s Arbitrage Theory

More information

INTERIM CORRELATED RATIONALIZABILITY IN INFINITE GAMES

INTERIM CORRELATED RATIONALIZABILITY IN INFINITE GAMES INTERIM CORRELATED RATIONALIZABILITY IN INFINITE GAMES JONATHAN WEINSTEIN AND MUHAMET YILDIZ A. We show that, under the usual continuity and compactness assumptions, interim correlated rationalizability

More information

Phys. Lett. A, 372/17, (2008),

Phys. Lett. A, 372/17, (2008), Phys. Lett. A, 372/17, (2008), 3064-3070. 1 Wave scattering by many small particles embedded in a medium. A. G. Ramm (Mathematics Department, Kansas State University, Manhattan, KS66506, USA and TU Darmstadt,

More information

Optimal robust bounds for variance options and asymptotically extreme models

Optimal robust bounds for variance options and asymptotically extreme models Optimal robust bounds for variance options and asymptotically extreme models Alexander Cox 1 Jiajie Wang 2 1 University of Bath 2 Università di Roma La Sapienza Advances in Financial Mathematics, 9th January,

More information

an example Giacomo Boffi

an example Giacomo Boffi an example http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano April 2, 27 vt P + x u(x, t L A uniform beam, (unit mass

More information

Financial Modelling Using Discrete Stochastic Calculus

Financial Modelling Using Discrete Stochastic Calculus Preprint typeset in JHEP style - HYPER VERSION Financial Modelling Using Discrete Stochastic Calculus Eric A. Forgy, Ph.D. E-mail: eforgy@yahoo.com Abstract: In the present report, a review of discrete

More information

Tangent 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. 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 information

STOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL

STOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL STOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL YOUNGGEUN YOO Abstract. Ito s lemma is often used in Ito calculus to find the differentials of a stochastic process that depends on time. This paper will introduce

More information

Large Deviations and Stochastic Volatility with Jumps: Asymptotic Implied Volatility for Affine Models

Large 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 information

Introduction to Probability Theory and Stochastic Processes for Finance Lecture Notes

Introduction to Probability Theory and Stochastic Processes for Finance Lecture Notes Introduction to Probability Theory and Stochastic Processes for Finance Lecture Notes Fabio Trojani Department of Economics, University of St. Gallen, Switzerland Correspondence address: Fabio Trojani,

More information

The Limiting Distribution for the Number of Symbol Comparisons Used by QuickSort is Nondegenerate (Extended Abstract)

The Limiting Distribution for the Number of Symbol Comparisons Used by QuickSort is Nondegenerate (Extended Abstract) The Limiting Distribution for the Number of Symbol Comparisons Used by QuickSort is Nondegenerate (Extended Abstract) Patrick Bindjeme 1 James Allen Fill 1 1 Department of Applied Mathematics Statistics,

More information

The rth moment of a real-valued random variable X with density f(x) is. x r f(x) dx

The rth moment of a real-valued random variable X with density f(x) is. x r f(x) dx 1 Cumulants 1.1 Definition The rth moment of a real-valued random variable X with density f(x) is µ r = E(X r ) = x r f(x) dx for integer r = 0, 1,.... The value is assumed to be finite. Provided that

More information

A Note on the No Arbitrage Condition for International Financial Markets

A Note on the No Arbitrage Condition for International Financial Markets A Note on the No Arbitrage Condition for International Financial Markets FREDDY DELBAEN 1 Department of Mathematics Vrije Universiteit Brussel and HIROSHI SHIRAKAWA 2 Department of Industrial and Systems

More information

Basic Concepts and Examples in Finance

Basic Concepts and Examples in Finance Basic Concepts and Examples in Finance Bernardo D Auria email: bernardo.dauria@uc3m.es web: www.est.uc3m.es/bdauria July 5, 2017 ICMAT / UC3M The Financial Market The Financial Market We assume there are

More information

A No-Arbitrage Theorem for Uncertain Stock Model

A No-Arbitrage Theorem for Uncertain Stock Model Fuzzy Optim Decis Making manuscript No (will be inserted by the editor) A No-Arbitrage Theorem for Uncertain Stock Model Kai Yao Received: date / Accepted: date Abstract Stock model is used to describe

More information

Lecture Quantitative Finance Spring Term 2015

Lecture Quantitative Finance Spring Term 2015 implied Lecture Quantitative Finance Spring Term 2015 : May 7, 2015 1 / 28 implied 1 implied 2 / 28 Motivation and setup implied the goal of this chapter is to treat the implied which requires an algorithm

More information

Math-Stat-491-Fall2014-Notes-V

Math-Stat-491-Fall2014-Notes-V Math-Stat-491-Fall2014-Notes-V Hariharan Narayanan December 7, 2014 Martingales 1 Introduction Martingales were originally introduced into probability theory as a model for fair betting games. Essentially

More information

Permutation Factorizations and Prime Parking Functions

Permutation Factorizations and Prime Parking Functions Permutation Factorizations and Prime Parking Functions Amarpreet Rattan Department of Combinatorics and Optimization University of Waterloo Waterloo, ON, Canada N2L 3G1 arattan@math.uwaterloo.ca June 10,

More information

Optimal Stopping in Infinite Horizon: an Eigenfunction Expansion Approach

Optimal Stopping in Infinite Horizon: an Eigenfunction Expansion Approach Optimal Stopping in Infinite Horizon: an Eigenfunction Expansion Approach Lingfei Li a, Vadim Linetsky b a Department of Systems Engineering and Engineering Management, The Chinese University of Hong Kong

More information

Order book resilience, price manipulations, and the positive portfolio problem

Order book resilience, price manipulations, and the positive portfolio problem Order book resilience, price manipulations, and the positive portfolio problem Alexander Schied Mannheim University PRisMa Workshop Vienna, September 28, 2009 Joint work with Aurélien Alfonsi and Alla

More information

LECTURE 2: MULTIPERIOD MODELS AND TREES

LECTURE 2: MULTIPERIOD MODELS AND TREES LECTURE 2: MULTIPERIOD MODELS AND TREES 1. Introduction One-period models, which were the subject of Lecture 1, are of limited usefulness in the pricing and hedging of derivative securities. In real-world

More information

Laurence Boxer and Ismet KARACA

Laurence Boxer and Ismet KARACA SOME PROPERTIES OF DIGITAL COVERING SPACES Laurence Boxer and Ismet KARACA Abstract. In this paper we study digital versions of some properties of covering spaces from algebraic topology. We correct and

More information

Local Volatility Dynamic Models

Local Volatility Dynamic Models René Carmona Bendheim Center for Finance Department of Operations Research & Financial Engineering Princeton University Columbia November 9, 27 Contents Joint work with Sergey Nadtochyi Motivation 1 Understanding

More information

Interpolation of κ-compactness and PCF

Interpolation of κ-compactness and PCF Comment.Math.Univ.Carolin. 50,2(2009) 315 320 315 Interpolation of κ-compactness and PCF István Juhász, Zoltán Szentmiklóssy Abstract. We call a topological space κ-compact if every subset of size κ has

More information

A note on the existence of unique equivalent martingale measures in a Markovian setting

A note on the existence of unique equivalent martingale measures in a Markovian setting Finance Stochast. 1, 251 257 1997 c Springer-Verlag 1997 A note on the existence of unique equivalent martingale measures in a Markovian setting Tina Hviid Rydberg University of Aarhus, Department of Theoretical

More information

The stochastic calculus

The stochastic calculus Gdansk A schedule of the lecture Stochastic differential equations Ito calculus, Ito process Ornstein - Uhlenbeck (OU) process Heston model Stopping time for OU process Stochastic differential equations

More information

Disaster risk and its implications for asset pricing Online appendix

Disaster risk and its implications for asset pricing Online appendix Disaster risk and its implications for asset pricing Online appendix Jerry Tsai University of Oxford Jessica A. Wachter University of Pennsylvania December 12, 2014 and NBER A The iid model This section

More information

DEPTH OF BOOLEAN ALGEBRAS SHIMON GARTI AND SAHARON SHELAH

DEPTH OF BOOLEAN ALGEBRAS SHIMON GARTI AND SAHARON SHELAH DEPTH OF BOOLEAN ALGEBRAS SHIMON GARTI AND SAHARON SHELAH Abstract. Suppose D is an ultrafilter on κ and λ κ = λ. We prove that if B i is a Boolean algebra for every i < κ and λ bounds the Depth of every

More information

Estimation of Value at Risk and ruin probability for diffusion processes with jumps

Estimation of Value at Risk and ruin probability for diffusion processes with jumps Estimation of Value at Risk and ruin probability for diffusion processes with jumps Begoña Fernández Universidad Nacional Autónoma de México joint work with Laurent Denis and Ana Meda PASI, May 21 Begoña

More information

Universität Regensburg Mathematik

Universität Regensburg Mathematik Universität Regensburg Mathematik Modeling financial markets with extreme risk Tobias Kusche Preprint Nr. 04/2008 Modeling financial markets with extreme risk Dr. Tobias Kusche 11. January 2008 1 Introduction

More information

A THREE-FACTOR CONVERGENCE MODEL OF INTEREST RATES

A 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 information

Optimum Thresholding for Semimartingales with Lévy Jumps under the mean-square error

Optimum Thresholding for Semimartingales with Lévy Jumps under the mean-square error Optimum Thresholding for Semimartingales with Lévy Jumps under the mean-square error José E. Figueroa-López Department of Mathematics Washington University in St. Louis Spring Central Sectional Meeting

More information

A NEW NOTION OF TRANSITIVE RELATIVE RETURN RATE AND ITS APPLICATIONS USING STOCHASTIC DIFFERENTIAL EQUATIONS. Burhaneddin İZGİ

A NEW NOTION OF TRANSITIVE RELATIVE RETURN RATE AND ITS APPLICATIONS USING STOCHASTIC DIFFERENTIAL EQUATIONS. Burhaneddin İZGİ A NEW NOTION OF TRANSITIVE RELATIVE RETURN RATE AND ITS APPLICATIONS USING STOCHASTIC DIFFERENTIAL EQUATIONS Burhaneddin İZGİ Department of Mathematics, Istanbul Technical University, Istanbul, Turkey

More information

In chapter 5, we approximated the Black-Scholes model

In chapter 5, we approximated the Black-Scholes model Chapter 7 The Black-Scholes Equation In chapter 5, we approximated the Black-Scholes model ds t /S t = µ dt + σ dx t 7.1) with a suitable Binomial model and were able to derive a pricing formula for option

More information

Econ 582 Nonlinear Regression

Econ 582 Nonlinear Regression Econ 582 Nonlinear Regression Eric Zivot June 3, 2013 Nonlinear Regression In linear regression models = x 0 β (1 )( 1) + [ x ]=0 [ x = x] =x 0 β = [ x = x] [ x = x] x = β it is assumed that the regression

More information

Girsanov s Theorem. Bernardo D Auria web: July 5, 2017 ICMAT / UC3M

Girsanov s Theorem. Bernardo D Auria   web:   July 5, 2017 ICMAT / UC3M Girsanov s Theorem Bernardo D Auria email: bernardo.dauria@uc3m.es web: www.est.uc3m.es/bdauria July 5, 2017 ICMAT / UC3M Girsanov s Theorem Decomposition of P-Martingales as Q-semi-martingales Theorem

More information

Homework Assignments

Homework Assignments Homework Assignments Week 1 (p. 57) #4.1, 4., 4.3 Week (pp 58 6) #4.5, 4.6, 4.8(a), 4.13, 4.0, 4.6(b), 4.8, 4.31, 4.34 Week 3 (pp 15 19) #1.9, 1.1, 1.13, 1.15, 1.18 (pp 9 31) #.,.6,.9 Week 4 (pp 36 37)

More information

IEOR E4703: Monte-Carlo Simulation

IEOR E4703: Monte-Carlo Simulation IEOR E4703: Monte-Carlo Simulation Simulating Stochastic Differential Equations Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com

More information

On generalized resolvents of symmetric operators of defect one with finitely many negative squares

On generalized resolvents of symmetric operators of defect one with finitely many negative squares 21 On generalized resolvents of symmetric operators of defect one with finitely many negative squares Jussi Behrndt and Carsten Trunk Abstract Behrndt, Jussi and Carsten Trunk (2005). On generalized resolvents

More information

Functional vs Banach space stochastic calculus & strong-viscosity solutions to semilinear parabolic path-dependent PDEs.

Functional vs Banach space stochastic calculus & strong-viscosity solutions to semilinear parabolic path-dependent PDEs. Functional vs Banach space stochastic calculus & strong-viscosity solutions to semilinear parabolic path-dependent PDEs Andrea Cosso LPMA, Université Paris Diderot joint work with Francesco Russo ENSTA,

More information

δ j 1 (S j S j 1 ) (2.3) j=1

δ j 1 (S j S j 1 ) (2.3) j=1 Chapter The Binomial Model Let S be some tradable asset with prices and let S k = St k ), k = 0, 1,,....1) H = HS 0, S 1,..., S N 1, S N ).) be some option payoff with start date t 0 and end date or maturity

More information

Global convergence rate analysis of unconstrained optimization methods based on probabilistic models

Global convergence rate analysis of unconstrained optimization methods based on probabilistic models Math. Program., Ser. A DOI 10.1007/s10107-017-1137-4 FULL LENGTH PAPER Global convergence rate analysis of unconstrained optimization methods based on probabilistic models C. Cartis 1 K. Scheinberg 2 Received:

More information

Option Pricing Formula for Fuzzy Financial Market

Option Pricing Formula for Fuzzy Financial Market Journal of Uncertain Systems Vol.2, No., pp.7-2, 28 Online at: www.jus.org.uk Option Pricing Formula for Fuzzy Financial Market Zhongfeng Qin, Xiang Li Department of Mathematical Sciences Tsinghua University,

More information

The Semi-Weak Square Principle

The Semi-Weak Square Principle The Semi-Weak Square Principle Maxwell Levine Universität Wien Kurt Gödel Research Center for Mathematical Logic Währinger Straße 25 1090 Wien Austria maxwell.levine@univie.ac.at Abstract Cummings, Foreman,

More information

Non replication of options

Non replication of options Non replication of options Christos Kountzakis, Ioannis A Polyrakis and Foivos Xanthos June 30, 2008 Abstract In this paper we study the scarcity of replication of options in the two period model of financial

More information

Project 1: Double Pendulum

Project 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 information

Yao s Minimax Principle

Yao s Minimax Principle Complexity of algorithms The complexity of an algorithm is usually measured with respect to the size of the input, where size may for example refer to the length of a binary word describing the input,

More information

Exam in TFY4275/FY8907 CLASSICAL TRANSPORT THEORY Feb 14, 2014

Exam in TFY4275/FY8907 CLASSICAL TRANSPORT THEORY Feb 14, 2014 NTNU Page 1 of 5 Institutt for fysikk Contact during the exam: Professor Ingve Simonsen Exam in TFY4275/FY8907 CLASSICAL TRANSPORT THEORY Feb 14, 2014 Allowed help: Alternativ D All written material This

More information