seemps.evolution.euler2#

seemps.evolution.euler2(H: MPO, time: float | tuple[float | int | floating | integer, float | int | floating | integer] | Sequence[float | int | floating | integer] | ndarray[tuple[Any, ...], dtype[_ScalarT]], state: MPS, steps: int = 1000, strategy: Strategy = DEFAULT_STRATEGY, callback: Callable[[float, MPS], Any] | None = None, itime: bool = False)[source]#

Solve a Schrodinger equation using the 2nd order Euler method.

Implements a two-step integration method. In imaginary time this is

\[\xi(t_{n}) = 2 \psi(t_{n}) - i \delta t H \psi(t_{n}) \psi(t_{n+1}) = \psi(t_{n}) - (i \delta t/2) H \xi(t_{n})\]

The Euler algorithm is a very bad integrator and is offered only for illustrative purposes. See euler() to understand the parameters and the function’s output.