A dynamical concept defining how quantum states or operators change with time under a specified Hamiltonian. It governs time propagation, phase accumulation, and the effect of driving or slowly varying parameters when present. It does not ensure accurate prediction without correct initial conditions, boundary conditions, and validated model assumptions. It provides the basis for computing transition probabilities, energy spectra, and time-dependent expectation values. The concept is generally stable, though approximation techniques and simulation tools evolve over time.

Principle

Discretize space (and possibly time) in a way that preserves key physical properties: unitarity or probability conservation, numerical stability, and sufficient resolution of relevant length and energy scales. Choose algorithms with controlled error, known stability regions, and, when possible, conservation of invariants.

Demonstration

Simulate a Gaussian wavepacket scattering off a finite Gaussian barrier using the split-operator Fourier method: apply kinetic evolution in momentum space and potential evolution in position space alternately, ensuring time-step small enough for accuracy while implementing absorbing boundaries to avoid reflections from simulation edges.

Misapplication

Using an explicit Euler time integrator or too-large time steps that break unitarity and cause norm growth/decay, employing too coarse spatial grids that alias high-frequency components, or failing to implement absorbing boundary conditions, producing spurious reflections that contaminate results.

Consequence

A properly implemented numerical simulation yields time-resolved probabilities, transmission/reflection coefficients, bound-state energies, and expectation values with quantifiable errors and enables exploration of parameter regimes inaccessible to analytic methods.

Reversal

The reversal is reliance solely on analytic solutions when unavailable: refusing numerical simulation constrains exploration to solvable idealized models and prevents quantitative predictions for complex or experimentally realistic setups.

Boundary

Applies to discretized approximations of Schrödinger problems; the fidelity to the continuous problem depends on grid resolution, time-step, boundary treatments and operator representation. It does not equate to exact solutions and may require extrapolation or convergence testing to validate results.

Semantic Tension

There is tension between high-order spectral methods (very accurate for smooth problems) and local finite-difference/finite-element methods (more flexible for complex geometries); trade-offs include accuracy, ease of implementing boundary conditions, and computational cost.

Synthesis

Numerical simulation of the Schrödinger equation is the disciplined choice and implementation of discretization and propagation algorithms that balance unitarity, stability, accuracy, and computational feasibility to produce reliable approximations of quantum dynamics and spectra, accompanied by convergence checks and uncertainty quantification.