GIG random matrices and a Yang-Baxter extension of the Matsumoto-Yor property
Gérard Letac, Mauro Piccioni, Jacek Wesołowski
TL;DR
This work extends the Matsumoto–Yor independence property to the matrix variate setting by employing a matrix MGIG framework on the SPD cone $Ω_+$. Through a functional equation approach linked to the Jacobian of a matrix Yang–Baxter map $φ^{(α,β)}$, the authors prove that the MGIG law is characterized by an independence property induced by $φ^{(α,β)}$, and that the corresponding matrix map is a parametric Yang–Baxter map. They explicitly identify the density parameters $a,b$ and the common scalar $λ$ by analyzing log-densities, derivatives, and scaling relations, culminating in a full characterization: $X∼MGIG(λ,α a,b)$, $Y∼MGIG(-λ,a,β b)$, $U∼MGIG(λ,α b,a)$, and $V∼MGIG(-λ,b,β a)$. The results bridge probabilistic independence properties with integrable systems in the matrix setting and open paths to symmetric-cone extensions and broader noncommutative generalizations.
Abstract
Sasada and Uozumi, \cite{SasUoz2024}, identified independence preserving $[2:2]$ quadrirational parametric Yang-Baxter maps, see \eqref{YBEQ}, on $(0,\infty)$. In particular, the map denoted there by $H_{III,B}^{(α,β)}$, see \eqref{CS}, was connected to the independence preserving property of the GIG distributions on $(0,\infty)$. Remarkably, the property appears also naturally in probabilistic integrable models of discrete Korteweg de Vries type, as observed by Croydon and Sasada, \cite{CroSas2020}. In the case of $(α,β)=(1,0)$ the independence reduces to the classical Matsumoto-Yor property, \cite{MatYor2001}. In \cite{LetWes2024} we proposed an extension of $H_{III,B}^{(α,β)}$ to a map on the cone of symmetric positive definite matrices of a fixed dimension, showing that such extended map preserves independence of GIG random matrices. In the present paper we prove two results: (i) the matrix GIG distributions are characterized by the independence property governed by this map; (ii) the matrix variate extension of $H_{III,B}^{(α,β)}$ we use, is a parametric Yang-Baxter map.
