# The Many Body problem

(This is a simple question, with likely a rather involved answer.)

What are the primary obstacles to solve the many-body problem in quantum mechanics?

Specifically, if we have a Hamiltonian for a number of interdependent particles, why is solving for the time-independent wavefunction so hard? Is the problem essentially just mathematical, or are there physical issues too? The many-body problem of Newtonian mechanics (for example gravitational bodies) seems to be very difficult, with no solution for $n > 3$. Is the quantum mechanical case easier or more difficult, or both in some respects?

In relation to this, what sort of approximations/approaches are typically used to solve a system composed of many bodies in arbitrary states? (We do of course have perturbation theory which is sometimes useful, though not in the case of high coupling/interaction. Density functional theory, for example, applies well to solids, but what about arbitrary systems?)

Finally, is it theoretically and/or practically impossible to simulate high-order phenomena such as chemical reactions and biological functions precisely using Schrodinger's quantum mechanics, over even QFT (quantum field theory)?

(Note: this question is largely intended for seeding, though I'm curious about answers beyond what I already know too!)

-
Why do you restrict it to quantum problems ? –  Cedric H. Nov 4 '10 at 23:19
You could say restrict, but in many ways it's generalising! In any case, the problem is rather different for quantum mechanics, and certainly more interesting I find. –  Noldorin Nov 4 '10 at 23:30

First let me start by saying that the $N$-body problem in classical mechanics is not computationally difficult to approximate a solution to. It is simply that in general there is not a closed form analytic solution, which is why we must rely on numerics.

For quantum mechanics, however, the problem is much harder. This is because in quantum mechanics, the state space required to represent the system must be able to represent all possible superpositions of particles. While the number of orthogonal states is exponential in the size of the system, each has an associated phase and amplitude, which even with the most coarse grain discretization will lead to a double exponential in the number of possible states required to represent it. Thus in quantum systems you need $O(2^{2^n})$ variables to reasonable approximate any possible state of the system, versus only $O(2^n)$ required to represent an analogous classical system. Since we can represent $2^m$ states with $m$ bits, to represent the classical state space we need only $O(n)$ bits, versus $O(2^n)$ bits required to directly represent the quantum system. This is why it is believed to be impossible to simulate a quantum computer in polynomial time, but Newtonian physics can be simulated in polynomial time.

Calculating ground states is even harder than simulating the systems. Indeed, in general finding the ground state of a classical Hamiltonian is NP-complete, while finding the ground state of a quantum Hamiltonian is QMA-complete.

-
An informative answer. Interestingly, I've heard that (universal) qunatum computers should be able to simulate other quantum computers in polynomial times, just that classical computers can't. –  Noldorin Nov 5 '10 at 16:14
Yes, that is true. Calculating ground states seems to be beyond their reach though. –  Joe Fitzsimons Nov 5 '10 at 16:39
Ah I see. I suppose calculating ground states and performing complete quantum simulations of systems typically cover different ranges of application, however, so it's not all bad news. Anyway, cheers for the detail, you seem to be very knowledgeable on the subject; the answer is yours. –  Noldorin Nov 9 '10 at 3:22
@Noldorin: Thanks. I only know this stuff because I have spent quite a while working in this exact field. By the way,ground states are to some extent less relevant because the systems for which is is computationally hard to calculate the ground state of (at least on a QC) don't cool efficiently either. –  Joe Fitzsimons Nov 9 '10 at 3:50
the N-body problem in classical mechanics is not computationally difficult to approximate a solution to. This is somewhat misleading. The systems are typically chaotic, so they're fundamentally impossible to predict, for the same reason that long-term weather forecasts are fundamentally impossible. Calculating ground states is even harder than simulating the systems. This may be true in some abstract sense, but it is not true in a practical sense. For example, it's relatively straightforward to get a good approximation to the ground state of a nucleus using the Hartree-Fock method. –  Ben Crowell Aug 16 '13 at 15:06

The answer is fairly simple -- classical N-body problem has its solution in $6N$ 1D functions of time, quantum N-body problem has its solution in one complex function, but $3N$-dimensional (not counting spin and similar stuff). Then, there is no wonder why one can find analytical solutions only for trivial problems or at least make $N$ huge and escape into statistical mechanics. And yes, this is only the problem of mathematical complexity here.
From modelling point of view exact solving also seems hopeless, with only memory complexity of $\mathcal{O}(K^{3N})$.

For the rest of the answer I will restrict myself to quantum chemistry/material science, since this is the most exploited region -- this means we are now talking about atoms. First of all, atoms have small and very heavy nuclei, which thus can be treated as almost stationary sources of electrostatic potential; this reduces the problem to electrons only (Born-Oppenheimer approx.). Now, there are two main routes to follow: Hartree-Fock or Density Functional Theory.
In HF, one roughly represents the many-body weavefunction as a combination of some standard base functions -- then one can optimize their contributions to get minimal energy, yet using extended Hamiltonian to adjust the effects of such approximation. In DFT, one encouraged by Hohenberg-Kohn theorems reduces the many body weavefunction to electron probability density field (3-dimensional), and accordingly Shroedinger equation terms into density functionals (and there approximations are applied). Next, it can be either solved as this 3D field or in Kohn-Sham way, which is pretty much Hartree-Fock for DFT (one represents density with base functions). People sometimes are making something analytical here, but those are mostly theories made to support computational approaches.

And finally your last question: those approximate methods (but still ab initio -- there are no experimental parameters there) do predict things like chemical reactions, various spectra and other measurable quantities; accuracy is problematic though. Biology is mostly out of reach because of a time scale; at least there are hybrid methods able to mix for instance the classical simulation of protein motion with quantum simulation of the binding site when it is squeezed enough so something quantum like enzymatic reaction can take place.

-
Looks like a pretty good answer, I'll read it properly tomorrow. In any case, it's important to make clearer that that although the "properties of the solution* are "fairly simple", the solutions themselves are certainly not! –  Noldorin Nov 5 '10 at 1:37
Note that H-F, DFT are the main approximation techniques to the quantum many-body problem, though neither are "well-controlled approximations" in the sense that they are used as the first term in a convergent expansion to the actual solution. And I'm not sure what level of computational complexity they reduce the problem to, though that's an important question. –  j.c. Nov 5 '10 at 14:42
@j.c. Those are approximated theories rather than approximated ways of solving equations. Reduction of complexity is obvious -- 3N-dim function to a vector of parameters in case of HF or to 3-dim field in case of DFT. –  mbq Nov 5 '10 at 17:23

In addition to what mbq said, it might be interesting to know that things get really funny in relativistic quantum mechanics, that is using the Klein-Gordon and the Dirac equation (but without the "second" quantization of Quantum Field Theory). There, there's one wave function per particle sort, so no matter how many particles of one kind you consider, the only thing that changes is the field itself. You only get more degrees of freedom by actually adding another kind of particles. Of course, since Fermions require Spinors, you may end up with other computational issues then...

-
The problem here, of course, is that the field modes are continuous variables. –  Joe Fitzsimons Nov 5 '10 at 6:47
By which I mean the problem with simulating the system, not a problem with your answer. –  Joe Fitzsimons Nov 5 '10 at 7:04
Yeah, I was curious as to whether QFT actually makes things easier in some respect. It's a tricky scenario. –  Noldorin Nov 5 '10 at 16:15
It definitely can only make things harder as you can encode a discrete system in the CV but not necessarily the other way around. –  Joe Fitzsimons Nov 5 '10 at 16:38
Noldorin: QFT would probably make things even more complicated, I was only wondering whether the unquantized relativistic QM equations would yield an advantage over the non-QFT Schrödinger equation, but as @Joe mentions, this may not be the case... –  Tobias Kienzler Nov 8 '10 at 8:00

The many-body equation is immensely difficult to study, both classically and quantum-mechanically. The late John Pople, of Northwestern University, won a Nobel Prize in 1998 for his numerical models of wave functions of atoms, developing a theoretical basis for their chemical properties. Here is a link:

http://nobelprize.org/nobel_prizes/chemistry/laureates/1998/

-
Thanks for the info. I may just have to read some of Pople's papers some day, out of curiosity. :) –  Noldorin Nov 7 '10 at 1:14

On a more abstract level, the problem is linearity versus non-linearity. It's straightforward to solve a number of linear equations, and they always yield an analytic answer. However, non-linear equations produce chaotic behaviour, which cannot be generalised in most cases.

As an example, the 3-body Newtonian problem involves 2C3 = 3 non-linear equations; the nonlinearity comes from the r2 relationship. And 3 non-linear relations are the minimum requirement for a chaotic system.

Similarly, quantum mechanics involves a large number of non-linear equations - given a set of 3 electrons, each will repel the others via a non-linear relation, and with even more complexity than the Newtonian problem where all things are known and determinable.

So, the simple answer is that the problem is mathematics that can't be solved for the general case, which result from the physics, and that the quantum case is indeed worse than the classical one.

-