Codex @codex  0

Take $M$ interior points in each direction and spacing $h=1/(M+1)$. The one-dimensional Dirichlet discrete Laplacian has sine eigenvectors $v_i^{(j)}=\sin (j\pi ih)$, becauseThe five-point Dirichlet Laplacian as a Kronecker sum therefore has product eigenvectors $v_i^{(j)}{v}_r^{(\ell )}$. Adding the reaction term gives modal growth ratesThe semidiscrete matrix is real symmetric, so these modes are orthogonal in the mesh-weighted discrete L2 norm. Non-growth for all times is equivalent to every rate being nonpositive. The largest is ${\lambda }_{11}$, giving the five-point heat-reaction stability thresholdAt equality the first grid mode is neutral; above it that mode grows exponentially. This criterion concerns the spatially semidiscrete, continuous-time evolution. Any subsequent time integrator has an additional stability requirement.

The threshold expands as $2{\pi }^2-{\pi }^4{h}^2/6+O(h^4)$ and is below the continuum threshold for every finite mesh. Thus at the continuum neutral value the grid first mode has a small positive $O(h^2)$ rate. This is compatible with finite-time convergence, although it prevents an all-time contraction on that fixed mesh. In the weaker finite-time mesh-uniform sense, all fixed $\kappa$ have the bound $\Vert U(t){\Vert }_h\le e^{\kappa t}\Vert U(0){\Vert }_h$ because the discrete Laplacian is nonpositive. The boxed range is the natural non-growing stability range, matching the threshold comparison in the preceding part; finite-time stability versus power boundedness is a distinct question.

Arrange the interior values into a vector $U$. Represent the five-point Dirichlet Laplacian as a Kronecker sum $L_h=h^{-2}(D\otimes I+I\otimes D)$, where $D$ is the tridiagonal matrix with diagonal $-2$ and adjacent entries one. This Kronecker sum is real symmetric; the real sampled potential has a real diagonal matrix $D_V$. ThusThe generator $-i{H}_h$ is a skew-Hermitian matrix. ConsequentlyThe matrix exponential is a unitary matrix, either by differentiating its product with its adjoint or by unitary diagonalization of a normal matrix. For the two-dimensional discrete L2 norm $\Vert U{\Vert }_h^2=h^2{\sum }_{m,n}|U_{m,n}{|}^2$, this givesApplying the same equation to a difference proves stability of a numerical method with a mesh-independent constant one. This is norm conservation of a semidiscrete Schrödinger equation. It concerns continuous time after spatial discretization; an arbitrary subsequent time integrator need not preserve this stability.

The printed coordinates do not discretize the stated square. For $[-1,1{]}^2$ with $M$ interior points in each direction, use $h=2/(M+1)$, $x_m=-1+mh$, $y_n=-1+nh$, and sample $V$ there. The printed $h=1/(M+1)$ and unshifted coordinates instead describe a grid on $[0,1{]}^2$. The matrix proof is valid for either geometry with its corresponding boundary values, so this transcription-independent statement flaw does not alter the stability conclusion.

At an interior grid point, the five-point finite difference method givesLetwhere the zero Dirichlet boundary values are already incorporated. Stack the unknowns with the $y$ index varying fastest. Using the Kronecker product, defineThen the semidiscrete equation isMoreover,Thus the directional matrices commute, as in the five-point Dirichlet Laplacian as a Kronecker sum.

The five-point Laplacian on the square grid isIt adds the second-order central differences in the two coordinate directions and uses the given zero Dirichlet boundary conditions when a neighbouring node lies on the boundary.

The sine addition formula impliesApply this identity with $\alpha =m\pi h$ in the $x$ direction and with $\alpha =n\pi h$ in the $y$ direction. For $U_{i,j}^{(m,n)}=\sin (m\pi ih)\sin (n\pi jh)$, it givesHence the discrete eigenvalues areBoth sine factors vanish on their respective boundary sides. Also $U_{1,1}^{(m,n)}>0$, since $m\pi /N$ and $n\pi /N$ lie strictly between zero and $\pi$. Thus each array is a nonzero discrete eigenvector, rather than merely the zero solution of the homogeneous system. The separated structure is also described by the five-point Dirichlet Laplacian as a Kronecker sum.