Solving a non-linear Schrödinger equation

Question

I am trying to use NDSolve to solve numerically the well-known Schrödinger-Newton equation for a certain initial value of the wave function.

ClearAll[r, t, ψ, ϕ]

σ = 1

sol = NDSolve[{I*D[ψ[r, t], t] == -1/(2)*D[ψ[r, t], {r, 2}] + ψ[r, t]*ϕ[r, t], D[ϕ[r, t], {r, 2}] == σ*ψ[x,t]*ψ[x, t]\[Conjugate], ψ[r, 0] == (2/Pi)^(3/4)*Exp[-r^2], ψ[-100,t] == (2/Pi)^(3/4)*Exp[-100^2], ψ[100, t] == (2/Pi)^(3/4)*Exp[-100^2], ϕ[-100, t] == 0, ϕ[100, t] == 0, ϕ[r, 0] == 0}, ψ, {r, -100, 100}, {t, 0, 1}, Method -> Automatic]

Mathematica says 'NDSolve::delpde: Delay partial differential equations are not currently supported by NDSolve'. Do any of you know how to re-express the problem in order to get the solution?

Answer 1

There is a typo:

D[ϕ[r, t], {r, 2}] == σ*ψ[x,t]*ψ[x, t]\[Conjugate]

x should be r in the above snippet from your code. It is interpreting the equation as a delayed equation because it thinks you are meaning to call ψ at the (undetermined) value x. Although, my machine still is not powerful enough to produce the solution. Maybe yours is.