Unlike in Maxwell’s wave Equation, where perturbations propagate at the speed of light, in Schrödinger dynamics, any local perturbation is instantaneously felt all over space, just as in the diffusion Equation . However, unlike the diffusion Equation where the nonlocal effect is exponentially small, in the Schrödinger Equation, it decays much slower—as a power law.
In both cases, i.e., in the diffusion and the Schrödinger cases, the causality is violated due to the asymmetry between space and time.
In the Klein-Gordon’s (KG), or similarly in the Dirac’s, Equation, due to the symmetry between space and time, causality reappears. In the KG case, the nonlinear dispersion relation distorts the wavefunction in high agreement with the Schrödinger Equation only as far as causality allows, i.e., as far as the distance from the local perturbation  . That is, any local perturbation affects over the entire domain. Clearly, in the non-relativistic regime (i.e., the Schrödinger case) this domain is the entire space. As a result, an initial discontinuous wavefunction can kindle currents all over space instantaneously  .
There are several numerical methods to calculate non-relativistic quantum mechanics dynamics. Every method emphasizes a different feature of this dynamics. The Schrödinger Equation is a differential Equation and therefore it emphasizes the locality nature of the dynamics (see, for example, Ref.  ). Consequently, it is not optimized to evaluate the dynamics of singular, discontinuous or even sharp functions. Eigenstates expansion, in general, and Fourier expansion, in particular, is a powerful technique when the eigenstates, as well as their boundary conditions, are well-known (see, for example, Ref.  ). It is inconvenient to use this method when the boundary conditions are ill-defined, such as in scattering problems. Path integrals method is an elegant method to formulate the wavefunction as a superposition of infinite paths . Due to its complexity, it is mainly used for pedagogical purposes to explain the nature of particle interferences (or with conjunction with Monte-Carlo simulation  ). When other options are possible, such as kernel methods, they are usually preferred. Kernel, or Green function, methods are very useful in cases where the initial wavefunction is given  . However, if the wavefunction is spread all over space, the kernel’s oscillations frequencies, which increase indefinitely may cause numerical difficulties.
In this paper we discuss the case, where the initial wavefunction is spectrally bounded, i.e., the results of any energy measurement are bounded by a certain value. In this case, all the information of the wavefunction can be encapsulated in a series of discrete values. Consequently, the quantum dynamics can be written as an exact infinite set of ordinary differential Equations. It will be shown that the couplings between these Equations consist of universal numbers, which inversely depend on the spatial distance between points on the wavefunction.
2. The Dynamics
The differential version of the free Schrödinger Equation
sometimes conceals the nonlocal properties of the Schrödinger dynamics. However, its integral presentation
with the free-space Schrödinger Kernel 
illustrates the nonlocality more vividly.
Nevertheless, the oscillations’ frequency increase so rapidly that their averages (which is equivalent to the integral operation) quickly converges to zero and the locality properties reappear. In particular, in some cases, the latter nonlocal Equation (2) was used to derive the former local one (1) .
A local analysis of the Schrödinger Equation is an excellent approximation in the quasi-classical regimes, which mathematically equivalent to the stationary phase approximation . However, the locality is questionable in the quantum regime. The problem is that Equation (3) is the impulse response of the Schrödinger Equation (1), i.e., it is the quantum system’s response to the initial state of a delta function. However, a delta function can never be a physical state (it is based on infinite energies and it is not normalizable). To take a more physical initial state, it is usually accustomed to replacing the impulse response with a more physical, finite-width pulse-response, i.e.,
In which case the pulse response (after a period t)
decays (in space) exponentially as well.
Therefore, it may seem as if the locality approximation is justified since in
principle one can choose to be arbitrarily small. However, this result is
based on the premises that the spatial spectrum of the wavefunction is unbounded.
When the spatial spectrum of the wavefunction is bounded, i.e., when the spectral coefficient beyond the spatial frequency domain are all zero (k represents the wavenumber), then, according to the Nyquist theorem, the wavefunction can be written as a superposition of sinc functions . That is, all the information in the wavefunction can be written as an infinite discrete series
of complex numbers for . In the spectral domain, the wavefunction occupies the spatial spectral bandwidth . Consequently, from the Nyquist theorem, the initial wavefunction can be
written as an infinite sequence of overlapping Nyquist-sinc functions (for applications in the optical communication sphere see Ref.  -  ) (see Figure 1), i.e.,
where is the “sinc” function, which corresponds to the Fourier
transform of the rectangular function (the spectral shape of the channel).
Figure 1. Illustration of the method, in which any spectrally bounded function can be written as an infinite series of sinc pulses. In the figure the sinc pulses, which are centered at , 4 and 5 are plotted by dashed curves, while the final function is presented by solid curves (real/imaginary part in the upper/lower panel). The dashed curves represent single sinc’s, while the solid curve represents the superposition function. The Upper panel represents the real part of the wavefunction, while the lower one represents its imaginary component.
Due to the linear nature of the system, Equation (6) can be solved directly
where “dsinc” is the dynamic-sync function 
To simplify the derivation we use the dimensionless variables
and . (9)
With notations (9) at hand, Equations (1) and (7) can be rewritten
3. Matrix Formulation
When is the spatial resolution of the problem, then the wavefunction at the
Figure 2. The relation between the real and imaginary parts of the dsinc function for the discrete values .
Figure 3. The temporal dependence of the absolute value of dsinc for the discrete values .
center of the mth point, i.e., at , is a simple discrete convolution
and . (13)
and , (14)
then Equation (10) can be written as a linear set of differential Equations
with the universal (and dimensionless) vector
and the asterisk stands for discrete convolution.
Note that due to Euler’s formula  .
This Equation is universal in the sense that the vector is time independent. This is a unique property of the sinc pulses, which does not exist in other sets of orthogonal pulses (like rectangular pulses).
Moreover, since the Schrödinger dynamics is a unitary operation, normalisation is kept and there is no change in the wavefunction spectrum. Therefore, Equation (15) is valid for any given time.
In a matrix form, Equation (15) can be written
where M is a matrix with the coefficients
The nonlocality of this form is clearly emphasized, when compared to the ordinary numerical form of the Schrödinger Equation with the ordinary 1D Cartesian local Laplacian
In the presence of a non-zero potential, whose maximum spatial frequencies is also lower than the Schrödinger Equation can be rewritten in a matrix form
Therefore, a pulse which is initially located at has an instantaneous effect over the entire space, and its effect on any other point (say ) is inversely proportional to the distance between them, i.e., .
On the other hand, when the local Laplacian is used, i.e., (20), then the pulse effect on a point afar, would be felt only after n consecutive steps. Therefore, if there is a barrier (i.e., non zero potential) between these points ( and ) then with a local Laplacian it may seem that the effect of the one on the other (and vice versa) must take into account the barrier in between the two points. However, in fact, as the nonlocal form teaches, in the short time its effect is negligible since the Schrödinger Equation can be approximated by
and its short-time solution
shows that in the short time the potential has only a local effect (provided it is a smooth function), i.e., is affected only by (and the effect is a simple phase change). However, the wavefunction has a nonlocal effect, i.e., is affected by any non zero (for any m). This result is consistent with Ref. , where it was demonstrated that in short time, singular wavefunction are unaffected by the barrier despite their nonlocal effect.
4. Inverse Square Law
By multiplying Equation (15) with the complex conjugate of the wavefunction and taking the real part of the Equation one finds a nonlocal Equation for the probability density
Using the notation then Equation (26) is simply
and the equivalent phase Equation reads
If is presented as a 2D vector in a 3D space
instead of a complex number in a complex plane, then the numerator in the summation can be presented as the cross product of two vectors, i.e.,
In this terminology, is the norm of the vector and the cross represents cross (3D vectors) product.
It is instructive to see the resemblance between this law and any other inverse square law.
Equation (30) can also be written in terms of the derivative of the vector’s norm
is the unit vector in the direction.
It is therefore clear, that maximum probability density transfer occurs when the relative phase between the two points is , i.e. when the “vectors” and are orthogonal.
In Figure 4 such a density transfer is illustrated. In this case the wavefunction
was taken as the initial state (N is the normalisation constant). As can be seen, in the short time regime ( in this case), the probability was transferred from the pulse at to the one at .
Moreover, it is clear from Equation (27) and Equation (30) that maximum probability transfer to a certain location (say ) occurs (provided the initial state is bounded) when the initial state oscillates in signs, i.e.,
In this case the rate in which the probability increases (or decreases in the opposite case) is exactly since (using Euler formula, see Ref.   )
and the probability density at can increase almost four fold before it start to decay (see Figure 5).
Figure 4. Illustration of probability transfer. The dashed curve represents the initial state , while the solid curve stands for . In both cases the real part is plotted in the upper panel, while the imaginary part is plotted in the lower one.
Figure 5. The short time dynamics of the wavefunction (32). The dashed curves represent the initial ( ) state, while the solid curves represent the state after a period of , where the probability density at increases by a factor of 4.
The source of this nonlocality is the fact that each one of the sinc is spread over the entire space. However, the important result is, that this nonlocal presentation of the Schrödinger Equation is universal and it is independent of , which can be as short as the spatial measurement’s accuracy.
It has been shown that when a given wavefunction is spectrally bounded, then the Schrödinger dynamics can be formulated in a universal nonlocal form. Instead of a local partial differential Equation, it can be formulated as an infinite set of ordinary differential Equation, where the coupling is pure numbers, which
are strongly related to Euler’s formula .
Therefore, the mutual effect of every two points on the wavefunction is instantaneous and can be formulated by an inverse square law.
 Garca-Calderon, G., Rubio, A. and Villavicencio, J. (1999) Low-Energy Relativistic Effects and Nonlocality in Time-Dependent Tunneling. Physical Review A, 59, 1758.
 Delgado, F., Muga, J.G., Ruschhaupt, A., Garcia-Calderon, G. and Villavicencio, J. (2003) Tunneling Dynamics in Relativistic and Nonrelativistic Wave Equations. Physical Review A, 68, Article ID: 032101.
 Soto, M.A., Alem, M., Shoaie, M.A., Vedadi, A., Brès, C.-S., Thévenaz, L. and Schneider, T. (2013) Optical Sinc-Shaped Nyquist Pulses of Exceptional Quality. Nature Communications, 4, 2898.
 Schmogrow, R., Bouziane, R., Meyer, M., Milder, P.A., Schindler, P.C., Killey, R.I., Bayvel, P., Koos, C., Freude, W. and Leuthold, J. (2012) Real-Time OFDM or Nyquist Pulse Generation—Which Performs Better with Limited Resources? Optics Express, 20, B543.
 Hirooka, T., Ruan, P., Guan, P. and Nakazawa, M. (2012) Highly Dispersion-Tolerant 160 Gbaud Optical Nyquist pulse TDM Transmission over 525km. Optics Express, 20, 15001-15007.
 Bosco, G., Carena, A., Curri, V., Poggiolini, P. and Forghieri, F. (2010) Performance Limits of Nyquist-WDM and CO-OFDM in High-Speed PM-QPSK Systems. IEEE Photonics Technology Letters, 22, 1129-1131.
 Granot, E. (2018) Information Loss in Quantum Dynamics, Advanced Technologies of Quantum Key Distribution Sergiy Gnatyuk. IntechOpen, London.
 Granot, E. (2019) Derivation of Euler’s Formula and ζ(2k) Using the Nyquist-Shannon Sampling Theorem. American Journal of Signal Processing, 9, 1-5.
 Granot, E. and Marchewka, A. (2007) Universal Potential-Barrier Penetration by Initially Confined Wave Packets. Physical Review A, 76, Article ID: 012708.