Matrix model with 3D BMS constraints

Our goal is to find a matrix model with BMS₃ constraints built in. These constraints are imposed through loop equations. We solve them using a free field realization of the algebra and write down the partition function in eigenvalue form. We comment on the nature of this partition function and its relations with integrable systems.