Integrable structure of products of finite complex Ginibre random matrices

We consider the squared singular values of the product of [Formula presented] standard complex Gaussian matrices. Since the squared singular values form a determinantal point process with a particular Meijer G-function kernel, the gap probabilities are given by a Fredholm determinant based on this kernel. It was shown by Strahov (2014) that a hard edge scaling limit of the gap probabilities is described by Hamiltonian differential equations which can be formulated as an isomonodromic deformation system similar to the theory of the Kyoto school. We generalize this result to the case of finite matrices by first finding a representation of the finite kernel in integrable form. As a result we obtain the Hamiltonian structure for finite size matrices and formulate it in terms of a [Formula presented] matrix Schlesinger system. The case [Formula presented] reproduces the Tracy and Widom theory which results in the Painlevé V equation for the [Formula presented] gap probability. Some integrals of motion for [Formula presented] are identified, and a coupled system of differential equations in two unknowns is presented which uniquely determines the corresponding [Formula presented] gap probability.