The production of $\Upsilon$ mesons in Pb-Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}}$ = 5 TeV is measured with the muon spectrometer of the ALICE detector at the LHC. The yields as well as the nuclear modification factors are determined in the forward rapidity region $2.5<y<4.0$, as a function of rapidity, transverse momentum and collision centrality. The results show that the production of the $\Upsilon$(1S) meson is suppressed by a factor of about three with respect to the production in proton-proton collisions. For the first time, a significant signal for the $\Upsilon$(2S) meson is observed at forward rapidity, indicating a suppression stronger by about a factor 2-3 with respect to the ground state. The measurements are compared with transport, hydrodynamic, comover and statistical hadronisation model calculations.
Rapidity-differential yield of $\Upsilon(1\mathrm{S}) \rightarrow \mu^{+}\mu^{-}$ divided by the average nuclear overlap function $\langle T_{\mathrm{AA}} \rangle$ for the 0–90% centrality interval ($\langle T_{\mathrm{AA}} \rangle$ = 6.28 $\pm$ 0.06 mb$^{-1}$).
Rapidity-differential yield of $\Upsilon(2\mathrm{S}) \rightarrow \mu^{+}\mu^{-}$ divided by the average nuclear overlap function $\langle T_{\mathrm{AA}} \rangle$ for the 0–90% centrality interval ($\langle T_{\mathrm{AA}} \rangle$ = 6.28 $\pm$ 0.06 mb$^{-1}$).
$p_{\mathrm{T}}$-differential yield of $\Upsilon(1\mathrm{S}) \rightarrow \mu^{+}\mu^{-}$ divided by the average nuclear overlap function $\langle T_{\mathrm{AA}} \rangle$ for the 0–90% centrality interval ($\langle T_{\mathrm{AA}} \rangle$ = 6.28 $\pm$ 0.06 mb$^{-1}$).