The expression for can be readily applied to the probabilities and positions as defined above, resulting in the first term given by
And the second term given by the square of Eq. (4). The second term of is given by the square of Eq. (5), but the first term is more nuanced. This is owing to the complex number returned when acting the derivative operator twice on the probability. To overcome this, Fourier transforms were used to change Eq. (5) into the momentum basis which then allowed for the efficient calculation of .
Denoting as the momentum-space probability obtained through a 2-dimensional, non-uniform Fourier transform of the position space pseudo-wavefunction, Eq. (5) can be rewritten as,
Leading to the first term in the expression to be written as,
The FINUFFT python wrapper was used to take the Fourier transform using a type 3, 2d non-uniform FFT[1][2], and the minimum value in time of the uncertainty relation was found. Points in momentum space were sampled on and along with the two additional points () and ().
Figure 4 shows the position and momentum probabilities respectively in their own basis. An animation showing how these evolve in time for the different conditions is presented in Supplementary Material 2.
To compute the values reported in Table 2, the corresponding value was found for each subject, and these were used to calculate the group average reported here.
- ↑ Barnett AH, Magland J, Klinteberg LAF. A parallel nonuniform fast Fourier transform library based on an “Exponential of semicircle” kernel. SIAM J. Sci. Comput. 2019;41:C479–C504. doi: 10.1137/18M120885X. [CrossRef] [Google Scholar]
- ↑ Barnett, A. H. Aliasing error of the kernel in the nonuniform fast Fourier transform. arXiv:2001.09405 [math.NA] (2020).