Correction: The article states that Thomson's problem has only been solved for 2, 3, 4, 6, and 12 charges. However, in 2013 the 5-electron case was solved by Richard Schwartz. Here's the paper:
Needless to say, there has been a lot of work on this subject from a number of different directions, and any particular new paper is unlikely a priori to be a revolutionary breakthrough. Look at Google Scholar or a similar citation graph and you can find thousands of papers about it. (Some key words: “spherical cap discrepancy”, “Fekete points”, “Riesz energy”, “logarithmic potential”, obviously also “Thomson’s problem”, etc.)
I actually was asked this problem in a programming interview. This is first I've read that it is an unsolved problem, which definitely makes me chuckle.
I guess if you give someone a puzzle, there's a possibility they might solve the puzzle from memory. If you give an unsolved problem, you know they won't have seen the solution before and can see their thought process better(since you know thinking will be required).
The exact phrasing matters too, I guess. There are lots of ways to find a good packing, and even good bounds on how far your good packing is from the global optimizer. So a reasonable question to your interviewer would have been "does it need to be the exact global optimizer, or can it be a good approximation ?". If they wanted a good approximation it might not be a bad question, and the question may have been designed to get you to ask for the distinction? Or maybe they were just dicks, who knows.
Not sure. For n points on a sphere, you're trying to prove a global minimum a very bumpy function over 2n-dimensional space. It may become computationally intractable.
And even if you can find a numerical solution to arbitrary precision, at best you only get arbitrarily close to the mathematically optimal solution (at best because a purely numerical approach cannot rule out that it misses a very small, but high peak in the function to be optimized)
As an extreme example, for n=2 and using binary floating point, you probably will find the mathematically correct solution, but there's no way to tell numerically whether your answer is off by a fraction of your floating point's epsilon value.
That's probably not important to physicists who want to know an answer, but it is for mathematicians.
"Numerical solutions for N=8 and 20 are not the regular convex polyhedral configurations of the remaining two Platonic solids, whose faces are square and pentagonal, respectively."
I would like to see the visualizations of the better solutions for n=20 (wikipedia links to the one for n=8) and/or hear a heuristic argument as to how that can happen.
The function here has a simple continuity, so there is no risk of narrow peaks beyond some computable threshold. The challenge is the high dimensionality.
The unsolved problem referred to in the title is finding a configuration close to the global minimum. Most optimization algorithms need such a configuration as a starting point.
http://www.tandfonline.com/doi/abs/10.1080/10586458.2013.766...