  We compute the suppression of the bottomonia states Upsilon(1S), Upsilon(2S),
Upsilon(3S), chi_b(1P), chi_b(2P), and chi_b(3P) states in Large Hadron
Collider (LHC) sqrt(s_NN)) = 5.023 TeV Pb-Pb collisions. For the background
evolution we use 3+1d anisotropic hydrodynamics with conditions extrapolated
from sqrt(s_NN) = 2.76 TeV and we self-consistently compute bottomonia decay
rates including non-equilibrium corrections to the interaction potential. For
our final results, we make predictions for R_AA as function of centrality,
rapidity, and p_T for the Upsilon(1S) and Upsilon(2S) states, including feed
down effects. In order to assess the dependence on some of the model
assumptions, we vary the shear viscosity-to-entropy density ratio, 4 pi eta/s =
(1, 2, 3), and the initial momentum-space anisotropy parameter, xi_0 = (0, 10,
50), while holding the total light hadron multiplicity fixed.
