Dirac Equation computational complexity

  • Thread starter Astudious
  • Start date
  • #1
61
0

Main Question or Discussion Point

How fast does the computational complexity of the Dirac equation, with regards to full* solution, grow with number of particles N? can we specify the order of time t(N) for this solution in terms of t(N=1)?

(I assume that number of protons, neutrons and electrons combined is N - i.e. that complexity grows similarly with particles of every kind.)

*Full solution meaning a numerical solution reached without making any approximations or neglecting any terms (beyond those neglected even for cases we can solve "exactly" like that of a hydrogen atom)
 

Answers and Replies

  • #2
An exact solution (or, at least, one that can be computed to arbitrarily good precision) means, in effect, simulating the time evolution of the system. With the best known algorithms for classical simulations of quantum systems, simulation complexity of arbitrary composite systems grows exponentially in the number of subsystems. It's not possible to do better than that in full generality. That's pretty much the main impetus behind quantum computing. A quantum computer can simulate other quantum systems (including, it seems, particles in quantum field theory—this a current research project by Preskill and others) in polynomial time.

Certain special cases can have efficient classical solutions, of course. In this case, it depends on what you mean by solutions to the Dirac equation. If you mean the equation for free fields, then, well, the solution is basically constant time. Once you've solved it, you can apply it to as many fermions as you like. But if you mean an interacting Dirac equation—mediated by, say, a Yukawa interaction or the EM field—then you're most likely looking at exponential classical complexity. If you can do better, you can probably become quite famous...
 
Last edited:
  • #3
61
0
An exact solution (or, at least, one that can be computed to arbitrarily good precision) means, in effect, simulating the time evolution of the system. With the best known algorithms for classical simulations of quantum systems, simulation complexity of arbitrary composite systems grows exponentially in the number of subsystems. It's not possible to do better than that in full generality.
Ah, so t(N) = t(N=1)*e^(kN) in full generality for some constant k, is the best we can guarantee?

Let us take the case of atomic/molecular systems, can we predict k for that case?

But if you mean an interacting Dirac equation—mediated by, say, a Yukawa interaction or the EM field—then you're most likely looking at exponential classical complexity. If you can do better, you can probably become quite famous...
Yes. We could consider simple chemical systems for now.
 

Related Threads on Dirac Equation computational complexity

Replies
1
Views
1K
  • Last Post
Replies
1
Views
3K
  • Last Post
2
Replies
37
Views
6K
  • Last Post
Replies
2
Views
945
  • Last Post
Replies
11
Views
678
  • Last Post
Replies
5
Views
752
  • Last Post
Replies
4
Views
569
  • Last Post
Replies
1
Views
577
  • Last Post
Replies
16
Views
4K
  • Last Post
Replies
3
Views
2K
Top