In a [recent publication](http://arxiv.org/abs/1603.04246) by the Ukrainian mathematician Maryna Viazovska the Kepler problem for dimension $8$ and $24$, namely the densest packing of spheres, was solved.


Admittedly it is very difficult to grasp the work for someone not involved in the field, but it would be incredibly valuable if, briefly, the methodology adopted, and the main breakthrough could be explained at a conceptual level. Such that one can have a rough idea how the solution came to be.


*This is merely an attempt to get more input in order to understand some of the main ideas involved, as this is very exciting news given that the problem of densest packing of spheres in different dimensions has a long history.*



Understanding sphere packing in higher dimensions
There are two things you need to understand. The first is how to prove sphere packing bounds via harmonic analysis ("linear programming bounds"). My [lecture notes from PCMI 2014](http://arxiv.org/abs/1603.05202) give an exposition of this theory, which covers the period up to, but not including, Viazovska's paper on eight dimensions. This reduces the problem to finding radial functions in eight and twenty-four dimensions with certain roots and sign conditions, and the lecture notes discuss why everybody believed such functions exist. The functions themselves are surprisingly subtle; see [this paper](http://arxiv.org/abs/1603.04759) with Steve Miller for numerical experiments and conjectures.


Then Viazovska's breakthrough lies in how to construct the right functions. It's convenient to split them up into eigenfunctions of the Fourier transform, which have to have roots at certain locations. One way to look at it is that she brute forces those roots by including a sine squared factor that produces the roots, and then she multiplies it by something that's essentially the Laplace transform of a (quasi)modular form. To make this work, one has to answer three questions:


1. What sort of modular forms play nice with the sine squared factor and produce radial eigenfunctions of the Fourier transform?
2. Do there exist such modular forms that actually give the desired functions?
3. How can we prove the inequalities these functions are supposed to satisfy?


Viazovska's methodology gives nice answers:


1. This amounts to identifying the right level, weight, and depth (for quasimodular forms).
2. Yes. One can reverse engineer them based on the desired properties.
3. In the best way one could hope for, namely at the level of the modular form itself (rather than having to worry about something subtle happening in the Laplace transform).


For details, see the [eight](http://arxiv.org/abs/1603.04246) and [twenty-four](http://arxiv.org/abs/1603.06518) dimensional papers. As Peter Sarnak points out in Erica Klarreich's [Quanta article](https://www.quantamagazine.org/20160330-sphere-packing-solved-in-higher-dimensions/), the actual construction is surprisingly simple. If you know about linear programming bounds and basic facts about modular forms, the proof itself is very down to earth.


