We derive the quantitative estimates of propagation of chaos for the large interacting particle systems in terms of the relative entropy between the joint law of the particles and the tensorized law of the mean field PDE. We resolve this problem for the first time for the viscous vortex model that approximating 2D Navier-Stokes equation in the vorticity formulation on the whole space. We obtain as key tools the Li-Yau-type estimates and Hamilton-type heat kernel estimates for 2D Navier-Stokes on the whole space. This is based on a joint work with Xuanrui Feng from Peking University.