We give a report of medium modification of Z+jet correlations in Pb+Pb collisions at the Large Hadron Collider using Sherpa to generate initial Z+jet at next-to-leading-order matrix element matched to the parton shower, and the Linear Boltzmann Transport Model for jet propagation in the expanding quark-gluon-plasma. Our numerical calculations show excellent agreement with all available observables of Z+jet simultaneously in both proton + proton and Pb+Pb collisions. Our results can well explain the shift of momentum asymmetry $x_{jZ}=p_T^{jet}/p_T^Z$ as well as its mean values, the suppression of the jet yields per Z trigger $R_{jZ}$ and the modification of azimuthal angle correlations $\Delta \phi_{jZ}$. We also demonstrate that it is the energy loss effect on multi-jet from higher order corrections that leads to the suppression of the Z+jet correlations in small azimuthal angle difference $\Delta \phi_{jZ}$ region and in small $x_{jZ}$ region. The jet shape reflecting transverse momentum distribution inside the jet is also calculated, which indicates that large fraction of jet energy is carried away from the jet axis in Pb+Pb collisions.