Here's an attempt at answering your Question A:
Currently, one of the most powerful methods for proving upper bounds on sphere packing densities is the linear programming bound of Cohn and Elkies (which is what Viazovska used).
Exact answers
According to the numerical computations of Cohn and Elkies, we do not know any dimensions other than 1, 2, 8, and 24 where their linear program has a chance of proving the exact bound (see their Figure 1). There are exceptionally good lattices in those four dimensions that match the bounds given by the Cohn-Elkies linear program. It could be the case (likely?) that these are the only dimensions where the bounds match, although we don't know how to rule out other possibilities. There could be some deep reason (perhaps related to finite simple groups) why these four dimensions are special.
Hales' proof of Kepler's conjecture (sphere packing in dimension 3) uses a completely different method.
Asymptotic bounds
An important open question, for which there is still an exponential amount of room for improvement, is the problem of sphere packing bounds in very high dimensions.
The current best asymptotic upper bound to the packing density $\Delta_n$ in $\mathbb{R}^n$ is by Kabatiansky and Levenshtein (1977)
$$
\Delta_n \le 2^{-(0.5990\dots + o(1))n}
$$
(see, e.g., my paper with Henry Cohn on the matter; in particular, this asymptotic bound lies within the scope of the Cohn-Elkies linear program). The current best lower bound is due to Venkatesh. For all $n$,
$$
\Delta_n \ge c n 2^{-n}
$$
for some explicit $c$ (that has been improved over time), and for a specific infinite set of $n$,
$$
\Delta_n \ge c n (\log\log n) 2^{-n}.
$$
Closing the gap seems to be a difficult problem.
Further reading