Evaluation of integrals for the emptiness formation probability in the square-ice model

We study the emptiness formation probability (EFP) in the six-vertex model with domain wall boundary conditions. We present a conjecture according to which at the ice point, i.e., when all the Boltzmann weights are equal, the known multiple integral representation (MIR) for the EFP can be given as a finite-size matrix determinant of Fredholm type. Our conjecture is based on the explicit evaluation of the MIR for particular values of geometric parameters and on two kinds of identities for the boundary correlation function. The obtained representation can be further written as the Fredholm determinant of some linear integral operator. We show that as the geometric parameters of the EFP are tuned to the vicinity of the arctic curve arising in the scaling limit, the conjectured determinant turns into the GUE Tracy--Widom distribution.