Two-periodic Aztec diamond

Similar documents
Central Limit Theorems for linear statistics for Biorthogonal Ensembles

Universality of distribution functions in random matrix theory Arno Kuijlaars Katholieke Universiteit Leuven, Belgium

Domino tilings, non-intersecting Random Motions and Critical Processes

Proof of a generalization of Aztec diamond theorem

Uniformly Random Lozenge Tilings of Polygons on the Triangular Lattice

OPSF, Random Matrices and Riemann-Hilbert problems

The Hermitian two matrix model with an even quartic potential. Maurice Duits Arno B.J. Kuijlaars Man Yue Mo

Rectangular Young tableaux and the Jacobi ensemble

Tiling expression of minors

Airy and Pearcey Processes

RANDOM MATRIX THEORY AND TOEPLITZ DETERMINANTS

Correlation function of Schur process with application to local geometry of a random 3-dimensional Young diagram

Double contour integral formulas for the sum of GUE and one matrix model

OPSF, Random Matrices and Riemann-Hilbert problems

Statistical Mechanics and Combinatorics : Lecture IV

Domino tilings, non-intersecting Random Motions and Integrable Systems

Statistical mechanics via asymptotics of symmetric functions

Multiple orthogonal polynomials. Bessel weights

Schur processes and dimer models

Determinantal Processes And The IID Gaussian Power Series

Tacnode GUE-minor Processes and Double Aztec Diamonds

Properties of orthogonal polynomials

Multiple Orthogonal Polynomials

Enumeration of domino tilings of a double Aztec rectangle

arxiv: v1 [math-ph] 30 Jun 2011

Fredholm determinant with the confluent hypergeometric kernel

arxiv: v2 [math-ph] 20 Dec 2015

CORRELATION FUNCTION OF SCHUR PROCESS WITH APPLICATION TO LOCAL GEOMETRY OF A RANDOM 3-DIMENSIONAL YOUNG DIAGRAM

Correlations for the Novak process

Definition We say that a topological manifold X is C p if there is an atlas such that the transition functions are C p.

Discrete Toeplitz Determinants and their Applications

1 Unitary representations of the Virasoro algebra

Near extreme eigenvalues and the first gap of Hermitian random matrices

ON GENERALIZED LAGUERRE POLYNOMIALS WITH REAL AND COMPLEX PARAMETER.

Maximal height of non-intersecting Brownian motions

Cointegrated VAR s. Eduardo Rossi University of Pavia. November Rossi Cointegrated VAR s Fin. Econometrics / 31

Riesz bases, Meyer s quasicrystals, and bounded remainder sets

limit shapes beyond dimers

The Arctic Circle re-revisited

Random tilings and Markov chains for interlacing particles

Painlevé VI and Hankel determinants for the generalized Jacobi weight

Markov operators, classical orthogonal polynomial ensembles, and random matrices

Complex Analysis Important Concepts

A b r i l l i a n t young chemist, T h u r e Wagelius of N e w Y o r k, ac. himself with eth

Central Limit Theorem for discrete log gases

Beyond the Gaussian universality class

Determinantal point processes and random matrix theory in a nutshell

Large Deviations, Linear Statistics, and Scaling Limits for Mahler Ensemble of Complex Random Polynomials

Differential Equations for Dyson Processes

Orthogonal polynomials with respect to generalized Jacobi measures. Tivadar Danka

Zeros and ratio asymptotics for matrix orthogonal polynomials

ENUMERATION OF DOMINO TILINGS OF AN AZTEC RECTANGLE WITH BOUNDARY DEFECTS

ORTHOGONAL POLYNOMIALS

Möbius Transformation

Basic Concepts in Linear Algebra

RHOMBUS TILINGS OF A HEXAGON WITH TWO TRIANGLES MISSING ON THE SYMMETRY AXIS

Bergman shift operators for Jordan domains

1 Determinants. 1.1 Determinant

MATRIX INTEGRALS AND MAP ENUMERATION 2

Random Matrices and Determinantal Process

LECTURE 16: LIE GROUPS AND THEIR LIE ALGEBRAS. 1. Lie groups

Refined Cauchy/Littlewood identities and partition functions of the six-vertex model

The Pfaan Schur Process

Lecture Notes Introduction to Cluster Algebra

Gaussian Free Field in (self-adjoint) random matrices and random surfaces. Alexei Borodin

CALCULUS ON MANIFOLDS. 1. Riemannian manifolds Recall that for any smooth manifold M, dim M = n, the union T M =

Hyperbolic Geometry on Geometric Surfaces

Department of Mathematics, University of California, Berkeley. GRADUATE PRELIMINARY EXAMINATION, Part A Fall Semester 2014

DOMINO TILING. Contents 1. Introduction 1 2. Rectangular Grids 2 Acknowledgments 10 References 10

Compact symetric bilinear forms

Dual Temperley Lieb basis, Quantum Weingarten and a conjecture of Jones

Spectra of Semidirect Products of Cyclic Groups

INTEGRATION WORKSHOP 2004 COMPLEX ANALYSIS EXERCISES

BERGMAN KERNEL ON COMPACT KÄHLER MANIFOLDS

Taylor Expansions with Finite Radius of Convergence Work of Robert Jentzsch, 1917

Random Matrices and Multivariate Statistical Analysis

HERMITE-PADÉ APPROXIMANTS FOR A PAIR OF CAUCHY TRANSFORMS WITH OVERLAPPING SYMMETRIC SUPPORTS

Kernel Methods. Machine Learning A W VO

Explicit Computations of the Frozen Boundaries of Rhombus Tilings of Polygonal Domains

ELE/MCE 503 Linear Algebra Facts Fall 2018

Random walks in Beta random environment

ETNA Kent State University

K. T.-R. McLaughlin and P. D. Miller. 1 Introduction. 1.1 Asymptotic analysis of Riemann-Hilbert problems

arxiv: v2 [math.pr] 28 May 2013

MATRIX KERNELS FOR THE GAUSSIAN ORTHOGONAL AND SYMPLECTIC ENSEMBLES

Determinantal point processes and random matrix theory in a nutshell

INTRODUCTION TO PADÉ APPROXIMANTS. E. B. Saff Center for Constructive Approximation

Commutants of Finite Blaschke Product. Multiplication Operators on Hilbert Spaces of Analytic Functions

Determinantal Probability Measures. by Russell Lyons (Indiana University)

Phys 201. Matrices and Determinants

Problem 1A. Use residues to compute. dx x

Centrum voor Wiskunde en Informatica

Cluster Variables and Perfect Matchings of Subgraphs of the dp 3 Lattice

Inverses of regular Hessenberg matrices

A factorisation theorem for the number of rhombus tilings of a hexagon with triangular holes

From Bernstein approximation to Zauner s conjecture

On rational approximation of algebraic functions. Julius Borcea. Rikard Bøgvad & Boris Shapiro

A Note on Tiling under Tomographic Constraints

Lecture Notes in Linear Algebra

Orthogonal Polynomial Ensembles

Transcription:

Two-periodic Aztec diamond Arno Kuijlaars (KU Leuven) joint work with Maurice Duits (KTH Stockholm) Optimal and Random Point Configurations ICERM, Providence, RI, U.S.A., 27 February 2018

Outline 1. Aztec diamond 2. The model and main result 3. Non-intersecting paths 4. Matrix Valued Orthogonal Polynomials (MVOP) 5. Analysis of RH problem 6. Saddle point analysis 7. Periodic tilings of a hexagon

1. Aztec diamond

Aztec diamond North West East South

Tiling of an Aztec diamond North West East South Tiling with 2 1 and 1 2 rectangles (dominos) Four types of dominos

Large random tiling Deterministic pattern near corners Solid region or Frozen region Disorder in the middle Liquid region Boundary curve Arctic circle

Recent development Two-periodic weighting Chhita, Johansson (2016) Beffara, Chhita, Johansson (2018 to appear)

Two-periodic weights A new phase within the liquid region: gas region

Phase diagram solid solid liquid gas solid solid

2. The model and main result

Two periodic weights a b b b a a a a b b b a a a a b b b a Aztec diamond of size 2N b Weight w(t ) of a tiling T is the product of the weights of dominos Partition function Z N = T w(t ) Probability for T Prob(T ) = w(t ) Z N

Equivalent weights α = a 2 and β = b 2 α β α North and East dominos have weight 1 β α α β Without loss of generality β β α αβ = 1 and α 1 Since North dominos have weight 1, we can transfer the weights to non-intersecting paths.

Particles in West and South dominos Particles along diagonal lines are interlacing Positions of particles are random in the two-periodic Aztec diamond. Structure of determinantal point process We found explicit formula for kernel K N using matrix valued orthogonal polynomials (MVOP).

Coordinates n runs from 0 to 2N 1 m runs from 0 to 2N

Formula for correlation kernel THEOREM 1 Assume N is even and m + n and m + n are even. ( ) KN (m, n; m, n ) K N (m, n + 1; m, n ) K N (m, n; m, n + 1) K N (m, n + 1; m, n + 1) = χ m>m A m m (z)z m m+n n dz 2 2πi γ 0,1 z + 1 dz dw z N m n 2 (z 1) N A N m (w)f (w)a N+m (z) (2πi) 2 γ 0,1 z γ 1 z w w N m n 2 (w 1) N where A(z) = 1 ( ) 2αz α(z + 1) z 1 βz(z + 1) 2βz F (z) = 1 ( ) 2 I 1 (α β)z α(z + 1) 2 + 2 z(z + α 2 )(z + β 2 ) βz(z + 1) (α β)z

3. Non-intersecting paths

Non-intersecting paths Line segments on West, East and South dominos North West East South

Double Aztec diamond 2N particles along each diagonal line

Non-intersecting paths on a graph Paths are transformed to fit on a graph 6 5 4 3 2 1 0-1 -2 0 0.5 1 1.5 2 2.5 3 3.5 4

Weights on the graph 6 β β β β β 5 β β β β β α α α α α 4 α α α α α β β β β β 3 β β β β β α α α α α 2 α α α α α β β β β β 1 β β β β β α α α α α 0 α α α α α β β β β β -1 β β β β β α α α α α -2 α α α α α 0 0.5 1 1.5 2 2.5 3 3.5 4

Weights on non-intersecting paths Any tiling of double Aztec diamond is equivalent to system (P 0,..., P 2N 1 ) of 2N non-intersecting paths P j is path on the graph from (0, j) to (2N, j), P i is vertex disjoint from P j if i j.

Transitions and LGV theorem There are 2N + 1 levels, 0, 1,..., 2N. Transition from level m to level m > m T m,m (x, y) = w(p), x, y Z P:(m,x) (m,y)

Transitions and LGV theorem There are 2N + 1 levels, 0, 1,..., 2N. Transition from level m to level m > m T m,m (x, y) = w(p), x, y Z P:(m,x) (m,y) Lindström-Gessel-Viennot theorem Probability that paths at level m are at positions x (m) 0 < x (m) 1 < < x (m) 2N 1 : 1 Z N [ ] 2N 1 [ ] 2N 1 det T 0,m (i, x (m) k ) det T m,2n (x (m) k, j) i,k=0 k,j=0 Lindström (1973) Gessel-Viennot (1985)

Determinantal point process Corollary: The positions at level m are determinantal with kernel K N,m (x, y) = 2N 1 i,j=0 T 0,m (i, x) [ G t] i,j T m,2n(y, j) where G = [T 0,2N (i, j)] 2N 1 i,j=0 Multi-level extension is known as Eynard-Mehta theorem.

Block Toeplitz matrices In our case: Transition matrices are 2 periodic T (x + 2, y + 2) = T (x, y) Block Toeplitz matrices, infinite in both directions, with block symbol A(z) = B j z j j=............. B0 B 1.. if T =... B 1 B 0 B 1....... B 1 B 0...........

Double contour integral formula THEOREM 2: Suppose transition matrices are 2-periodic. Then ( ) KN,m (2x, 2y) K N,m (2x + 1, 2y) K N,m (2x, 2y + 1) K N,m (2x + 1, 2y + 1) = 1 w y A (2πi) 2 m,2n (w)r N (w, z)a 0,m (z) z x+1 w dzdw N γ γ A m,2n and A 0,m are block symbols for the transition matrices T m,2n and T 0,m. R N (w, z) is a reproducing kernel for matrix valued polynomials.

4. Matrix Valued Orthogonal Polynomials (MVOP)

MVOP Matrix valued polynomial of degree j, P j (z) = j C i z i i=0 each C i is d d matrix, det C j 0 W (z) is d d matrix valued weight Orthogonality 1 P j (z)w (z)p t 2πi k(z) dz = H j δ j,k γ

Reproducing kernel N 1 R N (w, z) = Pj t (w)h 1 j P j (z) j=0 is reproducing kernel for matrix polynomials of degree N 1 If Q has degree N 1, then 1 Q(w)W (w)r N (w, z)dw = Q(z) 2πi γ There is a Christoffel-Darboux formula for R N and a Riemann Hilbert problem

Riemann-Hilbert problem Y : C \ γ C 2d 2d satisfies Y is analytic, ( ) Id W Y + = Y on γ, 0 d I d ( ) z Y (z) = (I 2d + O(z 1 )) N I d 0 d 0 d z N I d as z. Grünbaum, de la Iglesia, Martínez-Finkelshtein (2011)

Solution of RH problem Unique solution (provided P N uniquely exists) is 1 P N (s)w (s) P N (z) ds Y (z) = 2πi γ s z 1 Q N 1 (s)w (s) Q N 1 (z) ds 2πi s z where P N is monic MVOP of degree N and Q N 1 = H 1 N 1 P N 1 has degree N 1 γ

Solution of RH problem Unique solution (provided P N uniquely exists) is 1 P N (s)w (s) P N (z) ds Y (z) = 2πi γ s z 1 Q N 1 (s)w (s) Q N 1 (z) ds 2πi s z where P N is monic MVOP of degree N and Q N 1 = H 1 N 1 P N 1 has degree N 1 Christoffel Darboux formula R N (w, z) = 1 ( ) ( ) 0d I d Y 1 Id (w)y (z) z w 0 d γ Delvaux (2010)

Our case of interest Weight matrix in special case of two periodic Aztec diamond is W N (z), with ( ) 1 (z + 1) W (z) = 2 + 4α 2 z 2α(α + β)(z + 1) (z 1) 2 2β(α + β)z(z + 1) (z + 1) 2 + 4β 2 z No symmetry in W. Existence and uniqueness of MVOP are not immediate.

Our case of interest Weight matrix in special case of two periodic Aztec diamond is W N (z), with ( ) 1 (z + 1) W (z) = 2 + 4α 2 z 2α(α + β)(z + 1) (z 1) 2 2β(α + β)z(z + 1) (z + 1) 2 + 4β 2 z No symmetry in W. Existence and uniqueness of MVOP are not immediate. Scalar valued analogue Weight ( z+1 N z 1) on circle around z = 1 and OPs are Jacobi polynomials P ( N,N) j (z) with nonstandard parameters

5. Analysis of RH problem

Surprise Steepest descent analysis of RH problem leads to explicit formula RH problem is solved in terms of contour integrals. For example: MVOP is P N (z) = (z 1) N W N/2 W N/2 (z), if N is even.

Surprise Steepest descent analysis of RH problem leads to explicit formula RH problem is solved in terms of contour integrals. For example: MVOP is P N (z) = (z 1) N W N/2 W N/2 (z), if N is even. It leads to proof of THEOREM 1

6. Saddle point analysis

Asymptotic analysis Saddle point analysis on the double contour integral 1 dz (2πi) 2 γ 0,1 z γ 1 when N dw z w z N 2x 2 (z 1) N w N 2y 2 (w 1) N AN m (w)f (w)a N+m (z) m, x, y scale with N in such a way that m (1 + ξ 1 )N, x, y (1 + ξ 1+ξ 2 2 )N Saddle points are critical points of 2 log(z 1) (1 + ξ 2 ) log z + ξ 1 log λ(z) where λ(z) is an eigenvalue of W (z) = A2 (z). z

Saddle point analysis Let 1 < ξ 1, ξ 2 < 1. There are always four saddle points, depending on ξ 1, ξ 2, and they lie on the Riemann surface for y 2 = z(z + α 2 )(z + β 2 ) (genus one) with branch points α 2 < β 2 < 0 and infinity. At least two saddles are in z [ α 2, β 2 ].

Classification of phases Location of other two saddles determines the phase. Two saddles are in [0, ): solid phase Two saddles are in C \ ([ α 2, β 2 ] [0, )): liquid phase All four saddles are in [ α 2, β 2 ]: gas phase Transitions between phases occur when saddles coalesce.

Phase diagram solid solid liquid liquid gas solid liquid liquid solid

7. Periodic tilings of a hexagon

Tiling of a hexagon Lozenge tiling of a regular hexagon Also admits a non-intersecting path formulation

Large random tiling

Two periodic tiling of a hexagon Ongoing work with Charlier, Duits, and Lenells

Thanks Thank you for your attention