# Effect on methods like Crank-Nicolson of adding a potential term, changing heat equation to Schrodinger equation

I'm studying up on methods for numerically solving the Schrodinger equation. The Schrodinger equation with a zero potential is formally identical to the heat equation in the sense that we just make one of the coefficients in the equation imaginary.

A practitioner of finite-difference methods for the heat equation has to be aware of certain things. From a casual reading of some Wikipedia articles, the main things I've learned are that:

1. the forward method can be unstable;

2. Crank-Nicolson has better asymptotics; and

3. I need to satisfy the Courant condition in the Crank-Nicolson method.

Do any of these facts change if I switch to the Schrodinger equation with a nonzero potential?

In the Schrodinger equation, a kink in the wavefunction contains arbitrarily short wavelengths and therefore has components that propagate with arbitrarily high velocities. Does this mean that the Courant condition can never be satisfied if there is a kink? In a discretized representation, however, there is no clear distinction between a kink and a differentiable point, so is there some adaptive criterion that can be used to set an appropriate time step?

• I think we'll need to review some of those "facts" first. What do you mean by "the forward method"? What do you mean by "Crank-Nicolson has better asymptotics"? Where did you find on Wikipedia that Crank-Nicolson has the Courant condition? We should edit that article. Why do you think a wavefunction can have a kink? That would mean infinite kinetic energy, unless the wavefunction vanishes at that point. If there are no delta-functions in the potential, the wavefunction should be smooth. Jul 3, 2020 at 4:32
• What do you mean by "the forward method"? en.wikipedia.org/wiki/Finite_difference_method#Explicit_method What do you mean by "Crank-Nicolson has better asymptotics"? Its errors are proportional to the square of the time step, whereas the forward method's errors are proportional to the time step. Where did you find on Wikipedia that Crank-Nicolson has the Courant condition? Same article: "The von Neumann and Courant-Friedrichs-Lewy criteria are often evaluated to determine the numerical model stability." If I'm interpreting this wrong, it would be good to know that.
– user28109
Jul 3, 2020 at 17:24
• I think that you could add some clarification to your question. Comments are not intended to discuss but to clarify questions. For example, you could say "explicit Euler" or "forward Euler" instead of "forward method". Jul 3, 2020 at 18:38