NUMERICAL CALCULATION OF RANDOM MATRIX DISTRIBUTIONS AND ORTHOGONAL POLYNOMIALS. Sheehan Olver NA Group, Oxford
|
|
- Rosanna Blair
- 5 years ago
- Views:
Transcription
1 NUMERICAL CALCULATION OF RANDOM MATRIX DISTRIBUTIONS AND ORTHOGONAL POLYNOMIALS Sheehan Olver NA Group, Oxford
2 We are interested in numerically computing eigenvalue statistics of the GUE ensembles, i.e., Hermitian matrices with the distribution (given a particular V) e ntr V (M) dm Z n First question: can we automatically calculate the global mean distribution of eigenvalues, i.e., the equilibrium measure? Second question: how do statistics which satisfy universality laws differ when n is finite? In short, we do the Riemann Hilbert approach in a numerical way
3 Many other physical and mathematical objects have since been found to also be described by random matrix theory, including: Quantum billiards Random Ising model for glass Trees in the Scandinavian forest Resonance frequencies of structural materials Even the nontrivial zeros of the Riemann zeta function!
4 OUTLINE. Relationship of random matrix theory and orthogonal polynomials 2. Equilibrium measures supported on a single interval. Equivalence to a simple Newton iteration 2. Equilibrium measures supported on multiple intervals 3. Computation of large order orthogonal polynomials, through Riemann Hilbert problems 4. Calculation of gap eigenvalue statistics
5 Relationship of random matrix theory with orthogonal polynomials For the Gaussian unitary ensemble (i.e., large random Hermitian matrices), gap eigenvalue statistics can be reduced to the Fredholm determinant det [I + K n Ω ] where the kernel the weight K n is constructed from pn and pn+; the orthogonal polynomials with respect to e nv (x) dx Fredholm determinants can be easily computed numerically (Bornemann 200), as long as we can evaluate the kernel Therefore, computing distributions depends only on evaluating the orthogonal polynomials We will construct a numerical method for computing pn and pn+ whose computational cost is independent of n, which requires first calculating the equilibrium measure
6 EQUILIBRIUM MEASURES Suppose the real line is a conductor, on which n discrete charges are placed, with total unit charge Suppose further an external field V is present 4 3 The equilibrium measure dµ = ψ(x)dx is the limiting distribution (weak* limit) of the charges On the right is a sketch for V(x) = x (see eg. Saff and Totik 997)
7 Applications of the equilibrium measure Global mean distribution of the eigenvalues of GUE matrix with distribution Z n e ntr V (M) dm Distribution of near optimum interpolation points (Fekete points) Distribution of zeros of orthogonal polynomials with the weight e nv (x) dx For V(x) = x 2 these are scaled Hermite polynomials In the last slide, I cheated and plotted these roots Computing orthogonal polynomials Best rational approximation
8 FORMAL DEFINITION Given an external field V : R R, the equilibrium measure is the unique Borel measure dµ = ψ(x)dx such that log t s dµ(t)dµ(s)+ V (s)dµ(s) is minimal.
9 This definition can be reduced to the following Euler Lagrange formulation: 2 log 2 log dµ + V (z) = for z supp µ x z dµ + V (z) for all real z x z We let g(z) = log(z x)dµ so that (where ± imply limit from the left and right) Differentiating we get g + + g = V and g log z φ + + φ = V and φ z for φ(z) =g (z) = dµ(x) z x
10 Now suppose we manage to find an analytic function such that φ + + φ = V and φ z on some subset Γ of the real line In certain cases (for example, if V is convex), there is only one Γ such that this is possible, hence Γ must be supp μ Plemelj s lemma then tells us that we can find dµ = ψ(x)dx by: i 2π φ + (x) φ (x) = ψ(x)
11 For a given Γ, we can find all solutions to φ + + φ = V and φ(z) c Γ z The goal, then, is to choose Γ so that: cγ is precisely φ is bounded This is the inverse Cauchy transform, which we denote by P Γf
12 φ + (z)+φ (z) =f(z) and φ( ) =0 PROBLEM ON THE CIRCLE f(z) = k= ˆf k z k φ + φ
13 φ + (z)+φ (z) =f(z) and φ( ) =0 PROBLEM ON THE CIRCLE f(z) = k= ˆf k z k φ + = ˆf k z k φ = ˆf k z k k=0 k=
14 PROBLEM ON THE UNIT INTERVAL Consider the Joukowski map from the unit circle to the unit interval J(z) = 2 z + z Functions analytic inside and outside the unit circle are mapped to functions analytic off the unit interval.
15 We define four inverses to the Joukowski map: J + (x) =x x x + J (x) =x + x x + J (x) =x +i x +x J (x) =x i x +x
16 Suppose we compute the inverse Cauchy transform of the mapped function on the unit circle φ = P f(j( )) 2 Then the following function satisfies the necessary jump φ(j + (x)) + φ(j (x)) f(j(z)) φ + f(j(z)) φ 2 2 φ + (J φ + (J (x)) + φ (J (x)) (x)) + φ (J (x))
17 The problem: 2 φ(j + ( )) + φ(j ( )) = 2 [φ(0) + φ( )] = ˆf = ˆf 0 2 =0 Fortunately, κ(z) = z + z has no jump on the contour: κ + + κ = 0 and κ(z) z Therefore, we find that P ξ f = 2 φ(j + (x)) + φ(j (x)) ˆf 0 2 zκ(z)+ξκ(z) where there is now a free parameter
18 Now suppose f is sufficiently smooth, so that it can be expanded into Chebyshev series (in O(n log n) time using the DCT): f(x) = k=0 ˇf k T k (x) k Then f(j(z)) = ˇf k= ˇf k z k It follows that (a numerically stable, uniformly convergent expression) P ξ f(z) = 2 k=0 ˇf k J + (z) k ˇf 0 2 zκ(z)+ξκ(z)
19 OTHER INTERVALS M (a,b) (z) = 2z a b b a maps the interval (a, b) to the unit interval. Let ˇf (a,b),k be the Chebyshev coefficents over (a, b), so that f(z) = k=0 ˇf (a,b),k T k (M (a,b) (z)) k Then P (a,b),ξ f(z) = 2 k=0 ˇf (a,b),k J + (M (a,b) (z)) k ˇf (a,b),0 2 M (a,b) (z)κ(m (a,b) (z))+ξκ(m (a,b) (z))
20 P (a,b),ξ f(z) = 2 k=0 ˇf (a,b),k J + (M (a,b) (z)) k ˇf (a,b),0 2 M (a,b) (z)κ(m (a,b) (z)) + ξκ(m (a,b) (z)) Let f =. The only way the solution is bounded is if Since we require V ξ = 0 and ˇf(a,b),0 =0 J + (z) 2z + O z 2 and M (a,b) (z) 2z b a + O() so that b a ˇf (a,b), = 8 P (a,b) f(z) z (equivalent to standard conditions in terms of moments)
21 We thus have two unknowns, a and b, with two conditions: b a ˇf (a,b),0 = 0 and ˇf (a,b), = 8 The Chebyshev coefficients can be computed efficiently using the FFT Moreover, each of these equations are differentiable with respect to a and b Thus we can setup a trivial Newton iteration to determine a and b This iteration is guaranteed to converge to supp µ whenever V is convex or, if supp µ is a single interval and the initial guess of a and b are sufficiently accurate We then recover the equilibrium measure using the formula M(a,b) (x) dµ = π k= ˇf k U k (M (a,b) (x)) dx
22 x 00 e x x dµ dµ
23 5 x2 4 5 x x x sin 3x + 0x dµ dµ
24 g + + g = V µ g + + g V g = Pf Pf
25 æ κ Σ (z) z = a b J+ (M Σ (z)), 2 M Σ (z)κ Σ (z) z = b a MΣ (z) J+ (M Σ (z)), 2 J+ (M Σ (z)) z = b a 4 2 J + (M Σ (z)) 2 J+ (M Σ (z)), J+ (M Σ (z)) k z = b a 4 k + J + (M Σ (z)) k+ k J + (M Σ (z)) k (,b) (a, b) J (M Σ(x)) J (M Σ(x)) g(z) z = O z.
26 g + + g V 5 x2 4 5 x x x sin 3x + 0x
27 MULTIPLE INTERVALS P Γ Γ Γ = Γ Γ 2 P Γ f = P Γ r + P Γ2 r 2 r r 2
28 r r 2 P Γλ r κ (z) = ř λ,k P Γλ T k (M Γλ (z)) x Γ λ,...,xγ λ N r λ,k r (x Γ j )+[P+ Γ 2 + P Γ 2 ]r 2 (x Γ j )=f(xγ j ) P + Γ + P Γ r (x Γ 2 j )+r 2(x Γ 2 j )=f(xγ 2 j ) Γ Γ 2 f O(n n)
29 ř,0 =ř 2,0 =0 z Γ =(a,b ) Γ 2 =(a 2,b 2 ) b a 8 ř, + b 2 a 2 8 ř 2, g = Pf g + + g = V Γ Γ 2 g + (b )+g (b ) V (b )=g + (a 2 )+g (a 2 ) V (a 2 )
30 V (x) = ( 3+x)( 2+x)( + x)(2 + x)(3 + x)( +2x) α = α
31 ORTHOGONAL POLYNOMIALS
32 Once we have computed the equilibrium measure, we can set up the Riemann Hilbert problem for the orthogonal polynomials RH problems can be computed numerically (as described last week) We need to write the RH problem so that the solution is smooth (no singularities) along each contour, and localized around stationary points This will be accurate in the asymptotic regime: arbitrarily large n The key result: we can systematically compute arbitrarily large order orthogonal polynomials with respect to general weights We could also use this to compute the nodes and weights of Gaussian quadrature rules with respect to general weights in O(n) time
33 The unaltered Riemann Hilbert problem for orthogonal polynomials is Φ + = Φ e nv z n O z and Φ n O z n z n We know how to compute g such that g(z) log z and g + + g = V on supp dµ We thus let Φ = Ψ e ng e ng where Ψ I and Ψ + = Ψ e ng e nv e ng+ e ng = Ψ e n(g g + ) e n(g+ +g V ) e n(g+ g ) e ng+ (based on Deift 999)
34 From the construction of g, the matrix e n(g g + ) e n(g+ +g V ) has the following nice properties: e n(g+ g ) Off of supp µ, it has the form e n(g + +g V ) which decays exponentially (g decays at infinity, so V dominates) On supp µ, it has the form e n(g g + ) e n = L e n(g+ g ) for L = e n/2 e n/2 e n(v 2g ) e n(v 2g+ ) L (based on Deift 999)
35 We thus want to solve the RH problem (where all the jumps are computable numerically!) e n(v 2g) e n(g + +g V +) e n(2g V +) e n(v 2g) (based on Deift 999)
36 P (z) P + = P P (z) I Q Q 2
37 P e n(v 2g) P UP UP UP UQ UQ 2 PQ PQ 2 UP UP UP P e n(v 2g) P
38 RANDOM MATRIX DISTRIBUTIONS
39 We want to compute gap probabilities These are expressible as det [I + K n Ω ] for (where Φ is still the solution the RH problem) K n (x, y) = 2πi e n/2[v (x)+v (y)] Φ (x)φ 2 (y) Φ (y)φ 2 (x) x y
40 ALGORITHM Input: potential V, dimension of random matrix n, gap interval Ω Output: probability that there are no eigenvalues in Ω Step : Compute the equilibrium measure using Newton iteration Step 2: Construct and solve the orthogonal polynomial RH problem numerically Step 3: Use Bornemann s Fredholm determinant solver
41 UNIVERSALITY Let Ω = a + ( s, s) K n (a, a) Then it is known that, for a inside the support of the equilibrium measure and any polynomial potential, the Fredholm determinant converges to the Fredholm determinant over ( s,s) with the sine kernel K (x, y) = sin π(x y) x y However, for finite n, the distributions vary
42 n = 30, 50, x 2 x 4 sin 3x + 0x 20
43 Potential of numerical RH approach for random matrix theory Discovery of unknown universality laws? Inverse problem: determining information about the potential V from observed data? Can we determine which random matrix ensemble generates the zeros of the Riemann Zeta function? Modelling physical systems?
Computation of equilibrium measures
Report no. 0/8 Computation of equilibrium measures Sheehan Olver Abstract We present a new way of computing equilibrium measures, based on the Riemann Hilbert formulation. For equilibrium measures whose
More informationNumerical solution of Riemann Hilbert problems: random matrix theory and orthogonal polynomials
Numerical solution of Riemann Hilbert problems: random matrix theory and orthogonal polynomials Sheehan Olver and Thomas Trogdon February, 013 Abstract Recently, a general approach for solving Riemann
More informationOPSF, Random Matrices and Riemann-Hilbert problems
OPSF, Random Matrices and Riemann-Hilbert problems School on Orthogonal Polynomials in Approximation Theory and Mathematical Physics, ICMAT 23 27 October, 207 Plan of the course lecture : Orthogonal Polynomials
More informationOPSF, Random Matrices and Riemann-Hilbert problems
OPSF, Random Matrices and Riemann-Hilbert problems School on Orthogonal Polynomials in Approximation Theory and Mathematical Physics, ICMAT 23 27 October, 2017 Plan of the course lecture 1: Orthogonal
More informationUniversality of distribution functions in random matrix theory Arno Kuijlaars Katholieke Universiteit Leuven, Belgium
Universality of distribution functions in random matrix theory Arno Kuijlaars Katholieke Universiteit Leuven, Belgium SEA 06@MIT, Workshop on Stochastic Eigen-Analysis and its Applications, MIT, Cambridge,
More informationMarkov operators, classical orthogonal polynomial ensembles, and random matrices
Markov operators, classical orthogonal polynomial ensembles, and random matrices M. Ledoux, Institut de Mathématiques de Toulouse, France 5ecm Amsterdam, July 2008 recent study of random matrix and random
More informationA Riemann Hilbert approach to Jacobi operators and Gaussian quadrature
A Riemann Hilbert approach to Jacobi operators and Gaussian quadrature Thomas Trogdon trogdon@cims.nyu.edu Courant Institute of Mathematical Sciences New York University 251 Mercer St. New York, NY, 112,
More informationExponential tail inequalities for eigenvalues of random matrices
Exponential tail inequalities for eigenvalues of random matrices M. Ledoux Institut de Mathématiques de Toulouse, France exponential tail inequalities classical theme in probability and statistics quantify
More informationUniversality for random matrices and log-gases
Universality for random matrices and log-gases László Erdős IST, Austria Ludwig-Maximilians-Universität, Munich, Germany Encounters Between Discrete and Continuous Mathematics Eötvös Loránd University,
More informationUniversality of local spectral statistics of random matrices
Universality of local spectral statistics of random matrices László Erdős Ludwig-Maximilians-Universität, Munich, Germany CRM, Montreal, Mar 19, 2012 Joint with P. Bourgade, B. Schlein, H.T. Yau, and J.
More informationSimultaneous Gaussian quadrature for Angelesco systems
for Angelesco systems 1 KU Leuven, Belgium SANUM March 22, 2016 1 Joint work with Doron Lubinsky Introduced by C.F. Borges in 1994 Introduced by C.F. Borges in 1994 (goes back to Angelesco 1918). Introduced
More informationPolynomial Approximation: The Fourier System
Polynomial Approximation: The Fourier System Charles B. I. Chilaka CASA Seminar 17th October, 2007 Outline 1 Introduction and problem formulation 2 The continuous Fourier expansion 3 The discrete Fourier
More informationCOMPLEX HERMITE POLYNOMIALS: FROM THE SEMI-CIRCULAR LAW TO THE CIRCULAR LAW
Serials Publications www.serialspublications.com OMPLEX HERMITE POLYOMIALS: FROM THE SEMI-IRULAR LAW TO THE IRULAR LAW MIHEL LEDOUX Abstract. We study asymptotics of orthogonal polynomial measures of the
More informationFredholm determinant with the confluent hypergeometric kernel
Fredholm determinant with the confluent hypergeometric kernel J. Vasylevska joint work with I. Krasovsky Brunel University Dec 19, 2008 Problem Statement Let K s be the integral operator acting on L 2
More informationRandom Matrix: From Wigner to Quantum Chaos
Random Matrix: From Wigner to Quantum Chaos Horng-Tzer Yau Harvard University Joint work with P. Bourgade, L. Erdős, B. Schlein and J. Yin 1 Perhaps I am now too courageous when I try to guess the distribution
More informationLarge Deviations, Linear Statistics, and Scaling Limits for Mahler Ensemble of Complex Random Polynomials
Large Deviations, Linear Statistics, and Scaling Limits for Mahler Ensemble of Complex Random Polynomials Maxim L. Yattselev joint work with Christopher D. Sinclair International Conference on Approximation
More informationOrthogonal Polynomial Ensembles
Chater 11 Orthogonal Polynomial Ensembles 11.1 Orthogonal Polynomials of Scalar rgument Let wx) be a weight function on a real interval, or the unit circle, or generally on some curve in the comlex lane.
More informationFree probabilities and the large N limit, IV. Loop equations and all-order asymptotic expansions... Gaëtan Borot
Free probabilities and the large N limit, IV March 27th 2014 Loop equations and all-order asymptotic expansions... Gaëtan Borot MPIM Bonn & MIT based on joint works with Alice Guionnet, MIT Karol Kozlowski,
More information4 Uniform convergence
4 Uniform convergence In the last few sections we have seen several functions which have been defined via series or integrals. We now want to develop tools that will allow us to show that these functions
More informationOrthogonal polynomials with respect to generalized Jacobi measures. Tivadar Danka
Orthogonal polynomials with respect to generalized Jacobi measures Tivadar Danka A thesis submitted for the degree of Doctor of Philosophy Supervisor: Vilmos Totik Doctoral School in Mathematics and Computer
More informationKernel Method: Data Analysis with Positive Definite Kernels
Kernel Method: Data Analysis with Positive Definite Kernels 2. Positive Definite Kernel and Reproducing Kernel Hilbert Space Kenji Fukumizu The Institute of Statistical Mathematics. Graduate University
More informationORTHOGONAL POLYNOMIALS WITH EXPONENTIALLY DECAYING RECURSION COEFFICIENTS
ORTHOGONAL POLYNOMIALS WITH EXPONENTIALLY DECAYING RECURSION COEFFICIENTS BARRY SIMON* Dedicated to S. Molchanov on his 65th birthday Abstract. We review recent results on necessary and sufficient conditions
More informationA general framework for solving Riemann Hilbert problems numerically
A general framework for solving Riemann Hilbert problems numerically Sheehan Olver Oxford University Mathematical Institute 24 29 St Giles Oxford, UK Sheehan.Olver@sjc.ox.ac.uk Abstract A new, numerical
More informationA general framework for solving Riemann Hilbert problems numerically
Report no 0/5 A general framework for solving Riemann Hilbert problems numerically Sheehan Olver Abstract A new, numerical framework for the approximation of solutions to matrix-valued Riemann Hilbert
More informationChange of variable formulæ for regularizing slowly decaying and oscillatory Cauchy and Hilbert transforms
Change of variable formulæ for regularizing slowly decaying and oscillatory Cauchy and Hilbert transforms Sheehan Olver November 11, 2013 Abstract Formulæ are derived for expressing Cauchy and Hilbert
More informationRandom matrices and the Riemann zeros
Random matrices and the Riemann zeros Bruce Bartlett Talk given in Postgraduate Seminar on 5th March 2009 Abstract Random matrices and the Riemann zeta function came together during a chance meeting between
More informationCentral Limit Theorems for linear statistics for Biorthogonal Ensembles
Central Limit Theorems for linear statistics for Biorthogonal Ensembles Maurice Duits, Stockholm University Based on joint work with Jonathan Breuer (HUJI) Princeton, April 2, 2014 M. Duits (SU) CLT s
More informationSPRING 2006 PRELIMINARY EXAMINATION SOLUTIONS
SPRING 006 PRELIMINARY EXAMINATION SOLUTIONS 1A. Let G be the subgroup of the free abelian group Z 4 consisting of all integer vectors (x, y, z, w) such that x + 3y + 5z + 7w = 0. (a) Determine a linearly
More informationLectures 6 7 : Marchenko-Pastur Law
Fall 2009 MATH 833 Random Matrices B. Valkó Lectures 6 7 : Marchenko-Pastur Law Notes prepared by: A. Ganguly We will now turn our attention to rectangular matrices. Let X = (X 1, X 2,..., X n ) R p n
More informationQualifying Exams I, 2014 Spring
Qualifying Exams I, 2014 Spring 1. (Algebra) Let k = F q be a finite field with q elements. Count the number of monic irreducible polynomials of degree 12 over k. 2. (Algebraic Geometry) (a) Show that
More informationRANDOM MATRIX THEORY AND TOEPLITZ DETERMINANTS
RANDOM MATRIX THEORY AND TOEPLITZ DETERMINANTS David García-García May 13, 2016 Faculdade de Ciências da Universidade de Lisboa OVERVIEW Random Matrix Theory Introduction Matrix ensembles A sample computation:
More informationNUMERICAL SOLUTION OF DYSON BROWNIAN MOTION AND A SAMPLING SCHEME FOR INVARIANT MATRIX ENSEMBLES
NUMERICAL SOLUTION OF DYSON BROWNIAN MOTION AND A SAMPLING SCHEME FOR INVARIANT MATRIX ENSEMBLES XINGJIE HELEN LI AND GOVIND MENON Abstract. The Dyson Brownian Motion (DBM) describes the stochastic evolution
More information1 Math 241A-B Homework Problem List for F2015 and W2016
1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let
More informationNumerical Methods I Orthogonal Polynomials
Numerical Methods I Orthogonal Polynomials Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course G63.2010.001 / G22.2420-001, Fall 2010 Nov. 4th and 11th, 2010 A. Donev (Courant Institute)
More informationSeparation of Variables in Linear PDE: One-Dimensional Problems
Separation of Variables in Linear PDE: One-Dimensional Problems Now we apply the theory of Hilbert spaces to linear differential equations with partial derivatives (PDE). We start with a particular example,
More informationNon-perturbative effects in ABJM theory
October 9, 2015 Non-perturbative effects in ABJM theory 1 / 43 Outline 1. Non-perturbative effects 1.1 General aspects 1.2 Non-perturbative aspects of string theory 1.3 Non-perturbative effects in M-theory
More informationINTEGRATION WORKSHOP 2004 COMPLEX ANALYSIS EXERCISES
INTEGRATION WORKSHOP 2004 COMPLEX ANALYSIS EXERCISES PHILIP FOTH 1. Cauchy s Formula and Cauchy s Theorem 1. Suppose that γ is a piecewise smooth positively ( counterclockwise ) oriented simple closed
More informationDeterminantal point processes and random matrix theory in a nutshell
Determinantal point processes and random matrix theory in a nutshell part I Manuela Girotti based on M. Girotti s PhD thesis and A. Kuijlaars notes from Les Houches Winter School 202 Contents Point Processes
More informationMultiple Orthogonal Polynomials
Summer school on OPSF, University of Kent 26 30 June, 2017 Plan of the course Plan of the course lecture 1: Definitions + basic properties Plan of the course lecture 1: Definitions + basic properties lecture
More informationNotation. General. Notation Description See. Sets, Functions, and Spaces. a b & a b The minimum and the maximum of a and b
Notation General Notation Description See a b & a b The minimum and the maximum of a and b a + & a f S u The non-negative part, a 0, and non-positive part, (a 0) of a R The restriction of the function
More informationTriangular matrices and biorthogonal ensembles
/26 Triangular matrices and biorthogonal ensembles Dimitris Cheliotis Department of Mathematics University of Athens UK Easter Probability Meeting April 8, 206 2/26 Special densities on R n Example. n
More informationComplex Analysis, Stein and Shakarchi Meromorphic Functions and the Logarithm
Complex Analysis, Stein and Shakarchi Chapter 3 Meromorphic Functions and the Logarithm Yung-Hsiang Huang 217.11.5 Exercises 1. From the identity sin πz = eiπz e iπz 2i, it s easy to show its zeros are
More informationA quick derivation of the loop equations for random matrices
Probability, Geometry and Integrable Systems MSRI Publications Volume 55, 2007 A quick derivation of the loop equations for random matrices. M. ERCOLAI AD K. D. T-R MCLAUGHLI ABSTRACT. The loop equations
More informationPhysics 342 Lecture 17. Midterm I Recap. Lecture 17. Physics 342 Quantum Mechanics I
Physics 342 Lecture 17 Midterm I Recap Lecture 17 Physics 342 Quantum Mechanics I Monday, March 1th, 28 17.1 Introduction In the context of the first midterm, there are a few points I d like to make about
More informationUniversality in Numerical Computations with Random Data. Case Studies.
Universality in Numerical Computations with Random Data. Case Studies. Percy Deift and Thomas Trogdon Courant Institute of Mathematical Sciences Govind Menon Brown University Sheehan Olver The University
More informationarxiv: v1 [math.cv] 30 Jun 2008
EQUIDISTRIBUTION OF FEKETE POINTS ON COMPLEX MANIFOLDS arxiv:0807.0035v1 [math.cv] 30 Jun 2008 ROBERT BERMAN, SÉBASTIEN BOUCKSOM Abstract. We prove the several variable version of a classical equidistribution
More information1. The COMPLEX PLANE AND ELEMENTARY FUNCTIONS: Complex numbers; stereographic projection; simple and multiple connectivity, elementary functions.
Complex Analysis Qualifying Examination 1 The COMPLEX PLANE AND ELEMENTARY FUNCTIONS: Complex numbers; stereographic projection; simple and multiple connectivity, elementary functions 2 ANALYTIC FUNCTIONS:
More informationSection 9 Variational Method. Page 492
Section 9 Variational Method Page 492 Page 493 Lecture 27: The Variational Method Date Given: 2008/12/03 Date Revised: 2008/12/03 Derivation Section 9.1 Variational Method: Derivation Page 494 Motivation
More informationJim Lambers MAT 460/560 Fall Semester Practice Final Exam
Jim Lambers MAT 460/560 Fall Semester 2009-10 Practice Final Exam 1. Let f(x) = sin 2x + cos 2x. (a) Write down the 2nd Taylor polynomial P 2 (x) of f(x) centered around x 0 = 0. (b) Write down the corresponding
More informationDouble contour integral formulas for the sum of GUE and one matrix model
Double contour integral formulas for the sum of GUE and one matrix model Based on arxiv:1608.05870 with Tom Claeys, Arno Kuijlaars, and Karl Liechty Dong Wang National University of Singapore Workshop
More informationGENERALIZED STIELTJES POLYNOMIALS AND RATIONAL GAUSS-KRONROD QUADRATURE
GENERALIZED STIELTJES POLYNOMIALS AND RATIONAL GAUSS-KRONROD QUADRATURE M. BELLO HERNÁNDEZ, B. DE LA CALLE YSERN, AND G. LÓPEZ LAGOMASINO Abstract. Generalized Stieltjes polynomials are introduced and
More informationElectronic Transactions on Numerical Analysis Volume 50, 2018
Electronic Transactions on Numerical Analysis Volume 50, 2018 Contents 1 The Lanczos algorithm and complex Gauss quadrature. Stefano Pozza, Miroslav S. Pranić, and Zdeněk Strakoš. Gauss quadrature can
More informationCommutative Banach algebras 79
8. Commutative Banach algebras In this chapter, we analyze commutative Banach algebras in greater detail. So we always assume that xy = yx for all x, y A here. Definition 8.1. Let A be a (commutative)
More informationChapter 2 Orthogonal Polynomials and Weighted Polynomial Approximation
Chapter 2 Orthogonal Polynomials and Weighted Polynomial Approximation 2.1 Orthogonal Systems and Polynomials 2.1.1 Inner Product Space and Orthogonal Systems Suppose that X is a complex linear space of
More informationNumerical solution of Riemann Hilbert problems: Painlevé II
Report no 09/9 Numerical solution of Riemann Hilbert problems: Painlevé II Sheehan Olver Abstract We describe a new, spectrally accurate method for solving matrixvalued Riemann Hilbert problems numerically
More informationDeterminantal point processes and random matrix theory in a nutshell
Determinantal point processes and random matrix theory in a nutshell part II Manuela Girotti based on M. Girotti s PhD thesis, A. Kuijlaars notes from Les Houches Winter School 202 and B. Eynard s notes
More informationPhysics 221A Fall 1996 Notes 16 Bloch s Theorem and Band Structure in One Dimension
Physics 221A Fall 1996 Notes 16 Bloch s Theorem and Band Structure in One Dimension In these notes we examine Bloch s theorem and band structure in problems with periodic potentials, as a part of our survey
More informationMatrix Functions and their Approximation by. Polynomial methods
[ 1 / 48 ] University of Cyprus Matrix Functions and their Approximation by Polynomial Methods Stefan Güttel stefan@guettel.com Nicosia, 7th April 2006 Matrix functions Polynomial methods [ 2 / 48 ] University
More informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
More informationDeterminantal Processes And The IID Gaussian Power Series
Determinantal Processes And The IID Gaussian Power Series Yuval Peres U.C. Berkeley Talk based on work joint with: J. Ben Hough Manjunath Krishnapur Bálint Virág 1 Samples of translation invariant point
More informationF (z) =f(z). f(z) = a n (z z 0 ) n. F (z) = a n (z z 0 ) n
6 Chapter 2. CAUCHY S THEOREM AND ITS APPLICATIONS Theorem 5.6 (Schwarz reflection principle) Suppose that f is a holomorphic function in Ω + that extends continuously to I and such that f is real-valued
More informationNumerical Analysis Preliminary Exam 10.00am 1.00pm, January 19, 2018
Numerical Analysis Preliminary Exam 0.00am.00pm, January 9, 208 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationMeasures, orthogonal polynomials, and continued fractions. Michael Anshelevich
Measures, orthogonal polynomials, and continued fractions Michael Anshelevich November 7, 2008 MEASURES AND ORTHOGONAL POLYNOMIALS. µ a positive measure on R. A linear functional µ[p ] = P (x) dµ(x), µ
More informationMath 259: Introduction to Analytic Number Theory Functions of finite order: product formula and logarithmic derivative
Math 259: Introduction to Analytic Number Theory Functions of finite order: product formula and logarithmic derivative This chapter is another review of standard material in complex analysis. See for instance
More informationBulk scaling limits, open questions
Bulk scaling limits, open questions Based on: Continuum limits of random matrices and the Brownian carousel B. Valkó, B. Virág. Inventiones (2009). Eigenvalue statistics for CMV matrices: from Poisson
More informationRANDOM HERMITIAN MATRICES AND GAUSSIAN MULTIPLICATIVE CHAOS NATHANAËL BERESTYCKI, CHRISTIAN WEBB, AND MO DICK WONG
RANDOM HERMITIAN MATRICES AND GAUSSIAN MULTIPLICATIVE CHAOS NATHANAËL BERESTYCKI, CHRISTIAN WEBB, AND MO DICK WONG Abstract We prove that when suitably normalized, small enough powers of the absolute value
More informationFrom the mesoscopic to microscopic scale in random matrix theory
From the mesoscopic to microscopic scale in random matrix theory (fixed energy universality for random spectra) With L. Erdős, H.-T. Yau, J. Yin Introduction A spacially confined quantum mechanical system
More informationMeasures, orthogonal polynomials, and continued fractions. Michael Anshelevich
Measures, orthogonal polynomials, and continued fractions Michael Anshelevich November 7, 2008 MEASURES AND ORTHOGONAL POLYNOMIALS. MEASURES AND ORTHOGONAL POLYNOMIALS. µ a positive measure on R. 2 MEASURES
More informationAiry and Pearcey Processes
Airy and Pearcey Processes Craig A. Tracy UC Davis Probability, Geometry and Integrable Systems MSRI December 2005 1 Probability Space: (Ω, Pr, F): Random Matrix Models Gaussian Orthogonal Ensemble (GOE,
More informationA Remark on Hypercontractivity and Tail Inequalities for the Largest Eigenvalues of Random Matrices
A Remark on Hypercontractivity and Tail Inequalities for the Largest Eigenvalues of Random Matrices Michel Ledoux Institut de Mathématiques, Université Paul Sabatier, 31062 Toulouse, France E-mail: ledoux@math.ups-tlse.fr
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationContents. I Basic Methods 13
Preface xiii 1 Introduction 1 I Basic Methods 13 2 Convergent and Divergent Series 15 2.1 Introduction... 15 2.1.1 Power series: First steps... 15 2.1.2 Further practical aspects... 17 2.2 Differential
More informationQuantum Harmonic Oscillator
Quantum Harmonic Oscillator Chapter 13 P. J. Grandinetti Chem. 4300 Oct 20, 2017 P. J. Grandinetti (Chem. 4300) Quantum Harmonic Oscillator Oct 20, 2017 1 / 26 Kinetic and Potential Energy Operators Harmonic
More informationz b k P k p k (z), (z a) f (n 1) (a) 2 (n 1)! (z a)n 1 +f n (z)(z a) n, where f n (z) = 1 C
. Representations of Meromorphic Functions There are two natural ways to represent a rational function. One is to express it as a quotient of two polynomials, the other is to use partial fractions. The
More informationBernstein s inequality and Nikolsky s inequality for R d
Bernstein s inequality and Nikolsky s inequality for d Jordan Bell jordan.bell@gmail.com Department of athematics University of Toronto February 6 25 Complex Borel measures and the Fourier transform Let
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationSynopsis of Complex Analysis. Ryan D. Reece
Synopsis of Complex Analysis Ryan D. Reece December 7, 2006 Chapter Complex Numbers. The Parts of a Complex Number A complex number, z, is an ordered pair of real numbers similar to the points in the real
More informationarxiv:math/ v2 [math.pr] 10 Aug 2005
arxiv:math/0403090v2 [math.pr] 10 Aug 2005 ORTHOGOAL POLYOMIAL ESEMBLES I PROBABILITY THEORY By Wolfgang König 1 10 August, 2005 Abstract: We survey a number of models from physics, statistical mechanics,
More informationNonlinear steepest descent and the numerical solution of Riemann Hilbert problems
Nonlinear steepest descent and the numerical solution of Riemann Hilbert problems Sheehan Olver and Thomas Trogdon 2 School of Mathematics and Statistics The University of Sydney NSW 2006, Australia 2
More informationReferences 167. dx n x2 =2
References 1. G. Akemann, J. Baik, P. Di Francesco (editors), The Oxford Handbook of Random Matrix Theory, Oxford University Press, Oxford, 2011. 2. G. Anderson, A. Guionnet, O. Zeitouni, An Introduction
More informationSET-A U/2015/16/II/A
. n x dx 3 n n is equal to : 0 4. Let f n (x) = x n n, for each n. Then the series (A) 3 n n (B) n 3 n (n ) (C) n n (n ) (D) n 3 n (n ) f n is n (A) Uniformly convergent on [0, ] (B) Po intwise convergent
More informationPhysics 137A Quantum Mechanics Fall 2012 Midterm II - Solutions
Physics 37A Quantum Mechanics Fall 0 Midterm II - Solutions These are the solutions to the exam given to Lecture Problem [5 points] Consider a particle with mass m charge q in a simple harmonic oscillator
More informationNumerical Methods for Random Matrices
Numerical Methods for Random Matrices MIT 18.95 IAP Lecture Series Per-Olof Persson (persson@mit.edu) January 23, 26 1.9.8.7 Random Matrix Eigenvalue Distribution.7.6.5 β=1 β=2 β=4 Probability.6.5.4.4
More informationPart II Probability and Measure
Part II Probability and Measure Theorems Based on lectures by J. Miller Notes taken by Dexter Chua Michaelmas 2016 These notes are not endorsed by the lecturers, and I have modified them (often significantly)
More informationFast inverse transform sampling in one and two dimensions
Fast inverse transform sampling in one and two dimensions Sheehan Olver Alex Townsend Abstract We develop a computationally efficient and robust algorithm for generating pseudo-random samples from a broad
More informationComputation of Rational Szegő-Lobatto Quadrature Formulas
Computation of Rational Szegő-Lobatto Quadrature Formulas Joint work with: A. Bultheel (Belgium) E. Hendriksen (The Netherlands) O. Njastad (Norway) P. GONZÁLEZ-VERA Departament of Mathematical Analysis.
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationValerio Cappellini. References
CETER FOR THEORETICAL PHYSICS OF THE POLISH ACADEMY OF SCIECES WARSAW, POLAD RADOM DESITY MATRICES AD THEIR DETERMIATS 4 30 SEPTEMBER 5 TH SFB TR 1 MEETIG OF 006 I PRZEGORZAłY KRAKÓW Valerio Cappellini
More informationApplied/Numerical Analysis Qualifying Exam
Applied/Numerical Analysis Qualifying Exam August 9, 212 Cover Sheet Applied Analysis Part Policy on misprints: The qualifying exam committee tries to proofread exams as carefully as possible. Nevertheless,
More informationLucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche
Scuola di Dottorato THE WAVE EQUATION Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Lucio Demeio - DIISM wave equation 1 / 44 1 The Vibrating String Equation 2 Second
More informationMATH 311: COMPLEX ANALYSIS CONTOUR INTEGRALS LECTURE
MATH 3: COMPLEX ANALYSIS CONTOUR INTEGRALS LECTURE Recall the Residue Theorem: Let be a simple closed loop, traversed counterclockwise. Let f be a function that is analytic on and meromorphic inside. Then
More informationyou expect to encounter difficulties when trying to solve A x = b? 4. A composite quadrature rule has error associated with it in the following form
Qualifying exam for numerical analysis (Spring 2017) Show your work for full credit. If you are unable to solve some part, attempt the subsequent parts. 1. Consider the following finite difference: f (0)
More informationQualification Exam: Mathematical Methods
Qualification Exam: Mathematical Methods Name:, QEID#41534189: August, 218 Qualification Exam QEID#41534189 2 1 Mathematical Methods I Problem 1. ID:MM-1-2 Solve the differential equation dy + y = sin
More informationGlobal Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations
Global Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations C. David Levermore Department of Mathematics and Institute for Physical Science and Technology
More informationNumerical Evaluation of Standard Distributions in Random Matrix Theory
Numerical Evaluation of Standard Distributions in Random Matrix Theory A Review of Folkmar Bornemann s MATLAB Package and Paper Matt Redmond Department of Mathematics Massachusetts Institute of Technology
More informationElements of Positive Definite Kernel and Reproducing Kernel Hilbert Space
Elements of Positive Definite Kernel and Reproducing Kernel Hilbert Space Statistical Inference with Reproducing Kernel Hilbert Space Kenji Fukumizu Institute of Statistical Mathematics, ROIS Department
More informationLecture Note 3: Polynomial Interpolation. Xiaoqun Zhang Shanghai Jiao Tong University
Lecture Note 3: Polynomial Interpolation Xiaoqun Zhang Shanghai Jiao Tong University Last updated: October 24, 2013 1.1 Introduction We first look at some examples. Lookup table for f(x) = 2 π x 0 e x2
More informationHilbert Space Problems
Hilbert Space Problems Prescribed books for problems. ) Hilbert Spaces, Wavelets, Generalized Functions and Modern Quantum Mechanics by Willi-Hans Steeb Kluwer Academic Publishers, 998 ISBN -7923-523-9
More informationStatistical Interpretation
Physics 342 Lecture 15 Statistical Interpretation Lecture 15 Physics 342 Quantum Mechanics I Friday, February 29th, 2008 Quantum mechanics is a theory of probability densities given that we now have an
More informationRiemann Mapping Theorem (4/10-4/15)
Math 752 Spring 2015 Riemann Mapping Theorem (4/10-4/15) Definition 1. A class F of continuous functions defined on an open set G is called a normal family if every sequence of elements in F contains a
More information